Saved in:
Bibliographic Details
Main Author: Piccoli, Giovanni
Format: Preprint
Published: 2026
Subjects:
Online Access:https://arxiv.org/abs/2605.27231
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866917536636862464
author Piccoli, Giovanni
author_facet Piccoli, Giovanni
contents Scalar Induced Gravitational Waves (SIGW) are generated at second order in perturbation theory and to achieve observational relevance, inflationary dynamics must evade the standard slow-roll scenario at small scales, generating large curvature perturbations following strongly non-Gaussian statistics. We propose a method to efficiently compute the SIGW spectrum including arbitrary non-Gaussianities. First, we solve the wave equation adopting semi-analytic methods; this results in an expression involving integrals in Fourier space which are impossible to solve directly on a lattice. We overcome this bottleneck by recasting these integrals as a sum of about 50 convolutions, each of which can be computed efficiently with FFT methods. Finally, the power spectrum is measured directly from the lattice realization. We implement this in FLAN-SIGW, a GPU-accelerated code capable of computing fully non-perturbative, non-Gaussian SIGW spectra in seconds with an error within 10% with modest computational resources. The code is made public at https://github.com/giovannipiccoli99/FLAN-SIGW. In this first implementation, in order to assess the performance of the method, we adopt a standard radiation-dominated background with $w = 1/3$.
format Preprint
id arxiv_https___arxiv_org_abs_2605_27231
institution arXiv
publishDate 2026
record_format arxiv
spellingShingle A Fast Method to Compute Scalar Induced Gravitational Waves on a Lattice with Primordial Non-Gaussianities
Piccoli, Giovanni
Cosmology and Nongalactic Astrophysics
General Relativity and Quantum Cosmology
Computational Physics
Scalar Induced Gravitational Waves (SIGW) are generated at second order in perturbation theory and to achieve observational relevance, inflationary dynamics must evade the standard slow-roll scenario at small scales, generating large curvature perturbations following strongly non-Gaussian statistics. We propose a method to efficiently compute the SIGW spectrum including arbitrary non-Gaussianities. First, we solve the wave equation adopting semi-analytic methods; this results in an expression involving integrals in Fourier space which are impossible to solve directly on a lattice. We overcome this bottleneck by recasting these integrals as a sum of about 50 convolutions, each of which can be computed efficiently with FFT methods. Finally, the power spectrum is measured directly from the lattice realization. We implement this in FLAN-SIGW, a GPU-accelerated code capable of computing fully non-perturbative, non-Gaussian SIGW spectra in seconds with an error within 10% with modest computational resources. The code is made public at https://github.com/giovannipiccoli99/FLAN-SIGW. In this first implementation, in order to assess the performance of the method, we adopt a standard radiation-dominated background with $w = 1/3$.
title A Fast Method to Compute Scalar Induced Gravitational Waves on a Lattice with Primordial Non-Gaussianities
topic Cosmology and Nongalactic Astrophysics
General Relativity and Quantum Cosmology
Computational Physics
url https://arxiv.org/abs/2605.27231