Solitonic Boson Stars

• Solitonic boson stars are horizonless, self-gravitating configurations of a complex scalar field whose nonlinear self-interactions enable ultra-compact structures.
• They exhibit rich features such as thin-wall regimes, oscillatory scalar profiles, and multiple stability branches that can mimic black-hole phenomenology.
• Their dynamics include rotation, multipolar patterns, and distinct gravitational-wave and lensing signatures, offering promising probes of strong-field gravity.

Solitonic boson stars are horizonless, self-gravitating, localized solutions of the Einstein–Klein–Gordon system in which nonlinear scalar self-interactions are structurally essential. In the cited literature, the term is used for several closely related classes of models: sextic Q-ball-type potentials with false-vacuum structure, degenerate-vacuum scalar potentials that generate thin-wall configurations, and axion-inspired periodic potentials that yield very compact high-density branches (Collodel et al., 2022, Guerra et al., 2019). Across these realizations, the defining features are regularity, asymptotic flatness, harmonic time dependence of a complex scalar condensate, and compactness substantially larger than that of mini boson stars, with some families developing light rings and black-hole-mimicking strong-field phenomenology (Sukhov, 2024, Rosa et al., 3 Apr 2025).

1. Field-theoretic definition

A standard starting point is the Einstein–Klein–Gordon action

$$S=\int d^4x\,\sqrt{-g}\left[\frac{R}{2\kappa} - g^{ab}\partial_a\Phi^*\partial_b\Phi - V(|\Phi|^2)\right],$$

or equivalent normalizations thereof, with a complex scalar field and a $U(1)$-invariant potential. For static, spherically symmetric stars one typically uses

$$ds^2=-e^{v(r)}dt^2+e^{u(r)}dr^2+r^2d\Omega^2,\qquad \Phi(t,r)=\phi(r)e^{-i\omega t},$$

so that the scalar oscillates while the stress–energy tensor and geometry remain stationary (Guerra et al., 2019, Collodel et al., 2022).

The global $U(1)$ symmetry implies a conserved Noether current and charge, interpreted as particle number. This is one reason boson stars are classified as non-topological solitons rather than fluid stars: the matter source is a coherent condensate rather than a perfect-fluid equation of state, and the pressures are generically anisotropic (Visinelli, 2021, Collodel et al., 2022).

Relative to mini boson stars with $V=\mu^2|\Phi|^2$, the solitonic subclass is distinguished by self-interactions that can support localized configurations even in the weak-gravity or flat-space limit. In that sense, gravitating Q-balls, Q-stars, and several ultra-compact boson-star models lie in the same conceptual family (Bošković et al., 2021, Liebling et al., 2012).

2. Potentials and structural regimes

In the cited literature, “solitonic” does not denote a single unique potential. It denotes a family of strongly self-interacting scalar models that support compact, localized configurations.

Potential family	Representative form	Characteristic regime
Mini boson star	$V=\mu^2|\Phi|^2$	Centrally peaked, lower compactness
Sextic Q-ball type	$$V=\mu^2 |\Phi|^2 \left(1-\frac{2|\Phi|^2}{\sigma_0^2}\right)^2$$	Almost shell-like, thin-wall limit
Degenerate-vacuum sextic	$$U(|\Phi|)=\mu^2 \Phi^2 \left(1-\frac{\Phi^2}{\sigma^2}\right)^2$$	False-vacuum interior, sharp wall
Axion periodic	Periodic axion-inspired $V(\phi)$	Multiple branches, very compact ABS

For sextic and degenerate-vacuum models, the scalar field is nearly constant in the interior and then rapidly drops to the vacuum across a thin shell. The pseudo-spectral study of rotating and spherical solitonic boson stars describes these solutions as “almost shell-like”: the field is nearly constant inside up to a sharply located wall and becomes increasingly stiff as the shell thickness shrinks (Sukhov, 2024). The beyond-thin-wall analysis of static spherical solutions reaches the same conclusion and emphasizes that smaller $\sigma$ drives the system deeper into a thin-wall regime with strong numerical stiffness (Collodel et al., 2022).

The thin-wall approximation is useful but not exact. In the degenerate-vacuum sextic model, the idealized configuration with $U(1)$0 and $U(1)$1 is not an exact boson-star solution carrying nonzero Noether charge; the full Einstein–Klein–Gordon equations require a small but nonzero deviation from the false vacuum whenever $U(1)$2 (Collodel et al., 2022). This is one of the main technical corrections to older thin-shell intuition.

Axion boson stars extend the same theme to a periodic scalar potential. In the limit $U(1)$3, the periodic potential reduces to the free massive case, but as the decay constant decreases, strong self-interactions generate additional high-density branches and qualitatively new compact configurations (Guerra et al., 2019).

3. Compactness, stability, and ultracompact limits

Solitonic boson stars are studied through one-parameter equilibrium sequences labeled by a central field value or central amplitude. Stability changes are tracked by turning points in $U(1)$4 or analogous diagrams, as in neutron-star theory and ordinary boson-star theory. For mini boson stars there is typically one stable branch up to the maximum mass, whereas solitonic and axion-like models can exhibit several alternating stable and unstable branches (Guerra et al., 2019, Collodel et al., 2022).

A central result of the axion-boson-star analysis is the emergence of high-density stable branches with compactness $U(1)$5, approaching $U(1)$6 as $U(1)$7 decreases; within numerical reach, the stable solutions do not exceed that value (Guerra et al., 2019). In the broader self-gravitating Q-ball and soliton-boson-star analysis, the high-compactness limit is described as an effectively linear equation of state saturating the Buchdahl limit with the causality constraint (Bošković et al., 2021). These results place solitonic models among the few boson-star families that naturally enter the ultracompact-object regime.

Because scalar profiles have exponential tails, radius is necessarily operational. The cited literature uses enclosed-mass conventions such as $U(1)$8 or $U(1)$9, and compactness is then defined by $$ds^2=-e^{v(r)}dt^2+e^{u(r)}dr^2+r^2d\Omega^2,\qquad \Phi(t,r)=\phi(r)e^{-i\omega t},$$0 (Guerra et al., 2019, Rosa et al., 3 Apr 2025). Different radius definitions can shift quoted compactness values, which is why one dedicated study explicitly compared five radius prescriptions (Collodel et al., 2022).

The stability problem is richer than the usual bound-state heuristic. A 2025 nonlinear evolution study found perturbatively stable solitonic boson stars with positive binding energy that remain nonlinearly radially stable under explicit perturbations, challenging the conventional claim that negative binding energy is required for dynamical stability (Marks, 15 Aug 2025). Conversely, unstable models still follow the expected pattern: positive binding energy correlates with dispersal, while negative binding energy correlates with collapse or migration.

Photon spheres and light rings are another hallmark of the ultracompact regime. In axion boson stars, sufficiently small $$ds^2=-e^{v(r)}dt^2+e^{u(r)}dr^2+r^2d\Omega^2,\qquad \Phi(t,r)=\phi(r)e^{-i\omega t},$$1 produces two light rings—an inner stable one and an outer unstable one—despite the absence of an event horizon (Guerra et al., 2019). This feature is central to later lensing, imaging, and polarimetric studies.

4. Rotation, multipoles, and numerical construction

Rotating solitonic boson stars require stationary, axisymmetric solutions of a coupled elliptic system. A standard ansatz is

$$ds^2=-e^{v(r)}dt^2+e^{u(r)}dr^2+r^2d\Omega^2,\qquad \Phi(t,r)=\phi(r)e^{-i\omega t},$$2

with scalar field

$$ds^2=-e^{v(r)}dt^2+e^{u(r)}dr^2+r^2d\Omega^2,\qquad \Phi(t,r)=\phi(r)e^{-i\omega t},$$3

where $$ds^2=-e^{v(r)}dt^2+e^{u(r)}dr^2+r^2d\Omega^2,\qquad \Phi(t,r)=\phi(r)e^{-i\omega t},$$4 is the azimuthal quantum number. The angular momentum is quantized in the usual boson-star manner and the scalar profile becomes toroidal for rotating states (Sukhov, 2024, Gervalle, 2022).

The main technical obstacle is stiffness. In the ultra-compact or thin-shell regime, derivatives across the wall are very large over a narrow region, which makes shooting and conventional finite-difference schemes delicate. The multidomain pseudo-spectral solver introduced in 2024 addresses this by decomposing the computational region into Chebyshev collocation subdomains and concentrating resolution around the shell; the method was shown to display exponential convergence and to work especially well for almost shell-like solitonic solutions (Sukhov, 2024).

No comparably simple thin-shell analytic approximation exists for rotating solitonic boson stars. A survey of ultra-compact rotating SBS spacetimes therefore focused on direct numerical construction and on Geroch–Hansen multipole moments. It found that, for an intermediate range $$ds^2=-e^{v(r)}dt^2+e^{u(r)}dr^2+r^2d\Omega^2,\qquad \Phi(t,r)=\phi(r)e^{-i\omega t},$$5, the most compact stable rotating SBS can have reduced multipole moments close to those of Kerr, whereas for smaller $$ds^2=-e^{v(r)}dt^2+e^{u(r)}dr^2+r^2d\Omega^2,\qquad \Phi(t,r)=\phi(r)e^{-i\omega t},$$6 the deviations become larger (Sukhov, 2024). The same study tracked ergoregions and found that, in the most relevant range, stable ultra-compact configurations avoid or nearly avoid ergoregion formation.

Rotation also admits multi-solitonic generalizations. “Chains of rotating boson stars” are stationary multi-peak configurations whose even-constituent families show spiral-like dependence on frequency, mass, angular momentum, and Noether charge, while odd-constituent families develop loops that begin and end at the flat vacuum; this suggests that rotating chains are excitations of single rotating boson stars or of boson-star pairs (Gervalle, 2022).

5. Nonlinear dynamics, collisions, and mergers

Solitonic boson stars are dynamically active objects rather than purely equilibrium configurations. In the nonrelativistic limit, head-on collisions of self-interacting boson-star solitons can be understood with an effective potential built from kinetic, self-interaction, and gravitational terms. Relative phase, mass ratio, and self-coupling qualitatively control whether two solitons merge, rebound, pass through, or tidally disrupt (Cotner, 2016).

Full general-relativistic simulations of highly compact solitonic boson-star binaries showed that the post-merger gravitational-wave signal is mainly governed by the remnant’s fundamental frequency as it settles down to a non-rotating boson star, or to a black hole for sufficiently compact initial stars (Palenzuela et al., 2017). A notable outcome of those simulations is that mergers did not dynamically form rotating boson-star remnants; the post-merger state radiates significant gravitational energy while relaxing to a non-rotating configuration.

The scattering problem has recently been recast in an analytic effective-one-body language. For equal-mass solitonic boson stars, the conservative dynamics can be decomposed into point-mass gravitational, tidal, and short-range scalar-field interactions, with attractive behavior for in-phase binaries and repulsive behavior for out-of-phase binaries at small impact parameter; this yielded the first effective-one-body potential tailored to solitonic boson-star binaries (Damour et al., 30 Nov 2025).

A complementary strong-interaction problem is the piercing of a solitonic boson star by a black hole. In the hexic-potential study, the boson star was found to be easily captured by the black hole through an extreme tidal-capture process, and more than $$ds^2=-e^{v(r)}dt^2+e^{u(r)}dr^2+r^2d\Omega^2,\qquad \Phi(t,r)=\phi(r)e^{-i\omega t},$$7 of the boson-star material was accreted regardless of the black-hole initial mass and velocity (Zhong et al., 2023). In that setting, self-interaction increased compactness but did not change the qualitative outcome.

6. Astrophysical scales and gravitational-wave diagnostics

The physical mass and radius of a solitonic boson star depend on the underlying scalar mass and self-interaction scale. In axion boson stars, combining the dimensionless solutions with the QCD axion relation $$ds^2=-e^{v(r)}dt^2+e^{u(r)}dr^2+r^2d\Omega^2,\qquad \Phi(t,r)=\phi(r)e^{-i\omega t},$$8 yields stellar-mass or supermassive configurations depending on $$ds^2=-e^{v(r)}dt^2+e^{u(r)}dr^2+r^2d\Omega^2,\qquad \Phi(t,r)=\phi(r)e^{-i\omega t},$$9. For $U(1)$0, corresponding to $U(1)$1, the maximum axion-boson-star mass is approximately $U(1)$2, with compactness approaching the ultracompact regime (Guerra et al., 2019). The same analysis states that a boson star made of QCD axions with $U(1)$3 at the GUT scale can have a mass up to ten solar masses and be more compact than a neutron star.

Tidal deformability is one of the cleanest inspiral diagnostics. For spherically symmetric boson stars, the tidal-deformability study found a lower bound $U(1)$4 for quartic self-interaction models and $U(1)$5 for solitonic models (Sennett et al., 2017). The latter is much closer to the black-hole value $U(1)$6, which explains why solitonic boson stars are the most difficult boson-star class to separate from black holes using inspiral tides alone. The same study concluded that Advanced LIGO would struggle to distinguish solitonic boson stars from black holes using only inspiral tidal effects, whereas third-generation detectors should be able to separate binary black holes from binary solitonic boson stars (Sennett et al., 2017).

These gravitational-wave results fit the broader picture from merger simulations. Solitonic boson stars can look black-hole-like during inspiral and may collapse to black holes after merger, but the post-merger spectrum, tidal response, and absence of a horizon remain potential discriminants (Palenzuela et al., 2017).

7. Imaging, polarimetry, and black-hole-mimicker phenomenology

Because solitonic boson stars are horizonless, their lensing and imaging properties differ from those of black holes even when the exterior metric is very compact. One polarimetric hot-spot study found that the primary images of all analyzed solitonic boson-star models are comparable to those of Schwarzschild, and that distinguishable differences arise from secondary and plunge-through images rather than from the direct image itself (Rosa et al., 3 Apr 2025). In that framework, a low-inclination dilute solitonic model was found to be inconsistent with the observed QU-loop structure of Sgr A* flares, effectively excluding that dilute interpretation.

The absence of an event horizon implies that null geodesics and polarization vectors can probe the interior. In thin-disk ray-tracing calculations, polarization vectors were found to penetrate the stellar interior, unlike the black-hole case where no polarization signal exists within the horizon. The same study reported a positive correlation between polarization intensity and optical brightness, with the strongest polarization associated with the direct image, and emphasized that lensing images and photon rings become more prominent as the initial scalar field increases even though their polarization intensity remains weak (Zeng et al., 16 Aug 2025).

Near-infrared flare astrometry with GRAVITY provides a complementary constraint. A 2026 analysis of five spherically symmetric solitonic boson-star models found that fitted masses are systematically larger than the established $U(1)$7 value for Sgr A*, with more diffusive boson stars producing images closer to Schwarzschild and therefore mass estimates closer to the established value (Wang et al., 23 Apr 2026). The conclusion was restrictive but not absolute: GRAVITY flare astrometry places stringent constraints on solitonic-boson-star interpretations of Sgr A*, although it does not completely rule them out.

Taken together, these results clarify a common misconception. Solitonic boson stars can be exceptionally compact and can mimic several black-hole observables, but they are not simply horizonless copies of Schwarzschild or Kerr. Their thin walls, additional stable branches, possible light rings, finite tidal deformability, plunge-through images, and interior polarization signatures supply a set of discriminants that remain active targets for numerical relativity, strong-field imaging, and multimessenger inference (Sukhov, 2024, Rosa et al., 3 Apr 2025).