Spectrum of the product of independent random Gaussian matrices

We show that the eigenvalue density of a product X=X1X2...XM of M independent NxN Gaussian random matrices in the limit N-->infinity is rotationally symmetric in the complex plane and is given by a simple expression rho(z,z)=1/Mpisigma(-2/M)|z|(-2+(2/M)) for |z|sigma. The parameter sigma corresponds to the radius of the circular support and is related to the amplitude of the Gaussian fluctuations. This form of the eigenvalue density is highly universal. It is identical for products of Gaussian Hermitian, non-Hermitian, and real or complex random matrices. It does not change even if the matrices in the product are taken from different Gaussian ensembles. We present a self-contained derivation of this result using a planar diagrammatic technique. Additionally, we conjecture that this distribution also holds for any matrices whose elements are independent centered random variables with a finite variance or even more generally for matrices which fulfill Pastur-Lindeberg's condition. We provide a numerical evidence supporting this conjecture.