Probing Scalar-Tensor-Induced Gravitational Waves in the nHz Band: $\texttt{NANOGrav}$ and SKA

This paper computes the energy density of scalar-tensor-induced gravitational waves (STGWs) generated during early matter-dominated and radiation-dominated eras, demonstrating that while these signals dilute in a purely matter-dominated universe, they remain viable and potentially dominant contributors to the nanohertz stochastic gravitational wave background observed by NANOGrav and detectable by the future Square Kilometre Array.

The Gist

The Big Picture: Listening to the Universe's Hum

Imagine the universe isn't just a silent void, but a giant, vibrating drum. When massive objects crash into each other (like black holes), they create ripples in space-time called Gravitational Waves (GWs).

Recently, a group of astronomers called NANOGrav (using a network of pulsars, which are like cosmic lighthouses) detected a low-frequency "hum" in the universe. They call this the Stochastic Gravitational Wave Background (SGWB). It's like hearing the static on a radio, but this static is made of gravitational waves.

The big question is: What is making this noise?

• The Obvious Answer: Supermassive black holes orbiting each other (like two dancers spinning slowly).
• The New Idea: This paper explores a more exotic possibility: Scalar-Tensor-Induced Gravitational Waves (STGWs).

The Core Concept: The "Cosmic Mixer"

To understand STGWs, imagine the early universe was a kitchen with two main ingredients:

1. Scalar Perturbations: Think of these as density bumps. Some spots in the universe were slightly denser than others (like lumps in dough).
2. Tensor Perturbations: Think of these as ripples in the fabric itself (the gravitational waves).

Usually, scientists thought these two ingredients stayed separate. But this paper investigates what happens if you mix them together.

• The Analogy: Imagine you have a trampoline (space-time). If you drop a heavy bowling ball (a density bump) on it, it creates a dip. If you then shake the trampoline (a gravitational wave), the bowling ball doesn't just sit there; it interacts with the shaking.
• The Result: This interaction creates new ripples. The paper calculates exactly how loud these new ripples would be.

The Plot Twist: The "Early Matter" Era

The paper looks at two different eras in the universe's history:

1. The Standard Era (Radiation Dominated): The universe was hot and filled with light particles.
2. The "Early Matter" Era (eMD): A hypothetical time before the standard era where the universe was filled with slow-moving matter (like dust) instead of light.

The Discovery:

• In the Standard Era: Mixing the ingredients creates a signal, but it's relatively quiet.
• In the Early Matter Era: The authors found a "secret sauce." If the universe went through a short phase of being "matter-dominated" and then suddenly switched to being "radiation-dominated" (like a sudden reheating), the mixing creates a massive, loud signal.

They call this the "Poltergeist Mechanism."

• Why? Imagine a ghost (the gravitational wave) that seems to disappear when the house is quiet (the matter era), but then suddenly screams when the lights turn on (the radiation era). The signal is generated after the transition, even though the "ghost" was created earlier.

The Investigation: NANOGrav vs. The Future (SKA)

The authors used this theory to check two things:

1. Can it explain what NANOGrav sees today?
They ran the numbers using the current data from NANOGrav.

• The Problem: To make the signal loud enough to match NANOGrav's "hum," you need a lot of density bumps. But if you have that many bumps, they would collapse into too many Primordial Black Holes (PBHs).
• The Constraint: If you create too many black holes, they would eat up all the dark matter in the universe, which we know isn't happening. So, while the theory could fit the data, it's on a very tight leash. It's like trying to fill a cup with water without spilling over the edge; it's possible, but very difficult.

2. What will the Square Kilometre Array (SKA) see?
The SKA is a future, super-powerful radio telescope that will listen to the universe with much better ears.

• The Forecast: The authors simulated what SKA would see. They found that if this "Scalar-Tensor" signal is real, SKA will be able to pinpoint exactly what it is.
• The Result: Unlike current data, which is blurry and ambiguous, SKA will be able to tell the difference between "black holes dancing" and "this exotic cosmic mixing." It will measure the "shape" of the sound so precisely that we can confirm if the universe went through that "Early Matter" phase.

The "Ghost" in the Machine (Technical Nuance)

The paper also points out a tricky detail. In the "Early Matter" era, the math gets messy.

• Scalar waves (the density bumps) behave one way.
• Tensor waves (the ripples) behave differently.
• When the universe switches from "Matter" to "Radiation," the tensor waves don't match up perfectly with the standard math used for scalar waves. The authors had to write new equations to handle this "mismatch." They found that while the "Poltergeist" effect is real, it's not exactly the same as the scalar version; it requires a bit more fine-tuning.

Summary: Why Does This Matter?

1. New Physics: If we detect this specific signal, it proves the universe had a weird, "matter-dominated" childhood that we didn't know about.
2. Black Hole Limits: It helps us understand how many primordial black holes could exist without breaking the universe.
3. Future Proof: It tells the SKA team exactly what to look for. If the SKA hears this specific "hum," we will know that the early universe was a chaotic place where density bumps and space-time ripples were dancing together.

In a nutshell: This paper is a recipe book for a new type of cosmic noise. It says, "If the universe had a specific childhood (Early Matter Era), it would make a specific sound. NANOGrav might be hearing a faint echo of it, but the SKA will be able to hear the whole song clearly."

Technical Summary

1. Problem Statement and Motivation

The paper addresses the origin of the stochastic gravitational wave background (SGWB) signal recently reported by Pulsar Timing Arrays (PTAs), specifically the NANOGrav 15-year dataset, in the nanohertz (nHz) frequency band.

• Context: While supermassive black hole binaries (SMBHBs) are the standard astrophysical explanation, cosmological sources are strong alternatives. Scalar-induced gravitational waves (SIGWs), generated by large primordial scalar perturbations, are a leading candidate.
• Gap in Knowledge: Most studies focus on purely scalar-induced GWs. However, Scalar-Tensor-Induced Gravitational Waves (STGWs)—arising from the coupling between linear scalar perturbations and linear tensor (primordial gravitational wave) perturbations at second order—have been shown to potentially dominate the signal on small scales.
• Specific Objective: The authors aim to compute the STGW energy density generated during:
  1. A generic Matter-Dominated (MD) era.
  2. An Early Matter-Dominated (eMD) epoch followed by a sudden transition to the standard Radiation-Dominated (RD) era (a scenario relevant for Primordial Black Hole formation and "poltergeist" mechanisms).
  3. The standard RD era.
• Goal: To determine if STGWs can explain the NANOGrav signal and forecast the detection prospects and parameter constraints for the future Square Kilometre Array (SKA).

2. Methodology

Theoretical Framework

• Gauge and Formalism: The authors work in the Poisson gauge within General Relativity, assuming a flat FLRW metric with scalar ($\Phi$) and tensor ($\gamma_{ij}$) perturbations. They neglect vector modes and tensor-tensor couplings to isolate the scalar-tensor mixing effect.
• Evolution Equations: They derive the equation of motion for the second-order tensor modes sourced by the coupling of first-order scalar and tensor modes. The source term depends on the transfer functions of the background era (MD or RD).
• Green's Function Approach: The solution for the induced tensor modes is obtained using Green's functions, integrating the source term over conformal time.
• Energy Density Calculation: The dimensionless energy density spectrum $\Omega_{GW}(k)$ is computed by averaging the two-point correlation function of the induced tensor modes over oscillation periods.

Specific Scenarios Analyzed

1. Pure MD Era: The authors derive the kernel for STGWs in a standard matter-dominated universe. They find that, unlike the scalar-induced case, the STGW signal rapidly dilutes ($\propto 1/x^2$) because tensor modes decay while scalar potentials remain constant.
2. eMD to RD Transition: They extend the "poltergeist mechanism" (originally for scalar-induced GWs) to the scalar-tensor case.
   • They model an instantaneous reheating at conformal time $\eta_R$.
   • They derive a specific kernel ($I_{eMD}$) accounting for the matching of tensor modes across the eMD-RD transition.
   • Crucial Distinction: Unlike scalar perturbations, tensor modes do not exhibit the same "poltergeist" enhancement immediately after reheating due to their different decay behavior in the eMD phase. The authors explicitly derive the matching conditions for tensor modes, showing that the standard scalar matching procedure does not apply directly.
3. Regularization: To handle UV divergences in the integral of the energy density (caused by flat power spectra), the authors adopt peaked primordial power spectra (log-normal and monochromatic/Dirac delta) and introduce a damping function $\Upsilon(u,v)$.

Data Analysis and Forecasting

• Datasets:
  • NANOGrav 15-year (NG15): Used for current Bayesian inference.
  • Mock SKA10: A simulated dataset assuming 200 pulsars observed for 10.33 years.
• Statistical Tools:
  • PTArcade: Used for Bayesian inference on NG15 data.
  • fastPTA: Used for forecasting SKA constraints.
• Parameters: The analysis assumes log-normal scalar ($A_\zeta$) and tensor ($A_t$) power spectra with a peak frequency ($f_*$) and width ($\sigma$).
• Constraints: The parameter space is restricted by the requirement to avoid Primordial Black Hole (PBH) overproduction ($f_{PBH} \leq 1$).

3. Key Contributions

1. Derivation of STGWs in MD and eMD: The paper provides the first complete derivation of the STGW energy density kernel in a generic MD era and extends the eMD-RD transition formalism to include scalar-tensor mixing.
2. Dilution in Pure MD: It establishes that STGWs generated during a pure MD era are negligible due to rapid dilution, contrasting with the behavior of SIGWs.
3. Tensor Matching at Reheating: The authors provide a rigorous analytical treatment for matching tensor perturbations across an eMD-RD transition, demonstrating that the "poltergeist" enhancement mechanism for scalar-tensor mixing is distinct from the purely scalar case and requires careful handling of initial conditions.
4. Bayesian Forecasts: The paper presents the first Bayesian parameter estimation for STGWs using NANOGrav data and detailed forecasts for the SKA, including the interplay between second-order STGWs and first-order Primordial Gravitational Waves (PGWs).

4. Key Results

Theoretical Findings

• Pure MD: The STGW energy density vanishes as $\sim 1/x^2$ (where $x=k\eta$) in a pure MD era, making it an unviable source for the PTA signal if generated solely during this epoch.
• eMD Scenario: In an eMD scenario with a sudden transition to RD, the STGW signal remains non-vanishing. The spectral shape is highly sensitive to the width ($\sigma$) of the primordial power spectrum. Narrow widths produce visible oscillations in the spectrum due to the kernel structure, while broad widths smooth these features.

Observational Constraints (NANOGrav 15-year)

• Degeneracy: There is a strong negative correlation (degeneracy) between the scalar amplitude ($A_\zeta$) and tensor amplitude ($A_t$), as the STGW signal depends on their product.
• PBH Constraints: When applying the PBH overproduction bound, the posterior distribution for the scalar amplitude is heavily constrained. In many cases, the region allowed by NANOGrav data is excluded by the requirement that PBHs do not overclose the universe, unless the scalar amplitude is significantly lower than the median best-fit value.
• Shape Parameters: The peak frequency ($f_*$) and width ($\sigma$) are better constrained than the amplitudes, with the data favoring a narrow, curved region in the parameter space.

Future Prospects (SKA)

• Precision: The SKA is forecasted to constrain the spectral shape parameters ($f_*$, $\sigma$) and the tensor amplitude ($A_t$) with high precision (percent-level).
• Scalar Amplitude: The scalar amplitude ($A_\zeta$) remains less constrained (broader posterior) because the linear PGW signal (first-order) fixes the tensor sector, allowing the second-order STGW normalization to vary without spoiling the fit.
• PBH Viability: Unlike the NANOGrav case, the SKA forecasts show that there are viable regions in the parameter space where STGWs dominate the SGWB signal while satisfying PBH constraints, particularly when the linear PGW component is included.
• RD vs. eMD: While SKA can distinguish between RD and eMD origins, the constraints on shape parameters are slightly tighter for the eMD case due to the explicit dependence of the kernel on the reheating time.

5. Significance and Conclusion

• Viability of STGWs: The study confirms that Scalar-Tensor-Induced Gravitational Waves are a viable and potentially dominant source of the nHz SGWB, provided they are generated during an eMD epoch (or RD) rather than a pure MD epoch.
• Distinctive Signatures: The paper highlights that STGWs carry unique spectral signatures (oscillatory features dependent on $\sigma$) that can be resolved by future PTA experiments like SKA, distinguishing them from other cosmological sources.
• Interplay of Orders: The analysis underscores the importance of considering both linear (PGW) and second-order (STGW) contributions simultaneously. The presence of a linear PGW signal helps break degeneracies and relaxes constraints on the scalar amplitude, making the STGW scenario more consistent with PBH limits.
• Future Directions: The authors note that fully numerical treatments of the eMD-RD transition and the inclusion of third-order contributions are necessary for future refinements, particularly to fully resolve the "poltergeist" analogue for tensor modes.

In summary, this paper provides a comprehensive theoretical and statistical framework for testing the scalar-tensor origin of the PTA signal, demonstrating that while current data (NANOGrav) faces tension with PBH limits, future observations (SKA) will likely be able to detect and characterize this signal with high precision.