Emergent Turbulence in Nonlinear Gravity

Gravity in nonlinear and dynamical regimes underpins spectacular astrophysical phenomena and observable consequences, from the early universe to black hole collisions. In these extreme environments, inverse energy cascades - mediated by nonlinear interactions - may help explain the near scale-invariance of cosmic structure and the simplicity of gravitational waves from binary black hole mergers. Yet the presence, characteristics, and generality of such interactions in full General Relativity remain largely unexplored. Here we show that two types of nonlinear interactions - a four-mode and a three-mode interaction - emerge in the fully nonlinear regime, and can indeed channel inverse energy cascades by inducing resonant and anti-damping instabilities. This establishes what was previously only hinted at in highly specialized perturbative contexts. We further demonstrate a laminar'' to turbulent’’ transition for the largest-possible angular structure in General Relativity, whereas finer structures remain persistently turbulent. Our results reveal the impact and generality of these nonlinear interactions (instabilities), which can be key to understanding observations ranging from cosmological to kilometer scales. We anticipate that our work will shed new light on nonlinear gravitational phenomena and their consequences, such as constructing gravitational wave templates and testing General Relativity in the most extreme regime. Moreover, our work is a starting point for addressing nonlinear gravitational interactions using ideas and methods inspired by fluid dynamics.

💡 Research Summary

The paper investigates fully nonlinear dynamics in General Relativity (GR) by continuously “stirring” asymptotically flat spacetimes with quasi‑steady gravitational waves (GWs) injected at the outer boundary. Using the generalized harmonic formulation of Einstein’s equations, discretized with pseudo‑spectral methods and constraint‑preserving boundary conditions, the authors simulate two setups: a non‑spinning Schwarzschild black hole (BH) and pure flat spacetime. The injected waves are spin‑weighted spherical harmonic modes (ℓ,m) with ℓ = 2 (the largest possible angular structure in GR) or ℓ = 6 (allowing an inverse cascade to lower ℓ). The driver amplitude Aℓ,m(t) is ramped up linearly to a maximum value Ai (up to 600 Mi) and then held constant; the frequency ω is fixed and real.

The response of the spacetime is monitored via the tendicity scalar E, constructed from the Weyl tensor and orthonormal observers. Harmonic decomposition of E reveals how energy is transferred among modes and frequencies. Two distinct nonlinear interaction mechanisms emerge:

1. Four‑mode (quartic) coupling – a trilinear source term S^(3) involving three parent modes (e.g., ±ω from the ℓ = 2 driver) generates a daughter mode at 2 ω in the ℓ = 6 sector. This interaction acts as an anti‑damping term: the daughter amplitude grows exponentially with a rate proportional to Ai², saturating at a level ∝ Ai⁴. The process is visualized with a Feynman‑type diagram and explains the emergence of a 2 ω component after the initial laminar excitation of 3 ω and ω modes.

2. Three‑mode (cubic) coupling – a bilinear source term S^(2) couples two parent modes (e.g., 2 ω and – ω from an ℓ = 6 driver) to resonantly excite a daughter mode at ω in lower‑ℓ sectors (ℓ = 2, 4). The daughter amplitude grows linearly in time with a coefficient ∝ Ai³, reflecting a resonant “harmonic‑oscillator‑like” instability. This mechanism drives a spatial inverse cascade: energy moves from higher angular indices (ℓ = 6) toward the largest‑scale mode (ℓ = 2).

When the driver is ℓ = 2, early‑time laminar responses (3 ω and ω) appear, phase‑locked to the driver and scaling as Ai³ and Ai¹ respectively. For sufficiently large Ai (≈2 × 10⁻⁴), the four‑mode coupling triggers exponential growth of the 2 ω mode, leading to a clear laminar‑to‑turbulent transition. Time‑frequency plots display a cascade 3 ω → 2 ω → ω, i.e., an inverse cascade toward lower frequencies.

With an ℓ = 6 driver, each angular harmonic contains both 2 ω (quadratic coupling, amplitude ∝ Ai²) and ω components. The ω component in ℓ = 2 and ℓ = 4 grows continuously, redistributing its angular spectrum toward lower ℓ, evidencing a spatial inverse cascade independent of the presence of a BH. In flat spacetime the same pattern appears but with growth rates three to four orders of magnitude smaller, indicating that the phenomenon is generic to GR nonlinearity but amplified by BH dynamics.

The BH acts as a nonlinear amplifier: low‑ω waves are reflected, high‑ω waves are absorbed, and near the quasinormal‑mode frequency the absorption rate rises sharply, causing the BH mass to increase by up to a factor of six. This mass growth further enhances the nonlinear couplings.

Quantitative scaling laws are established: the saturation amplitude of the 2 ω mode (four‑mode channel) scales as Ai⁴, its exponential growth rate as Ai²; the ω mode (three‑mode channel) grows with a rate ∝ Ai³. These scalings match the perturbative expectations from the wave‑operator equations O(ψ)=S(ψ) with trilinear and bilinear source terms.

Overall, the study demonstrates that fully nonlinear GR exhibits both temporal and spatial inverse energy cascades analogous to fluid turbulence. The laminar‑to‑turbulent transition occurs only for the largest angular structure (ℓ = 2), while finer angular structures remain persistently turbulent. The findings have implications for cosmological structure formation, the simplicity of gravitational‑wave signals from binary black‑hole mergers, and the construction of more accurate GW templates. Moreover, the work opens a pathway for applying fluid‑dynamic concepts—such as cascades, instabilities, and steady‑state driving—to the extreme‑gravity regime, offering a new theoretical and computational framework for probing GR beyond the perturbative frontier.