Poisson statistics in the high temperature QCD Dirac spectrum

We analyze the eigenvalue statistics of the staggered Dirac operator above $T_{c}$ in QCD with 2+1 flavors of dynamical quarks. We use physical quark masses in our simulations. We compare the eigenvalue statistics from several parts of the Dirac spectrum with the predictions of Random Matrix Theory for this universality class and with Poisson statistics. We show that at the low end of the spectrum the eigenmodes are localized and obey Poisson statistics. Above a boundary region the eigenmodes become delocalized and obey Random Matrix statistics. Thus the QCD Dirac spectrum with physical dynamical quarks also has the Poisson to Random Matrix transition previously seen in the quenched SU(2) theory.