Gravitational Waves from Matter Perturbations of Spectator Scalar Fields

Abstract: We compute the stochastic gravitational wave background sourced at second order by a spectator scalar field $χ$ coupled to the inflaton $φ$ through a portal interaction $σφ^2χ^2$ and with quartic self-interaction $λχχ^4/4!$. In the large portal coupling regime ($σ/λ\gg 1$, with $λ$ the inflaton normalization), parametric resonance during reheating amplifies the spectator power spectrum by many orders of magnitude near the resonance band until Hartree backreaction from the quartic coupling detunes the instability, while the large inflationary effective mass suppresses superhorizon power and ensures compatibility with CMB isocurvature bounds. We focus on the direct field-gradient source $\partial_aχ\,\partial_bχ$ in the second-order Einstein equations and derive a master formula that factorizes into a spectral integral over the frozen, vacuum-subtracted spectator spectrum and a time integral encoding the post-inflationary expansion history. For our benchmark reheating history we obtain analytic scaling relations, including a peak amplitude $Ω{\rm GW}\propto T_{\rm reh}^{8/3}$, strong dependence on the portal strength, and weak sensitivity to $m_χ$. We validate the framework against nonlinear lattice simulations, demonstrating complementarity: the Hartree treatment captures superhorizon evolution inaccessible to the lattice, while the lattice resolves rescattering and fragmentation near the spectral peak. For $σ/λ\simeq 10^4$ and $T_{\rm reh}=2 \times 10^{14}\,\mathrm{GeV}$, the signal reaches $Ω{\rm GW}h^2\sim 10^{-11}$ at $f\sim10^{7}$-$10^{8}\,\mathrm{Hz}$. Increasing $λχ$ at fixed $σ$ has a non-monotonic effect: small values enhance the signal via rescattering, whereas larger values suppress it by detuning the resonance.

• The paper presents a novel analytic-numerical framework to compute the SGWB from amplified spectator scalar field fluctuations during reheating.
• It employs both Hartree approximations and lattice simulations to capture nonlinear dynamics and parametric resonance effects.
• The gravitational wave signal exhibits a steep blue tilt with amplitude ΩGWh² around 10⁻¹¹ at frequencies of ~10⁷–10⁸ Hz, probing reheating microphysics.

Gravitational Wave Production from Spectator Scalar Field Matter Perturbations

Introduction and Motivation

This work ["Gravitational Waves from Matter Perturbations of Spectator Scalar Fields" (2604.05078)] provides a systematic computation of the stochastic gravitational wave background (SGWB) generated at second order during reheating by spatial gradients in the fluctuations of a spectator scalar field $\chi$. These fields, minimally coupled and energetically subdominant during inflation, are coupled to the inflaton $\phi$ via a portal $\sigma \phi^2 \chi^2$ and may possess a quartic self-interaction $\lambda_\chi \chi^4 / 4!$. The analysis targets the nontrivial, high-frequency SGWB sourced not by the usual linearly-generated tensor metrics, but by field-gradient terms arising from parametric resonance and subsequent non-equilibrium dynamics in the early universe.

Parametric resonance, regulated by Hartree backreaction, leads to many orders of magnitude amplification of $\chi$ fluctuations near the resonance band. Importantly, the large inflationary mass of $\chi$ suppresses superhorizon power, keeping CMB isocurvature under control—this allows strong small-scale features without violating large-scale constraints. By focusing on the direct field-gradient source in the Einstein equations, the study develops a factorized (spectral × temporal) analytic framework complemented and validated by nonlinear lattice simulations.

Dynamics of Energy Transfer and Spectator Field Evolution

The background cosmology assumes a standard inflationary sector transitioning to reheating via oscillating inflaton dynamics. The benchmark inflaton model is a T-model $\alpha$-attractor with quadratic potential, yielding post-inflation matter domination.

The energy densities of the primary and spectator components are shown in the following:

Figure 1: Energy densities of the inflaton $\rho_\phi$, visible radiation $\rho_R$, and spectator field $\rho_\chi$ as functions of $\phi$0-folds, highlighting dynamics with and without $\phi$1 self-interaction.

The onset of reheating triggers substantial energy transfer via parametric resonance. The efficiency and temporal structure of this process depend strongly on the $\phi$2-$\phi$3 portal ($\phi$4) and the $\phi$5 self-coupling ($\phi$6):

Figure 2: Early reheating evolution of $\phi$7 for varying $\phi$8 (left) and $\phi$9 (right), compared with $\sigma \phi^2 \chi^2$0 (black). Both Hartree and lattice methods shown.

Key features demonstrated include non-monotonic dependence on $\sigma \phi^2 \chi^2$1 due to resonance band structure and fragmentation, and strong suppression of amplification for large $\sigma \phi^2 \chi^2$2 due to Hartree detuning. Lattice simulations confirm Hartree for IR dynamics but reveal significant rescattering beyond mean-field at high $\sigma \phi^2 \chi^2$3.

Spectator Power Spectrum and Non-Equilibrium Enhancement

The post-inflation amplification yields a frozen, vacuum-subtracted spectator spectrum sharply peaked near the resonance scale, with a universal $\sigma \phi^2 \chi^2$4 white-noise IR scaling inherited from the large effective mass regime. This behavior is corroborated numerically:

Figure 3: Time evolution of the rescaled $\sigma \phi^2 \chi^2$5 power spectrum during resonance, and transfer function for IR growth.

The exponential enhancement—up to 15 orders of magnitude over a few $\sigma \phi^2 \chi^2$6-folds—is shut off either by Hubble redshifting of the inflaton envelope or by backreaction from $\sigma \phi^2 \chi^2$7 self-interactions. Lattice and Hartree results are congruent for IR, with UV deviations arising from turbulent energy transfer to short modes.

Gravitational Wave Production Formalism

The direct $\sigma \phi^2 \chi^2$8 source in the Einstein equations is isolated and treated with a vacuum-subtracted, four-point Wick contraction. The resulting SGWB spectrum factorizes into a momentum integral over the frozen spectrum and a temporal integral governed by the expansion history. The spectral integral, $\sigma \phi^2 \chi^2$9, admits a $\lambda_\chi \chi^4 / 4!$0 IR scaling for white noise input, which is validated numerically:

Figure 4: Exact and analytic approximation for the GW spectral integral $\lambda_\chi \chi^4 / 4!$1. The IR scaling $\lambda_\chi \chi^4 / 4!$2 dominates until resonance cutoff.

The time-dependent build-up of gravitational wave energy is decomposed into contributions from: early-time (parametric amplification), matter-dominated reheating, and the subsequent radiation-dominated era.

Figure 5: Relative contributions of each phase to the time-dependent factor $\lambda_\chi \chi^4 / 4!$3 for two $\lambda_\chi \chi^4 / 4!$4 benchmarks.

The dominant source is set by the largest available $\lambda_\chi \chi^4 / 4!$5 and amplifies at higher frequencies, with the spectrum sharply rising as $\lambda_\chi \chi^4 / 4!$6 near the resonance scale.

Gravitational-Wave Signal Predictions and Parameter Dependence

The resulting gravitational wave energy density exhibits a steep blue tilt in the IR ($\lambda_\chi \chi^4 / 4!$7), peaking at ultra-high frequencies $\lambda_\chi \chi^4 / 4!$8--$\lambda_\chi \chi^4 / 4!$9 Hz. The total amplitude can reach $\chi$0 for strong coupling ($\chi$1) and $\chi$2 GeV, as shown:

Figure 6: GW spectra for varying $\chi$3 ($\chi$4), with lattice (dashed) and analytic (dot-dashed) results.

The dependence on $\chi$5 is highly non-monotonic, mirroring resonance band structure:

Figure 7: GW spectra for varying portal coupling ($\chi$6). Black dashed: lattice spectra.

The GW signal is largely insensitive to $\chi$7 as long as parametric resonance is efficiently driven by the portal.

Figure 8: GW spectrum as a function of $\chi$8 for fixed $\chi$9 and $\chi$0.

When $\chi$1 self-interaction is increased, the GW amplitude is suppressed for large $\chi$2 due to earlier resonance shutdown; at intermediate values, mode-mode rescattering can boost UV power, an effect fully captured by lattice simulations:

Figure 9: GW spectra for varying $\chi$3, showing enhancement at intermediate and suppression at high $\chi$4. Black dashed: lattice computations.

Implications and Theoretical Significance

The analysis demonstrates that second-order SGWB from spectator sector matter gradients can achieve significant amplitudes at ultra-high frequencies—orders of magnitude above those attainable by the linear inflationary background. The sensitivity to $\chi$5, $\chi$6, and $\chi$7 highlights the potential for GW detection to probe reheating microphysics and dark sector interactions inaccessible to electromagnetic or even large-scale structure probes. The signal's frequency localization at $\chi$8 LIGO/LISA bands places it beyond current detector capabilities, but amplitudes $\chi$9 are sufficient to motivate future resonant-cavity or other ultra-high-frequency GW concepts.

This framework is robust against CMB isocurvature constraints due to the suppression of superhorizon $\alpha$0 power by heavy effective mass during inflation. Lattice simulations validate the IR predictions of the analytic Hartree approach and reveal the role of nonlinear fragmentation and rescattering—essential for correct UV behavior, but subdominant for IR-dominated GW production.

Extension to alternative inflationary potentials (e.g., quartic) and to curvaton-like or stable relic scenarios are possible, with enhancements generally anticipated if post-inflation expansion is more radiation-like.

Conclusion

This study establishes an analytic and numerical approach for predicting gravitational wave backgrounds from spectator scalar matter perturbations amplified during reheating via inflaton-induced parametric resonance. The main findings include:

• Parametric resonance, strongly regulated by $\alpha$1 and $\alpha$2, can enhance the GW background by $\alpha$3 orders of magnitude relative to the gravitational production baseline.
• Maximum GW amplitude scales as $\alpha$4, reaching $\alpha$5 at $\alpha$6--$\alpha$7 Hz, well above BBN/CMB bounds but outside current detector reach.
• The signal is highly sensitive to inflaton-spectator portal structure and reheating dynamics, and only weakly on $\alpha$8 mass in the strong-coupling regime.
• Lattice and Hartree approaches agree in the IR, with lattice revealing important additional features at high $\alpha$9 and in the UV.
• The framework developed provides a target and motivation for continued experimental innovation in ultra-high-frequency GW detection and offers a cosmological probe of reheating and dark-sector microphysics otherwise inaccessible to conventional means.

The methodology invites future generalization to nonminimal couplings, full metric-sourced (isocurvature) second-order tensor channels, and broader exploration of parameter space in contemporary and next-generation cosmological studies.