Mesoscopic spectral CLT for block correlated random matrices

For random matrices with block correlation structure we show that the fluctuations of linear eigenvalue statistics are Gaussian on all mesoscopic scales with universal variance which coincides with that of the Gaussian unitary or Gaussian orthogonal ensemble, depending on the symmetry class of the model. The main tool used for determining this variance is a two-point version of the matrix-valued Dyson equation, that encodes the asymptotic behavior of the product of resolvents at different spectral parameters.