Spectral correlations of the massive QCD Dirac operator at finite temperature

We use the graded eigenvalue method, a variant of the supersymmetry technique, to compute the universal spectral correlations of the QCD Dirac operator in the presence of massive dynamical quarks. The calculation is done for the chiral Gaussian unitary ensemble of random matrix theory with an arbitrary Hermitian matrix added to the Dirac matrix. This case is of interest for schematic models of OCD at finite temperature.