In this article we obtain a superexponential rate of convergence in total variation between the traces first m powers n×n random unitary matrices and 2m-dimensional Gaussian variable. This generalizes previous results scalar case to multivariate setting, also give precise dependence on dimensions n estimates with explicit constants. We are especially interested regime where grows our main resul...