Correlated Dirac eigenvalues around the transition temperature on $N_{\tau}=8$ lattices

TL;DR

This study examines how the Dirac eigenvalue spectrum's derivatives with respect to light quark mass reveal the nature of the chiral phase transition near the critical temperature in lattice QCD, challenging the dilute instanton gas approximation.

Contribution

It introduces a method to analyze the temperature dependence of Dirac eigenvalue spectrum derivatives and estimates the power-law behavior of the spectrum near the transition.

Findings

Deviations from proportionality between derivatives at lower temperatures.

Negative second derivatives indicating breakdown of dilute instanton gas approximation.

A proposed method to estimate the power-law exponent of the spectrum.

Abstract

We investigate the criticality of chiral phase transition manifested in the first and second order derivatives of Dirac eigenvalue spectrum with respect to light quark mass in (2+1)-flavor lattice QCD. Simulations are performed at temperatures from about 137 MeV to 176 MeV on $N_{\tau}=8$ lattices using the highly improved staggered quarks and the tree-level improved Symanzik gauge action. The strange quark mass is fixed to its physical value $m_s^{\text{phy}}$ and the light quark mass is set to $m_s^{\text{phy}}/40$ which corresponds to a Goldstone pion mass $m_{\pi}=110$ MeV. We find that in contrast to the case at $T\simeq 205$ MeV $m_l^{-1} \partial \rho(\lambda, m_l)/\partial m_l$ is no longer equal to $\partial ^2\rho(\lambda, m_l)/\partial m_l^2$ and $\partial ^2\rho(\lambda, m_l)/\partial m_l^2$ even becomes negative at certain low temperatures. This means that as temperature getting closer to $T_c$ $\rho(\lambda, m_l)$ is no longer proportional to $m_l^2$ and thus dilute instanton gas approximation is not valid for these temperatures. We demonstrate the temperature dependence can be factored out in $\partial \rho(\lambda, m_l)/ \partial m_l$ and $\partial^2 \rho(\lambda, m_l)/ \partial m_l^2$ at $T \in [137, 153]$ MeV, and then we propose a feasible method to estimate the power $c$ given $\rho \propto m_l^{c}$.