Scalar-induced gravitational waves including isocurvature perturbations with lattice simulations

This paper introduces a lattice simulation framework to compute scalar-induced gravitational waves from isocurvature and mixed perturbations, demonstrating strong agreement with semi-analytical predictions and revealing how spectral features in early matter-dominated eras depend on microphysical properties like primordial black hole characteristics.

The Gist

Imagine the universe as a giant, invisible ocean. For a long time, scientists have been studying the "waves" on the surface of this ocean—specifically, ripples in space-time called Gravitational Waves.

Most of what we know about these waves comes from studying the "calm" parts of the ocean, where everything flows smoothly together. This is called adiabatic behavior (think of a crowd of people walking in perfect unison).

However, this new paper by Xiang-Xi Zeng suggests we've been ignoring the "chaos" in the ocean. There are times when different parts of the universe move at different speeds or have different densities, creating a kind of internal friction. This is called isocurvature behavior (think of a crowd where some people are running while others are standing still, creating a jostle).

Here is a simple breakdown of what this paper does, using some everyday analogies:

1. The Problem: We Only Watched the Calm Waves

For years, scientists have tried to predict how the early universe created gravitational waves. They mostly looked at the "perfect unison" scenario. But the early universe was messy. It had different types of matter (like radiation and dark matter) interacting in complex ways. The old math (semi-analytical formulas) works well for simple cases but gets messy and inaccurate when things get complicated or when different types of matter are mixed together.

2. The Solution: A Digital "Sandbox" (Lattice Simulations)

Instead of trying to solve a giant, impossible math equation, the author built a digital sandbox.

• The Analogy: Imagine trying to predict how a drop of ink spreads in a glass of water. You could try to write a complex formula for every molecule, or you could fill a computer with a grid of pixels and watch the ink move pixel by pixel.
• The Paper's Method: The author used a "lattice simulation." They created a virtual box representing the early universe, filled it with different types of "fluids" (matter and radiation), and let the computer simulate how they interacted over time. This allowed them to see the gravitational waves that naturally emerged from the chaos, including the messy "isocurvature" parts that old math missed.

3. What They Found: The "Echoes" of the Early Universe

The simulation revealed some fascinating things:

• The "Double-Track" Effect: When you mix the "calm" (adiabatic) and "chaotic" (isocurvature) movements, the resulting gravitational waves don't just look like a simple sum of the two. They create a unique multi-peak structure.

  • Analogy: If you hit a drum (calm) and shake a bucket of marbles (chaos) at the same time, the sound isn't just a drumbeat plus a rattle. It creates a complex chord with specific high and low notes. The paper shows that these "notes" (peaks in the wave spectrum) appear in very specific places, acting like a fingerprint for the early universe.

• The "Black Hole Battery": The paper also looked at a scenario where the early universe was dominated by tiny Primordial Black Holes (PBHs). Imagine these black holes as tiny batteries that slowly leak energy (via Hawking radiation) and eventually disappear.

  • The Discovery: How fast these "batteries" leak energy changes the sound of the gravitational waves.
  • Analogy: If you have a leaky bucket, the sound of the water dripping changes depending on how big the hole is and how full the bucket is. The paper found that the size of the black holes and how fast they decay dictate the "pitch" and "volume" of the gravitational waves. A faster decay creates a louder, sharper signal.

4. Why This Matters: Listening to the Universe's History

Why do we care about these digital simulations?

• The New Microphone: We are entering a new era of "Gravitational Wave Astronomy." Detectors like LISA (a space-based telescope) and Pulsar Timing Arrays are about to start "listening" to the universe's background noise.
• Decoding the Message: If we hear a specific pattern of waves, we need to know what caused it. Was it a smooth inflation? Or was it a chaotic mix of black holes and radiation?
• The Takeaway: This paper provides the "decoder ring." It tells us that if we see these specific "multi-peak" patterns or specific "decay slopes" in the gravitational wave data, it might be the signature of isocurvature perturbations or decaying primordial black holes.

Summary

In short, this paper builds a super-accurate computer model to simulate the messy, chaotic early universe. It shows us that the "noise" left over from that chaos (gravitational waves) carries a unique signature. By understanding this signature, future telescopes might be able to tell us exactly what the universe was made of billions of years ago, potentially revealing the existence of tiny black holes or exotic particles that we can't see any other way.

It's like realizing that the static on an old radio isn't just noise; it's a recorded message from the Big Bang, and this paper gives us the manual to translate it.

Technical Summary

1. Problem Statement

Scalar-induced gravitational waves (SIGWs) are a promising probe of early-universe physics, generated when primordial scalar perturbations re-enter the horizon and source tensor fluctuations via second-order gravitational interactions. While SIGWs from adiabatic perturbations (consistent with single-field inflation) have been extensively studied, the contribution from isocurvature perturbations remains largely unexplored.

Key challenges addressed in this work include:

• Complexity of Mixed Conditions: Semi-analytical formalisms for isocurvature-induced SIGWs are strictly valid only when modes enter the horizon well before matter-radiation equality. They become inaccurate or intractable for mixed initial conditions (adiabatic + isocurvature) or when modes enter near equality.
• Non-Gaussianity: Isocurvature perturbations, particularly those of Poissonian origin (e.g., from Primordial Black Holes, PBHs), are inherently large ($O(1)$) and non-Gaussian, requiring non-perturbative methods for accurate modeling.
• Early Matter-Dominated (eMD) Eras: Scenarios involving PBHs or non-topological solitons (like Q-balls) create eMD eras. The transition from eMD to radiation domination significantly affects SIGW spectra, but numerical studies for isocurvature modes in these eras are scarce.

2. Methodology

The authors developed a lattice simulation framework to compute the stochastic gravitational wave background (SGWB) from both pure isocurvature and mixed initial conditions.

• Theoretical Formalism:
  • The universe is modeled as a two-fluid system (radiation and non-relativistic matter).
  • The evolution is governed by perturbed Einstein equations and energy-momentum conservation equations in the Newtonian gauge.
  • Crucially, the framework includes energy transfer ($Q$) between components, essential for modeling PBH evaporation or soliton decay.
  • The source term for tensor modes ($h_{ij}$) includes contributions from scalar potentials ($\Phi$), fluid velocities ($v$), and a relative velocity term ($V_{rel} = v_m - v_r$), which is often neglected in literature but explicitly tested here.
• Numerical Implementation:
  • Simulations are performed in a comoving cubic box ($N=256$) with periodic boundary conditions.
  • Spatial derivatives are computed using a Fourier pseudospectral method.
  • Time integration uses a fourth-order Runge-Kutta scheme.
  • Initial conditions are set using Gaussian random fields for curvature ($\zeta$) and isocurvature ($S$) perturbations, with specific power spectra (Gaussian bumps or Dirac delta functions).
• Scenarios Investigated:
  1. No Energy Transfer: Pure isocurvature and mixed initial conditions in standard radiation/matter eras.
  2. With Energy Transfer: PBH-dominated eras where PBHs decay via Hawking radiation, and Q-ball scenarios with varying decay rates.

3. Key Contributions

• Extension to Isocurvature Modes: This is the first systematic application of lattice simulations to SIGWs generated specifically by isocurvature perturbations, extending previous work limited to adiabatic fluctuations.
• Validation of Semi-Analytical Methods: The authors benchmarked their lattice results against existing semi-analytical formulas for pure isocurvature cases, confirming excellent agreement when modes enter the horizon early. They identified where semi-analytical methods fail (near matter-radiation equality).
• Analysis of Multi-Peak Structures: The study extends the understanding of multi-peak GW spectra to mixed initial conditions, showing that while peak locations follow momentum conservation rules similar to adiabatic cases, isocurvature contributions introduce scale-dependent suppression factors.
• Microphysical Sensitivity in eMD: The work establishes a direct link between the microphysical properties of the dominant field (PBH mass, abundance, soliton decay rate) and the resulting SIGW spectral shape (peak amplitude and slope).

4. Key Results

A. Pure Isocurvature and Mixed Conditions

• Agreement with Theory: In the absence of energy transfer, lattice simulations match semi-analytical predictions perfectly for pure isocurvature cases.
• Peak Structures: For mixed initial conditions (adiabatic + isocurvature), the GW spectrum exhibits multi-peak structures. The number and location of peaks are governed by momentum conservation ($|k_i - k_j| \le k \le k_i + k_j$), similar to adiabatic cases.
• Distinction: While the peak locations are similar, the amplitudes differ. Isocurvature peaks are suppressed by factors involving the equality scale ($k_{eq}$), allowing potential discrimination between initial conditions via precise spectral shape analysis.

B. Early Matter-Dominated (eMD) Scenarios

• PBH-Induced SIGWs:
  • The peak amplitude of the GW spectrum is dictated solely by the PBH mass.
  • The location of the spectral break (transition in power-law scaling) is governed solely by the initial PBH abundance ($\beta$).
  • A larger PBH mass or higher abundance leads to a steeper power-law index near the peak.
• Q-Ball Decay:
  • The decay rate of non-topological solitons (parameterized by $n$ in $\Gamma \propto n/(t_{eva}-t)$) significantly impacts the signal.
  • Slower decay rates (smaller $n$) result in shallower power-law indices near the peak and reduced GW amplitudes.
• Relative Velocity: The inclusion of the relative velocity term ($V_{rel}$) in the source term was found to be negligible in the PBH-dominated scenario, consistent with previous analytical estimates, though it is prominent in early stages.

C. General Initial Conditions in eMD

• When simulating mixed conditions where PBHs form from a peak-like adiabatic perturbation, the GW signal evolves dynamically: initially sourced by adiabatic perturbations, then dominated by isocurvature modes as PBHs form and decay.
• The rapid transition from eMD to radiation domination amplifies the adiabatic part of the signal, indicating that simple linear superposition of adiabatic and isocurvature components is insufficient; full dynamical simulation is required.

5. Significance and Future Outlook

• Robust Tool: The study validates lattice simulations as a robust tool for predicting SIGW spectra from complex, non-Gaussian primordial perturbations, filling a gap left by semi-analytical methods.
• Observational Implications: The findings provide specific spectral signatures (peak heights, slopes, and break locations) that future gravitational wave detectors (LISA, Taiji, Einstein Telescope, DECIGO) can use to distinguish between different early-universe scenarios (e.g., PBH formation vs. soliton decay, adiabatic vs. isocurvature dominance).
• Future Work: The authors note that while perturbative equations are solved, future work must address non-linear cutoff scales, PBH accretion, clustering effects, and potentially full numerical relativity simulations to capture non-perturbative regimes.

In summary, this paper provides a critical numerical framework for understanding how isocurvature perturbations and early matter-dominated eras shape the stochastic gravitational wave background, offering new avenues for interpreting upcoming GW data.