Streamlined data analysis in Python

Python is a particularly appealing language to carry out data analysis, owing in part to its user-friendly character as well as its access to well maintained and powerful libraries like NumPy and SciPy. Still, for the purpose of analyzing data in a lattice QCD context, some desirable functionality is missing from these libraries. Moreover, scripting languages tend to be slower than compiled ones. To help address these points we present the LatticeToolbox, a collection of Python modules to facilitate lattice QCD data analysis. Some highlighted features include general-purpose jackknife and bootstrap routines; modules for reading in and storing gauge configurations; a module to carry out hadron resonance gas model calculations; and convenience wrappers for SciPy integration, curve fitting, and splines. These features are sped up behind the scenes using parallelization and just-in-time compilation.Comment: 9 pages, 1 figure, presented at the 40th International Symposium on Lattice Field Theor

Heavy Quark Diffusion from 2+1 Flavor Lattice QCD

We present the first calculations of the heavy flavor diffusion coefficient using lattice QCD with light dynamical quarks. For temperatures $195\,\mathrm{MeV}<T<352\,\mathrm{MeV}$, the heavy quark spatial diffusion coefficient is found to be significantly smaller than previous quenched lattice QCD and recent phenomenological estimates. The result implies very fast hydrodynamization of heavy quarks in the quark-gluon plasma created during ultrarelativistic heavy-ion collision experiments

Continuum extrapolation of the gradient-flowed color-magnetic correlator at 1.5 Tc1.5\,T_c

In a recently published work we employ gradient flow on the lattice to extract the leading contribution of the heavy quark momentum diffusion coefficient in the heavy quark limit from calculations of a well-known two-point function of color-electric field operators. In this article we want to report the progress of calculating the recently derived color-magnetic correlator that encodes a finite mass correction to this transport coefficient. The calculations we present here are based on the same ensemble of quenched gauge configurations at $1.5\, T_c$ that we previously used for the color-electric correlator.Comment: 7 pages, 4 figures, presented at the 38th International Symposium on Lattice Field Theor

HotQCD on Multi-GPU Systems

We present $\texttt{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.Comment: 7 pages, 3 figures, presented at the 38th International Symposium on Lattice Field Theor

Quark Mass Dependence of Heavy Quark Diffusion Coefficient from Lattice QCD

We present the first study of the quark mass dependence of the heavy quark momentum and spatial diffusion coefficients using lattice QCD with light dynamical quarks corresponding to a pion mass of 320 MeV. We find that, for the temperature range 195 MeV $<T<$ 293 MeV, the spatial diffusion coefficients of the charm and bottom quarks are smaller than those obtained in phenomenological models that describe the $p_T$ spectra and elliptic flow of open heavy flavor hadrons.Comment: 15 pages, 12 figure

Lattice QCD noise reduction for bosonic correlators through blocking

We propose a method to substantially improve the signal-to-noise ratio of lattice correlation functions for bosonic operators or other operator combinations with disconnected contributions. The technique is applicable for correlations between operators on two planes (zero momentum correlators) when the dimension of the plane is larger than the separation between the two planes which are correlated. In this case, the correlation arises primarily from points whose in-plane coordinates are close, but noise arises from all pairs of points. By breaking each plane into bins and computing bin-bin correlations, it is possible to capture these short-distance correlators exactly while replacing (small) correlators at large spatial extent with a fit, with smaller uncertainty than the data. The cost is only marginally larger than averaging each plane before correlating, but the improvement in signal-to-noise can be substantial. We test the method on correlators of the gradient-flowed topological charge density and squared field strength, finding noise reductions by a factor of ∼3–7 compared to the conventional approach on the same ensemble of configurations

SIMULATeQCD: A simple multi-GPU lattice code for QCD calculations

The rise of exascale supercomputers has fueled competition among GPU vendors, driving lattice QCD developers to write code that supports multiple APIs. Moreover, new developments in algorithms and physics research require frequent updates to existing software. These challenges have to be balanced against constantly changing personnel. At the same time, there is a wide range of applications for HISQ fermions in QCD studies. This situation encourages the development of software featuring a HISQ action that is flexible, high-performing, open source, easy to use, and easy to adapt. In this technical paper, we explain the design strategy, provide implementation details, list available algorithms and modules, and show key performance indicators for SIMULATeQCD, a simple multi-GPU lattice code for large-scale QCD calculations, mainly developed and used by the HotQCD collaboration. The code is publicly available on GitHub.Comment: 17 pages, 7 figure