Jun-Peng Li, Sai Wang, Zhi-Chao Zhao, Kazunori Kohri

The scalar-induced gravitational waves (SIGWs) produced nonlinearly by the enhanced cosmological curvature perturbations can serve as a potentially powerful probe of primordial non-Gaussianity (PNG) in the early Universe. In this work, we comprehensively investigate the imprints of local-type PNG on the SIGW background beyond the widely used quadratic and cubic approximations. We develop a diagrammatic approach capable of analyzing SIGWs for PNG up to any order. Utilizing this approach, we derive semi-analytic formulas for the energy-density fraction spectrum of the isotropic background, the angular power spectrum of the energy-density anisotropies, as well as the angular bispectrum and trispectrum of the energy-density non-Gaussianity. Specializing to PNG up to quartic approximation (parameterized by $f_\mathrm{NL}$, $g_\mathrm{NL}$, and $h_\mathrm{NL}$), we numerically compute all contributions to these SIGW spectra. We find that PNG can significantly alter the amplitude of the energy-density spectrum and generate substantial anisotropies through the initial inhomogeneities. Furthermore, we observe that the angular bispectrum and trispectrum always vanish when the primordial curvature perturbations are Gaussian; otherwise, they do not, indicating their potential utility as probes of PNG. Therefore, we anticipate that both the isotropic background and the anisotropies of the SIGW background will provide essential information about the early Universe.