Scalar Tsunamis from Black Hole Formation

This paper refines estimates of energy released as "scalar tsunamis" when light scalar fields surrounding collapsing stars are liberated upon black hole formation, demonstrating that while the total energy remains comparable to flat-space predictions, general relativity and improved initial modeling significantly alter the resulting energy spectrum.

The Gist

The Big Picture: The Cosmic "Pop"

Imagine a star is like a giant, heavy sponge soaked in a special, invisible liquid (a "scalar field"). This liquid is held in place because the star's matter is constantly "gluing" it there.

Now, imagine that star suddenly collapses into a black hole (like a supernova explosion or two neutron stars crashing together). When the star collapses, the "glue" disappears instantly. The invisible liquid, which was previously stuck to the star, is suddenly released. It rushes outward in a massive wave, like a tsunami.

The authors of this paper wanted to know: What happens to this "liquid wave" when it tries to escape a black hole?

The Old Idea vs. The New Idea

• The Old Idea (Flat Space): Previous scientists imagined the universe was empty and flat, like a calm pond. They thought that when the star vanished, the wave would split perfectly in half: 50% would rush inward and get sucked into the black hole, and 50% would rush outward and travel to Earth.
• The New Idea (Curved Space): This paper says, "Wait, the universe isn't flat near a black hole; it's curved and warped." The black hole acts like a giant, invisible hill or a bumpy wall. The authors used complex math and computer simulations to see how this "bumpy wall" changes the wave.

The Key Findings

1. The "Split" is Still Roughly 50/50

Even with the black hole's gravity warping space, the total amount of energy that escapes is surprisingly close to the old guess.

• The Analogy: Imagine throwing a ball at a trampoline with a hole in the middle. You might think the ball would either fall in or bounce out. The authors found that, generally, about half the energy falls in, and half escapes.
• The Twist: If the "sponge" (the star) was very large compared to the black hole, actually more than half might escape. This is because the black hole's "bumpy wall" (gravity) acts like a mirror for slow-moving waves, bouncing them back out instead of letting them fall in.

2. The Wave Changes Shape (The "Redshift")

While the amount of energy is similar, the type of wave changes significantly.

• The Analogy: Think of a siren on a passing ambulance. As it moves away, the pitch drops (it sounds deeper). This is the "Doppler effect."
• The Paper's Claim: The black hole's gravity does something similar. It stretches the waves, making them "lower pitched" (lower frequency) than scientists previously thought.
• Why it matters: If we are building detectors on Earth to listen for these waves, we need to know what "note" to listen for. If we listen for a high-pitched squeak, we might miss the signal because the black hole turned it into a low rumble.

3. The "Hair" Problem

There is a famous rule in physics called the "No-Hair Theorem," which says black holes are simple: they only have mass, spin, and charge. They shouldn't have any "hair" (extra messy fields sticking out).

• The Paper's Explanation: The authors show that while the field seems to stay stuck near the black hole for a long time, it's actually slowly leaking away or getting "swallowed" by the black hole growing slightly larger. Eventually, the black hole "eats" its own hair, and the field disappears, keeping the "No-Hair" rule intact.

The "Tsunami" Scenarios

The authors tested different shapes for the initial "sponge" to see how the wave behaves:

• The Uniform Sponge: If the field was spread out evenly, the wave behaves predictably.
• The Clumped Sponge: If the field was bunched up tightly near the star, the wave behaves differently, with more energy getting reflected back out by the gravity "wall."
• The Collapsing Sponge: They also simulated a star that was shrinking before it became a black hole. They found that even if the star was moving while it collapsed, the final result (the wave escaping) wasn't too different from the static case. The main change was a small "dip" in the wave pattern, but the overall tsunami still happened.

The Conclusion

The paper concludes that while the total energy released is roughly what we expected (about half escapes), the signal we would detect on Earth is different. The gravity of the black hole acts like a filter and a lens:

1. It changes the frequency (pitch) of the wave, making it lower.
2. It changes the shape of the wave, sometimes reflecting more energy back out than we thought.

So, if we want to find these "Scalar Tsunamis" from exploding stars, we need to tune our detectors to listen for lower-pitched, slightly different waves than we previously thought.

Technical Summary

Technical Summary: Scalar Tsunamis from Black Hole Formation

Problem Statement
Extreme astrophysical events, such as supernovae or neutron star mergers leading to black hole (BH) formation, offer a unique environment for searching for new physics, specifically feebly interacting, nearly massless scalar fields. Previous work [1] proposed that if a macroscopic object is surrounded by a large field configuration of such a scalar, the sudden disappearance of the matter source (as it collapses into a BH) releases the field, generating a propagating "scalar tsunami." However, the initial estimates in [1] relied on a simplified assumption: the evolution of the scalar field was modeled as a free field in flat space, neglecting the strong gravitational background of the newly formed black hole. This paper addresses the validity of that approximation by investigating how the Schwarzschild geometry affects the field's evolution, the total energy released, and the frequency spectrum of the signal arriving at large distances.

Methodology
The authors model a real, essentially massless scalar field $\phi$ minimally coupled to gravity. The mass is assumed to be sufficiently small ($m_\phi \lesssim 10^{-25}$ eV) to avoid signal dispersion over astrophysical distances. The study proceeds in two main stages:

1. Initial Conditions: The authors define several spherically symmetric initial field configurations representing the state of the scalar field just before the black hole forms. These include:

   • Static configurations: Yukawa-like profiles (generated by a coupling to matter density), homogeneous charged spheres, "compact" profiles (where non-linear interactions cluster the field), and spherical shells.
   • Dynamical configurations: Scenarios mimicking a collapsing source, including purely ingoing waves and a contracting homogeneous sphere with a finite velocity ($v \approx 0.4$), where the source vanishes once the Schwarzschild radius is reached.

2. Numerical Evolution: The field evolution is solved numerically in a static Schwarzschild black hole background. To handle the coordinate singularity at the horizon and the infinite domain, the authors employ tortoise coordinates ($x$), transforming the equation of motion into a 1D wave equation with an effective potential barrier $V_s(x)$:
   $$\ddot{\psi} = \frac{d^2\psi}{dx^2} - V_s(x)\psi$$
   where $\phi = \psi/r$. The potential $V_s(x)$ peaks near the photon sphere ($r \approx 1.33 r_H$) and decays exponentially near the horizon and as $1/x^3$ at infinity. The authors use the method of lines with second-order finite-difference operators to solve the partial differential equation.

3. Energy and Spectrum Analysis: The total energy released is calculated by integrating the energy density flux across a large sphere at late times. The authors compare the curved spacetime energy ($E_c$) against the flat space estimate ($E_f$) and calculate the fraction of energy escaping ($\epsilon$). The frequency spectrum is analyzed via the Power Spectral Density (PSD) of the time derivative of the field, $\dot{\phi}$, to determine the distribution of frequencies in the transient signal.

Key Contributions and Results

• Energy Release: The study confirms that the naive flat-space estimate—that approximately 50% of the field energy escapes while 50% is absorbed by the black hole—is generally a robust order-of-magnitude estimate even when general relativity is included.

  • For static Yukawa-like and homogeneous sphere configurations, the escape fraction $\epsilon \approx 0.52$.
  • For compact initial configurations, the escape fraction depends on the size of the source relative to the horizon. For large sources ($R \gg r_H$), the escape fraction approaches 1. This is attributed to the potential barrier reflecting low-frequency ingoing modes (which dominate large sources) back to infinity, rather than absorbing them.
  • Even for "purely ingoing" initial conditions, a small but non-zero fraction of energy escapes due to scattering off the potential barrier.

• Spectral Modifications: While the total energy is similar to flat-space estimates, the frequency spectrum receives significant corrections due to the black hole background:

  • Redshift: The spectrum is generally redshifted compared to flat space, with peak frequencies occurring at lower values. This is particularly pronounced for sources with radii comparable to the Schwarzschild radius.
  • Quasi-Normal Modes (QNMs): The signal exhibits peaks corresponding to the black hole's quasi-normal modes. For the scalar field, the lowest s-wave mode is identified at a dimensionless frequency $\hat{f} \approx 0.035$.
  • Suppression of High Frequencies: High-frequency ingoing modes are largely transmitted through the potential barrier into the black hole, leading to a suppression of high-frequency power in the escaping signal compared to flat space. Conversely, low-frequency modes are reflected, enhancing their relative power.
  • Late-Time Tails: The numerical results validate "Price's law," showing that the field decays asymptotically as $t^{-3}$ for Yukawa-like profiles and $t^{-4}$ for compact profiles, consistent with backscattering off the long-range curvature potential.

• Dynamical Effects: When modeling a contracting source (simulating a supernova collapse), the authors find that the resulting field evolution and energy release are qualitatively similar to the static case. The primary difference is a small dip in the field distribution caused by the source dragging the field during contraction, but the total escape fraction remains $\epsilon \approx 0.5$.

Significance and Claims
The paper claims that while the inclusion of general relativity does not drastically alter the total energy budget of the scalar tsunami (supporting the utility of simple flat-space estimates for order-of-magnitude energy constraints), it introduces notable and crucial corrections to the frequency spectrum.

The authors emphasize that experimental sensitivity depends heavily on the frequency content of the transient signal. The gravitational redshift and the filtering effects of the black hole potential (reflection of low frequencies, absorption of high frequencies) shift the signal away from the naive flat-space predictions. Consequently, searches for such signals in terrestrial experiments must account for these spectral distortions, particularly for sources where the initial field configuration is compact or where the source size is comparable to the resulting black hole's horizon. The work provides a more faithful numerical framework for predicting the observable signatures of scalar fields released during black hole formation.