Non-central limit of densities of some functionals of Gaussian processes

We establish the convergence of the densities of a sequence of nonlinear functionals of an underlying Gaussian process to the density of a Gamma distribution. We provide precise estimates on the distance between densities while developing the techniques of Malliavin calculus and Stein's method suitable to Gamma approximation. We first focus our study on the case of random variables living in a fixed Wiener chaos of an even order for which the bound for the difference of the densities can be dominated by a linear combination of moments up to order four. We then study the case of general Gaussian functionals with possibly infinite chaos expansion. Finally, we provide two applications of our results, namely the limit of random variables living in the second Wiener chaos, and the limit of sums of independent Pareto distributed random variables.