Spacetime of rotating black holes surrounded by massive scalar charges

This paper presents accurate spectral methods to construct the spacetime of rotating black holes surrounded by massive scalar fields with nonminimal couplings, enabling the calculation of horizon properties and paving the way for testing fundamental scalar degrees of freedom through electromagnetic and gravitational-wave observations.

The Gist

Imagine the universe as a giant, invisible trampoline. In our standard understanding of physics (General Relativity), if you place a heavy bowling ball (a black hole) in the center, the trampoline curves smoothly around it. If you spin that ball, the fabric twists and drags along with it. This is the "Kerr" black hole, the standard model we use today.

However, this paper explores a more complex scenario: What if the trampoline isn't just empty space, but is covered in a thick, invisible "fog" or "cloud" of heavy particles? And what if the rules of how the trampoline bends are slightly different from the standard rules?

Here is a simple breakdown of what the author, Adrian Ka-Wai Chung, actually did and found:

1. The Setup: Spinning Black Holes in a "Fog"

The paper looks at spinning black holes surrounded by a specific type of "fog" called massive scalar fields.

• The Fog: Think of this as a cloud of invisible particles that have weight (mass). In some theories of physics, these particles could be the "dark matter" that holds galaxies together, or they could be a side effect of a deeper theory of gravity.
• The Twist: These particles don't just sit there; they interact with the curvature of space itself. The paper studies three specific ways they might interact (called axi-dilaton, dynamical Chern–Simons, and scalar Gauss–Bonnet couplings).
• The Goal: The author wanted to build a precise mathematical map (a "spacetime") of what this spinning black hole looks like when it is wrapped in this heavy fog.

2. The Challenge: The "Stiff" Problem

Building this map is incredibly hard.

• The Analogy: Imagine trying to draw a picture of a cloud that is both swirling around a spinning top and shrinking exponentially fast as you move away from it.
• The Problem: Because these particles have mass, they die out very quickly as you move away from the black hole (like a flashlight beam that gets dimmer and dimmer). Standard mathematical tools (spectral methods) usually struggle with things that change this rapidly. It's like trying to take a high-resolution photo of a fast-moving object with a slow camera; the image gets blurry or "unstable."

3. The Solution: A New Mathematical "Lens"

The author developed a clever new way to use spectral methods (a type of high-precision math tool) to solve this.

• The Trick: Instead of trying to draw the whole cloud directly, the author mathematically "peeled off" the part that shrinks so fast (the exponential decay). They then focused on drawing the remaining "core" of the cloud, which is much smoother and easier to map.
• The Result: This allowed them to create a highly accurate map of the spacetime around the black hole, even when the "fog" is very heavy and shrinks very quickly. They tested this on black holes spinning at up to 80% of the maximum speed allowed by physics.

4. What They Found: The Shape of the Fog

When they looked at the maps they built, they discovered some interesting things about the "fog":

• The Shape Doesn't Change Much: Even though the particles are heavy, the overall shape of the cloud (whether it looks like a dipole or a quadrupole) stays very similar to what we see with massless particles. The mass mainly just makes the cloud shrink faster and become smaller in size.
• The Black Hole Changes: The presence of this heavy fog does change the black hole itself, but only slightly.
  • Spin: The fog makes the black hole spin a tiny bit slower (in some theories) or changes its spin speed in a specific pattern (in others).
  • Surface Heat: The "surface gravity" (which relates to the heat or temperature of the black hole's edge) changes slightly. In some theories, the black hole gets a tiny bit "hotter" or "colder" depending on how fast it spins.

5. Why This Matters (According to the Paper)

The paper claims these results are a "blueprint" for future detective work.

• The Blueprint: By having an accurate map of what the spacetime looks like with this "fog," scientists can now predict exactly how these black holes would behave if we could see them.
• The Tools: The author mentions two specific ways this map will be used:
  1. Gravitational Waves: When black holes crash into each other, they send out ripples in space (gravitational waves). If a black hole has this "fog" around it, the ripples will sound slightly different. This map helps scientists listen for those specific sounds.
  2. Black Hole "Ringdown": After a black hole is hit, it "rings" like a bell. The pitch of that ring depends on the black hole's spin and surface gravity. The author is currently using their map to calculate exactly what that "ring" sounds like for these heavy-fog black holes.

Summary

In short, the author built a high-precision mathematical model of a spinning black hole surrounded by a heavy, invisible cloud of particles. They found a clever math trick to handle the cloud's rapid shrinking, proved that the cloud changes the black hole's spin and "temperature" slightly, and provided the necessary data to help future telescopes and gravitational wave detectors search for these mysterious particles in the real universe.

Technical Summary

Technical Summary: Spacetime of Rotating Black Holes Surrounded by Massive Scalar Charges

Problem Statement
General Relativity (GR) faces challenges in explaining phenomena such as galactic rotation curves and cosmic acceleration without invoking dark matter or dark energy. Extensions to GR, including axi-dilaton, dynamical Chern–Simons (dCS), and scalar Gauss–Bonnet (sGB) gravity, introduce additional scalar degrees of freedom that couple nonminimally to spacetime curvature. Furthermore, massive scalar fields are well-motivated dark matter candidates capable of forming macroscopic condensates around black holes via superradiance.

While the spacetime of rotating black holes with massless scalar charges has been constructed using spectral methods, the presence of a non-zero scalar mass introduces significant theoretical and numerical challenges. The massive scalar field obeys a Klein–Gordon equation with an exponentially decaying asymptotic behavior ($e^{-\mu r}$), which complicates the construction of accurate solutions for the coupled, nonlinear system of field equations. Existing methods struggle to resolve these fields accurately, particularly for scalar masses where the Compton wavelength is comparable to the black hole mass. This paper addresses the need to construct the spacetime of rotating black holes surrounded by massive scalar charges within the small-coupling regime of dCS, sGB, and axi-dilaton gravity.

Methodology
The authors employ a perturbative approach in the dimensionless coupling parameter $\zeta$. The metric and scalar field are expanded as $g_{\mu\nu} = g^{(0)}_{\mu\nu} + \zeta g^{(1)}_{\mu\nu}$ and $\vartheta = \vartheta^{(0)} + \zeta \vartheta^{(1)}$, respectively. The background metric $g^{(0)}_{\mu\nu}$ is fixed as the Kerr spacetime.

1. Spectral Construction of Scalar Fields:

   • To handle the exponential decay of the massive scalar field, the authors define an auxiliary field $\phi$ such that $\vartheta = e^{-\mu r}\phi$. This isolates the stiff radial dependence analytically, leaving $\phi$ as a smoother function suitable for spectral expansion.
   • The radial coordinate $r \in (r_+, \infty)$ is compactified to $z \in [-1, 1]$ via $z = 2r_+/r - 1$.
   • The Klein–Gordon equation is transformed into a system of linear algebraic equations by expanding $\phi$ in Chebyshev polynomials (radial) and Legendre polynomials (angular).
   • Crucially, the exponential factor $e^{-\mu r}$ is retained explicitly on the derivative side of the differential operator during the spectral projection to ensure numerical stability and convergence.
   • Boundary conditions are enforced: the field vanishes at spatial infinity, and regularity is maintained at the event horizon. "Bare terms" (unphysical constant terms at infinity) are subtracted at each spectral order to enforce correct asymptotic behavior.

2. Spectral Construction of Metric Modifications:

   • Using the optimal scalar field solutions, the modified Einstein equations are solved for the metric deformations $H_i(r, \chi)$.
   • The source terms in these equations depend on the scalar field and its derivatives. The authors retain exponential factors ($e^{-2\mu r}$ and $e^{-4\mu r}$) in the source terms to maintain stability.
   • A similar spectral expansion is applied to the metric functions, incorporating boundary conditions that preserve asymptotic flatness and the definitions of ADM mass and angular momentum.

3. Diagnostics and Optimization:

   • Convergence is assessed using the "backward modulus difference" (BMD) and the absolute error of the field equations (residuals).
   • The authors identify a "twist spectral order" ($N_{twist}$) and an optimal spectral order ($N_{opt}$) where the error reaches a minimum. Beyond $N_{opt}$, increasing the spectral order leads to over-resolution and potential numerical instability, particularly for larger scalar masses.
   • Solutions are selected at $N_{opt}$ to minimize error propagation in subsequent physical calculations.

Key Contributions and Results

• Accurate Spacetime Construction: The authors successfully construct spacetimes for rotating black holes with dimensionless spins up to $a \le 0.8$ surrounded by massive scalar fields in dCS, sGB, and axi-dilaton gravity.
• Resolution of Massive Fields: The method accurately resolves scalar fields with Compton wavelengths as short as 5 times the black hole mass ($\mu \approx 0.2$ in geometric units).
  • Scalar field residuals are $\lesssim 10^{-5}$.
  • Metric modification residuals are $\lesssim 10^{-3}$ (specifically $\sim 10^{-3}$ at $a=0.8$).
• Impact of Scalar Mass:
  • Increasing the scalar mass $\mu$ leads to more rapid radial decay of the field but does not significantly alter the multipolar structure (dipolar for dCS, quadrupolar for sGB) of the scalar field or the geometric structure of the metric deformations.
  • However, higher masses increase the numerical difficulty, reducing the overall accuracy of the spectral method due to the sharper radial features introduced by the exponential factor.
• Physical Observables:
  • Horizon Angular Velocity ($\Omega_H$): The massive scalar field generally decreases $\Omega_H$ in dCS gravity. In sGB gravity, it initially increases $\Omega_H$ for low spins ($a \lesssim 0.7$) before decreasing it at higher spins.
  • Surface Gravity ($\kappa$): The presence of massive scalar charges generally increases the surface gravity relative to the GR value. In the non-rotating limit ($a \to 0$), $\kappa^{(1)} \approx 0$ for dCS but $\kappa^{(1)} \sim 10^{-1}$ for sGB.
  • Consistency Checks: The $L_2$ norms of the derivatives of $\Omega_H$ and $\kappa$ with respect to the polar angle $\chi$ are computed. These values remain small ($\lesssim 10^{-2}$ for $\Omega_H$ and $\lesssim 10^{-1}$ for $\kappa$), confirming that the constructed horizons remain approximately Killing horizons and validating the internal consistency of the solutions up to $a=0.8$.

Significance
The paper claims that these results pave the way for incorporating massive scalar charges into electromagnetic observations and gravitational-wave detections. Specifically:

• The constructed spacetimes provide the necessary background for computing quasinormal mode spectra, which could be used to search for massive scalar charges via black-hole ringdown spectroscopy (using tools like METRICS).
• The metric modifications can be integrated into waveform models for extreme-mass-ratio inspirals (EMRIs), potentially enabling future space-based detectors like LISA to probe fundamental scalar degrees of freedom.
• The work extends existing spectral approaches to gravitational systems with massive degrees of freedom, offering a systematic procedure for selecting the most accurate solutions based on error diagnostics.

The authors note that while the method is robust for $a \le 0.8$ and $\mu \le 0.2$, further refinements (such as alternative spectral bases or analytic projections of source terms) may be required for higher spins and larger scalar masses. The scalar field and metric solutions are provided as supplementary material to facilitate further research.