Existence of nonlinearly scalarized black holes in Einstein-scalar-Gauss-Bonnet theory with polynomial couplings

This paper investigates the existence and stability of nonlinearly scalarized black holes in Einstein-scalar-Gauss-Bonnet theory with polynomial couplings, demonstrating that while quartic-only couplings stabilize Schwarzschild black holes, specific higher-order polynomial couplings induce instability and transition to scalarized states above certain pulse thresholds, with the resulting solution branches exhibiting universal features in the probe limit but distinct dependencies on coupling strength when backreaction is included.

The Gist

Imagine the universe as a giant, quiet ocean. In this ocean, black holes are like massive, deep whirlpools. For a long time, physicists believed these whirlpools were perfectly smooth and featureless—like a simple, dark sphere. This was the "No-Hair Theorem": a black hole has no "hair" (no extra features like a magnetic field or a scalar field) other than its mass, spin, and electric charge.

However, this paper explores a fascinating new idea: What if we can make these black holes grow "hair"?

The authors are studying a specific theory called Einstein-scalar-Gauss-Bonnet (EsGB) theory. Think of this theory as a new set of rules for how gravity works, where gravity is coupled to a mysterious "scalar field" (let's call it the "Ghost Field").

Here is the story of their discovery, broken down into simple concepts:

1. The Setup: The "Ghost" and the "Whirlpool"

In this theory, the "Ghost Field" usually sits quietly at zero. But under certain conditions, it can wake up and wrap around a black hole, giving it "hair."

The researchers wanted to see if they could force a black hole to grow this hair by hitting it with a "pulse" of energy (like a wave crashing into the whirlpool). They used different mathematical recipes (called coupling functions) to describe how the Ghost Field interacts with the black hole's gravity.

2. The Three Recipes (Coupling Functions)

The team tested three different "recipes" for this interaction:

• Recipe A (The Runaway): A simple formula where the interaction gets stronger and stronger the more hair you try to grow.
  • The Result: Disaster. If you push the black hole too hard, the hair doesn't just grow; it explodes. The math breaks down, and the black hole becomes unstable. It's like trying to inflate a balloon that has no limit to how big it can get—it just pops.
• Recipe B (The Balanced): A complex formula with two parts: one part that wants to grow the hair (like a spring pushing out), and a second part that acts as a brake (like a rubber band pulling back).
  • The Result: Success. When they hit the black hole with a strong enough pulse, the hair grows, but the "brake" kicks in before it explodes. The hair settles into a stable, comfortable size. The black hole now has a permanent "hairstyle."
• Recipe C (The Middle Ground): Similar to Recipe B but with a slightly different shape.
  • The Result: Success. It also works, creating stable hairy black holes, though the exact shape of the hair depends on the strength of the "brake."

3. The "Plateau" and the "Valley"

The authors used a clever trick to explain why Recipe B works. They imagined the interaction as a landscape (a topographical map).

• For the Runaway Recipe: The landscape is a cliff that goes down forever. If the scalar field (the hair) starts moving, it just slides down the cliff forever, accelerating until it crashes. This explains the "divergence" or explosion.
• For the Balanced Recipe: The landscape is a W-shaped valley.
  • The black hole starts at the top of the middle hill (zero hair).
  • When you hit it with a pulse, it slides down into the valley.
  • But the valley has steep walls on the sides. The hair slides down, hits the bottom, and gets "trapped" there. It can't go back up the hill easily, and it can't slide off the edge because the walls are too high.
  • This "trapping" is what creates the stable, hairy black hole. The "plateau" they saw in their computer simulations is just the hair sitting comfortably at the bottom of this valley.

4. The Threshold: How Hard to Hit?

One of the key findings is that you can't just give the black hole a tiny tap.

• If you tap it gently, the hair fades away, and the black hole returns to being smooth (no hair).
• You need to hit it with a specific threshold of energy (a strong enough pulse). Once you cross that line, the hair grows and stays. It's like pushing a heavy boulder over a small hill; if you don't push hard enough, it rolls back down. If you push hard enough, it rolls into the valley and stays there.

5. The "Probe" vs. The "Real Deal"

The researchers did two types of calculations:

1. The Probe Limit: They pretended the hair was so light it didn't affect the black hole's shape at all. This gave them a universal pattern that looked the same for all recipes.
2. The Backreaction: They let the hair actually weigh something and change the shape of the black hole. This is the "real deal."
   • They found that when the hair gets heavy, the number of possible "hairstyles" changes. Depending on how strong the "brake" (the $\beta$ term) is, the black hole can have two stable hairstyles or three different ones. It's like a chameleon that can change into different patterns depending on the environment.

The Big Picture

This paper is significant because it shows that nonlinear effects (where things interact in complex, non-straightforward ways) can stabilize black holes that would otherwise be unstable.

• Analogy: Imagine a tightrope walker.
  • In the old view, the tightrope was perfectly balanced, but any wind would knock them off.
  • In this new view, the tightrope walker has a stabilizing pole (the nonlinear term). If they wobble too much, the pole pushes them back to the center. They can now walk the rope even in a storm, as long as the wind isn't too strong.

In summary: The authors discovered that by using specific mathematical rules, they can force black holes to grow stable "hair" (scalar fields) if they are hit with a strong enough energy pulse. This hair settles into a stable state because the laws of physics create a "trap" that prevents the hair from growing out of control. This opens up new possibilities for understanding how black holes might look and behave in our universe.

Technical Summary

Here is a detailed technical summary of the paper "Existence of nonlinearly scalarized black holes in Einstein-scalar-Gauss-Bonnet theory with polynomial couplings."

1. Problem Statement

The paper investigates the phenomenon of nonlinear scalarization of Schwarzschild Black Holes (SBHs) within the framework of Einstein-scalar-Gauss-Bonnet (EsGB) theory.

• Context: While spontaneous scalarization is well-known in EsGB theory via tachyonic instability (requiring $\zeta''(0) \neq 0$), this study focuses on coupling functions where the second derivative vanishes at the origin ($\zeta''(0) = 0$). In such cases, SBHs are stable against linear perturbations.
• The Gap: Previous studies utilized exponential coupling functions (e.g., $\zeta \sim 1 - e^{-\kappa \phi^4}$) to demonstrate that SBHs can become unstable against nonlinear perturbations if the perturbation amplitude exceeds a threshold. This paper aims to explore whether polynomial coupling functions can support similar nonlinear scalarization and stable scalarized black hole (SBH) solutions, specifically looking for solutions with constant scalar fields $\phi_s \neq 0$.

2. Methodology

The authors employ a combination of time-domain numerical evolution and static boundary value problem solving:

• Theoretical Framework:

  • The action includes the Ricci scalar, a scalar field $\phi$, and a coupling to the Gauss-Bonnet term $R^2_{GB}$ via a function $\zeta(\phi)$.
  • Coupling Functions Tested:
    1. $\zeta_1(\phi) = \alpha\phi^4 - \beta\phi^8$
    2. $\zeta_2(\phi) = \alpha\phi^4 - \beta\phi^6$
    3. $\zeta_3(\phi) = \alpha\phi^4$ (Pure quartic)
  • These functions are chosen such that $\zeta'(0)=0$ and $\zeta''(0)=0$, ensuring no linear tachyonic instability.

• Time Evolution (Nonlinear Instability):

  • The Klein-Gordon equation is evolved on a fixed SBH background using a (1+1)-dimensional formulation with a tortoise coordinate.
  • Initial Conditions: Gaussian pulses with varying amplitudes $A$ are introduced.
  • Numerical Scheme: 7th-order finite differences for spatial derivatives and 4th-order Runge-Kutta for time integration.
  • Goal: Determine the threshold amplitude $A_{th}$ above which the SBH transitions to a scalarized state.

• Static Solutions (Construction of Branches):

  • The full coupled Einstein-scalar equations are solved for spherically symmetric metrics.
  • Limits: Solutions are found in the probe limit (neglecting backreaction on the metric) and the exact limit (including backreaction).
  • Regularity: Solutions must satisfy a regularity condition at the horizon ($\Delta > 0$) to avoid curvature singularities.

• Effective Potential Analysis:

  • The term $\zeta(\phi)R^2_{GB}$ is treated as an effective potential $V_{eff}$ to explain the dynamical behavior (stabilization vs. divergence).

3. Key Contributions and Results

A. Threshold Amplitudes and Stability

• Polynomial Couplings ($\zeta_1, \zeta_2$): The authors identify a critical threshold amplitude $A_{th}$.
  • For $\zeta_1$ ($\alpha\phi^4 - \beta\phi^8$), $A_{th} \approx 0.05$.
  • For $\zeta_2$ ($\alpha\phi^4 - \beta\phi^6$), $A_{th} \approx 0.1$.
  • Outcome: If $A > A_{th}$, the scalar field grows and stabilizes into a time-independent configuration (scalarized phase). If $A < A_{th}$, the perturbation decays.
• Pure Quartic Coupling ($\zeta_3 = \alpha\phi^4$):
  • No stable scalarized phase is found.
  • For small amplitudes, the field decays.
  • For large amplitudes, the field undergoes a "runaway instability," diverging in a short timescale.

B. Mechanism of Stabilization (Effective Potential)

• The authors introduce an effective potential $V_{eff} \propto -\zeta(\phi)$.
• For $\zeta_3$ (Pure Quartic): The potential is unbounded from below (inverted bowl). The scalar field cascades down the slope, causing divergence.
• For $\zeta_1$ and $\zeta_2$: The higher-order negative terms ($-\beta\phi^8$ or $-\beta\phi^6$) create a W-shaped potential well.
  • The potential has a local minimum at a finite value $\phi_s = (\alpha/2\beta)^{1/4}$ (for $\zeta_1$).
  • This "traps" the scalar field, balancing the attractive $\alpha\phi^4$ term with the repulsive higher-order term. This explains the "plateau" observed in time evolution and the existence of stable solutions.

C. Solution Branches

The structure of the solution branches depends heavily on the coupling strength $\beta$:

1. Probe Limit (No Backreaction):

   • $\zeta_3$: A single linear branch where $\phi_H \propto M$.
   • $\zeta_1, \zeta_2$: Two branches emerge: a lower linear branch (starting from $\phi=0$) and an upper branch (starting near $\phi_s$). They form a triangular domain of existence.

2. Exact Solutions (With Backreaction):

   • Small/Medium $\beta$ (e.g., $\beta = 25/8, 100/8$): The solution space splits into three branches:
     1. Lower Branch: Similar to the $\zeta_3$ case; unstable and terminates when the regularity parameter $\Delta \to 0$.
     2. Primary Branch: Corresponds to the time-evolved stable solutions. $\phi_H$ remains close to $\phi_s$.
     3. Intermediate Branch: Connects the primary branch to the lower branch but terminates due to a curvature singularity (divergence of the Kretschmann scalar) outside the horizon.
   • Large $\beta$ (e.g., $\beta = 1000/8$): The intermediate branch disappears. Only two branches remain (Lower and Primary), and the domain of existence becomes more restricted.

D. Comparison with Exponential Couplings

• The results for polynomial couplings are qualitatively similar to those using exponential couplings ($\zeta_{ex} \sim 1 - e^{-\kappa \phi^4}$).
• Difference: In exponential models, the primary branch is not bounded by a finite $\phi_s$ (the field can grow indefinitely but is suppressed exponentially in the equations). In polynomial models, the primary branch is strictly bounded by the constant $\phi_s$ determined by the coefficients $\alpha$ and $\beta$.

4. Significance

• New Mechanism for Scalarization: The paper demonstrates that nonlinear scalarization is not exclusive to exponential couplings; polynomial couplings can also support stable scalarized black holes provided they include higher-order terms that create a confining potential well.
• Role of Higher-Order Terms: It highlights the critical role of the $\beta$ term (e.g., $\phi^8$) in preventing runaway instabilities and enabling the formation of stable "hairy" black holes.
• Universality vs. Specificity: While the lower (linear) branch of solutions is universal across different coupling types, the structure of the non-linear branches (number of branches, existence of singularities) is highly sensitive to the specific form of the coupling function and the strength of the non-linear coefficient $\beta$.
• Physical Insight: The effective potential analysis provides a clear geometric explanation for the "plateau" in time evolution and the divergence in unstable cases, linking dynamical behavior directly to the shape of the scalar potential.

In conclusion, the study confirms that EsGB theory with polynomial couplings supports a rich landscape of nonlinearly scalarized black holes, governed by the interplay between the driving quartic term and the stabilizing higher-order terms.