We apply the gradient flow on a color-electric two-point function that encodes the heavy quark momentum diffusion coefficient. The simulations are done on fine isotropic lattices in the quenched approximation at 1.5 Tc. The continuum extrapolation is performed at fixed flow time followed by a second extrapolation to zero flow time. Perturbative calculations of this correlation function under Wilson flow are used to enhance the extrapolations of the non-perturbative lattice correlator. The final estimate for the continuum correlator at zero flow time largely agrees with one obtained from a previous study using the multi-level algorithm. We perform a spectral reconstruction based on perturbative model fits to estimate the heavy quark momentum diffusion coefficient. The approach we present here yields high-precision data for the correlator and is also applicable for actions with dynamical fermions.

We present a detailed study of the applications of two stochastic approaches, stochastic optimization method (SOM) and stochastic analytical inference (SAI), to extract spectral functions from Euclidean correlation functions. SOM has the advantage that it does not require prior information. On the other hand, SAI is a more generalized method based on Bayesian inference. Under mean field approximation SAI reduces to the often-used maximum entropy method (MEM), and for a specific choice of the prior SAI becomes equivalent to SOM. To test the applicability of these two stochastic methods to lattice QCD, firstly, we apply these methods to various reasonably chosen model correlation functions, and present detailed comparisons of the reconstructed spectral functions obtained from SOM, SAI and MEM. Next, we present similar studies for charmonia correlation functions obtained from lattice QCD computations using clover-improved Wilson fermions on large, fine, isotropic lattices at 0.75 and 1.5Tc, Tc being the deconfinement transition temperature of a pure gluon plasma. We find that SAI and SOM give consistent results to MEM at these two temperatures.

No abstract

We present SIMULATeQCD, HotQCD's software for performing lattice QCD calculations on GPUs. Started in late 2017 and intended as a full replacement of the previous single GPU lattice QCD code used by the HotQCD collaboration, our software has been developed into an extensive framework for lattice QCD calculations distributed on multiple GPUs over many compute nodes. The code is built on C++, CUDA, and MPI and leverages modern C++ language features to provide high-level data structures, objects, and algorithms that allow users to express lattice QCD calculations in an intuitive way without sacrificing performance. Implemented algorithms range from gradient flow, correlator measurements, and mixed precision conjugate gradient solvers all the way to full HISQ gauge field configuration generation using RHMC. After successful deployment in large-scale computing projects, we want to share the result of our efforts with the lattice QCD community by making it publicly available. In these proceedings, we will present some of the key features of our code, demonstrate its ease of use, and show benchmarks of performance critical kernels on state-of-the-art supercomputers.

In this paper we report our results on quarkonium spectral functions in the vector channel obtained from quenched lattice QCD simulations at T ∈ [0.75, 2.25] T c . The calculations have been performed on very large and fine isotropic lattices where both charm and bottom quarks can be treated relativistically. The spectral functions are reconstructed using the Maximum Entropy Method. We study the dissociation of quarkonium states from the temperature dependence of the spectral functions and estimate heavy quark diffusion coefficients using the low-frequency behavior of the vector spectral functions.