Let $X^N = (X_1^N,\dots, X^N_d)$ be a d-tuple of $N\times N$ independent GUE random matrices and $Z^{NM}$ be any family of deterministic matrices in $\mathbb{M}_N(\mathbb{C})\otimes \mathbb{M}_M(\mathbb{C})$. Let $P$ be a self-adjoint non-commutative polynomial. A seminal work of Voiculescu shows that the empirical measure of the eigenvalues of $P(X^N)$ converges towards a deterministic measure defined thanks to free probability theory. Let now $f$ be a smooth function, the main technical result of this paper is a precise bound of the difference between the expectation of $$\frac{1}{MN}\text{Tr}\left( f(P(X^N\otimes I_M,Z^{NM})) \right)$$ and its limit when $N$ goes to infinity. If $f$ is six times differentiable, we show that it is bounded by $M^2\left\Vert f\right\Vert_{\mathcal{C}^6}N^{-2}$. As a corollary we obtain a new proof of a result of Haagerup and Thorbj\o rnsen, later developed by Male, which gives sufficient conditions for the operator norm of a polynomial evaluated in $(X^N,Z^{NM},{Z^{NM}}^*)$ to converge almost surely towards its free limit. Restricting ourselves to polynomials in independent GUE matrices, we give concentration estimates on the largest eingenvalue of these polynomials around their free limit. A direct consequence of these inequalities is that there exists some $\beta>0$ such that for any $\varepsilon_1<3+\beta)^{-1}$ and $\varepsilon_2<1/4$, almost surely for $N$ large enough, $$-\frac{1}{N^{\varepsilon_1}} \| P(X^N)\| - \left\Vert P(x)\right\Vert \leq \frac{1}{N^{\varepsilon_2}}.$$ Finally if $X^N$ and $Y^{M_N}$ are independent and $M_N = o(N^{1/3})$, then almost surely, the norm of any polynomial in $(X^N\otimes I_{M_N},I_N\otimes Y^{M_N})$ converges almost surely towards its free limit. This result is an improvement of a Theorem of Pisier, who was himself using estimates from Haagerup and Thorbj\o rnsen, where $M_N$ had size $o(N^{1/4})$.