Spectral functions at small energies and the electrical conductivity in hot, quenched lattice QCD

TL;DR

It is shown that at finite temperature the most commonly used algorithm, employing Bryan's method, is inherently unstable at small energies and a modification is given that avoids this.

Abstract

In lattice QCD, the Maximum Entropy Method can be used to reconstruct spectral functions from euclidean correlators obtained in numerical simulations. We show that at finite temperature the most commonly used algorithm, employing Bryan's method, is inherently unstable at small energies and give a modification that avoids this. We demonstrate this approach using the vector current-current correlator obtained in quenched QCD at finite temperature. Our first results indicate a small electrical conductivity above the deconfinement transition.

Figures (4)

• Figure 1: First four basis functions $u_i(\omega)$ as a function of $a\omega$ for $a\omega_{\rm max}=5$, $N_\omega=1000$, $N_\tau=24$, using the standard kernel. The inset shows a blow-up of the small energy region.
• Figure 2: As in Fig. \ref{['fig:svd']}, using the redefined kernel $\overline K(\omega,\tau)$.
• Figure 3: Vector spectral functions $\rho(\omega)/\omega T$ as a function of $\omega/T$. We used $N_\omega=1000$, $a\omega_{\rm max}=5$, $b=1$.
• Figure 4: Default model dependence of $\rho(\omega)/\omega T$ for $N_\tau=24$ (hot) and $16$ (very hot) in the low-energy region. We show results for $N_\omega=1000, 2000$ and $b=1.0, 0.5, 0.1$ at fixed $a\omega_{\rm max}=5$.