Multiplying unitary random matrices - universality and spectral properties

In this paper we calculate, in the large N limit, the eigenvalue density of an infinite product of random unitary matrices, each of them generated by a random hermitian matrix. This is equivalent to solving unitary diffusion generated by a hamiltonian random in time. We find that the result is universal and depends only on the second moment of the generator of the stochastic evolution. We find indications of critical behavior (eigenvalue spacing scaling like $1/N^{3/4}$) close to $\theta=\pi$ for a specific critical evolution time $t_c$.