Gauge dependence of scalar-induced gravitational waves from isocurvature perturbations: Analytical results

This paper analytically demonstrates that scalar-induced gravitational waves from isocurvature perturbations exhibit gauge-dependent divergences in most gauges during radiation domination, but by projecting out spurious contributions to isolate luminal modes, the authors derive a finite, gauge-independent late-time energy spectrum.

The Gist

Imagine the early universe as a giant, bubbling pot of soup. In this soup, there are two main ingredients: radiation (like hot photons) and matter (like cold dark matter). Usually, these ingredients mix perfectly, creating a smooth, uniform flavor. This is called an adiabatic state.

But sometimes, the mixing isn't perfect. You might have a spoonful of extra salt in one corner and a spoonful of extra pepper in another, even though the total amount of soup remains the same. This uneven distribution of ingredients is called an isocurvature perturbation.

According to Einstein's theory of gravity, when these ingredients slosh around, they don't just change the flavor; they create ripples in the fabric of space and time itself. These ripples are Gravitational Waves.

This paper is a detective story about a specific problem: How do we measure these ripples without getting tricked by the way we look at them?

The Problem: The "Camera Angle" Issue

Imagine you are filming a dancer.

• If you film them from the front, they look like they are moving forward.
• If you film them from the side, they look like they are moving sideways.
• If you zoom in and out, their speed looks different.

The dancer is the same, but the viewpoint (or "gauge" in physics) changes how the motion appears.

In the universe, scientists use different mathematical "camera angles" (called gauges) to describe these gravitational waves. The problem is that for a long time, when scientists looked at ripples caused by isocurvature (the uneven mixing), the results looked completely different depending on which camera angle they used.

• In some angles, the ripples seemed to grow infinitely large over time (like a sound getting louder and louder until it breaks the speakers).
• In other angles, the ripples settled down and behaved normally.

This was confusing! If the universe is real, the gravitational waves shouldn't depend on which math formula we choose to describe them. It's like saying the dancer is actually growing taller just because you changed the camera lens.

The Investigation: Testing Nine Different Lenses

The authors of this paper decided to test nine different camera angles (gauges) to see what was happening. They treated the universe during its "Radiation Domination" era (when the universe was very hot and filled with light).

They found a clear split in the results:

1. The "Broken Camera" Angles: In five of the angles (like the Uniform Density or Synchronous gauges), the math predicted that the gravitational waves would grow wildly out of control as time went on. It looked like the waves were getting stronger and stronger, eventually becoming infinite.

   • Analogy: It's like looking at a calm lake through a funhouse mirror that makes the waves look like they are turning into a tsunami. The lake is calm, but the mirror is lying.

2. The "Clear Lens" Angles: In the other four angles (like the Longitudinal or Newtonian gauges), the waves behaved normally. They settled down and acted like standard radiation, just fading away gently.

   • Analogy: This is looking at the lake through a clear, flat window. You see the gentle ripples exactly as they are.

The Solution: Filtering Out the "Static"

The authors realized that the "wild growth" in the broken angles wasn't real. It was mathematical noise—artifacts created by the way the coordinates were set up.

Think of it like listening to a radio station. Sometimes you hear the music (the real gravitational wave), but sometimes you also hear static and buzzing (the gauge artifacts). In the "broken" camera angles, the static was so loud it drowned out the music and made it sound like the volume was turning up to infinity.

The Fix:
The authors developed a "filter" (called a kernel projection). This filter isolates only the parts of the signal that are actually traveling waves—the parts that move at the speed of light, oscillating back and forth like a true wave (sin and cos functions).

When they applied this filter to all nine camera angles:

• The "wild growth" disappeared.
• The "static" vanished.
• All nine angles suddenly agreed! They all showed the exact same, calm, finite gravitational wave signal.

The Big Takeaway

1. The Waves are Real: Scalar-induced gravitational waves from isocurvature perturbations do exist and have a specific, predictable shape.
2. The Math Can Lie: If you use the wrong mathematical "camera angle," you might think the waves are exploding or growing infinitely. This is an illusion caused by the math, not the physics.
3. The Truth is Universal: Once you filter out the mathematical illusions and look only at the waves that actually travel through space, the result is the same no matter how you look at it.

Why Does This Matter?

Future telescopes (like LISA or pulsar timing arrays) are going to try to detect these faint ripples from the early universe. If scientists use the wrong "camera angle" or don't filter out the static, they might predict a signal that is way too strong or completely wrong.

This paper gives them the instruction manual on how to clean up the data. It ensures that when we finally "hear" the echo of the Big Bang, we are listening to the real music of the universe, not the static of our own math.

Technical Summary

1. Problem Statement

Scalar-induced gravitational waves (SIGWs), also known as secondary gravitational waves, are generated by the nonlinear coupling of first-order scalar perturbations to second-order tensor modes. While the gauge invariance of linear tensor perturbations is well-established, second-order SIGWs are not inherently gauge-invariant because the nonlinear structure of Einstein's equations couples different perturbation modes.

Previous literature has extensively analyzed the gauge dependence of SIGWs sourced by adiabatic (curvature) perturbations. However, the gauge dependence of SIGWs sourced by primordial isocurvature perturbations (variations in relative number densities of cosmic species) had not been systematically studied. Isocurvature modes evolve differently than adiabatic modes, potentially leading to distinct and more severe gauge artifacts. The authors aim to:

1. Analyze the gauge dependence of isocurvature-induced SIGWs across nine different gauges during the radiation-dominated (RD) era.
2. Identify the origin of divergences in certain gauges.
3. Establish a method to extract the physical, gauge-independent energy spectrum.

2. Methodology

The authors employ a comprehensive analytical framework to study SIGWs in the radiation-dominated era.

• Gauge Selection: The study covers nine distinct gauges, categorized by their slicing (time) and threading (space) conditions:
  • Divergent Gauges: Uniform-Density (UD), Total-Matter (TM), Uniform-Curvature (UC), Comoving-Orthogonal (CO), and Transverse-Traceless (TT).
  • Convergent Gauges: Longitudinal (Long.), Uniform-Expansion (UE), Newtonian-motion (Nm), and N-body (Nb).
• Formalism:
  • They utilize the perturbed Friedmann-Lemaître-Robertson-Walker (FLRW) metric with scalar ($\phi, \psi, B, E$) and tensor ($h_{ij}^{TT}$) components.
  • The source term for the tensor modes is derived from the quadratic combination of scalar perturbations.
  • The authors introduce dimensionless variables $x = k\eta$ (conformal time) and a transformation to the $(d, s)$ domain (where $d = |u-v|/\sqrt{3}$ and $s = (u+v)/\sqrt{3}$) to simplify the convolution integrals for the energy density spectrum.
• Analytical Integration:
  • Instead of relying solely on numerical integration, the authors derive exact analytical expressions for the source functions and the kernel functions ($I_c$ and $I_s$) for all nine gauges.
  • They utilize the xPand Mathematica package to handle the complex tensor algebra and gauge transformations.
  • They apply gauge transformation rules to relate results in arbitrary gauges to the longitudinal gauge, specifically isolating the gauge transformation term $\chi_{ij}^{TT}$.

3. Key Contributions

• First Comprehensive Analytic Study: This is the first work to provide a unified analytical treatment of isocurvature-induced SIGWs across nine different gauges, revealing a stark dichotomy in their late-time behavior.
• Identification of Divergence Mechanisms: The paper demonstrates that isocurvature perturbations induce stronger gauge divergences than adiabatic perturbations. For instance, gauges like TT and CO, which are often well-behaved for adiabatic sources, exhibit severe polynomial growth for isocurvature sources.
• Radiative Projection Technique: The authors propose and implement a "kernel projection" method. By isolating only the luminal (freely propagating) components of the tensor solution (terms oscillating as $\sin(k\eta)$ and $\cos(k\eta)$) and discarding non-radiative, gauge-dependent terms, they recover a unique physical spectrum.

4. Key Results

A. Dichotomy of Gauge Behaviors

The authors found a clear split in how the energy density spectrum $\Omega_{GW}$ evolves with conformal time $\eta$ (or $x=k\eta$) in the late-time limit ($x \gg 1$):

1. Divergent Gauges (Polynomial Growth):
   In these gauges, the energy density grows as a power of time, $\Omega_{GW} \propto \eta^n$, indicating unphysical behavior before projection:

   • Uniform-Density (UD): $\Omega_{GW} \propto \eta^2$
   • Total-Matter (TM) & Uniform-Curvature (UC): $\Omega_{GW} \propto \eta^4$
   • Comoving-Orthogonal (CO): $\Omega_{GW} \propto \eta^6$
   • Transverse-Traceless (TT): $\Omega_{GW} \propto \eta^8$ (The most severe divergence).

2. Convergent Gauges (Radiation-like Behavior):
   In these gauges, the late-time spectrum converges to a constant, behaving as standard radiation ($\Omega_{GW} \propto \text{const.}$):

   • Longitudinal, Uniform-Expansion (UE), Newtonian-motion (Nm), and N-body (Nb).

B. Subhorizon Divergence

For subhorizon modes ($k\eta \gg 1$), the divergence in the "bad" gauges becomes increasingly severe, confirming that SIGWs are gauge-dependent observables in this regime if the full tensor solution is considered.

C. Resolution via Radiative Projection

The authors demonstrate that the apparent gauge dependence arises from non-radiative tensor modes (spurious gauge artifacts) that contaminate the metric in specific slicings.

• By applying a projection that retains only the oscillatory $\sin(k\eta)$ and $\cos(k\eta)$ terms (the free gravitational wave components), the kernel decays as $(k\eta)^{-1}$.
• Result: After this projection, all nine gauges yield the exact same finite, gauge-independent energy spectrum, which matches the longitudinal gauge benchmark.

D. Spectral Features

For a Dirac-delta isocurvature peak at $k_p$, the physical spectrum exhibits:

• An infrared rise scaling as $k^2 \ln^2 k$ for $k \ll k_p$.
• A resonant peak at $k \approx 2c_s k_p$ (where $c_s = 1/\sqrt{3}$ is the sound speed).
• A sharp UV cutoff at $k = 2k_p$ due to momentum conservation.

5. Significance

• Theoretical Consistency: The paper resolves the long-standing issue of gauge ambiguity in second-order cosmological perturbation theory for isocurvature sources. It confirms that while intermediate calculations are gauge-dependent, the physical observable (the energy density of free gravitational waves) is gauge-invariant.
• Methodological Framework: The "radiative projection" method provides a robust prescription for extracting physical SIGW signals from any gauge, ensuring that theoretical predictions are not contaminated by coordinate artifacts.
• Observational Implications: By establishing a consistent, gauge-independent spectrum, the work provides reliable theoretical templates for future gravitational wave experiments (e.g., LISA, TianQin, PTA) searching for signals from primordial isocurvature perturbations or Primordial Black Hole (PBH) formation scenarios.
• Insight into Perturbation Theory: The severe divergences found in gauges like TT and CO for isocurvature modes suggest that linear perturbation theory may break down earlier in these specific slicings compared to adiabatic cases, highlighting the need for careful gauge selection in numerical simulations and analytical studies.

In conclusion, the paper establishes that only the luminal, freely propagating modes represent physical SIGWs, and by isolating these modes, one obtains a unique, gauge-independent prediction for the gravitational wave background sourced by isocurvature perturbations.