Micro-Boson Stars: Compact Bosonic Objects

• Micro-boson stars are horizonless, self-gravitating bosonic condensates defined as the low-mass limit of solutions to the Einstein–Klein–Gordon system using a free massive scalar field.
• They exhibit distinct equilibrium structures, stability branches, and scaling laws, with variations arising from self-interactions and alternative formation channels.
• Observational signatures include unique gravitational lensing, electromagnetic imaging, and gravitational-wave patterns that differentiate them from black holes.

Searching arXiv for recent and foundational papers on micro-boson stars, mini-boson stars, and closely related compact boson-star models. Micro-boson stars are horizonless, self-gravitating bosonic condensates whose defining content depends on the literature in which the term is used. In the most common usage, they are the low-mass or “mini” limit of boson stars, namely solutions of the Einstein–Klein–Gordon system built from a free massive complex scalar field with no self-interaction beyond the mass term. In other contexts, the same expression denotes very small or microscopic bosonic compact objects, including ultralight-field condensates relevant to dark matter, sub-kilogram stars generated by attractive self-interactions, or picometer-scale Bose/Proca solitons produced from primordial black holes. Across these usages, the common structural elements are a localized bosonic field profile, harmonic time dependence, asymptotic flatness or a controlled nonrelativistic limit, and the absence of an event horizon (Shnir, 2022, Visinelli, 2021, March-Russell et al., 2022).

1. Terminology and conceptual scope

The term mini-boson star has a precise technical meaning: it is the solution obtained when the scalar potential is reduced to a pure mass term, $U(|\Phi|^2)=\mu^2|\Phi|^2$, with no quartic or higher self-interaction. Reviews of boson stars use this as the canonical free-field model and explicitly identify “micro-boson stars” with the smallest or low-mass members of that class (0801.0307, Shnir, 2022).

A second usage is scale-based rather than model-based. A review of boson stars and oscillatons defines micro-boson stars as gravitationally bound compact objects composed of a large number of ultralight bosons, with total mass much less than that of typical astrophysical objects, and notes that this class is also termed mini-boson stars in the literature (Visinelli, 2021). By contrast, a primordial-black-hole production scenario studies “micro-Bose/Proca” stars with radius $\lesssim {\rm pm}$ and mass $\sim 10^{\rm few}$ kg, emphasizing microscopic size and relic dark-matter phenomenology rather than the free-field Kaup limit (March-Russell et al., 2022).

A third usage appears in recent self-interacting models. In the Einstein–Friedberg–Lee–Sirlin framework, self-interaction is described as especially relevant for “micro-boson stars (smaller, denser stars),” where positive quartic terms enlarge the stable and compact sector (Sá et al., 24 Nov 2025). This suggests that “micro-boson star” functions less as a unique Lagrangian label than as a shorthand for the small-mass, high-density, or microscopic end of boson-star parameter space.

2. Field-theoretic formulation

The standard relativistic description is the Einstein–Klein–Gordon system for a complex scalar field minimally coupled to gravity,

$$\mathcal{S} = \int d^4x \sqrt{-g} \left[ \frac{R}{16\pi G} - \frac{1}{2} g^{\mu\nu} \left( \Phi_{,\mu}^* \Phi_{,\nu} + \Phi_{,\nu}^* \Phi_{,\mu} \right) - U(|\Phi|^2) \right] ,$$

with Einstein equations for the metric and a Klein–Gordon equation for the matter field (Shnir, 2022). The harmonic ansatz

$$\Phi(t,\mathbf{x}) = f(r,\theta,\varphi)e^{-i(\omega t+n\varphi)}$$

or, in spherical symmetry,

$\Phi({\bf r},t)=\phi(r)e^{i\omega t},$

renders the stress tensor stationary even though the scalar carries a time-dependent phase (Shnir, 2022, Choi et al., 2019).

Localization requires a bound-state regime. For a massive free field, the asymptotic scalar profile decays exponentially,

$$f(r) \sim \frac{1}{\sqrt{r}} \exp\!\left[- r \sqrt{\mu^2-\omega^2}\right],$$

and the frequency lies in the interval $0<\omega<\mu$ (Shnir, 2022). The global $U(1)$ symmetry yields a conserved Noether charge,

$Q = \int_{\Sigma} j^\mu n_\mu dV ,$

which is interpreted as particle number in many constructions (Shnir, 2022). Mass may be read from Komar-type expressions or asymptotically from the metric; for spherically symmetric stars this is often encoded through a mass function such as

$\lesssim {\rm pm}$0

or its equivalent coordinate realizations (Guerra et al., 2019, Sá et al., 24 Nov 2025).

In the nonrelativistic ultralight regime, the dynamics reduce to the Schrödinger–Poisson system,

$\lesssim {\rm pm}$1

which provides the solitonic background for perturbation theory, halo-core applications, and evaporation-rate calculations (Chan et al., 2023).

3. Equilibrium structure, scaling, and stability

For the free mini-boson-star model, the characteristic maximum mass is the Kaup limit,

$\lesssim {\rm pm}$2

and the radius obeys the inverse mass relation

$\lesssim {\rm pm}$3

with $\lesssim {\rm pm}$4 from numerical solutions (Visinelli, 2021). Reviews of the Einstein–Klein–Gordon model summarize the same scaling as $\lesssim {\rm pm}$5 and $\lesssim {\rm pm}$6, with the stable branch ending at the first maximum of the mass–frequency curve (Shnir, 2022). The compactness of mini-boson stars remains modest; at the Kaup limit one review quotes $\lesssim {\rm pm}$7 (Visinelli, 2021).

The global structure of solution families is organized by the usual spiral in $\lesssim {\rm pm}$8 or $\lesssim {\rm pm}$9. The first branch, extending from $\sim 10^{\rm few}$0 to the first critical point, is stable; beyond the maximum mass, the family enters an unstable spiral (Shnir, 2022). Catastrophe-theory analysis sharpens this statement: mini-boson stars possess only a single stable phase, whereas self-interacting boson stars can exhibit two stable regions separated by an unstable one, with an additional unstable spiralling sector near the black-hole limit (Kleihaus et al., 2011). Binding-energy criteria appear in equivalent forms,

$\sim 10^{\rm few}$1

or

$\sim 10^{\rm few}$2

with positive $\sim 10^{\rm few}$3 or negative $\sim 10^{\rm few}$4 identifying bound configurations on the stable branch (Kleihaus et al., 2011, Sá et al., 24 Nov 2025).

The simplest micro/mini sector is not universally stable once the spacetime framework is altered. In non-compact spacetime dimensions $\sim 10^{\rm few}$5, free mini boson stars are dynamically unstable: the perturbation eigenvalue $\sim 10^{\rm few}$6 is always negative, and the binding energy remains positive for all physically relevant backgrounds (Franzin, 2024). Radial and nonradial excitations add further structure. For ultralight solitons in the Schrödinger–Poisson regime, gravity mixes positive and negative frequencies, producing zero-energy monopole and dipole modes, non-standard orthogonality relations, and low-lying excitations that differ significantly from those of a Schrödinger equation in a fixed potential (Chan et al., 2023).

4. Self-interactions and alternative micro-boson-star regimes

Self-interaction changes both the scale and the internal taxonomy of micro-boson stars. A quartic repulsive interaction can raise the maximum mass from the Kaup scaling to

$\sim 10^{\rm few}$7

making astrophysical masses possible for boson masses that would otherwise generate only very small objects (Visinelli, 2021). Earlier reviews summarize the same transition through the Colpi–Shapiro–Wasserman scaling

$\sim 10^{\rm few}$8

and emphasize that self-interaction is what allows boson-star phenomenology to parallel neutron-star-like stable phases (0801.0307).

The Einstein–Friedberg–Lee–Sirlin model provides a two-scalar realization. Its matter sector contains a complex scalar with quartic self-interaction and a real scalar field, with potential

$\sim 10^{\rm few}$9

Positive $$\mathcal{S} = \int d^4x \sqrt{-g} \left[ \frac{R}{16\pi G} - \frac{1}{2} g^{\mu\nu} \left( \Phi_{,\mu}^* \Phi_{,\nu} + \Phi_{,\nu}^* \Phi_{,\mu} \right) - U(|\Phi|^2) \right] ,$$0 increases both the maximum mass and compactness, allows masses comparable to the Chandrasekhar limit without ultralight bosonic masses, and enlarges the bound sector. In the limit $$\mathcal{S} = \int d^4x \sqrt{-g} \left[ \frac{R}{16\pi G} - \frac{1}{2} g^{\mu\nu} \left( \Phi_{,\mu}^* \Phi_{,\nu} + \Phi_{,\nu}^* \Phi_{,\mu} \right) - U(|\Phi|^2) \right] ,$$1, where the real scalar becomes massless, the stars develop larger effective radii and a broader range of stable solutions (Sá et al., 24 Nov 2025).

Periodic axion potentials generate a different departure from the mini limit. As the axion decay constant $$\mathcal{S} = \int d^4x \sqrt{-g} \left[ \frac{R}{16\pi G} - \frac{1}{2} g^{\mu\nu} \left( \Phi_{,\mu}^* \Phi_{,\nu} + \Phi_{,\nu}^* \Phi_{,\mu} \right) - U(|\Phi|^2) \right] ,$$2 decreases, the mass–central-density curve develops multiple extrema, new stability branches emerge at high density, and some of the most compact configurations acquire a photon sphere. For GUT-scale QCD-axion parameters, the maximum mass can be up to ten solar masses, with compactness exceeding that of neutron stars (Guerra et al., 2019).

The sign of the self-interaction can also force the micro-boson-star regime toward genuinely microscopic masses. In self-interacting dark-matter models with attractive quartic coupling, the maximum mass scales as

$$\mathcal{S} = \int d^4x \sqrt{-g} \left[ \frac{R}{16\pi G} - \frac{1}{2} g^{\mu\nu} \left( \Phi_{,\mu}^* \Phi_{,\nu} + \Phi_{,\nu}^* \Phi_{,\mu} \right) - U(|\Phi|^2) \right] ,$$3

and the paper emphasizes that the resulting objects are “micro-boson stars,” with maximum masses below $$\mathcal{S} = \int d^4x \sqrt{-g} \left[ \frac{R}{16\pi G} - \frac{1}{2} g^{\mu\nu} \left( \Phi_{,\mu}^* \Phi_{,\nu} + \Phi_{,\nu}^* \Phi_{,\mu} \right) - U(|\Phi|^2) \right] ,$$4 kg in the relevant parameter space. For repulsive coupling, by contrast, masses up to $$\mathcal{S} = \int d^4x \sqrt{-g} \left[ \frac{R}{16\pi G} - \frac{1}{2} g^{\mu\nu} \left( \Phi_{,\mu}^* \Phi_{,\nu} + \Phi_{,\nu}^* \Phi_{,\mu} \right) - U(|\Phi|^2) \right] ,$$5 are possible for MeV dark matter (Eby et al., 2015).

5. Formation channels and dynamical extensions

Several distinct formation mechanisms have been studied. Reviews identify gravitational relaxation or gravitational cooling, primordial formation from early-universe overdensities, and condensation during halo evolution as standard pathways for diffuse bosonic condensates (Visinelli, 2021). These mechanisms readily produce nonrelativistic or weakly bound objects, but compact micro-boson-star formation is more model-dependent.

One microscopic production channel uses small primordial black holes with initial mass $$\mathcal{S} = \int d^4x \sqrt{-g} \left[ \frac{R}{16\pi G} - \frac{1}{2} g^{\mu\nu} \left( \Phi_{,\mu}^* \Phi_{,\nu} + \Phi_{,\nu}^* \Phi_{,\mu} \right) - U(|\Phi|^2) \right] ,$$6 kg. Simultaneous Hawking evaporation and black-hole superradiance can generate extremely dense gravitationally-bound dark-matter Bose or Proca soliton stars with radius $$\mathcal{S} = \int d^4x \sqrt{-g} \left[ \frac{R}{16\pi G} - \frac{1}{2} g^{\mu\nu} \left( \Phi_{,\mu}^* \Phi_{,\nu} + \Phi_{,\nu}^* \Phi_{,\mu} \right) - U(|\Phi|^2) \right] ,$$7 and mass $$\mathcal{S} = \int d^4x \sqrt{-g} \left[ \frac{R}{16\pi G} - \frac{1}{2} g^{\mu\nu} \left( \Phi_{,\mu}^* \Phi_{,\nu} + \Phi_{,\nu}^* \Phi_{,\mu} \right) - U(|\Phi|^2) \right] ,$$8 kg that survive after black-hole decay. In the vector case these relics can constitute $$\mathcal{S} = \int d^4x \sqrt{-g} \left[ \frac{R}{16\pi G} - \frac{1}{2} g^{\mu\nu} \left( \Phi_{,\mu}^* \Phi_{,\nu} + \Phi_{,\nu}^* \Phi_{,\mu} \right) - U(|\Phi|^2) \right] ,$$9 of the dark-matter density (March-Russell et al., 2022).

A separate cosmological mechanism employs two coupled scalar fields: a complex field that forms the star and a spatially homogeneous background field whose slow evolution increases the effective mass of the star-forming boson. Within an adiabatic approximation, a non-relativistic boson cloud can evolve into a compact boson star on cosmological timescales if the background variation is as large as the Planck scale. The same work stresses a limitation: the required initial states are not configurations that can be described by the well-studied Schrödinger–Poisson system (Miyauchi et al., 14 May 2025).

Micro-boson-star dynamics also include multi-centre and non-spherical equilibria. Fully nonlinear general relativity admits two non-spinning mini-boson stars in equilibrium when the scalar fields have a relative $$\Phi(t,\mathbf{x}) = f(r,\theta,\varphi)e^{-i(\omega t+n\varphi)}$$0 phase, so that short-range scalar repulsion balances long-range gravity; the same paper shows that no analogous flat-spacetime equilibrium exists (Cunha et al., 2022). Chains of mini-boson stars, with even or odd numbers of constituents along the symmetry axis, have likewise been constructed without introducing any self-interacting potential, broadening the free-field sector beyond isolated spherical stars (Sun et al., 2022).

6. Electromagnetic, lensing, and gravitational-wave signatures

Because micro-boson stars are compact yet horizonless, their observational signatures are organized around transparency, lensing, and the absence of horizon absorption. Early reviews emphasized gravitational redshift, micro-lensing, and the possibility that boson stars could contribute to the MACHO sector of dark matter, while later reviews noted that microlensing surveys such as EROS-2 and OGLE-IV constrain the dark-matter fraction in point-like micro-boson stars over part of the subsolar-mass range (0801.0307, Visinelli, 2021).

Fast-radio-burst lensing and compact-binary gravitational waves probe self-interacting micro-boson-star models more directly. For Liouville and logarithmic benchmark potentials, numerical analysis shows that the maximum mass and compactness depart from those of quartic boson stars, and that FRB lensing and GW measurements can discriminate between the two potentials. In the preferred parameter window, Liouville stars are more accessible to FRB lensing because their repulsive self-interactions yield higher compactness, whereas logarithmic stars are harder to lens but may still produce detectable GW signals in mergers (Choi et al., 2019).

Electromagnetic imaging can approach black-hole mimicry without eliminating distinguishing features. In self-interacting Einstein–Friedberg–Lee–Sirlin stars, backward ray-tracing shows strong gravitational lensing, extended shadows, bright photon rings, and Doppler-brightened disk emission for the most compact configurations. Yet the absence of an event horizon permits photons to cross the star and re-emerge, allowing central bright spots and multiple ring structures in optically thin disks; the disk can also extend nearly to the origin instead of terminating at a black-hole ISCO (Sá et al., 24 Nov 2025). Equilibrium binaries of two mini-boson stars produce a different pattern: lensing develops pronounced dipolar features as the system becomes relativistic, but no photon spheres or light rings are found on or outside the equatorial plane (Cunha et al., 2022).

Gravitational-wave structure is sensitive to multipoles. For rotating scalar and vector boson stars, dimensionless observables such as $$\Phi(t,\mathbf{x}) = f(r,\theta,\varphi)e^{-i(\omega t+n\varphi)}$$1, $$\Phi(t,\mathbf{x}) = f(r,\theta,\varphi)e^{-i(\omega t+n\varphi)}$$2, $$\Phi(t,\mathbf{x}) = f(r,\theta,\varphi)e^{-i(\omega t+n\varphi)}$$3, $$\Phi(t,\mathbf{x}) = f(r,\theta,\varphi)e^{-i(\omega t+n\varphi)}$$4, and compactness fall on family-dependent approximate universal surfaces. The analysis explicitly notes that these relations are scale-free once rendered dimensionless, so the same framework applies to ultralight micro-boson stars and can distinguish scalar from vectorial stars if multipoles are measured accurately enough (Adam et al., 2024).

Accretion provides a further discriminant in binaries. A study of Gaia BH1 compares Bondi–Michel accretion onto a Schwarzschild black hole and a non-rotating boson star, emphasizing that matter can accumulate and radiate from the deep interior region of the horizonless object. For Gaia BH1, black-hole accretion yields luminosities of $$\Phi(t,\mathbf{x}) = f(r,\theta,\varphi)e^{-i(\omega t+n\varphi)}$$5, whereas boson-star accretion can produce observable luminosities in the order of $$\Phi(t,\mathbf{x}) = f(r,\theta,\varphi)e^{-i(\omega t+n\varphi)}$$6 to $$\Phi(t,\mathbf{x}) = f(r,\theta,\varphi)e^{-i(\omega t+n\varphi)}$$7 (Passos et al., 13 Apr 2025).

Micro-boson stars therefore occupy a broad but coherent research domain. The canonical free-field mini-boson star remains the reference model for mass–radius scaling, stability theory, and perturbative dynamics, while self-interactions, coupled fields, primordial-black-hole production, and multi-centre equilibria extend the concept toward compact astrophysical objects and genuinely microscopic dark-matter relics. The principal misconception is that all such objects are interchangeable with black holes or with one another. The current literature instead shows a stratified taxonomy: free and self-interacting sectors differ sharply in maximum mass and stable branches, microscopic relic Bose/Proca stars need not share the same formation channel as ultralight-field solitons, and horizonless lensing or accretion signatures remain the decisive phenomenological separator across the class (Kleihaus et al., 2011, March-Russell et al., 2022, Sá et al., 24 Nov 2025).