Title: Dirac Spectral Density in $N_f = 2+1$ QCD at $T = 230$ MeV

Abstract: We compute the renormalized Dirac spectral density in $N_f = 2+1$ QCD at physical quark masses, temperature $T = 230$ MeV and system size $L_s = 3.4$ fm. To that end, we perform a point-wise continuum limit of the staggered density in lattice QCD with staggered quarks. We find, for the first time, that a clear infrared structure (IR peak) emerges in the density of Dirac operator describing dynamical quarks. Features of this structure are consistent with those previously attributed to the recently-proposed IR phase of thermal QCD. Our results (i) provide solid evidence that these IR features are stable and physical; (ii) improve the upper bound for IR-phase transition temperature $T_{\mathrm{IR}}$ so that the new window is $200 < T_{\mathrm{IR}} < 230\,$MeV; (iii) support non-restoration of anomalous U$_{\mathrm A}$(1) symmetry (chiral limit) below $T \!=\! 230\,$MeV.