Exotic Compact Objects (ECOs)

• Exotic Compact Objects (ECOs) are horizonless, ultracompact gravitational structures arising from modified or quantum gravity theories, mimicking black hole observables.
• They exhibit measurable departures such as nonzero tidal Love numbers, echo signatures in gravitational waves, and distinctive spectral features.
• Observational data from gravitational-wave detectors and electromagnetic imaging offer practical approaches to differentiate ECOs from classical black holes.

Exotic Compact Objects (ECOs) are ultracompact, horizonless solutions of (possibly extended or quantum-modified) gravity that mimic the macroscopic observables of black holes while exhibiting departures in their near-horizon structure. Such objects arise across a broad landscape of theoretical settings, including self-gravitating fundamental fields, modifications of the Standard Model or General Relativity, quantum gravitational completions, and are increasingly relevant in the interpretation of gravitational-wave and electromagnetic observations. ECOs admit Schwarzschild or Kerr-like exteriors down to surfaces or effective “walls” located at radii $r_0 = r_+ (1+\epsilon)$, with $\epsilon \ll 1$, and are characterized by nonzero tidal Love numbers, reflective or partially absorbing boundary conditions, nontrivial matter content, and often dynamically testable stability properties.

1. Principal ECO Classes and Theoretical Motivations

ECOs are grouped according to the new physics or exotic matter content responsible for their structure:

• Boson/Proca Stars: Composed of self-gravitating complex scalar (boson stars) or vector (Proca stars) fields, governed by the Einstein–Klein–Gordon/Proca equations. These solutions can accommodate dark-matter candidates, arise in axion or hidden sector extensions, and exhibit equilibrium by balancing gravitational collapse with field gradients or repulsive self-interactions. The scalar ansatz is $$\Phi(t,\vec{x}) = \Phi_0(r,\theta) e^{i(m\varphi-\omega t)}$$, parametrizing rotating configurations (Bezares et al., 2024).
• Hybrid Fermion–Boson Stars: Incorporate baryonic (fluid) and bosonic sectors, allowing configurations with coexistence of standard and dark sector matter, and governed by coupled Einstein–Klein–Gordon–Euler dynamics.
• Strange (quark) Stars: Supported by equations of state with deconfined $u,d,s$ quark matter (e.g., MIT bag model), possibly relevant in high-density QCD phases (Bezares et al., 2024).
• Ultracompact Anisotropic Stars, Gravastars, and Fuzzballs: Feature pressure anisotropy ($P_r \neq P_t$) to circumvent Buchdahl’s compactness limit, or replace the black hole interior with a de Sitter core matched at a thin shell (gravastars) or a string-theory-derived microstate structure (fuzzballs). Extensions of GR, such as Palatini $f(R)$ or higher-curvature corrections, facilitate further equilibrium scenarios.
• Hairy Black Holes and Superradiant Clouds: Comprise BHs with additional macroscopic scalar/vector "hair" enabled via synchronization at superradiant thresholds (e.g., $\omega = m\Omega_H$), evading standard no-hair theorems (Bezares et al., 2024).

Horizonless ECOs are also realized as wormholes, 2-2-holes of higher-derivative gravity, collapsed polymers, or alternatives invoking Planck-scale structure at the would-be horizon (Raidal et al., 2018).

2. Equilibrium Structures and Stability Analysis

The equilibrium of ECOs is defined by the solution of Einstein’s equations (or suitable extensions), possibly coupled to exotic matter fields. Representative frameworks include:

• Einstein–Klein–Gordon (EKG) System: Governs boson stars, leading to stationary solutions with compactness $C\equiv M/R$ spanning the range $0 < C < C_{\rm max}$ (the latter set by field mass and interaction strengths) (Bezares et al., 2024). Harmonic field ansätze reduce the coupled PDEs to ODEs for radial metric and field profiles.
• Modified TOV Equations: For self-bound matter (quark, anisotropic stars), equilibrium is modified by effective pressures and energy densities, possibly including field self-interactions or curvature effects.
• Boundary Conditions: For horizonless ECOs, the "surface" is placed at $r_0$, a proper distance $\delta$ outside the would-be horizon, with a boundary condition on the perturbation variable: $$u_{\omega l}(x\to x_0) \propto e^{-i\omega x} + K e^{-2i\omega x_0} e^{+i\omega x}$$, where $K$ encodes reflectivity/dissipation (Macedo et al., 2018).

Dynamical Stability:

• Linear radial perturbations determine stable and unstable solution branches; nonlinear numerical relativity simulations demonstrate that stable boson stars persist under large perturbations, while unstable branches disperse or collapse to black holes or migrate to less compact configurations (Bezares et al., 2024).
• Ultracompact Stars and Instability: Light rings (photon spheres) and ergoregions can trigger dynamical instabilities—ergoregion instabilities in spinning, horizonless objects limit allowed spins (see Sec. 5).

3. Gravitational-Wave Signatures and Observational Discriminants

ECOs generically feature departures from black-hole predictions in gravitational-wave and electromagnetic observables:

• Echoes: The post-merger ringdown of an ECO exhibits, after an initial BH-like decay (photon sphere QNMs), a train of "echoes" spaced by the round-trip time between the surface and the photon sphere ($\Delta t_{\rm echo} \sim 2R_s \ln(R_s/\epsilon)$) (Mark et al., 2017, Cardoso et al., 2016, Micchi et al., 2019). For Kerr-like objects, rotation further introduces spectral asymmetries in the echoes: rotation breaks the positive/negative frequency symmetry and introduces subdominant frequencies (Micchi et al., 2019).
• Tidal Love Numbers: Horizonless ECOs possess nonzero tidal deformability, parameterized by $k_2$ and $\Lambda$. Causality imposes a lower bound: physically motivated ECOs cannot reach the vanishing $\Lambda$ of black holes, leading to a "tidal gap" in the compactness–deformability plane (Russo et al., 22 Dec 2025). For a linear EoS at the causal limit, $k_{2,\min} \simeq 0.018$ and $\Lambda_{\min} \simeq 2.19$; solitonic boson stars further reduce $\Lambda_\text{min} \simeq 1.2$ via strong-energy-condition violation, but always yielding $\Lambda>0$.
• Resonant Absorption and Spectral Lines: Trapped modes between the ECO surface and the photon-sphere barrier generate Breit–Wigner "spectral lines" in the absorption cross section, with spacing and width set by the cavity length and reflectivity, absent in true BHs (Macedo et al., 2018).
• Resonant Excitations in Inspiral: As the inspiral frequency sweeps through ECO cavity resonances, sharp dephasing imprints ("kicks") arise in the GW phase evolution. While current datasets show no statistically significant resonances (Asali et al., 2020), next-generation detectors may reach subrad dephasing sensitivity and probe the underlying ECO structure.

A summary of key observable properties is provided in the following table:

Signature	BHs	ECOs (horizonless)	Observable Consequence
Echoes	No	Yes, if $\epsilon \ll 1$	Late-time GW pulse train
Tidal Love #'s	0	$\Lambda_\text{min} > 0$	Inspiral phase deviation
Spectral Lines	Absent	Present ("cavity" modes)	Absorption resonances
Superradiance	Present (absorbing)	Modified (reflective)	Altered instability regions

[See (Cardoso et al., 2016, Russo et al., 22 Dec 2025, Macedo et al., 2018, Zhou et al., 2023) for analytic and numerical demonstrations.]

4. Numerical Relativity and Waveform Modeling

Numerical-relativity simulations are pivotal for establishing the dynamical realizability and observational imprints of ECOs. Techniques include:

• Evolution Formalisms: BSSN or CCZ4 formulations for gravity; EKG or GRHydro modules for matter; conformal thin sandwich initial data for binary boson stars.
• Wave Extraction: The Newman–Penrose scalar $\Psi_4$ is decomposed into spin-weighted harmonics to extract GW strain, enabling detailed analysis of inspiral, merger, post-merger, and possible echo signals (Bezares et al., 2024).
• Phenomenological Models: Frequency-domain waveform families parameterized by compactness or deformability. For example, PhenomDECO models encode an early inspiral cutoff at $f_{\rm ECO,ISCO} \sim (2C)^{3/2} f_{\rm BBH,ISCO}$, with compactness $C<0.5$ signifying ECOs (Ghosh et al., 22 May 2025).
• Template Construction: Echo templates (Gaussian sum, Green's function) incorporate cavity delays and damping, with amplitude, width, and frequency drift tuned to numerical results (Wang et al., 2018, Mark et al., 2017).
• Stability Criteria: Dynamical NR studies identify the threshold for prompt collapse in head-on boson-star collisions and probe the excitation of nonaxisymmetric instabilities (Bezares et al., 2024).

5. Instabilities and Population Constraints

ECO models face theoretical and empirical constraints arising from intrinsic instabilities and population-level gravitational-wave data:

• Ergoregion Instability: Spinning horizonless ECOs with ergoregions are unstable to gravitational-wave emission unless the surface reflectivity $\mathcal{R}$ is less than unity. The critical spin above which the ergoregion instability develops is $$\chi_{\rm crit}(\epsilon) \simeq \pi(1+q)/(m |\log_{10}\epsilon|)$$ for azimuthal number $m$ and sector $q$ (Mastrogiovanni et al., 11 Feb 2025). LIGO–Virgo–KAGRA spin statistics exclude a pure-ultracompact ECO population at 90% credibility for $\epsilon<10^{-30}$, with $f_{\rm ECO}<28\%$ (polar) and $f_{\rm ECO}<25\%$ (axial). Next-generation detectors will push $f_{\rm ECO}$ constraints to a few percent with a day of observation.
• Superradiant Instabilities: Reflectivity at the ECO wall modifies the mass and spin parameter regions excluded by superradiant instabilities in the presence of ultralight bosons. The growth rate is rescaled by a factor $g_{\mathcal{K}}$ encoding the energy flux through the surface and interference effects (Zhou et al., 2023). For moderate reflectivity, exclusion regions shift minimally, but for $\mathcal{K}\to1$, they can shrink substantially.
• GW-Induced Collapse: Energy input from gravitational waves can cause ECOs to collapse through horizon formation unless the compactness parameter $\epsilon$ exceeds a threshold set by the GW fluence, requiring fine tuning or non-causal expansion to maintain echo observability (Chen et al., 2019).

6. Environmental Effects and Astrophysical Context

• Dark Matter Halos: ECOs embedded in realistic galactic or astrophysical environments experience backreaction effects. Surrounding dark matter alters the effective potential, shifting the light ring, echo time delays ($\Delta t_{\rm echo}$), and, to a lesser extent, tidal Love numbers (Chakravarti et al., 3 Sep 2025). In extended (nonrelativistic) DM halos, echo delays can lengthen by an order of magnitude; these effects are accessible to future detectors targeting millisecond timing.
• Accretion Flows and Imaging: General relativistic magnetohydrodynamic (GRMHD) simulations reveal distinctive signatures in plasma flow morphology, jet production, and radio imaging. For example, horizonless boson stars can stall MRI-driven accretion flows, yielding “shadow-like” depressions smaller than Kerr shadows. Imaging and timing observations (EHT, GRAVITY, LISA) may help distinguish ECOs from true black holes in certain scenarios (Olivares-Sánchez et al., 2024).
• Modified Gravity Extensions: Models such as Rastall gravity impact the echo time delay, absorption cross section, and QNM spectra. Observable echo periods provide a probe of underlying gravitational coupling constants (Sakti et al., 2021).

7. Quantum Gravity, Primordial ECOs, and Dark Matter

• Quantum Gravity Motivations: ECOs arise as low-energy manifestations of quantum gravitational effects (e.g., fuzzballs, firewalls) or string-theory microstates, and some subclasses (e.g., wormholes, 2-2-holes) regularize singularities via Planck-scale structure (Raidal et al., 2018).
• Primordial Dark Matter Candidates: Horizonless ECOs evade Hawking evaporation and offer new cosmologically viable windows for sublunar-mass dark matter ($10\,{\rm TeV} \lesssim M \lesssim 10^{-16} M_\odot$), with suppressed radiation rates compared to PBHs. Proposed searches include femto/picolensing, γ-ray background, and GW echoes (Raidal et al., 2018).

---

Extensive ongoing research is quantifying the stability, formation, merger dynamics, and observational diagnostics of ECOs across the gravitational and electromagnetic spectra. Gravitational-wave observatories (current and upcoming) are increasingly sensitive to late-time waveforms, tidal deformability, echoes, and other discriminants, thereby progressively narrowing the allowed parameter space for horizonless ultracompact objects vis-à-vis classical black holes.