We extend key results from free probability on polynomials to arbitrary Sobolev functions. Applied to time Heisenberg evolution operators this proves the thermalization effect of the unitary group generated by Wigner matrices.

We determine the approximate location of the extreme eigenvalues for a large class of random matrix models. These two-sided bounds are fundamentally beyond the reach of classical matrix concentration inequalities.