# The Bayesian reconstruction of the in-medium heavy quark potential from lattice QCD and its stability

###### Abstract

We report recent results of a non-perturbative determination of the static heavy-quark potential in quenched and dynamical lattice QCD at finite temperature. The real and imaginary part of this complex quantity are extracted from the spectral function of Wilson line correlators in Coulomb gauge. To obtain spectral information from Euclidean time numerical data, our study relies on a novel Bayesian prescription that differs from the Maximum Entropy Method. We perform simulations on quenched lattices with , which cover . To investigate the potential in a quark-gluon plasma with light u,d and s quarks we utilize ASQTAD lattices with by the HotQCD collaboration, giving access to temperatures between . The real part of the potential exhibits a clean transition from a linear, confining behavior in the hadronic phase to a Debye screened form above deconfinement. Interestingly its values lie close to the color singlet free energies in Coulomb gauge at all temperatures. We estimate the imaginary part on quenched lattices and find that it is of the same order of magnitude as in hard-thermal loop perturbation theory. From among all the systematic checks carried out in our study, we discuss explicitly the dependence of the result on the default model and the number of datapoints.

address=Institut de Théorie des Phénomènes Physiques, Ecole Polytechnique Fédérale de Lausanne, CH-1015, Lausanne, Switzerland

address=Fakultät für Physik, Universität Bielefeld, D-33615 Bielefeld, Germany

address=Institute for Theoretical Physics, Heidelberg University, Philosophenweg 16, D-69120 Heidelberg, Germany

The Bayesian reconstruction of the in-medium heavy quark potential from lattice QCD and its stability

^{}^{,}

missing

Keywords: Heavy Quarkonium, Static Potential, Quark-Gluon Plasma, Lattice QCD

PACS: 12.38.Gc,03.65.Ge,12.38.Mh,14.40.Pq

## Introduction

Nuclear matter under extreme conditions provides ample challenges to our understanding of the mechanism of confinement based on microscopic quantum-chromo dynamics. Strong fluctuations, be it at the high temperatures that existed shortly after the big bang or at the high densities in the interior of neutron stars, lead to an emergence of partonic degrees of freedom, confined in hadrons at everyday energy scales. With the advent of modern collider facilities, such as RHIC and LHC, as well as future facilities, such as FAIR, it is possible to recreate these scenarios directly in the laboratory.

To gain quantitative insight into the dynamical aspects of deconfinement, the bound states of a heavy quark and antiquark are ideal probes Brambilla:2010cs (). The ground states of these two-body systems in vacuum possess lifetimes of the order of inverse keV Agashe:2014kda (), which makes them essentially stable from the point of view of QCD. If immersed in a thermal bath they are thus able to accumulate effects from the medium over time, while still remaining well defined particles that can be cleanly measured in experiment. Eventually, as temperature is increased, they will melt into their constituents. How and over which timescales this melting occurs, has been and remains the focus of a wealth of theoretical studies (for a selection of lattice QCD based studies to in-medium modification of heavy quarkonium see Umeda:2002vr (); Asakawa:2003re ()).

A first principles understanding of the melting process in turn opens up the possibility to extract from the measured yields in relativistic heavy-ion collisions Adare:2006ns () the properties of the medium created in the collision center. This essentially would turn heavy quarkonium into a dynamical thermometer for RHIC and LHC. While the classic idea of a suppression of yields in the presence of a QGP Matsui:1986dk () invigorated this area of research, we now know that an intricate interplay between melting and regeneration Liu:2009nb (), needs to be understood before quantitative predictions can be made.

What distinguishes heavy quarkonium, such as e.g. bottomonium (the vector channel), from lighter mesons is the inherent separation of scales between the constituent quark mass and the medium temperature, as well as the characteristic scale of QCD (). From the rough estimates for the initial temperature based on direct photons at RHIC (MeV) Adare:2008ab () and LHC (MeV) Wilde:2012wc () we see that indeed . Intuitively this tells us that under the conditions of current heavy-ion experiments, neither quantum nor thermal fluctuations can spontaneously create a heavy quark pair, so that a description in terms of non-relativistic physics should be feasible. Note that these temperatures are still close to the deconfinement temperature MeV as determined from lattice QCD Aoki:2006we (); Bazavov:2011nk () and hence the medium partons remain strongly interacting.

Using effective field theories (EFT) Brambilla:2004jw () it has been shown that the effects of the interactions between a static quark-antiquark pair and a thermal medium can be summarized in a Schrödinger equation with an effective in-medium potential. The derivation of the potential relies on a matching procedure to fundamental QCD, and relates its values to the late real-time behavior of the rectangular Wilson loop

(1) |

A first evaluation of this quantity in hard-thermal loop perturbation theory Laine:2007qy () showed that the potential at finite temperature is complex valued. As expected at temperatures far above , the presence of deconfined partons leads to Debye screening of the real part. The imaginary part on the other hand has been attributed to Landau damping (scattering with medium gluons) Beraudo:2007ky () and a subsequent systematic EFT analysis revealed Brambilla:2008cx () additional contributions from the transition of a singlet to an octet state (absorption of medium gluons)^{1}^{1}1The presence of the imaginary part has been given a dynamical interpretation in the language of open-quantum systems Akamatsu:2011se (), where it is related to the remnants of thermal medium fluctuations, which were integrated out when establishing the effective potential description.

Before the advent of a modern effective field theory for in-medium heavy quarkonium, the temperature dependence of the real-part of the static inter-quark potential has often been modeled in phenomenological studies. Some studies used the color singlet free energies others the internal energies or even linear combinations thereof Nadkarni:1986as (). Interestingly, resummed HTL perturbation theory to leading order tells us that and agree but this relation already appears to break down at next to leading order Burnier:2009bk (). To provide phenomenology with dependable input, the question of how to evaluate Eq.(1) in a non-perturbative manner becomes a pressing issue.

Simulations of lattice QCD provide us with non-perturbative access to the static properties of QCD at finite temperature, e.g. the equation of state Borsanyi:2013bia (). These calculations are carried out in Euclidean time, so that real-time quantities, such as the Wilson loop , are not directly accessible. This challenge has spurred remarkable progress over the last two years, both in terms of concepts and methods, moving us towards a reliable evaluation of Eq.(1) on the lattice. One of the central ingredients, acting as bridge between the Euclidean and Minkowski domain, is the following spectral decomposition Rothkopf:2009pk (); Rothkopf:2011db ()

Substituted into the definition Eq.(1), it relates the values of the potential to the real and positive definite spectral function

(2) |

We face two challenges when trying to obtain the values of the potential in practice. One is related to the fact that extracting the continuous function from a finite and noisy set of lattice datapoints is an inherently ill-defined problem. Indeed a simple fit would lead to an infinite set of degenerate spectra that all reproduce the datapoints within their error bars. In order to give meaning to such a task, we turn to the concept of Bayesian inference, an approach well established in mathematical statistics. It provides the formal framework to incorporate additional, so called prior information, such as the positivity of or smoothness properties, in order to regularize an otherwise under-determined fit. For our study we deploy a recently developed Bayesian prescription to spectral function reconstruction Burnier:2013nla (), which has been shown in mock data studies to give better accuracy and precision than both the conventional Maximum Entropy Method Asakawa:2000tr (), as well as the implementation with extended search space Rothkopf:2011ef ().

Carrying out the late time limit in Eq.(2) represents the second challenge. Simply inserting the reconstructed spectrum is impractical, as the numerical limit requires a precision of currently not obtainable. Instead Rothkopf:2011db () proposed that, since the dynamics of the Wilson loop at late times will eventually be dominated by the lowest lying peak in , its functional form should be fitted from the lattice spectra. The Fourier transform of Eq.(2) can subsequently be carried out analytically without difficulty. The question of which functional form to use in such a fit however remained and was only answered later in Burnier:2012az (), where the effects of different time-scales contributing to were disentangled. Based on the symmetries of the Wilson loop it was shown in a general argument that the lowest lying peak will take the form of a skewed Lorentzian:

(3) |

The fitted position thus corresponds to the real part of the potential, the fitted width to the imaginary part. Taking the skewness of the Lorentzian peak into account is crucial, as otherwise e.g. the peak position is systematically overestimated.

The viability of this strategy, sketched in Fig.1, to access the static potential from Euclidean correlator data has been successfully tested in Burnier:2013fca () using HTL resummed perturbation theory. At the same time indications were found that for the sole purpose of determining the potential, one can replace the Wilson loop with the correlator of Wilson lines in Coulomb gauge. The latter is devoid of the cusp divergences of the Wilson loop (see Berwein:2012mw () and references therein) and hence is much easier to measure on the lattice.

These conceptual and technical improvements combined, allow us to reliably extract the temperature dependence of the static inter-quark potential from quenched and dynamical lattice QCD. Besides presenting our main results Burnier:2014ssa () we show explicitly the default model dependence of the reconstructions. (For plots of the effects of a removal of datapoints along see Burnier:2014tbd ()).

## Quenched Lattice QCD

SU(3): | 24 | 32 | 40 | 48 | 56 | 64 | 72 | 80 | 96 |
---|---|---|---|---|---|---|---|---|---|

[MeV] | 839 | 629 | 503 | 419 | 360 | 315 | 280 | 252 | 210 |

3270 | 2030 | 1940 | 1110 | 1410 | 1520 | 860 | 1190 | 1800 |

The first part of our study investigates the real and imaginary part of the static potential in a purely gluonic medium. The reason is practicability, since the relatively low cost of quenched lattice QCD simulations allows us to use lattices with a temporal extend larger than . This is a prerequisite for the reliable determination of spectral peak positions (Re[V]), as well as for a robust order of magnitude estimate of the spectral width, i.e. Im[V], in the presence of our fine lattice spacing.

Based on the naive anisotropic Wilson action with bare parameters and Asakawa:2003re () we performed Monte Carlo simulations on lattices. In this fixed scale approach temperature is varied between by changing the temporal extend between , as tabulated in Tab.1. The physical lattice spacing of fm allows us to resolve the Coulombic region below fm and gives access to distances where the potential is already screened at MeV. Each configuration is fixed to Coulomb gauge iteratively before the Wilson line correlators are measured along each spatial axis, on the square- as well as on the cubic diagonals. We correct for some of the lattice spacing artifacts in the spatial distances by comparing the behavior of free propagators on the lattice to that in the continuum, as described in Necco:2001xg ().

With the Euclidean Wilson lines at hand, we carry out the spectral reconstruction using all datapoints except those at , to avoid possible overlap divergences Berwein:2012mw (). We discretize the frequency interval between GeV by points. of those are used to resolve the lowest lying peak in more detail. A flat default model is selected, not to introduce additional bias in the functional form of the end result. Our numerical implementation uses the LBFGS minimizer with high precision arithmetic ( bits) for the task of finding the unique Bayesian spectrum related to the measured data. The solution is accepted if a relative step size of is reached.

The values of the real and imaginary part shown in Fig.2 are obtained from a subsequent fit of the full width at half maximum of the reconstructed lowest lying peak with the skewed Lorentzian form of Eq.3. and the ’s are treated as fit parameters. The error bars are obtained from the Jackknife variance between ten reconstructions, each excluding a consecutive block of 10% of the measurements .

We find that the real part of the potential weakens gradually. I.e. the confining linear rise present at low temperatures transitions smoothly to a Debye screened behavior in the deconfined phase. Even though we excluded the datapoint in the reconstruction of the potential, Fig.2 shows that the values of Re[V] (open colored symbols) lie close to the color singlet free energies in Coulomb gauge (gray filled circles). Note that perturbation theory indicates Burnier:2009bk () that at intermediate distances Re[V] will lie slightly below . The spectral determination of Re[V] at temperatures around the phase transition shows much less variation compared to the free energies, since it does not only rely on information at but takes into account multiple datapoints along . From the fit of the spectral width we find that Im[V] takes on values that are of the same order of magnitude as those predicted by HTL resummed perturbation theory. This first principles result for the imaginary part is consistent with a recent study Bazavov:2014kva () which fitted the Wilson line correlator using HTL spectral functions.

Simply quoting the statistical errors on the potential is not sufficient to capture the uncertainties, which enter from the underlying spectral reconstruction. Hence our study attempts to quantify also the systematic effects, which are given as error bands in Fig.2. They are estimated using the maximum variation arising from three different changes. First is a reduction of the number of datapoints along by four or eight, second the change of the default model normalization (, ) or functional form () and last the removal of 10%, 20% or 30% of the statistics of the input data. Let us discuss the first two factors in more detail.

In Fig.3 we show the dependence of the potential reconstruction on changes to the default model at the lowest (MeV, ) and highest (MeV, ) temperature. At we find that, due to the abundance of available datapoints and the relatively high statistics, both Re[V] and Im[V] are very robust and do not show variation beyond the small statistical error bars. In the case of much fewer datapoints (), the real part still remains stable, while the spectral width shows a dependence to the form of for fm. Note however that this reduction in reconstruction reliability is accompanied with an increase in the Jackknife variance, so that these error bars overlap for all tested scenarios.

The explicit form of the reconstructed spectra at small and intermediate spatial separation distances and their variation with is shown in Fig.4. As the form of the lowest lying peak structure constitutes the basis for the potential extraction, the behavior of Re[V] and Im[V] observed in Fig.3 is reflected in these plots. The influence of on the position and width of the peak at fm is insignificant at both temperatures. The larger number of available datapoints for keeps the reconstruction of the lowest peak stable even at 5fm, while for , variation in the peak position and width can already be identified by eye.

We have also checked what happens to the result of the potential extraction if the number of datapoints is diminished. While this systematic effect dominates the error bands at high temperatures, it is reassuring that at , taking into account only the first points still leads to the same result for Re[V] within statistical errors. I.e. simply removing the late data does not mimic the effects an actual increase in temperature (to , MeV) would have. In the case of on the other hand, reducing the number of datapoints, produces significant changes in the values we obtain for Re[V]. The decrease in reliability of the reconstruction however is captured well in the simultaneous increase of the Jackknife variance, so that the best result still lies within the error bars. These observations give us confidence that the changes in Re[V] can indeed be attributed to the effects of a thermal medium and not simply to a degradation of the spectral reconstruction at higher T.

Im[V] reacts to the removal of datapoints even at . Its values appears to decrease, as can be expected. The width of the lowest lying peak dominates the correlator exactly at the late values we discard in this example. At intermediate imaginary time on the other hand its contribution is exponentially suppressed and is not easily extracted. Note that we are still able to obtain a finite results for Im[V] at if the correlator values close to are retained.

## Dynamical Lattice QCD

QGP: | 6.8 | 6.9 | 7 | 7.125 | 7.25 | 7.3 | 7.48 |
---|---|---|---|---|---|---|---|

[MeV] | 148 | 164 | 182 | 205 | 232 | 243 | 286 |

a [fm] | 0.111 | 0.1 | 0.09 | 0.08 | 0.071 | 0.068 | 0.057 |

1295 | 1340 | 1015 | 840 | 1220 | 1150 | 1130 |

The second part of our study is aimed at the extraction of the in-medium potential in a true quark-gluon plasma with light u, d and s quarks. Taking into account the effect of dynamical quarks is essential for phenomenological applications, as it e.g. shifts the (pseudo-) critical temperature by more than MeV. To this end we repeat the potential extraction on lattices generated by the HotQCD collaboration Bazavov:2011nk () for the study of the equation of state. Their isotropic configurations are based on the ASQTAD action () and span the temperature range . A fixed box approach is used, so that by varying the lattice spacing , at each T we have at our disposal the same number of points in temporal direction (see Tab.2). The significantly higher computational cost for dynamical quarks currently limits us to the use of , which unfortunately does not allow a reliable determination of spectral widths even with the high precision of the available data.

Since the number of datapoints does not change between temperatures, we carry out the spectral reconstruction with a common choice for and discretize the frequencies between in these units. From the overall points we use to resolve the lowest lying peak. Fitting this structure at the full-width at half maximum with the functional form of Eq.(3) allows us to extract Re[V] shown in Fig.5.

As already seen in the quenched case, the determination of Re[V] based on a spectral reconstruction is very robust at . This allows us to reach distances up to fm at which the signal for the free energies is already lost in statistical noise. The real part does not show any signs of string breaking at the separations we investigated, probably due to the still relatively large pion masses of MeV on our lattices. Nevertheless in the presence of dynamical quarks, the Debye screened behavior in Re[V] is already well pronounced at 286MeV. Interestingly the values of the real part again lie close to those of the color singlet free energies in Coulomb gauge.

Statistical uncertainty is represented by the error bars in Fig.5, which are obtained from the Jackknife variance between ten different reconstructions, from which consecutive blocks of 10% of the correlator measurements had been removed. The error bands contain also systematics and represent the maximum variance from either changing the number of available correlator datapoints by one or two, varying the normalization (, ) or the functional form () of the default model or decreasing the statistics of the input data by 10%, 20% or 30%.

One of the systematic uncertainties we need to understand is the dependence of the result on the default model. The two panels on the right of Fig.5 show explicitly how the extracted values of Re[V] change with . While at low temperature (MeV) the real part is virtually unaffected at all distances probed, we find that at MeV the reconstruction becomes more and more unreliable beyond fm. The Jackknife variance however captures this degradation of the extraction, in that the size of the error bars grows accordingly.

The influence of different choices of on the spectral reconstruction underlying the potential extraction is shown Fig.6 at the lowest (MeV, ) and highest (MeV, ) temperature. Even with the little number of datapoints available, the reconstructions at the small values of fm show a stable position of the lowest lying peak, with some variation at higher frequencies visible. While a separation distance of fm corresponds to only at it already amounts to at . At the same time a smaller lattice spacing leads to a higher UV cutoff, so that the reconstruction has to capture more structure overall. We thus expect and find in Fig.6 that for the larger separation distances, the lowest lying peak at is less reliably reconstructed than at .

In the case of dynamical lattices with a small number of datapoints to begin with, it is also important to quantify the effects of removing further points along . Note that discarding a single temporal step at and amounts to a different change in the actual time extend due to the difference in physical lattice spacing. With the relatively high statistics available, we find that the extraction of Re[V] is stable against small changes in . This again gives us confidence to attribute the change, from the confining behavior at low T to a Debye screened form at , to the effects of the thermal medium and not to a degradation of the spectral reconstruction.

## Conclusion

With the help of recent conceptual and technical progress, a reliable determination of the static inter-quark potential from first principles thermal lattice QCD has become possible. A major role is played by a novel Bayesian approach to spectral function reconstruction, which for the first time allows us to faithfully reconstruct the sharp peaked features of the Wilson line spectra from Euclidean time Monte Carlo data. Fitting the position and width of these spectra, we determined the real and imaginary part in quenched QCD and Re[V] on realistic dynamical lattices with light u,d and s quarks. We find a gradual transition of the linear confining behavior to a Debye screened form in Re[V] and observe that its values lie close to those of the color singlet free energies in Coulomb gauge. Our estimate for the imaginary part on quenched lattices takes on the same order of magnitude as the predictions from HTL resummed perturbation theory. We furthermore checked the systematic dependencies of the Bayesian reconstruction including various default models. The spectra at are very stable, while the decrease in physical extend or in the number of points in Euclidean time at higher temperatures degrades the outcome at separation distances fm, as expected.

## Acknowledgments

The authors thank H. B. Meyer, M. P. Lombardo, P. Petreczky and J.-I. Skullerud for fruitful discussions. YB is supported by SNF grant PZ00P2-142524.

## References

- (1) N. Brambilla et al., Eur. Phys. J. C 71, 1534 (2011).
- (2) K. A. Olive et al. [Particle Data Group Collaboration], Chin. Phys. C 38, 090001 (2014).
- (3) T. Umeda, K. Nomura and H. Matsufuru, Eur. Phys. J. C 39S1, 9 (2005); S. Datta, F. Karsch, P. Petreczky and I. Wetzorke, Phys. Rev. D 69, 094507 (2004); A. Jakovac, P. Petreczky, K. Petrov and A. Velytsky, Phys. Rev. D 75, 014506 (2007); G. Aarts, C. Allton, M. B. Oktay, M. Peardon and J. -I. Skullerud, Phys. Rev. D 76, 094513 (2007); H. Ohno et al., Phys. Rev. D 84, 094504 (2011); G. Aarts et al., JHEP 1111, 103 (2011); H. T. Ding, A. Francis, O. Kaczmarek, F. Karsch, H. Satz and W. Soeldner, Phys. Rev. D 86, 014509 (2012); G. Aarts, C. Allton, S. Kim, M. P. Lombardo, S. M. Ryan and J. -I. Skullerud, arXiv:1310.5467 [hep-lat]; G. Aarts, et al., JHEP 1407, 097 (2014); S. Kim, P. Petreczky and A. Rothkopf, arXiv:1409.3630.
- (4) M. Asakawa and T. Hatsuda, Phys. Rev. Lett. 92, 012001 (2004) [hep-lat/0308034].
- (5) A. Adare et al. [PHENIX Collaboration], Phys. Rev. Lett. 98 (2007) 232301; Z. Tang [STAR Collaboration], J. Phys. G 38 (2011) 124107; S. Chatrchyan et al. [CMS Collaboration], JHEP 1205 (2012) 063, Phys. Rev. Lett. 107 (2011) 052302; B. Abelev et al. [ALICE Collaboration], Phys. Rev. Lett. 109 (2012) 072301.
- (6) T. Matsui and H. Satz, Phys. Lett. B 178 (1986) 416; F. Karsch, M. T. Mehr and H. Satz, “Color Screening and Deconfinement for Bound States of Heavy Quarks,” Z. Phys. C 37 (1988) 617. H. Satz, “Quarkonium Binding and Dissociation: The Spectral Analysis of the QGP,” Nucl. Phys. A 783 (2007) 249 [arXiv:hep-ph/0609197].
- (7) Y. p. Liu, Z. Qu, N. Xu and P. f. Zhuang, Phys. Lett. B 678, 72 (2009) [arXiv:0901.2757 [nucl-th]]. R. Rapp, D. Blaschke and P. Crochet, “Charmonium and bottomonium production in heavy-ion collisions,” Prog. Part. Nucl. Phys. 65 (2010) 209 [arXiv:0807.2470 [hep-ph]].
- (8) A. Adare et al. [PHENIX Collaboration], Phys. Rev. Lett. 104, 132301 (2010) [arXiv:0804.4168 [nucl-ex]].
- (9) M. Wilde [ALICE Collaboration], Nucl. Phys. A 904-905, 573c (2013) [arXiv:1210.5958 [hep-ex]].
- (10) Y. Aoki, G. Endrodi, Z. Fodor, S. D. Katz and K. K. Szabo, Nature 443, 675 (2006) [hep-lat/0611014]; S. Borsanyi et al. [Wuppertal-Budapest Collaboration], JHEP 1009, 073 (2010) [arXiv:1005.3508 [hep-lat]];
- (11) A. Bazavov, T. Bhattacharya, M. Cheng, C. DeTar, H. T. Ding, S. Gottlieb, R. Gupta and P. Hegde et al., Phys. Rev. D 85, 054503 (2012) [arXiv:1111.1710 [hep-lat]].
- (12) N. Brambilla, A. Pineda, J. Soto and A. Vairo, Rev. Mod. Phys. 77, 1423 (2005).
- (13) M. Laine, O. Philipsen and M. Tassler, JHEP 0709, 066 (2007);
- (14) A. Beraudo, J. -P. Blaizot, C. Ratti, Nucl. Phys. A 806, 312 (2008);
- (15) N. Brambilla, J. Ghiglieri, A. Vairo, P. Petreczky, Phys. Rev. D 78, 014017 (2008).
- (16) Y. Akamatsu and A. Rothkopf, Phys. Rev. D 85, 105011 (2012); Y. Akamatsu, Phys. Rev. D 87, no. 4, 045016 (2013) [arXiv:1209.5068 [hep-ph]]; A. Rothkopf, JHEP 1404, 085 (2014) [arXiv:1312.3246 [hep-ph]]; Y. Akamatsu, arXiv:1403.5783 [hep-ph].
- (17) S. Nadkarni, Phys. Rev. D 34, 3904 (1986); C. Y. Wong, Phys. Rev. C 72, 034906 (2005); O. Kaczmarek, PoS CPOD 07, 043 (2007); H. Satz, J. Phys. G 36, 064011 (2009).
- (18) Y. Burnier, M. Laine and M. Vepsalainen, JHEP 1001, 054 (2010).
- (19) S. Borsanyi, Z. Fodor, C. Hoelbling, S. D. Katz, S. Krieg and K. K. Szabo, Phys. Lett. B 730, 99 (2014) [arXiv:1309.5258 [hep-lat]]; A. Bazavov et al. [HotQCD Collaboration], arXiv:1407.6387 [hep-lat].
- (20) A. Rothkopf, T. Hatsuda and S. Sasaki, PoS LAT 2009, 162 (2009);
- (21) A. Rothkopf, T. Hatsuda and S. Sasaki, Phys. Rev. Lett. 108, 162001 (2012)
- (22) Y. Burnier and A. Rothkopf, Phys. Rev. Lett. 111, 182003 (2013); ibid., PoS LATTICE 2013 (2013) 490.
- (23) M. Asakawa, T. Hatsuda and Y. Nakahara, Prog. Part. Nucl. Phys. 46 (2001) 459; ibid., Nucl. Phys. A 715, 863 (2003) [Nucl. Phys. Proc. Suppl. 119, 481 (2003)] [hep-lat/0208059];
- (24) A. Rothkopf, J. Comput. Phys. 238 (2013) 106.
- (25) Y. Burnier and A. Rothkopf, Phys. Rev. D 86 (2012) 051503.
- (26) Y. Burnier and A. Rothkopf, Phys. Rev. D 87 (2013) 11, 114019; ibid, PoS LATTICE 2013 (2013) 491.
- (27) M. Berwein, N. Brambilla, J. Ghiglieri and A. Vairo, JHEP 1303 (2013) 069.
- (28) A. Rothkopf, Mod. Phys. Lett. A 28, 1330005 (2013) [arXiv:1302.6195 [hep-ph]].
- (29) Y. Burnier, O. Kaczmarek and A. Rothkopf, arXiv:1410.2546 [hep-lat].
- (30) Y. Burnier, O. Kaczmarek and A. Rothkopf, arXiv:1410.7311 [hep-lat].
- (31) M. Luscher and P. Weisz, Nucl. Phys. B 445 (1995) 429; S. Necco and R. Sommer, Nucl. Phys. B 622 (2002) 328.
- (32) A. Bazavov, Y. Burnier and P. Petreczky, arXiv:1404.4267 [hep-lat].