Intrinsic non-Gaussianity of ultra slow-roll inflation

We study the non-Gaussian tail of the curvature fluctuation, $\zeta$, in an inflationary scenario with a transient ultra slow-roll phase that generates a localized large enhancement of the spectrum of $\zeta$. To do so, we implement a numerical procedure that provides the probability distribution of $\zeta$ order by order in perturbation theory. The non-Gaussianities of $\zeta$ can be shown to arise from its non-linear relation to the inflaton fluctuations and from the intrinsic non-Gaussianities of the latter, which stem from its self interactions. We find that intrinsic non-Gaussianities, which have often been ignored to estimate the abundance of primordial black holes in this kind of scenario, are important. The relevance of the intrinsic contribution depends on the rapidity with which the transient ultra slow-roll phase occurs, as well as on its duration. Our method cannot be used accurately when the perturbative in-in formalism fails to apply, highlighting the relevance of developing fully non-perturbative approaches to the problem.