Spacetime of rotating black holes surrounded by massive scalar charges

Abstract: Massive scalar charges are ubiquitous in extensions to General Relativity and the Standard Model in particle physics. We describe spectral methods which can accurately construct the spacetime of rotating black holes with dimensionless spin up to $a \leq 0.8$ surrounded by massive scalar fields nonminimally coupled to spacetime curvature. We consider axi dilaton, dynamical Chern Simons, and scalar Gauss Bonnet couplings, and obtain leading order solutions for both the scalar field and the associated metric modifications. Our method accurately resolves massive scalar fields with Compton wavelengths as short as 5 times the black hole mass, achieving residual errors $\lesssim 10<sup>{-5}$, and yields the corresponding leading order spacetime modifications with residual errors $\lesssim 10<sup>{-3}$. Using the constructed spacetimes, we computes the leading-order shifts in the surface gravity and the angular velocity of the event horizon, important information for computing the quasinormal modes. These results pave the way to incorporate massive scalar charges into electromagnetic observations and gravitational-wave detections of black holes, potentially enabling new probes of fundamental scalar degrees of freedom.

• The paper develops leading-order spectral solutions for massive scalar fields and Kerr-spacetime deformations in dynamical Chern–Simons, scalar Gauss–Bonnet, and axi-dilaton gravity for spins up to 0.8 and masses up to 0.2/M.
• Massive fields preserve the characteristic angular structure of the scalar hair—dipolar in dCS and quadrupolar in sGB—while increasing mass shortens radial decay and reduces deformation amplitudes.
• The constructed horizon quantities support quasinormal-mode and inspiral modeling, with scalar errors below about 10⁻⁵ and metric errors around 10⁻³–10⁻⁴ within the explored parameter range.

This paper develops spectral methods to construct, at leading order in a small coupling parameter $\zeta$, the spacetime of rotating Kerr black holes surrounded by massive scalar fields nonminimally coupled to curvature. The framework covers axi-dilaton gravity (the sum of dynamical Chern–Simons and scalar Gauss–Bonnet contributions), with dimensionless spins up to $a \leq 0.8$ and scalar masses up to $\mu \leq 0.2/M$, corresponding to Compton wavelengths as short as five times the black-hole mass (2602.10462).

Theoretical setup

The author works from a unified Lagrangian density containing two quadratic-gravity couplings $\ell_{1,2}$, a scalar $\vartheta_1$ coupled to the Gauss–Bonnet invariant $\mathscr{G}$ and a pseudoscalar $\vartheta_2$ coupled to the Pontryagin invariant $\mathscr{P}$, each with mass parameters $\mu_{1,2}$. Setting $\ell_2=0$ recovers scalar Gauss–Bonnet (sGB) gravity; $a \leq 0.8$0, $a \leq 0.8$1 gives dynamical Chern–Simons (dCS) gravity; $a \leq 0.8$2, $a \leq 0.8$3 gives axi-dilaton gravity. The rescaled fields $a \leq 0.8$4 and the small-coupling expansion $a \leq 0.8$5 reduce the problem to solving a Klein–Gordon equation on a fixed Kerr background, followed by linearized modified Einstein equations sourced by curvature-scalar terms ($a \leq 0.8$6) and the trace-reversed stress-energy tensor ($a \leq 0.8$7). The background is taken to be exactly Kerr; this is justified only in the small-coupling regime, which existing observational constraints support.

Asymptotic analysis fixes the boundary structure: at spatial infinity the source decays as $a \leq 0.8$8 and the field behaves as $a \leq 0.8$9, while near the horizon regularity excludes the divergent Bessel-$\mu \leq 0.2/M$0 branch of the homogeneous solution. This motivates factoring out the exponential decay via $\mu \leq 0.2/M$1, leaving an auxiliary field $\mu \leq 0.2/M$2 that is finite and differentiable on the whole domain and hence amenable to spectral representation.

Spectral construction of massive scalar fields

The key numerical device is retaining the exponential factor explicitly in the differential operator rather than absorbing it into the unknown. After compactifying the radial coordinate via $\mu \leq 0.2/M$3, the equation takes the schematic form $\mu \leq 0.2/M$4. The field is expanded in Chebyshev polynomials $\mu \leq 0.2/M$5 times even Legendre polynomials for sGB (even parity) or odd Legendre polynomials for dCS (odd parity), projected onto the same basis, and solved as a least-squares system after appending the vanishing condition at spatial infinity. The projection integrals $\mu \leq 0.2/M$6, involving $\mu \leq 0.2/M$7 weights, are evaluated numerically at very high precision (working precision 700), and spurious "bare" non-decaying terms are subtracted at each spectral order.

Two diagnostics characterize convergence: the backward modulus difference (BMD) between successive spectral orders, and an absolute error measuring the residual of the Klein–Gordon equation under a regularization weight chosen so the residual integral can be evaluated analytically via exponential-integral functions. Both decrease exponentially with spectral order $\mu \leq 0.2/M$8 up to a characteristic order $\mu \leq 0.2/M$9, beyond which convergence slows; for $\ell_{1,2}$0 they reach a minimum at $\ell_{1,2}$1 before degrading. The paper explains this through the Compton wavelength partitioning the domain into an inner region where the exponential varies slowly and an outer region dominated by mass screening, so that once the Chebyshev resolution scale drops below the effective distance to infinity in $\ell_{1,2}$2-space, uniform improvement ceases. The optimal solution is selected by minimizing the error, and the least error shows little dependence on spin but increases noticeably from $\ell_{1,2}$3 to $\ell_{1,2}$4, reflecting the stiffness introduced by $\ell_{1,2}$5.

A central physical finding is that increasing the scalar mass does not significantly alter the multipolar geometry of the field: dCS solutions retain their dipolar structure antisymmetric about the equatorial nodal plane, sGB solutions remain quadrupolar and equatorially symmetric, consistent with the massless results. The mass changes only the radial decay rate and overall magnitude. This implies that mass-dependent signatures enter observables primarily through amplitude and horizon quantities rather than through the angular structure of the hair.

Spacetime modifications

The metric ansatz writes the deformation as four functions $\ell_{1,2}$6 multiplying corrections to the $\ell_{1,2}$7, $\ell_{1,2}$8, and spatial components, with boundary conditions enforcing asymptotic flatness while preserving the ADM interpretation of $\ell_{1,2}$9 and $\vartheta_1$0. Six independent components of the linearized modified Einstein equations are solved simultaneously. Two features distinguish this implementation from prior work: both even and odd Legendre polynomials are retained (required by the odd-parity $\vartheta_1$1 component), and the boundary conditions are incorporated directly into an augmented overdetermined system solved by least squares rather than used to eliminate coefficients beforehand.

Convergence behavior mirrors the scalar case, with one notable difference: the transition spectral order for the metric residual is roughly half that of the scalar residual, because the stress-energy source scales as $\vartheta_1$2 and is therefore twice as stiff radially. Nonphysical oscillations appear near the poles for $\vartheta_1$3, signaling over-resolution, so the analysis is restricted accordingly. The least error grows systematically with spin, reaching $\vartheta_1$4 at $\vartheta_1$5 — the stated reason for capping the spin there — while remaining approximately constant across $\vartheta_1$6. As with the scalar field, the geometric structure of $\vartheta_1$7 is essentially unchanged by the mass; heavier fields simply produce smaller deformations because their energy density decays faster.

Horizon quantities

From the constructed spacetimes, the leading-order corrections to the horizon angular velocity $\vartheta_1$8 and surface gravity $\vartheta_1$9 are computed. In dCS gravity the massive scalar decreases $\mathscr{G}$0 over the full range $\mathscr{G}$1; in sGB it first increases $\mathscr{G}$2 for $\mathscr{G}$3 and then decreases it. For the surface gravity, both theories show a general increase relative to GR, with $\mathscr{G}$4 in the nonrotating limit for dCS — consistent with Schwarzschild remaining a solution there — versus $\mathscr{G}$5 for sGB, where Schwarzschild is not a solution. These trends match the massless limits.

Consistency checks exploit black-hole mechanics: both $\mathscr{G}$6 and $\mathscr{G}$7 must be constant over the Killing horizon, so the $\mathscr{G}$8 norms of their $\mathscr{G}$9-derivatives should vanish identically. The computed norms increase monotonically with spin but stay below $\vartheta_2$0 for $\vartheta_2$1 and $\vartheta_2$2 for $\vartheta_2$3, i.e., at most about 10% relative variation even in the worst cases, validating the spacetimes up to $\vartheta_2$4. Since these quantities feed directly into quasinormal-mode computations via the METRICS framework, the results provide the necessary inputs for ringdown searches for massive scalar charges.

Limitations and open questions

The paper states its accuracy limits plainly. The achievable precision degrades at larger scalar masses because the exponential factor remains attached to the differential operator even after being factored out of the unknown, limiting spectral convergence; the reported errors of $\vartheta_2$5 for the scalar field and $\vartheta_2$6–$\vartheta_2$7 for the metric apply only within the explored parameter space. Higher spins are not treated: the implementation does not achieve sufficient accuracy beyond $\vartheta_2$8. Two refinements are proposed but left open — projecting the source onto associated Legendre harmonics to obtain a more diagonal radial-mode system, and adopting generalized Laguerre bases whose weight function matches the exponential stiffness. Whether these extend the reach to larger $\vartheta_2$9 and $\mathscr{P}$0 remains unanswered. Additionally, all results are strictly leading order in $\mathscr{P}$1 on an exact Kerr background, so second-order backreaction is outside the scope.

Conclusion

The paper extends established spectral techniques for massless scalar hair to massive scalar charges in dCS, sGB, and axi-dilaton gravity, delivering accurate leading-order scalar profiles, metric deformations, and horizon quantities for spins up to 0.8 and Compton wavelengths down to $\mathscr{P}$2. Its main structural conclusion — that finite mass leaves the multipolar geometry intact while modifying magnitudes and horizon properties — combined with the published solutions, enables incorporation of massive scalar charges into extreme-mass-ratio-inspiral waveform models and quasinormal-mode spectra, opening concrete routes to constrain fundamental scalars with LISA and ground-based ringdown observations.