Researcher profile

Olaf Kaczmarek

Olaf Kaczmarek contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
20works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2022arXiv

Spectral reconstruction details of a gradient-flowed color-electric correlator

In a recently published work we provide a proof-of-concept of a novel method to extract the heavy quark momentum diffusion coefficient from color-electric correlators on the lattice using gradient flow. The transport coefficient can be found in the infrared limit of the corresponding spectral function which is reconstructed through perturbative model fits of the correlator data. In this proceedings report we want to give more detailed insights into the systematic uncertainties of this procedure and compare our results with other studies.

preprint2022arXiv

Static quark anti-quark interactions at non-zero temperature from lattice QCD

We study the interactions of a static quark antiquark pair at non-zero temperature using realistic 2+1 flavor lattice QCD calculations. The study consists of two parts: the first investigates the properties of Wilson line correlators in Coulomb gauge and compares to predictions of hard-thermal loop perturbation theory. As a second step we extract the spectral functions underlying the correlators using four conceptually different methods: spectral function fits, a HTL inspired fit for the correlation function, Padé rational approximation and the Bayesian BR spectral reconstruction. We find that our high statistics Euclidean lattice data are amenable to different hypotheses for the shapes of the spectral function and we compare the implications of each analysis method for the existence and properties of a well defined ground state spectral peak.

preprint2022arXiv

Static quark anti-quark interactions at non-zero temperature from lattice QCD

We present results on the in-medium interactions of static quark anti-quark pairs using realistic 2+1 HISQ flavor lattice QCD. Focus is put on the extraction of spectral information from Wilson line correlators in Coulomb gauge using four complementary methods. Our results indicate that on HISQ lattices, the position of the dominant spectral peak associated with the real-part of the interquark potential remains unaffected by temperature. This is in contrast to prior work in quenched QCD and we present follow up comparisons to newly generated quenched ensembles.

preprint2021arXiv

Charm and beauty in the deconfined plasma from quenched lattice QCD

We present continuum extrapolated results of charmonium and bottomonium correlators in the vector channel at several temperatures below and above $T_c$. The continuum extrapolation jointly performed with the interpolations to have physical values of $J/ψ$ and $Υ$ masses in the confined phase is based on calculations on several large quenched isotropic lattices using clover-improved Wilson valence fermions carrying different quark masses. The extrapolated lattice correlators are confronted with perturbation theory results incorporating resummed thermal effects around the threshold from pNRQCD and vacuum asymptotics above the threshold. An additional transport peak is modelled below the threshold allowing for an estimate of the diffusion coefficients for charm and bottom quarks. We find that charmonium correlators in the vector channel can be well reproduced by perturbative spectral functions above $T_c$ where no resonance peaks for $J/ψ$ are needed at and above 1.1 $T_c$, while for bottomonium correlators a resonance peak for $Υ$ is still needed up to 1.5 $T_c$. By analyzing the transport contribution to the correlators we find that the drag coefficient of a charm quark is larger than that of a bottom quark.

preprint2021arXiv

Critical behavior towards the chiral limit at vanishing and non-vanishing chemical potentials

We study the scaling behavior of the (2+1)-flavor QCD crossover region towards the chiral limit with smaller-than-physical light quark mass gauge ensembles, generated using the HISQ fermion discretization. At zero chemical potential, we study the fluctuations of conserved charges and their correlations with the chiral condensate, towards the chiral limit. We analyse the role of universal and regular contributions to the above quantities. We find a preliminary estimate of the leading curvature coefficient of the chiral phase transition line using scaling arguments.

preprint2021arXiv

Heavy quark momentum diffusion from the lattice using gradient flow

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\,T_c$. 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.

preprint2021arXiv

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.

preprint2021arXiv

Sensitivity of the Polyakov loop and related observables to chiral symmetry restoration

While the Polyakov loop is an order parameter of the deconfinement transition in the heavy quark mass regime of QCD, its sensitivity to the deconfinement of light, dynamical quarks in QCD is not apparent. On the other hand, the quark mass dependence of the Polyakov loop is sensitive to the appearance of a chiral phase transition. Using lattice QCD calculations in the staggered fermion discretization scheme at finite values of the lattice spacing, $aT = 1/8$, we show here, for the first time, that the Polyakov loop expectation value, and the heavy quark free energy extracted from it, behave like energy-like observables in the vicinity of the chiral phase transition temperature $T_c$. Consistent with scaling behavior of energy-like observables in the 3-$d$, O(2) universality class, the quark mass derivatives diverge in the chiral limit at $T_c$ while the temperature derivatives stay finite. The latter will develop a characteristic spike at $T_c$. This, however, may be resolved only in calculations with quark masses being two orders of magnitude smaller than those currently accessible in lattice QCD calculations.

preprint2020arXiv

Dirac Eigenvalue spectrum of $N_f$=2+1 QCD towards the chiral limit using HISQ fermions

We utilize the eigenvalue filtering technique combined with the stochastic estimate of the mode number to determine the eigenvalue spectrum. Simulations of (2 + 1)-flavor QCD are performed using the Highly Improved Staggered Quarks (HISQ/tree) action on $N_τ$ = 8 lattices with aspect ratios $N_σ/N_τ$ ranging from 5 to 7. The strange quark mass is fixed to its physical value $m_{s}^{\rm phy}$, and the light quark masses $m_{l}$ are varied from $m_{s}^{\rm phy}/40$ to $m_{s}^{\rm phy}/160$ which correspond to pion mass $m_π$ ranging from 110 MeV to 55 MeV in the continuum limit. We compute the chiral condensate and $χ_π - χ_δ$ through the eigenvalue spectrum obtained from the the eigenvalue filtering method. We compare these results with those obtained from a direct calculation of the observables which involves inversions of the fermion matrix using the stochastic "noise vector" method. We find that these approaches yield consistent results. Furthermore, we also investigate the quark mass and temperature dependences of the Dirac eigenvalue density at zero eigenvalues to gain more insights about the $U_A(1)$ symmetry breaking in QCD.

preprint2020arXiv

Euclidean correlation functions of the topological charge density

We present first results of our study on the Euclidean topological charge density correlation function. In order to get a well defined topological charge density and to improve the signal of the correlation function at large separations we make use of the gradient flow. We investigate the flow-time dependence on fine quenched lattices. The final goal of this study is to perform a continuum extrapolation for the pure SU(3) plasma and to extract the related transport coefficient, the sphaleron rate.

preprint2020arXiv

Open-Boundary Conditions in the Deconfined Phase

In this work, we consider open-boundary conditions at high temperatures, as they can potentially be of help to measure the topological susceptibility. In particular, we measure the extent of the boundary effects at $T=1.5T_c$ and $T=2.7T_c$. In the first case, it is larger than at $T=0$ while we find it to be smaller in the second case. The length of this "boundary zone" is controlled by the screening masses. We use this fact to measure the scalar and pseudo-scalar screening masses at these two temperatures. We observe a mass gap at $T=1.5T_c$ but not at $T=2.7T_c$. Finally, we use our pseudo-scalar channel analysis to estimate the topological susceptibility. The results at $T=1.5T_c$ are in good agreement with the literature. At $T=2.7T_c$, they appear to suffer from topological freezing, impeding us from providing a precise determination of the topological susceptibility. It still provides us with a lower bound, which is already in mild tension with some of the existing results.

preprint2020arXiv

Polyakov Loop Susceptibility and Correlators in the Chiral Limit

In quenched QCD the Polyakov loop is an order parameter of the deconfinement transition, but with decreasing quark mass, the peak in the Polyakov loop susceptibility becomes less pronounced, and it loses its interpretation as an indicator for deconfinement. For this $N_f=2+1$ HISQ study, we fix the strange quark mass $m_s$ at its physical value and investigate the dependence of the Polyakov loop on the light quark mass $m_l$ in the range $m_s/m_l=27-160$, following $m_l$ toward the chiral limit. In particular we will look how the inflection point and susceptibility behave as we decrease $m_l$, to see whether one finds any indication of a crossover, and therefore whether the Polyakov loop is sensitive to the chiral phase transition. Preliminary results show no signal of a crossover from the real part of the Polyakov loop in the vicinity of the chiral crossover. Closely related is an investigation of Polyakov loop correlations and the Debye mass in this limit. Preliminary results suggest little or no dependence on $m_l$.

preprint2020arXiv

QCD phase transition in the chiral limit

We present a lattice QCD based determination of the chiral phase transition temperature in QCD with two massless (up and down) and one strange quark having its physical mass. We propose and calculate two novel estimators for the chiral transition temperature for several values of the light quark masses, corresponding to Goldstone pion masses in the range of $58~{\rm MeV}\lesssim m_π\lesssim 163~{\rm MeV}$. The chiral phase transition temperature is determined by extrapolating to vanishing pion mass using universal scaling analysis. After thermodynamic, continuum and chiral extrapolations we find the chiral phase transition temperature $T_c^0=132^{+3}_{-6}$ MeV. We also present some preliminary calculations on interplay of effective $U_A(1)$ restoration and chiral phase transition towards chiral limit.

preprint2020arXiv

Thermal modifications of quarkonia and heavy quark diffusion from a comparison of continuum-extrapolated lattice results to perturbative QCD

We investigate the in-medium modifications of heavy quarkonia in the vector channel and the heavy quark diffusion coefficient by comparing Euclidean correlators from the lattice to a perturbative spectral function. On the lattice side, we work with continuum extrapolated data from four different large and fine lattices with Clover-improved Wilson fermions in the quenched approximation at five temperatures (0.75, 1.1, 1.3, 1.5 and 2.25$T_c$). On the perturbative side, we use a combination of pNRQCD and vacuum asymptotics to describe the spectral function. After accounting for systematic errors, we obtain a spectral function that is suited to describe the bound state region. This spectral function describes charmonium well without a resonance peak at any of our analyzed temperatures above $T_c$, while we observe a thermally broadened resonance peak for bottomonium that is only melted at our largest temperature, $2.25T_c$. For the transport contribution we assume a Breit-Wigner shaped peak and find that the drag coefficient of charm quarks is larger than that of bottom quarks.

preprint2019arXiv

Chiral phase transition in (2 + 1)-flavor QCD

The chiral phase transition temperature $T_{c}^{0}$ is a fundamental quantity of QCD. To determine this quantity we have performed simulations of (2 + 1)-flavor QCD using the Highly Improved Staggered Quarks (HISQ/tree) action on $N_τ=6, 8$ and 12 lattices with aspect ratios $N_σ/N_τ$ ranging from 4 to 8. In our simulations the strange quark mass is fixed to its physical value $m_{s}^{\rm{phy}}$, and the values of two degenerate light quark masses $m_{l}$ are varied from $m_{s}^{\rm{phy}}/20$ to $m_{s}^{\rm{phy}}/160$ which correspond to a Goldstone pion mass $m_π$ ranging from 160 MeV to 55 MeV in the continuum limit. By investigating the light quark mass dependence and the volume dependence of various chiral observables, e.g. chiral susceptibilities and Binder cumulants, no evidence for a first order phase transition in our current quark mass window is found. Two estimators $T_{60}$ and $T_δ$ are proposed to extract the chiral phase transition temperature $T_{c}^{0}$ in the chiral and continuum limit and our current estimate for $T_{c}^{0}$ is $132_{-6}^{+3}$ MeV.

preprint2019arXiv

Color-electric correlation functions under gradient flow

We report on the progress of our study on the color-electric correlation functions under gradient flow on the lattice. This calculation is the first step of our long-term project to estimate a series of important transport coefficients, of which the heavy quark momentum diffusion coefficient is our first attempt, as it can be extracted from a color-electric correlation function that has been calculated non-perturbatively using noise reduction technique in the quenched approximation [1]. By comparing the flowed correlation function with those obtained by other signal-improving techniques, for instance the multi-level algorithm used in [1], one can gain insight into the applicability of the gradient flow and the renormalization of the correlation functions. We start with quenched, isotropic lattices. Currently we have finished measuring the color-electric correlation functions on 4 different lattices with $β$-values corresponding to a temperature of $T\approx 1.5T_c$. We perform a continuum extrapolation on the flowed correlators at fixed physical flow times followed by an extrapolation of the continuum estimate back to zero flow time. The next step is to extend the study to different temperatures and to extend the study to dynamical QCD.

preprint2019arXiv

Meson Screening Masses in (2+1)-Flavor QCD

We present lattice QCD results for mesonic screening masses in the temperature range 140 MeV $\lesssim T \lesssim$ 2500 MeV. Our calculations were carried out using (2+1)-flavors of the Highly Improved Staggered Quark (HISQ) action, with a physical value for the strange quark mass and two values of the light quark mass corresponding to pion masses of 160 MeV and 140 MeV. Continuum-extrapolated results were obtained using calculations with a variety of lattice spacings corresponding to temporal lattice extents $N_τ= 6 - 16$. We discuss the implications of these results for the effective restoration of various symmetries in the high temperature phase of QCD, as well as the approach toward the perturbative limit.

preprint2018arXiv

Continuum extrapolation of quarkonium correlators at non-zero temperature

In the investigation of in-medium modifications of quarkonia and for determining heavy quark diffusion coefficients, correlation functions play a crucial role. For the first time we perform a continuum extrapolation of charmonium and bottomonium correlators in the vector channel based on non-perturbatively clover-improved Wilson fermions in quenched lattice QCD. Calculations were done at 4 different lattice spacings, spatial extents between 96 and 192, aspect ratio from 1/6 to 1/2, for 5 temperatures between $T/T_c = 0.75$ and $T/T_c = 2.2$. We interpolate between different quark masses to match to the same vector meson mass over different lattice setups. Afterwards we extrapolate the renormalized correlators to the continuum. While we find a strong temperature dependence for charmonium, bottomonium states are only slightly affected at higher temperatures.

preprint2017arXiv

Thermal modifications of charmonia and bottomonia from spatial correlation functions

We present our study on the thermal modifications of charmonia and bottomonia from spatial correlation functions at zero and nonzero momenta in quenched QCD. To accommodate the heavy quarks on the lattice we performed simulations on very fine lattices at a fixed beta value corresponding to a lattice spacing $a^{-1}=22.8$ GeV on $192^3\times32$, $192^3\times48$, $192^3\times56$, $192^3\times64$ and $192^3\times96$ lattices using clover-improved Wilson fermions. These lattices correspond to temperatures of $2.25 T_c$, $1.50 T_c$, $1.25 T_c$, $1.10 T_c$ and $0.75 T_c$. To increase the signal to noise ratio in the axial-vector and scalar channels we used multi-sources for the measurement of spatial correlation functions. By investigating on the differences between spatial and temporal correlators as well as the temperature dependence of screening masses we will discuss the thermal effects in different channels of quarkonium states. Besides this the dispersion relation of the screening mass at different momenta is also discussed.

preprint2017arXiv

Thermal Simulations, Open Boundary Conditions and Switches

$SU(N)$ gauge theories on compact spaces have a non-trivial vacuum structure characterized by a countable set of topological sectors and their topological charge. In lattice simulations, every topological sector needs to be explored a number of times which reflects its weight in the path integral. Current lattice simulations are impeded by the so-called freezing of the topological charge problem. As the continuum is approached, energy barriers between topological sectors become well defined and the simulations get trapped in a given sector. A possible way out was introduced by Lüscher and Schaefer using open boundary condition in the time extent. However, this solution cannot be used for thermal simulations, where the time direction is required to be periodic. In this proceedings, we present results obtained using open boundary conditions in space, at non-zero temperature. With these conditions, the topological charge is not quantized and the topological barriers are lifted. A downside of this method are the strong finite-size effects introduced by the boundary conditions. We also present some exploratory results which show how these conditions could be used on an algorithmic level to reshuffle the system and generate periodic configurations with non-zero topological charge.