Gravitational Waves from Isocurvature Perturbations of Spectator Scalar Fields

TL;DR

This work demonstrates that blue-tilted isocurvature perturbations from a subdominant spectator scalar field during inflation can induce a stochastic GW background via second-order scalar perturbations. By performing full numerical evolution from inflation through reheating and incorporating gravitational particle production, reheating dynamics, and scalar-induced GW generation, the authors derive the GW spectrum as a sum of time-dependent and momentum-dependent factors, with $oxed{\Omega_{\text{GW}} h^2}$ reaching $10^{-20}-10^{-12}$ across $10^{-20}-1$ Hz for $0.5\lesssim m_{\chi,\mathrm{eff}}/H_I\lesssim 0.7$. They show that GW-induced constraints can be stronger than isocurvature bounds, and that the signal depends sensitively on reheating temperature and inflaton-spectator couplings, yielding distinct signatures in curvaton-like decays and DM scenarios. The results connect inflationary dynamics, dark matter production, and reheating through multi-frequency GW observations, offering a powerful, testable window into the early universe with prospects across PTAs, space-based interferometers, and next-generation CMB experiments.

Abstract

We present a mechanism for gravitational wave (GW) production from isocurvature perturbations in spectator scalar fields during inflation. These energetically subdominant fields develop blue-tilted power spectra through inflationary dynamics, generating second-order scalar perturbations that source a stochastic GW background. The mechanism naturally satisfies CMB constraints at large scales while producing enhanced signals at smaller scales across a broad frequency range $10^{-20} - 1$ Hz. We perform comprehensive numerical and analytical calculations of the complete isocurvature spectrum evolution, including gravitational particle production, reheating dynamics, and scalar-induced GW generation. For spectator fields with effective masses $0.5 \lesssim m_{χ,\mathrm{eff}}/H_I$, the resulting GW energy density reaches $Ω_{\text{GW}} h^2 \sim 10^{-20}$-$10^{-12}$, accessible to pulsar timing arrays, space-based interferometers, and next-generation CMB experiments. Our analysis reveals that GW-induced constraints exceed current isocurvature bounds. We examine both unstable (curvaton-like) and stable (dark matter) spectator fields, demonstrating strong sensitivity to reheating temperature, inflaton-spectator coupling, and decay dynamics. This framework establishes isocurvature-sourced GWs as a powerful probe of early universe physics, enabling simultaneous constraints on inflationary dynamics, dark matter production, and reheating through coordinated multi-frequency GW observations.

Figures (17)

• Figure 1: Energy densities of the inflaton, $\rho_{\phi}$, the visible sector radiation, $\rho_R$, and the spectator scalar, $\rho_{\chi}$, as functions of the number of $e$-folds, for a high reheating temperature, $T_{\rm reh}=10^{13}\ {\rm GeV}$. The continuous black curve follows $\rho_{\chi}$ until its decay, which is assumed to occur when $f_{\chi}(t_d)=10^{-2}$. The dashed line follows $\rho_{\chi}$ in the absence of a decay.
• Figure 2: Integration domain of Eq. (\ref{['eq:PSp4']}), which corresponds to the case $2k_{\rm IR} < k < k_{\rm UV} - k_{\rm IR}$. This condition ensures that the integration limits cover the regions of interest for both low and high momenta avoiding potential divergences in both IR and UV regimes. The domain is divided into three parts for numerical evaluation: (1) $k_{\rm IR} \leq p \leq k/2$ with $k-p \leq q \leq k+p$, (2) $k/2 \leq p \leq k_{\rm UV}-k$ with $p \leq q \leq k+p$, and (3) $k_{\rm UV}-k \leq p \leq k_{\rm UV}$ with $p \leq q \leq k_{\rm UV}$.
• Figure 3: Isocurvature power spectrum for the spectator field $\chi$ for pure minimal gravitational production (left), and for a fixed mass with varying inflaton coupling (right). Each mass parameter (left panel) and coupling strength (right panel) is represented by a different color. Here $N_{\rm tot}=76.5$$e$-folds of inflation. The vertical dashed lines show the positions of the present horizon scale ($k_0 = 1.4 \times 10^{-26} k_{\rm end}$) and the Planck pivot scale ($k_* = 3.2\times 10^{-24}k_{\rm end}$). The horizontal dotted line represents the current upper bound on isocurvature perturbations from Planck, $\Delta_{\mathcal{S}}^2 \simeq 8.3 \times 10^{-11}$ at the Planck pivot scale (\ref{['eq:isocurvatureupperlimit']}).
• Figure 4: Numerical computation of the exact isocurvature power spectrum (solid lines), compared with two analytical approximations given by Eq. (\ref{['eq:PSan2']}) (black dashed lines) and Eq. (\ref{['eq:PSan1']}) (red dotted lines), for representative values of spectator field masses and inflaton couplings. For the coupled case shown in the right panel, we implement a phenomenological re-scaling $m_{\chi,{\rm eff}}\rightarrow0.9 m_{\chi,{\rm eff}}$ to improve the accuracy of the analytical approximation.
• Figure 5: Left: Time ($e$-fold) evolution of the isocurvature power spectrum for the spectator field $\chi$, with $m_{\chi}/H_{\rm end} = 10^{-2}$ and $\sigma/\lambda = 0.032$, shown at different number of $e$-folds during and after inflation. Continuous lines represent super-horizon modes ($k < a H$), which maintain their spectral shape while undergoing amplitude modulation. Dashed lines indicate sub-horizon modes ($k > a H$) that exhibit the characteristic $k^3$ ultraviolet behavior arising from normal-ordering regularization (see the discussion around Eq. (\ref{['eq:UVk']})). The spectrum decreases during inflation but rapidly grows during the early reheating phase as the field becomes effectively non-relativistic. Right: Time evolution of the function $\mathcal{T}_{\mathcal{S}}$ during reheating (see Eq. (\ref{['eq:TS']})).