Gravitational wave interactions with a viscous fluid: Core collapse supernova, binary neutron star merger, and accretion around a black hole merger

This paper extends the theory of gravitational wave interactions with viscous fluids to general static, spherically symmetric spacetimes, revealing that damping and heating effects in astrophysical scenarios like core-collapse supernovae and black hole mergers can be orders of magnitude stronger than in flat space, potentially leading to complete signal damping and gamma-ray bursts.

The Gist

The Big Idea: Ripples in a Thick Soup

Imagine you are standing by a calm pond. If you throw a stone in, you see ripples (waves) travel across the surface. If the water is thin and clear, those ripples travel far without losing much energy.

Now, imagine that same pond, but instead of water, it is filled with thick, sticky honey. If you throw a stone in, the ripples still happen, but the honey resists the motion. The ripples slow down, lose their energy, and eventually stop. Where did that energy go? It turned into heat. The honey gets warmer because the ripples were rubbing against it.

This paper is about what happens when "ripples" in the fabric of space (Gravitational Waves) travel through "thick honey" (viscous fluid) in the universe.

For a long time, scientists thought this "rubbing" effect was so tiny it didn't matter. They assumed space was mostly empty, like a vacuum, so waves could travel forever without getting tired. But this team of researchers asked: "What if the waves are traveling through a very dense, sticky environment right next to where they were born?"

They found that in these specific, crowded cosmic neighborhoods, the waves don't just pass through; they get stuck, damped (weakened), and heat up the surrounding matter significantly.

---

The Three Cosmic Scenarios

The researchers tested this idea in three extreme situations where space is crowded with matter:

1. The Core-Collapse Supernova (The Exploding Star)

• The Scene: A massive star runs out of fuel and collapses inward, then bounces back with a massive explosion.
• The Analogy: Imagine a giant, dense ball of dough (the star's core) being squeezed and then suddenly released. It vibrates violently, sending out shockwaves.
• The Discovery: The researchers found that the gravitational waves generated by this explosion have to fight their way through the star's own thick, hot outer layers (the "mantle").
• The Result: The waves lose almost all their energy trying to escape. Instead of traveling far into space to be detected by our telescopes, the energy is dumped into the star's outer layers, heating them up to incredible temperatures.
  • Simple takeaway: The star might be "screaming" in gravitational waves, but the sound is being swallowed by its own body, turning into heat instead of a signal we can hear.

2. The Binary Neutron Star Merger (The Cosmic Dance)

• The Scene: Two incredibly dense stars (neutron stars) spiral into each other and crash together.
• The Analogy: Think of two ice skaters spinning faster and faster until they collide. The crash creates a new, super-dense object surrounded by a swirling disk of debris.
• The Discovery: After the crash, the new object is surrounded by a hot, sticky disk of matter. The gravitational waves trying to escape this disk get heavily damped.
• The Result: The friction from the waves heating up the debris could be a major source of energy. It might explain why these events are so hot and energetic, potentially contributing to the formation of powerful jets of energy (like gamma-ray bursts).

3. The Black Hole Merger with an Accretion Disk (The Black Hole with a Ring)

• The Scene: Two black holes merge. Usually, black holes are in empty space, but sometimes they are surrounded by a ring of gas and dust (an accretion disk).
• The Analogy: Imagine two whirlpools merging in a river. Usually, the river is empty, but here, there is a thick ring of mud swirling around them.
• The Discovery: When the black holes merge, they send out a massive burst of gravitational waves. If there is a ring of gas nearby, the waves hit it.
• The Result: The paper found that the waves heat up the gas in the ring massively.
  • The "Wow" Factor: The heating was so intense (trillions of degrees) that it could explain a mysterious flash of light (a Gamma-Ray Burst) that was seen when the black holes merged in 2015 (GW150914). It suggests the black hole merger didn't just make a "thud" in space; it also "fried" the surrounding gas, creating a burst of light.

---

Why This Matters: The "Sticky Bead" Experiment

The paper references a famous thought experiment from 1957 by Richard Feynman called the "sticky bead."

• The Experiment: Imagine a wire with a bead on it. If you shake the wire, friction between the bead and the wire creates heat. Feynman said, "If gravitational waves can shake a bead and create heat, then gravitational waves carry energy."
• The Paper's Contribution: For decades, we knew the bead could get warm, but we thought it would take a billion years to get even a little warm. This paper says, "Wait a minute! If the bead is in a very thick, sticky fluid and the shaking happens very close by, the bead gets hot almost instantly."

The "Minkowski" vs. "Real World" Difference

The researchers compared two ways of doing the math:

1. The "Empty Room" Model (Minkowski): This assumes space is flat and empty, like a vacuum. In this model, the waves barely lose any energy.
2. The "Real World" Model (Curved Spacetime): This accounts for the fact that massive objects (like stars and black holes) bend space, making it "curved" and "crowded."

The Surprise: When they used the "Real World" model, the waves lost energy much faster (sometimes by millions of times more) than the "Empty Room" model predicted.

Summary in a Nutshell

• Old View: Gravitational waves travel through the universe like ghosts, passing through matter without touching it.
• New View: If gravitational waves are born inside or near thick, sticky matter (like a dying star or a merging black hole's debris), they interact strongly with it.
• The Consequence: The waves get "stuck" and die out quickly, but in doing so, they dump a huge amount of energy into the matter, heating it up to extreme temperatures.
• Why it's cool: This might explain why some cosmic explosions are so hot and why we sometimes see flashes of light (Gamma-Ray Bursts) right when we detect gravitational waves. It changes how we understand the "aftermath" of the universe's most violent events.

In short: Gravitational waves aren't just silent messengers; in the right conditions, they are cosmic heaters.

Technical Summary

1. Problem Statement

While it is well-established that gravitational waves (GWs) can transfer energy to matter (the "sticky bead" argument), standard astrophysical models typically treat this interaction as negligible. This assumption relies on the premise that the damping factor is extremely small in cosmological or astrophysical contexts because the shear viscosity ($\eta$) and propagation distance ($\Delta$) product rarely approaches the scale required to make the damping term significant ($G\eta\Delta/c^3 \ll 1$).

However, recent work suggested that when the distance between the GW source and the matter is comparable to or smaller than the GW wavelength ($r \lesssim \lambda$), the damping and consequent heating effects can become astrophysically significant. Previous studies were limited to Minkowski backgrounds (flat spacetime) or Schwarzschild backgrounds (vacuum).

The Core Problem: This paper addresses the gap in modeling GW-matter interactions within general, non-vacuum, static, spherically symmetric spacetimes. Realistic astrophysical scenarios (like supernovae or merger remnants) involve dense fluid shells surrounding a core, creating a background geometry that is neither flat nor a simple vacuum Schwarzschild metric. The authors aim to quantify how this realistic background geometry alters GW damping and fluid heating compared to previous simplified models.

2. Methodology

The authors extend the formalism of Bondi-Sachs metric perturbations to a general static, spherically symmetric background containing matter.

• Background Geometry:

  • The model consists of an inner core of mass $M_c$ and radius $r_c$, surrounded by a static fluid shell ($r_c < r < r_t$) with density $\rho^{[B]}$ and pressure $p^{[B]}$, and an exterior vacuum region ($r > r_t$) described by the Schwarzschild metric.
  • The background metric functions ($\beta^{[B]}, W^{[B]}$) are determined by solving a system of three first-order Ordinary Differential Equations (ODEs) derived from the Einstein equations and matter conservation.

• Perturbation Analysis:

  • Small perturbations are introduced to the metric and fluid velocity, oscillating with frequency $\nu$.
  • The perturbations are decomposed using spin-weighted spherical harmonics, focusing on the dominant quadrupolar ($\ell=2, m=2$) mode.
  • In the fluid region, the perturbed Einstein equations yield a system of five coupled first-order ODEs for the metric perturbation variables ($\beta^{[2,2]}, W^{[2,2]}, U^{[2,2]}, J^{[2,2]}$) and fluid velocity perturbations.
  • Unlike the Minkowski case (where solutions are analytic) or the pure Schwarzschild case (which reduces to a single master ODE), the general case requires full numerical integration.

• Numerical Implementation:

  • The authors developed computer code (using Maple and Matlab/Octave) to solve the background ODEs and the perturbation system.
  • Shear and Heating: The fluid shear tensor $\sigma_{ab}$ is calculated from the velocity perturbations. The rate of energy transfer (heating) is derived from the time-averaged shear scalar $\langle \sigma^2 \rangle$, proportional to the shear viscosity $\eta$.
  • Damping Factor: The GW damping is calculated by integrating the energy transfer rate over the fluid shell, yielding a damping factor $f_E(r)$ that modifies the wave amplitude.

3. Key Contributions

1. Generalization of Background Spacetime: The paper successfully extends GW-matter interaction theory from Minkowski and Schwarzschild backgrounds to a general static, spherically symmetric spacetime with matter. This allows for the modeling of fluid shells with arbitrary density profiles.
2. Numerical Framework: The authors provide a robust numerical framework to solve the coupled ODEs for metric and fluid perturbations in non-vacuum backgrounds, where analytic solutions are unavailable.
3. Quantification of Background Effects: The study demonstrates that the presence of a massive fluid shell (non-vacuum background) significantly enhances GW damping and heating effects compared to Minkowski background calculations, particularly when $r \lesssim \lambda$.

4. Key Results

The authors applied their code to three specific astrophysical scenarios:

A. Core-Collapse Supernovae (CCSNe)

• Parameters: Modeled a proto-neutron star (PNS) with a dense core and a hot mantle.
• Findings:
  • GW damping is nearly complete for typical viscosity ranges ($10^{23} - 10^{25}$ kg m$^{-1}$ s$^{-1}$) and frequencies (100–1000 Hz).
  • The Minkowski code underestimates damping and heating by several orders of magnitude compared to the matter code.
  • Heating: The predicted temperature rise ($\Delta T$) in the mantle can reach $10^{14} - 10^{20}$ K (depending on viscosity and energy release), far exceeding the ambient temperature ($\sim 10^{12}$ K).
  • Implication: GW signals from CCSNe might be completely damped before escaping the star, making them unobservable unless the frequency is very high ($>1$ kHz).

B. Binary Neutron Star (BNS) Merger Post-Merger

• Parameters: Modeled the hypermassive neutron star remnant and accretion disk.
• Findings:
  • Similar to CCSNe, higher viscosity and lower frequencies lead to stronger damping.
  • For high viscosity ($\eta = 10^{30}$), the damping factor at the outer boundary drops to $<10^{-3}$.
  • Heating: GW-induced viscous heating can raise temperatures to $10^{11} - 10^{14}$ K, comparable to or exceeding the thermal state expected from shock heating alone. This suggests GW dissipation is a non-negligible contributor to the remnant's thermal budget.

C. Accretion around a Binary Black Hole (BBH) Merger

• Parameters: Modeled an accretion disk around the remnant of a BBH merger (e.g., GW150914) using a Schwarzschild background.
• Findings:
  • Damping: GW damping is negligible ($<10^{-14}$) because the disk mass is insignificant compared to the black hole mass.
  • Heating: Despite negligible damping, the heating effect is massive.
  • At the Innermost Stable Circular Orbit (ISCO), the Schwarzschild background code predicts a temperature increase of $\sim 10^{12}$ K at the equator, compared to $\sim 10^8$ K in the Minkowski approximation.
  • Significance: This temperature spike is sufficient to potentially generate a Gamma-Ray Burst (GRB), offering a theoretical explanation for the GRB associated with GW150914.

5. Significance and Conclusion

• Astrophysical Relevance: The study proves that ignoring the background curvature and matter distribution (using Minkowski approximations) leads to severe underestimations of GW damping and heating in dense environments.
• Observational Impact:
  • CCSNe: GW signals may be undetectable from Earth if the source is embedded in a dense, viscous fluid shell, altering multi-messenger search strategies.
  • BNS/BBH: GW heating could be a primary driver for high-energy electromagnetic counterparts (like GRBs) in merger events.
• Limitations: The authors note that the models assume static backgrounds and small perturbations, whereas real events are highly dynamic. Future work should utilize full Numerical Relativity to capture these dynamics.
• Conclusion: The interaction between GWs and viscous fluids in realistic astrophysical backgrounds is not just a theoretical curiosity but a potentially dominant physical process that can suppress GW signals and drive extreme heating events.