Finite Temperature Landau Gauge Lattice Quark Propagator

TL;DR

This study examines the quark propagator at finite temperature using lattice QCD, revealing how it varies across the deconfinement transition and linking gluon dynamics to chiral symmetry breaking.

Contribution

It provides new insights into the temperature dependence of the quark propagator and its relation to deconfinement and chiral symmetry restoration in lattice QCD.

Findings

Significant differences in propagator form factors below and above $T_c$.

Strong correlation between gluon dynamics and chiral symmetry breaking.

Support for quarks as free quasi-particles above $T_c$ at low momenta.

Abstract

The quark propagator at finite temperature is investigated using quenched gauge configurations. The propagator form factors are investigated for temperatures above and below the gluon deconfinement temperature $T_c$ and for the various Matsubara frequencies. Significant differences between the functional behaviour below and above $T_c$ are observed both for the quark wave function and the running quark mass. The results for the running quark mass indicate a strong link between gluon dynamics, the mechanism for chiral symmetry breaking and the deconfinement mechanism. For temperatures above $T_c$ and for low momenta, our results support also a description of quarks as free quasi-particles.

Figures

Tables

Table 1. Table 1: Lattice setup for the computation of the quark propagator in the quenched approximation. The coefficient c s ​ w subscript 𝑐 𝑠 𝑤 c_{sw} associated with the non-perturbative improvement of the simulations and κ c subscript 𝜅 𝑐 \kappa_{c} are taken from  Luscher:1996sc . The bare quark mass m b ​ a ​ r ​ e subscript 𝑚 𝑏 𝑎 𝑟 𝑒 m_{bare} is computed using Eq. ( 9 ). T 𝑇 T and m b ​ a ​ r ​ e subscript 𝑚 𝑏 𝑎 𝑟 𝑒 m_{bare} are given in MeV, the lattice spacing a 𝑎 a is given in fm and the inverse of the lattice spacing 1 / a 1 𝑎 1/a in GeV.

T	$\beta$	$L_s^3\times L_t$	$\kappa$	${\kappa }_c$	$a$	$1/a$	$m_{bare}$	$c_{sw}$
243	6.0000	${64}^3\times 8$	0.1350	0.13520	0.1016	1.9426	10	1.769
			0.1342				53
260	6.0347	${68}^3\times 8$	0.1351	0.13530	0.09502	2.0767	11	1.734
			0.1344				51
275	6.0684	${72}^3\times 8$	0.1352	0.13540	0.08974	2.1989	12	1.704
			0.1345				54
290	6.1009	${76}^3\times 8$	0.1347	0.13550	0.08502	2.3211	51	1.678
305	6.1326	${80}^3\times 8$	0.1354	0.13559	0.08077	2.4432	13	1.655
			0.1348				53
324	6.0000	${64}^3\times 6$	0.1342	0.13520	0.1016	1.9426	53	1.769

Equations

$$S^{-1}(p_4,\vec{p})=i{\gamma }_4{p}_4\omega (p_4,\vec{p})+i\vec{\gamma }\cdot \vec{p}Z(p_4,\vec{p})+\sigma (p_4,\vec{p})$$

$$S(p_4,\vec{p})=\frac{-i{\gamma }_4{p}_4\omega (p_4,\vec{p})-i\vec{\gamma }\cdot \vec{p}Z(p_4,\vec{p})+\sigma (p_4,\vec{p})}{p_4^2{\omega }^2(p_4,\vec{p})+\left(\vec{p}\cdot \vec{p}\right)Z^2(p_4,\vec{p})+{\sigma }^2(p_4,\vec{p})}.$$

$$S(x,y)=\bigg{(}1+2\,b_{q}\,a\,m\bigg{)}~{}\left\langle~{}L(x)\,M^{-1}_{SW}(x,y)\,R(y)~{}\right\rangle\ ,$$

$$L(x)$$

$$R(x)$$

$$\overrightarrow{/D}(x)\psi (x)$$

$$\overline{\psi }(x)\overleftarrow{/D}(x)$$

$$p_4=\frac{\pi }{L_t}\left(2n+1\right)\text{mbox}and{p}_i=\frac{2\pi }{L_s}n,$$

$$am=\frac{1}{2}\left(\frac{1}{\kappa }-\frac{1}{{\kappa }_c}\right),$$

$$\frac{\omega (p_4,\vec{p})}{p_4^2{\omega }^2(p_4,\vec{p})+\left(\vec{p}\cdot \vec{p}\right)Z^2(p_4,\vec{p})+{\sigma }^2(p_4,\vec{p})},$$

$$\frac{Z(p_4,\vec{p})}{p_4^2{\omega }^2(p_4,\vec{p})+\left(\vec{p}\cdot \vec{p}\right)Z^2(p_4,\vec{p})+{\sigma }^2(p_4,\vec{p})},$$

$$\frac{\sigma (p_4,\vec{p})}{p_4^2{\omega }^2(p_4,\vec{p})+\left(\vec{p}\cdot \vec{p}\right)Z^2(p_4,\vec{p})+{\sigma }^2(p_4,\vec{p})},$$

$$Z_c(p_4,\vec{p})=\frac{Z(p_4,\vec{p})}{\omega (p_4,\vec{p})},\text{mbox}andM(p_4,\vec{p})=\frac{\sigma (p_4,\vec{p})}{\omega (p_4,\vec{p})},$$

$$S^{(0)}(p)=\frac{-ia/\overline{p}A(m,p)+B(m,p)}{a^2{\overline{p}}^2{A}^2(m,p)+B^2(m,p)}.$$

Full text

††journal: Eur. Phys. J. C

1 Introduction and Motivation

The study of strong interactions including temperature and density effects has been driven both by the tentative to understand the QCD dynamics and by a strong experimental program that involves various facilities. From the theory side, this research includes the computation of two-point correlation functions that provide information on spectra, transport properties and other fundamental properties such as the confinement mechanism. Herein, we will not consider the dependence on the density of hadronic matter.

At low temperatures quarks and gluons are confined particles and appear only as constituents of mesons and baryons. At sufficiently high temperatures or densities quark and gluons are expected to become deconfined and behave as free quasiparticles in a new state of matter, the strong coupled quark gluon plasma. At extreme temperatures, due to asymptotic freedom, quarks and gluons can be considered as free particles.

This view of the quark and gluon dynamics for temperatures above the critical temperature $T_{c}$, where the deconfined phase transition takes place, is motivated by asymptotic freedom and lies at the heart of the quasiparticle description used to investigate the thermodynamical properties of hadronic matter for $T>T_{c}$. For temperatures above $T_{c}$ the perturbative approach of the hard thermal loop (HTL) expansion, see e.g. LeBellac:1991cq ; Kapusta:2006pm and Su:2015esa for a recent review, seems to be the framework to understand the dynamics of QCD. Furthermore, some of the HTL predictions have been confirmed by first principles non-perturbative lattice QCD simulations; an example being the behaviour of thermal masses with the temperature of the heat bath.

The standard approach to identify the confined and the deconfined phases relies on the Polyakov loop $L$ which is related to the free energy of quarks $F_{T}$. At $T<T_{c}$, the renormalised $L$ is small, $F_{T}$ is large, and quarks are confined particles. For temperatures above $T_{c}$, $L$ approaches unit, $F_{T}\approx 0$, and quarks become free quasiparticles. The curve $L(T)$ shows that the nature of the confined-deconfined transition is a cross-over for full QCD and first order for the pure Yang-Mills (quenched QCD) McLerran:1980pk ; McLerran:1981pb . If for full QCD the critical temperature is $T_{c}\sim 150-160$ MeV, for its quenched version one gets $T_{c}\sim 270$ MeV. In what concerns gluon dynamics, for the quenched version of the theory, the consideration of other quantities besides the Polyakov loop reproduce both the nature of the transition and the quoted value for the critical temperature, see e.g. Silva:2016onh .

One of the interests to study the quark and gluon two-point correlation functions with temperature comes from the connection between the propagators and confinement that can be linked, for example, by computing the corresponding spectral functions. If above $T_{c}$ quarks and gluon can be viewed as quasiparticles, their propagators should behave differently from the corresponding functions at temperatures below $T_{c}$.

The Landau gauge gluon propagator has been studied using non-perturbative methods for a wide range of temperatures both for full QCD and for the pure gauge theory Huber:2018ned ; Aouane:2011fv ; Silva:2013maa ; Aouane:2012bk ; Silva:2016msq ; Silva:2016onh . These studies show that in the quenched theory the gluon propagator is sensitive to the breaking of center symmetry, that electric and magnetic mass scales are generated dynamically and the gluon thermal mass associated with the electric propagator scales with the temperature according to the prediction of hard thermal loop approach to QCD, namely $m_{D}\propto T$ above the critical temperature. Furthermore, the calculations performed using the continuum formulation of QCD or lattice QCD simulations for the pure gauge sector are in good agreement; see Maas:2011se and references therein. Moreover, the non-perturbative approaches suggest that HTL QCD is a good framework to describe the dynamics of QCD for $T\sim T_{c}$ and above. From these studies one can claim to have now a good picture for the gluon dynamics.

The Landau gauge quark propagator at finite temperature was also studied within the continuum non-perturbative approaches to QCD and using first principles lattice simulations.

The quark gap equation in the Landau gauge was solved relying on gluon propagators obtained from lattice simulations. The lattice data was fitted to a functional form that reproduces the perturbative tail at sufficiently high momenta Ikeda:2001vc ; Mueller:2010ah ; Fischer:2010fx ; Contant:2017gtz . In order to solve the Dyson-Schwinger equation for the quarks, the quark-gluon vertex was parameterised keeping keeping only its tree level tensor structure, i.e. assuming $\Gamma_{\mu}\propto\gamma_{\mu}$ for the vertex, and, once more, ensuring the right perturbative tail at high momenta. The continuum studies have been focused on dynamical mass generation, through the computation of either the quark condensate Fischer:2009wc , the running quark mass, the spectral function Qin:2010pc at the chiral limit and order parameters for the deconfined phase transition. The form factors appearing on the quark propagator as function of the Matsubara frequencies, momentum and temperature have been only briefly reported but they show quite different behaviours above and below the critical temperature. How the various functions that appear on the quark propagator depend on the parameterisations introduced to solve the quark gap equation is not known and, therefore, the results obtained care for further study or confirmation from independent calculations. The interest on the spectral functions at the chiral limit comes from knowing that at high temperature the quark propagator has two sets of poles, corresponding to an usual mass term and a collective plasmino mode, i.e. two different types of dispersion relations can be associated with this propagator Weldon:1999th . Some authors also speculate on the presence of a third ultrasoft fermion mode that should appear on the quark propagator Qin:2010pc ; Nakkagawa:2011ci ; Nakkagawa:2012ip . Evidence for a third mode for $T>T_{c}$ was also observed on Yukawa like models Harada:2008vk , on QED like models where it becomes an ultrasoft mode Satow:2010ia ; Hidaka:2011rz and in low energy effective models of QCD Kitazawa:2005mp ; Kitazawa:2006zi ; Kitazawa:2006vh ; Kitazawa:2007ep where, again, the third mode is ultrasoft.

In what concerns continuum methods, the quark spectral function together with the quark self energy was computed within the framework of functional renormalization group Wang:2018osm . The outcome reveals a spectral function with multiple peaks, whose details and number of maxima depend on the truncation used in the calculation.

The lattice quark propagator at finite temperature, in the Landau gauge, was computed for $\mathcal{O}(a)$-improved Wilson fermions using quenched gauge configurations by a number of authors Hamada:2006ra ; Karsch:2007wc ; Karsch:2007bf ; Hamada:2008zz ; Karsch:2009tp ; Kitazawa:2009uw ; Hamada:2010zz ; Kaczmarek:2012mb . The lattice studies have been focused either on the mass function and on the same spectral function as investigated by the continuum methods. For the quantities studied, at the qualitative level, lattice and continuum results are in good agreement. On the lattice the curve $M/T$ as a function of the temperature, where $M$ is the quark mass, has a discontinuity at $T_{c}$ but its behaviour for $T>T_{c}$ has not been resolved Hamada:2010zz . For Wilson fermions, the authors Hamada:2006ra found a linear grow of $M/T$ above the transition temperature. However, given that there is some ambiguity on the definition of the quark mass and the poor agreement between the values of the Wilson fermion mass and the $\mathcal{O}(a)$-improved Wilson fermion mass, this result demands for confirmation. Furthermore, in what concerns the computation of the quark mass, a plateaux associated to an effective quark mass was not always identified Hamada:2006ra ; Hamada:2008zz suggesting that instead one should consider a running quark mass. The lattice spectral function was computed and investigated assuming a multiple pole ansätze Karsch:2007wc ; Karsch:2007bf ; Karsch:2009tp ; Kitazawa:2009uw ; Kaczmarek:2012mb that favours a pole structure associated with a thermal and a plasmino mode as predicted in Weldon:1999th . The lattice calculations for the propagator used typically spatial physical volumes under ($\sim$2 fm)3 with two simulations performed on physical volumes ($\sim$3.5 fm)3 Karsch:2009tp ; Kaczmarek:2012mb for $T/T_{c}=0.55$ and 1.5 at the chiral limit111The authors call the reader’s attention that in these simulations they used $T_{c}=300$ MeV for the transition temperature (quenched theory) and the current used value is $T_{c}=270$ MeV. A rescaling of the ratios just quoted results in $T/T_{c}=0.61$ and 1.72..

In the current work we compute the lattice Landau gauge quark propagator at finite temperature with non-perturbatively $\mathcal{O}(a)$-improved Wilson fermions with gauge configurations generated with the Wilson action for the pure gauge theory. These configurations were used to investigate the lattice Landau gauge gluon propagator in Silva:2013maa . As in this work, we take $T_{c}=270$ MeV and consider temperatures above and below the critical temperature to study how different the propagator is in the confined and deconfined phase. Our main focus is on the calculation of the three form factors that determine the quark propagator at finite temperature and the running quark mass as a function of momentum and for the various Matsubara frequencies. By providing the form factors as a function of $T$, we aim to understand the difference between the two phases and provide information than can be useful also for the continuum approaches to QCD. No attempt is made to investigate any of the spectral functions that can be associated with the quark propagator.

The computations of the quark propagator reported below use lattices whose spatial physical volume is ($\sim$ 6.5 fm)3 and various lattice spacings around $a\sim 0.1$ fm. The lattice investigations of the Landau gauge gluon and ghost propagators at zero temperature Oliveira:2012eh ; Duarte:2016iko ; Dudal:2018cli , suggest that finite volume effects are under control in the sense that they are small and below the statistical precision of the simulations. We take as definition of the temperature the inverse length on the time direction and typically consider ratios of the spatial and time direction $L_{s}/L_{t}=8$. Due to the anti-periodic boundary conditions along the time direction $p_{4}=\pi\,T\,(2\,n_{t}+1)$ where $n_{t}=0,\,1,\,\dots$ and the spatial momenta are $p_{i}=2\,\pi\,T\,(L_{t}/L_{s})\,n_{i}\approx 0.785\,T\,n_{i}$ for $n_{i}=0,\,1,\dots$ Furthermore, in order to experience the quark propagator at the chiral limit, we report on simulations using two values of the bare quark mass, namely $m_{0}\approx 10$ MeV and 50 MeV.

Our results show clearly that the nature of the quark propagators form factors changes for temperature above $T_{c}$, compared with temperatures in the confined phase. This is reported with detail for both the quark wave function and the running quark mass. In particular, for the running quark mass we find that it is highly suppressed above $T_{c}$, with typical infrared values being about half of the corresponding values for the temperatures below $T_{c}$. Our computation if performed with quenched configurations and this results indicates that the gluon dynamics plays an important role on the mechanism of chiral symmetry breaking. A similar change on the functional form of the finite temperature gluon propagator form factors for temperatures above and below $T_{c}$ was observed also in lattice simulations Aouane:2011fv ; Silva:2013maa . Indeed, these simulations show that the relative importance of its electric and magnetic components are inverted for $T\gtrsim T_{c}$, relative to confined phase, i.e. for small temperatures the electric form factor is larger than the magnetic one, while above $T_{c}$ the dominant form factor is associated with the gluon magnetic component.

The paper is organised as follows. In Sec. 2 we set the notation and definitions used through out the current work and how the various form factors are measured. Further, we detailed the setup of the simulations analysed herein. In Sec. 3 we report our results for the various form factors prior to the estimation of the lattice artefacts for the two quark masses considered in the simulations. The study of the lattice artefacts and the definition of the running quark mass are detailed in Sec. 4.

2 The Quark Propagator and Lattice Setup

On the continuum, the quark propagator is diagonal in color space. At finite temperature, the presence of a thermal bath breaks rotational invariance, and the inverse of the space-spin quark propagator in momentum space reads

[TABLE]

or, equivalently,

[TABLE]

Our analysis of the lattice propagator will assume that we are close to the continuum and, therefore, expressions (1) and (2) can be applied. The form factors $\omega$, $Z$ and $\sigma$ can be accessed by computing traces of the propagator times gamma matrices.

Our calculation of the lattice quark propagator relies on the use of the non-perturbative $\mathcal{O}(a)$ improved clover action Sheikholeslami:1985ij ; Luscher:1996sc with tree-level O(a)-improved sources Heatlie:1990kg ; Skullerud:2000un ; Skullerud:2001aw . The lattice quark propagator in real space reads

[TABLE]

where $M^{-1}_{SW}$ stands for the inverse of the improved fermionic matrix, the rotated sources are given by

[TABLE]

the left and right derivatives are

[TABLE]

$U_{\mu}(x)$ stands for the gauge configuration, $\hat{\mu}$ is the unit vector associated with direction $\mu$ and $a$ is the lattice spacing. For the coefficients associated with the improvement coming from the rotated sources we use their tree level value $b_{q}=c_{q}=1/4$.

The momentum space propagator is obtained from Eq. (3) after a Fourier transformation. Recall that on the lattice the fermionic boundary conditions are periodic in the spatial directions and anti-periodic in time and the available momenta $\vec{p}$ and $p_{4}$ assume the following discrete values

[TABLE]

where $n=0,\,1,\,2,\,\dots$

On the lattice the bare quark mass is given by

[TABLE]

where the critical hopping parameter $\kappa_{c}$ depends on $\beta$, i.e. on the lattice spacing, and is defined as the value of $\kappa$ corresponding to a vanishing mass for the lightest pseudo-scalar meson. In the following we will use $m$ as an indicator of how close the simulation is from the chiral limit. The $\kappa_{c}$ are taken from Luscher:1996ug , interpolating their reported values when necessary.

The traces of propagator give us the form factors

[TABLE]

and by taking ratios of these functions one gets a continuum-like quark wave function and running quark mass defined by

[TABLE]

respectively. Note that we use $\omega$ and not $Z$ to define the ratios $Z_{c}$ and $M$ because the simulation does not allow the computation of $Z(p_{4},\vec{p}=0)$.

For the computation of the Landau gauge quark propagator we use a subset of the quenched gauge configurations, rotated to the Landau gauge, generated for the work Silva:2013maa . We refer the read to this paper for the details on the generation of the Wilson action pure gauge configurations, the rotation to the Landau gauge and on the scale setting for the conversion into physical units222As discussed in Boucaud:2017ksi ; Duarte:2017wte there is a small uncertainty in the definition of the lattice spacing that can havw a small impact on the outcome specially for large statistical ensembles of configurations. In the current work, we do not take into account this uncertainty.. Furthermore, in order to minimise effects due to the breaking of the rotational invariance, we only report functions for momenta that verifies the cuts considered in Silva:2013maa and introduced in Aouane:2011fv ; see also Leinweber:1998uu . In order to improve the signal to noise ratio we perform always a $Z_{3}$ average of the various lattice quantities; for example when reporting a quantity associated $F$ with momenta $(p_{4},p,0,0)$ it refers to $(F(p_{4},p,0,0)+(p_{4},0,p,0)+(p_{4},0,0,p))/3$.

The statistical errors reported for the various quantities were computed using the bootstrap method with a confidence level of 67.5%. When reporting the ‘continuum” functions $Z_{c}$ and $M$, that are defined as ratios of traces of quark propagator multiplied by gamma matrices, an additional cut on the data surviving the momenta cuts is introduced to have a clear picture on the various functions and only the data whose relative error is below 50% is kept.

The lattice setup of the simulations reported below are given in Tab. 1. In all simulations, the results reported are for computations considering 100 gauge configurations rotated to the Landau gauge. Moreover, to increase the signal to noise ratio for all temperatures, with the exception of the highest, the quark propagator was computed for two point sources located at $(0,0,0,0)$ and $(L_{t}/2,0,0,0)$ and the results averaged before performing any analysis. For the highest temperature the quark propagator was computed for a single source located at the origin of the lattice.

For pure gauge theory the deconfinement temperature happens at $T_{c}=270$ MeV, see e.g. Silva:2016onh and references therein. Our quark propagator simulations consider two temperatures below $T_{c}$, a temperature just above the critical temperature and three temperatures clearly above $T_{c}$. In all cases, with the exception of two temperatures, the quark propagator is computed for two bare quark masses $m_{bare}\approx 10$ and $50$ MeV. In this way, one expects to reveal the fundamental properties of quark properties at finite temperature close to the chiral limit.

3 The Lattice Propagator

The lattice form factors (10), (11) and (12) are measured by taking traces of the lattice propagator times gamma matrices. Their behaviour as a function of $p=\sqrt{p^{2}_{4}+\vec{p}^{2}}$ is illustrated in Figs. 1,  2,  3 and  4 for temperatures below and above the critical temperature $T_{c}=270$ MeV and for all the Matsubara frequencies.

For the smallest bare quark mass and for the smallest temperature, the data shows large fluctuations but seems to follow the same pattern as the data for the heaviest bare quark mass. For the smallest $m_{bare}$ and $T$, due to the large fluctuations observed in particular at the lower momenta, some of data does not appear in the Fig. due our choice for the $y$ scale. Around and above $T_{c}$ the behaviour of the form factors for the two $m_{bare}$ is quite similar. The Figs. also show a clear violation of the rotational symmetry that occurs at all temperatures and that is particularly dramatic for (12).

The bare lattice quark wave function $Z_{c}(p_{4},|\vec{p}|)$ for the full range of momenta, and the running quark mass $M(p_{4},|\vec{p}|)$ for momenta up to 1.5 GeV can be seen on Figs. 5 and 6 for the first Matsubara frequency. Note that the data in the figures is not corrected for lattice spacing artefacts and that is the reason why for $M(p_{4},|\vec{p}|)$ we only report data for small $p=\sqrt{p^{2}_{4}+\vec{p}^{2}}$. Indeed, the bare running quark mass at finite temperature looks like that observed at zero temperature as seen, for example, in Fig. 10 in Oliveira:2018lln . This same Fig. also shows that up to momenta $ap\sim 0.5$ the lattice artefacts corrections are negligible or quite small. For the simulations reported here $ap\sim 0.5$ corresponds to a $p\sim 1$ GeV, see Tab. 1, and that is the reason why we plotted the running mass only up to $p=1.5$ GeV. In this section, our analysis of the running quark mass refers only to its low momentum behaviour. The function $M(p_{4},|\vec{p}|)$ for the full range of momenta is discussed in Sec. 4.

As Figs. 5 and 6 show, the quark propagator looks rather different below and above $T_{c}$. For high momenta, the quark wave function approaches a constant from above for $T<T_{c}$ and from below at $T>T_{c}$. Around the critical temperature the results favour a constant $Z_{c}$ for the full range of momenta accessed in our simulation. Fig 5 also shows that, for sufficiently high momenta, $\omega(p_{4},\vec{p})\cong Z(p_{4},\vec{p})$, up to a constant factor close to unity. At low momenta $\omega(p_{4},\vec{p})$ exceeds $Z(p_{4},\vec{p})$ above the critical temperature, while $Z(p_{4},\vec{p})$ exceeds $\omega(p_{4},\vec{p})$ for temperatures below $T_{c}$. Note also that for temperatures above $T_{c}$ the quark wave function approaches a constant value from below, reproducing the same type of behaviour observed for zero temperature and in agreement with the predictions of perturbation theory.

In order to compare our results with the predictions of the hard thermal loop (HTL) for the quark propagator, in Fig. 7 we show the lattice data together with the one loop HTL prediction for $Z_{c}(p_{4},\vec{p})$ and for $T=305$ MeV, i.e. $T/T_{c}=1.13$, for the first Matsubara frequency. We take the one loop HTL quark propagator from Andersen:1999va , and convert it to Euclidean space – see their Eqs. (5), (6) and (7). The exact definition of the quark mass that should appear in the HTL expressions is difficult to read directly from the lattice data, see Figs. 6 and 8, and we have included in the plot for the HTL quark wave function the corresponding numbers for a number of quark masses. Furthermore, we also show twice the same lattice data rescaled to reproduce the HTL function at $p=3$ GeV for $m_{q}=200$ MeV and for $m_{q}=400$ MeV. As can be observed, at higher momenta, i.e. for $p=(~p^{2}_{4}+\vec{p}^{2})^{1/2}\gtrsim 2$ GeV, the lattice data is in good agreement with HTL. In particular for $m_{q}=400$ MeV the lattice quark wave function and the HTL function are in good agreement for the full range of momenta. Our conclusion being that for $T/T_{c}=1.13$, the lattice data is in good qualitative agreement with HTL for the quark wave function.

Similarly as for the quark wave function, the running quark mass shows different behaviours above and below the deconfinement temperature (see Fig. 6). For temperatures below the critical temperature $M(p_{4},\vec{p})$ is a decreasing function of $p=\sqrt{p^{2}_{4}+\vec{p}^{2}}$, while above the deconfinement phase transition it becomes essentially constant. The observed behaviour of $M$ for $T>T_{c}$ supports the interpretation of quarks as quasiparticles with a constant mass. Note, however, that for the largest temperature $M(p_{4},\vec{p})$ seems to increase with $p$; see also the results reported in Sec. 4. For $T>T_{c}$ the constant quark mass is a function of the temperature and, for the range of $T$ considered here, $M$ decreases when $T$ increases. This can be better viewed in Fig. 8, where $M(p_{4},\vec{p}=0)$ for the first Matsubara frequency is shown for the full range of temperatures considered in the current simulations. The values of $M(p_{4},\vec{p}=0)$ reported in Fig. 8 were computed using the data from the lowest momentum that fulfill the cuts mentioned before and taking the data of the lattice form factors (10), (11), (12) with the errors computed assuming Gaussian error propagation. Note that the values reported in Fig. 8 do not refer exactly to the same $p_{4}$, that ranges from just above $\sim 0.83$ GeV to just below $\sim 1.1$ GeV. It follows that for $T>T_{c}$ typical values for the mass of the quasiparticle are about $\sim 100$ MeV.

Although, our simulations are not at the chiral limit, the data for the running mass clearly shows a strong mass suppression as the temperature crosses $T_{c}$, with the values of $M$ given in Fig. 8 above $T_{c}$ being about half of the values reported when the quarks are in the confined phase.

4 On the Correction of the Lattice Artefacts

As discussed in Skullerud:2000un ; Skullerud:2001aw for the zero temperature case, the lattice form factors computed directly from the inversion of the fermion matrix are contaminated by lattice artefacts. The solutions suggested to remove the lattice artefacts rely on the tree level lattice quark propagator Skullerud:2000un ; Skullerud:2001aw , on an expansion on the invariants of the H4 group Becirevic:1999uc ; deSoto:2007ht or on a combination of both methods Oliveira:2018lln .

In what concerns the quark wave function at finite temperature, i.e. the results reported in Sec. 3 and summarized in Fig. 5, the function $Z_{c}(p_{4},\vec{p})$ looks rather flat for momenta above $\sim 2$ GeV for the temperatures investigated. This is precisely the type of functional behaviour predicted by perturbation theory. However, these findings for $Z_{c}(p_{4},\vec{p})$ contrast with the results of previous simulations at zero temperature, where for sufficiently high momenta the lattice quark wave function is a decreasing function of $p$, even after the partial removing of the lattice artefacts based on the rotated tree level quark propagator Skullerud:2000un ; Skullerud:2001aw ; Oliveira:2018lln . These results follow because $Z_{c}(p_{4},\vec{p})$ is measured as a ratio of functions that, according to the procedure mentioned above, have exactly the same type of lattice artefacts corrections and, therefore, by taking ratios of these type of functions the lattice artefacts cancel exactly or, at least, are strongly suppressed. Given the good agreement between the lattice $Z_{c}(p_{4},\vec{p})$ reported in Fig. 5 and the results of perturbation theory for $p\gtrsim 2$ GeV, we assume that the computed $Z_{c}(p_{4},\vec{p})$ are essentially free of lattice artefacts.

On the other hand, the running quark mass at finite temperature is computed as in the zero temperature case. Not surprisingly, the lattice data for $M(p_{4},\vec{p})$ shows a similar pattern as observed at zero temperature and this function increases at higher momenta as observed in Figs. 1 and 10 of Oliveira:2018lln for the uncorrected lattice data. It follows that the computation of the running quark mass requires an estimation and subtraction of the lattice artefacts.

For the estimation of the lattice artefacts for the running mass we use the procedure outlined in Skullerud:2000un ; Skullerud:2001aw . The space-spin components of the rotated tree-level quark propagator is given by

[TABLE]

Expressions for $A(m,p)$ and $B(m,p)$ can be found in Skullerud:2000un ; Skullerud:2001aw . There is no unambiguous way to subtract the lattice artefacts. Typically, two main definitions of the “continuum” running mass are considered where one naïvely subtracts part of the lattice artefacts, defining $m_{s}(p_{4},\vec{p})$, and the so-called hybrid scheme that defines $m_{h}(p_{4},\vec{p})$, where in the subtraction of the lattice artefacts it is taking into consideration if the corrections give either positive or negative contributions to the pole mass. In the definition of $m_{h}(p_{4},\vec{p})$ the negative contributions are subtracted before rescaling the result to take into account the remaining corrections. Details of the procedure and definitions can be found in the above cited works. To illustrate the running quark mass definition, on Fig. 9 we report results obtained with the rotated tree level quark propagator, where the naïve pole mass grows with momenta and $m_{h}(p_{4},\vec{p})$ is constant for all $ap$.

The hybrid running quark mass for the simulations with a $m_{bare}\approx 50$ MeV and all temperatures are given in Fig. 10 for the first Matsubara frequency. The functional form at small momenta for $m_{h}(p_{4},\vec{p})$ changes dramatically when the deconfinement transition is approached. At the infrared scales $m_{h}$ is a decreasing function of $p$ for $T\lesssim T_{c}$, is essentially constant for $T$ above $T_{c}$ and in the range 290 - 305 MeV and for the highest temperature considered here $m_{h}$ is an increasing function of the momenta. We call the reader’s attention that in the simulation for $T=324$ MeV the number of lattice points in the time direction, used to define the temperature, is shorter $L_{t}=6$ compared with all other cases that use $L_{t}=8$. It is not clear if such a short time direction is responsible for the different observed behaviour on the hybrid running quark mass. At sufficiently high $p$, the data shows that the running masses seem to collapse into a single curve. The data reported on Fig. 10 suggests that at high $p$ and within the range of momenta accessed in our simulation $m_{h}(p_{4},\vec{p})$ is a decreasing function of $p$. We have tried to fit a constant value to the highest momenta and have verified that for $p\gtrsim 10$ GeV the lattice data for $m_{h}(p_{4},\vec{p})$ is compatible with a constant just above 50 MeV; for example, for $T=290$ MeV, we get $m_{h}=51.984\pm 0.030$ MeV with a $\chi^{2}/d.o.f=1.1$ for $p>9.8$ GeV. This fit value for the mass is within the same range of values for $m_{bare}$ considered here – see Tab. 1.

For completeness, in Fig. 11 we show $m_{h}(p_{4},\vec{p})$ for all temperatures and all the Matsubara frequencies. The data shows, once more, a clear violation of rotational symmetry.

The computed quark wave function $Z_{c}(p_{4},\vec{p})$ and running quark mass $m_{h}(p_{4},\vec{p})$ imply changes on the functional behaviour of the quark spectral functions below and above the critical temperature. However, given the small number of Matsubara frequencies accessed in our study, we do not attempt to compute any of the components of the quark spectral function.

5 Summary and Conclusions

In the current work we report on the computation of the finite temperature Landau gauge lattice QCD quark propagator in the quenched approximation at small quark masses $\sim 10$ MeV and $\sim 50$ MeV for temperatures below and above the deconfinement phase transition. The various quark form factors are investigated as a function of the temperature for the various Matsubara frequencies.

Our results show that both the quark wave function and the running quark mass have different functional forms for $T<T_{c}$ and for $T>T_{c}$. Above $T_{c}$ the quark wave function seems to reproduce the predictions of HTL expansion for $T/T_{c}\gtrsim 1.1$ being constant for $p\gtrsim 2$ GeV and slightly suppressed at smaller momenta. The suppression of the quark wave function at infrared scales is also observed for the zero temperature case; see Oliveira:2018lln and references therein. On the other hand, at low momenta and for $T<T_{c}$ the quark wave function is enhanced. These results show that $Z_{c}$ has a non-trivial dependence on the temperature and its suppression or enhancement at low momenta is not enough to identify if the quark is on the confined or deconfined phase.

The running quark mass also shows a non-trivial dependence with the temperature of the heat bath. Indeed, the running quark mass is highly suppressed when $T$ crosses $T_{c}$ from below, with typical values just above $T_{c}$ being about half of the corresponding values just below $T_{c}$. At low momenta $p\lesssim 1.5$ GeV for the temperatures considered here, the running quark mass is nearly constant above $T_{c}$. This can be viewed as favouring the description of quarks as quasi free particles with a constant mass, as is used e.g. by the HTL approach. Our simulations also show that at high momenta the running quark mass reproduces the values of the bare quark mass reported in Tab. 1. For temperatures below $T_{c}$ the running quark mass increase when one approaches the zero momentum limit, reproducing the observed behaviour for the running quark mass at zero temperature.

Our simulations are based on pure gauge configurations, where the pure gauge sector is solved exactly but without taking into account the dynamics of quarks. The computed quark propagators refer to the propagation of the quark fields on a background of gluons that is solved exactly for the bosonic sector of the theory. The observed significant suppression of the running quark mass for temperatures above $T_{c}$ is an indication of a major contribution of the gluons for the mechanism of chiral symmetry breaking. Furthermore, our results suggest an important link between (gluon) confinement and the restoration of chiral symmetry that needs to be further investigated.

Acknowledgements.

The authors acknowledge M. Strickland and P. Costa for helpful discussions. The authors acknowledge financial support from FCT under contract with reference UID/FIS/04564/2016. P.J.S. also acknowledges partial support by Fundação para a Ciência e a Tecnologia (FCT) under contracts SFRH/BPD/40998/2007 and SFRH/BPD/109971/2015. The authors also acknowledge the Laboratory for Advanced Computing at University of Coimbra (http://www.uc.pt/lca) for providing access to the HPC resource Navigator. The SU(3) lattice simulations were done using Chroma Edwards2005 and PFFT Pippig2013 libraries.