Source author record

O. Kaczmarek

O. Kaczmarek appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

44works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

44 published item(s)

preprint2022arXiv

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 $\sim$ 3$-$7 compared to the conventional approach on the same ensemble of configurations.

preprint2022arXiv

Taylor expansions and Padé approximants for cumulants of conserved charge fluctuations at non-vanishing chemical potentials

Using high statistics datasets generated in (2+1)-flavor QCD calculations at finite temperature we present results for low order cumulants of net baryon-number fluctuations at non-zero values of the baryon chemical potential. We calculate Taylor expansions for the pressure (zeroth order cumulant), net baryon-number density (first order cumulant) and the variance of the distribution on net-baryon number fluctuations (second order cumulant). We obtain series expansions from an eighth order expansion of the pressure and compare these to diagonal Padé approximants. This allows us to estimate the range of values for the baryon chemical potential in which these expansions are reliable. We find $μ_B/T\le 2.5$, $2.0$ and $1.5$ for the zeroth, first and second order cumulants, respectively. We furthermore, construct estimators for the radius of convergence of the Taylor series of the pressure. In the vicinity of the pseudo-critical temperature, $T_{pc}\simeq 156.5$ MeV, we find $μ_B/T \gtrsim\ 2.9$ at vanishing strangeness chemical potential and somewhat larger values for strangeness neutral matter. These estimates are temperature dependent and range from $μ_B/T \gtrsim\ 2.2$ at $T=135$ MeV to $μ_B/T\ \gtrsim\ 3.2$ at $T=165$ MeV. The estimated radius of convergences is the same for any higher order cumulant.

preprint2020arXiv

Skewness, kurtosis and the 5th and 6th order cumulants of net baryon-number distributions from lattice QCD confront high-statistics STAR data

We present new results on up to $6^{th}$ order cumulants of net baryon-number fluctuations at small values of the baryon chemical potential, $μ_B$, obtained in lattice QCD calculations with physical values of light and strange quark masses. Representation of the Taylor expansions of higher order cumulants in terms of the ratio of the two lowest order cumulants, $M_B/σ_B^2=χ_1^B(T,μ_B)/χ_2^B(T,μ_B)$, allows for a parameter free comparison with data on net proton-number cumulants obtained by the STAR Collaboration in the Beam Energy Scan at RHIC. We show that recent high statistics data on skewness and kurtosis ratios of net proton-number distributions, obtained at beam energy $\sqrt{s_{_{NN}}}=54.4$ GeV, agree well with lattice QCD results on cumulants of net baryon-number fluctuations close to the pseudo-critical temperature, $T_{pc}(μ_B)$, for the chiral transition in QCD. We also present first results from a next-to-leading order expansion of $5^{th}$ and $6^{th}$ order cumulants on the line of pseudo-critical temperatures.

preprint2019arXiv

Chiral crossover in QCD at zero and non-zero chemical potentials

We present results for pseudo-critical temperatures of QCD chiral crossovers at zero and non-zero values of baryon ($B$), strangeness ($S$), electric charge ($Q$), and isospin ($I$) chemical potentials $μ_{X=B,Q,S,I}$. The results were obtained using lattice QCD calculations carried out with two degenerate up and down dynamical quarks and a dynamical strange quark, with quark masses corresponding to physical values of pion and kaon masses in the continuum limit. By parameterizing pseudo-critical temperatures as $ T_c(μ_X) = T_c(0) \left[ 1 -κ_2^{X}(μ_{X}/T_c(0))^2 -κ_4^{X}(μ_{X}/T_c(0))^4 \right] $, we determined $κ_2^X$ and $κ_4^X$ from Taylor expansions of chiral observables in $μ_X$. We obtained a precise result for $T_c(0)=(156.5\pm1.5)\;\mathrm{MeV}$. For analogous thermal conditions at the chemical freeze-out of relativistic heavy-ion collisions, i.e., $μ_{S}(T,μ_{B})$ and $μ_{Q}(T,μ_{B})$ fixed from strangeness-neutrality and isospin-imbalance, we found $κ_2^B=0.012(4)$ and $κ_4^B=0.000(4)$. For $μ_{B}\lesssim300\;\mathrm{MeV}$, the chemical freeze-out takes place in the vicinity of the QCD phase boundary, which coincides with the lines of constant energy density of $0.42(6)\;\mathrm{GeV/fm}^3$ and constant entropy density of $3.7(5)\;\mathrm{fm}^{-3}$.

preprint2019arXiv

Insight into thermal modifications of quarkonia from a comparison of continuum-extrapolated lattice results to perturbative QCD$^{\dagger}$

In this work, we strive to gain insight into thermal modifications of charmonium and bottomonium bound states as well as the heavy quark diffusion coefficient. The desired information is contained in the spectral function which can not be calculated on the lattice directly. Instead, the correlator given by an integration over the spectral function times an integration kernel is obtained. Extracting the spectral function is an ill-posed inversion problem and various different solutions have been proposed. We focus on a comparison to a spectral function obtained from combining perturbative and pNRQCD calculations. In order to get precise results, continuum extrapolated correlators originating from large and fine lattices are used. We first analyze the pseudoscalar channel since the absence of a transport peak simplifies the analysis. The knowledge gained from this is then used to extend the analysis to the vector channel, where information on heavy quark transport is encoded in the low frequency regime of the spectral function. The comparison shows a qualitatively good agreement between perturbative and lattice correlators. Quantitative differences can be explained by systematic uncertainties.

preprint2016arXiv

Lattice constraints on the thermal photon rate

We estimate the photon production rate from an SU(3) plasma at temperatures of about 1.1Tc and 1.3Tc. Lattice results for the vector current correlator at spatial momenta k ~ (2-6)T are extrapolated to the continuum limit and analyzed with the help of a polynomial interpolation for the corresponding spectral function, which vanishes at zero frequency and matches to high-precision perturbative results at large invariant masses. For small invariant masses the interpolation is compared with the NLO weak-coupling result, hydrodynamics, and a holographic model. At vanishing invariant mass we extract the photon rate which for k \gsim 3T is found to be close to the NLO weak-coupling prediction. For k \lsim 2T uncertainties remain large but the photon rate is likely to fall below the NLO prediction, in accordance with the onset of a strongly interacting behaviour characteristic of the hydrodynamic regime.

preprint2016arXiv

The curvature of the freeze-out line in heavy ion collisions

We calculate the mean and variance of net-baryon number and net-electric charge distributions from Quantum Chromodynamics (QCD) using a next-to-leading order Taylor expansion in terms of temperature and chemical potentials. We compare these expansions with experimental data from STAR and PHENIX, determine the freeze-out temperature in the limit of vanishing baryon chemical potential, and, for the first time, constrain the curvature of the freeze-out line through a direct comparison between experimental data on net-charge fluctuations and a QCD calculation. We obtain a bound on the curvature coefficient, kappa_2^f < 0.011, that is compatible with lattice QCD results on the curvature of the QCD transition line.

preprint2015arXiv

A non-perturbative estimate of the heavy quark momentum diffusion coefficient

We estimate the momentum diffusion coefficient of a heavy quark within a pure SU(3) plasma at a temperature of about 1.5Tc. Large-scale Monte Carlo simulations on a series of lattices extending up to 192^3*48 permit us to carry out a continuum extrapolation of the so-called colour-electric imaginary-time correlator. The extrapolated correlator is analyzed with the help of theoretically motivated models for the corresponding spectral function. Evidence for a non-zero transport coefficient is found and, incorporating systematic uncertainties reflecting model assumptions, we obtain kappa = (1.8 - 3.4)T^3. This implies that the "drag coefficient", characterizing the time scale at which heavy quarks adjust to hydrodynamic flow, is (1.8 - 3.4) (Tc/T)^2 (M/1.5GeV) fm/c, where M is the heavy quark kinetic mass. The results apply to bottom and, with somewhat larger systematic uncertainties, to charm quarks.

preprint2015arXiv

Conserved Charge Fluctuations from Lattice QCD and the Beam Energy Scan

We discuss the next-to-leading order Taylor expansion of ratios of cumulants of net-baryon number fluctuations. We focus on the relation between the skewness ratio, $S_Bσ_B = χ_3^B/χ_1^B$, and the kurtosis ratio, $κ_Bσ_B^2 =χ_4^B/χ_2^B$. We show that differences in these two cumulant ratios are small for small values of the baryon chemical potential. The next-to-leading order correction to $κ_Bσ_B^2$ however is approximately three times larger than that for $S_Bσ_B$. The former thus drops much more rapidly with increasing beam energy, $\sqrt{s_{NN}}$. We argue that these generic patterns are consistent with current data on cumulants of net-proton number fluctuations measured by the STAR Collaboration at $\sqrt{s_{NN}}\ge 19.6$~GeV.

preprint2015arXiv

Critical point and scale setting in SU(3) plasma: An update

We explore a method developed in statistical physics which has been argued to have exponentially small finite-volume effects, in order to determine the critical temperature Tc of pure SU(3) gauge theory close to the continuum limit. The method allows us to estimate the critical coupling betac of the Wilson action for temporal extents up to Nt ~ 20 with < 0.1% uncertainties. Making use of the scale setting parameters r0 and sqrt{t0} in the same range of beta-values, these results lead to the independent continuum extrapolations Tc r0 = 0.7457(45) and Tc sqrt{t0} = 0.2489(14), with the latter originating from a more convincing fit. Inserting a conversion of r0 from literature (unfortunately with much larger errors) yields Tc / LambdaMSbar = 1.24(10).

preprint2015arXiv

Debye mass at the QCD transition in the PNJL model

We consider colour-electric screening as expressed by the quark contribution to the Debye mass calculated in a PNJL model with emphasis on confining and chiral symmetry breaking effects. We observe that the screening mass is entirely determined by the nonperturbative quark distribution function and temperature dependent QCD running coupling. The role of the gluon background (Polyakov loop) is to provide strong suppression of the number of charge carriers below the transition temperature, as an effect of confinement, while the temperature dependent dynamical quark mass contributes additional suppression, as an effect of chiral symmetry breaking. An alternative derivation of this result from a modified kinetic theory is given, which allows for a slight generalization and explicit contact with perturbative QCD. This gives the possibility to gain insights into the colour screening mechanism in the region near the QCD pseudocritical temperature and to provide a guideline for the interpretation of lattice QCD data.

preprint2014arXiv

Additional Strange Hadrons from QCD Thermodynamics and Strangeness Freeze-out in Heavy Ion Collisions

We compare lattice QCD results for appropriate combinations of net strangeness fluctuations and their correlations with net baryon number fluctuations with predictions from two hadron resonance gas (HRG) models having different strange hadron content. The conventionally used HRG model based on experimentally established strange hadrons fails to describe the lattice QCD results in the hadronic phase close to the QCD crossover. Supplementing the conventional HRG with additional, experimentally uncharted strange hadrons predicted by quark model calculations and observed in lattice QCD spectrum calculations leads to good descriptions of strange hadron thermodynamics below the QCD crossover. We show that the thermodynamic presence of these additional states gets imprinted in the yields of the ground-state strange hadrons leading to a systematic 5--8 MeV decrease of the chemical freeze-out temperatures of ground-state strange baryons.

preprint2014arXiv

Conjugate gradient solvers on Intel Xeon Phi and NVIDIA GPUs

Lattice Quantum Chromodynamics simulations typically spend most of the runtime in inversions of the Fermion Matrix. This part is therefore frequently optimized for various HPC architectures. Here we compare the performance of the Intel Xeon Phi to current Kepler-based NVIDIA Tesla GPUs running a conjugate gradient solver. By exposing more parallelism to the accelerator through inverting multiple vectors at the same time, we obtain a performance greater than 300 GFlop/s on both architectures. This more than doubles the performance of the inversions. We also give a short overview of the Knights Corner architecture, discuss some details of the implementation and the effort required to obtain the achieved performance.

preprint2014arXiv

HISQ inverter on Intel Xeon Phi and NVIDIA GPUs

The runtime of a Lattice QCD simulation is dominated by a small kernel, which calculates the product of a vector by a sparse matrix known as the "Dslash" operator. Therefore, this kernel is frequently optimized for various HPC architectures. In this contribution we compare the performance of the Intel Xeon Phi to current Kepler-based NVIDIA Tesla GPUs running a conjugate gradient solver. By exposing more parallelism to the accelerator through inverting multiple vectors at the same time we obtain a performance 250 GFlop/s on both architectures. This more than doubles the performance of the inversions. We give a short overview of both architectures, discuss some details of the implementation and the effort required to obtain the achieved performance.

preprint2014arXiv

Quark mass dependence of quarkonium properties at finite temperature

Quarkonium properties at finite temperature have been studied with quark masses of the charm and bottom quarks. Our simulations have been performed in quenched QCD with the $O(a)$-improved Wilson quarks on large and fine isotropic lattices with the spatial lattice extents $N_σ=$ 96, 192 and the corresponding lattice spacings $a =$ 0.0190, 0.00967 fm, respectively, at temperatures in a range between about 0.7$T_c$ and 1.4$T_c$. We show temperature and quark mass dependence of quarkonium correlation functions and related physical quantities: the quark number susceptibility and the heavy quark diffusion constant.

preprint2014arXiv

The melting and abundance of open charm hadrons

Ratios of cumulants of conserved net charge fluctuations are sensitive to the degrees of freedom that are carriers of the corresponding quantum numbers in different phases of strong interaction matter. Using lattice QCD with 2+1 dynamical flavors and quenched charm quarks we calculate second and fourth order cumulants of net charm fluctuations and their correlations with other conserved charges such as net baryon number, electric charge and strangeness. Analyzing appropriate ratios of these cumulants we probe the nature of charmed degrees of freedom in the vicinity of the QCD chiral crossover region. We show that for temperatures above the chiral crossover transition temperature, charmed degrees of freedom can no longer be described by an uncorrelated gas of hadrons. This suggests that the dissociation of open charm hadrons and the emergence of deconfined charm states sets in just near the chiral crossover transition. Till the crossover region we compare these lattice QCD results with two hadron resonance gas models --including only the experimentally established charmed resonances and also including additional states predicted by quark model and lattice QCD calculations. This comparison provides evidence for so far unobserved charmed hadrons that contribute to the thermodynamics in the crossover region.

preprint2013arXiv

Strangeness at high temperatures: from hadrons to quarks

Appropriate combinations of up to fourth order cumulants of net strangeness fluctuations and their correlations with net baryon number and electric charge fluctuations, obtained from lattice QCD calculations, have been used to probe the strangeness carrying degrees of freedom at high temperatures. For temperatures up to the chiral crossover separate contributions of strange mesons and baryons can be well described by an uncorrelated gas of hadrons. Such a description breaks down in the chiral crossover region, suggesting that the deconfinement of strangeness takes place at the chiral crossover. On the other hand, the strangeness carrying degrees of freedom inside the quark gluon plasma can be described by a weakly interacting gas of quarks only for temperatures larger than twice the chiral crossover temperature. In the intermediate temperature window these observables show considerably richer structures, indicative of the strongly interacting nature of the quark gluon plasma.

preprint2013arXiv

Thermal dilepton rates from quenched lattice QCD

We present new lattice results on the continuum extrapolation of the vector current correlation function. Lattice calculations have been carried out in the deconfined phase at a temperature of 1.1 Tc, extending our previous results at 1.45 Tc, utilizing quenched non-perturbatively clover-improved Wilson fermions and light quark masses. A systematic analysis on multiple lattice spacings allows to perform the continuum limit of the correlation function and to extract spectral properties in the continuum limit. Our current analysis suggests the results for the electrical conductivity are proportional to the temperature and the thermal dilepton rates in the quark gluon plasma are comparable for both temperatures. Preliminary results of the continuum extrapolated correlation function at finite momenta, which relates to thermal photon rates, are also presented.

preprint2013arXiv

Towards the continuum limit in transport coefficient computations

The analytic continuation needed for the extraction of transport coefficients necessitates in principle a continuous function of the Euclidean time variable. We report on progress towards achieving the continuum limit for 2-point correlator measurements in thermal SU(3) gauge theory, with specific attention paid to scale setting. In particular, we improve upon the determination of the critical lattice coupling and the critical temperature of pure SU(3) gauge theory, estimating r0*Tc ~ 0.7470(7) after a continuum extrapolation. As an application the determination of the heavy quark momentum diffusion coefficient from a correlator of colour-electric fields attached to a Polyakov loop is discussed.

preprint2012arXiv

Charmonium dissociation and heavy quark transport in hot quenched lattice QCD

We study the properties of charmonium states at finite temperature in quenched lattice QCD on large and fine isotropic lattices. We perform a detailed analysis of charmonium correlation and spectral functions both below and above Tc. Our analysis suggests that the S wave states disappear at about 1.5 Tc. The charm diffusion coefficient is estimated and found to be approximately 1/πT at 1.5Tc {\leq} T {\leq} 3Tc.

preprint2012arXiv

Charmonium properties in hot quenched lattice QCD

We study the properties of charmonium states at finite temperature in quenched QCD on large and fine isotropic lattices. We perform a detailed analysis of charmonium correlation and spectral functions both below and above $T_c$. Our analysis suggests that both S wave states ($J/ψ$ and $η_c$) and P wave states ($χ_{c0}$ and $χ_{c1}$) disappear already at about $1.5 T_c$. The charm diffusion coefficient is estimated through the Kubo formula and found to be compatible with zero below $T_c$ and approximately $1/πT$ at $1.5 T_c\lesssim T\lesssim 3 T_c$.

preprint2012arXiv

Freeze-out Conditions in Heavy Ion Collisions from QCD Thermodynamics

We present a determination of chemical freeze-out conditions in heavy ion collisions based on ratios of cumulants of net electric charge fluctuations. These ratios can reliably be calculated in lattice QCD for a wide range of chemical potential values by using a next-to-leading order Taylor series expansion around the limit of vanishing baryon, electric charge and strangeness chemical potentials. From a computation of up to fourth order cumulants and charge correlations we first determine the strangeness and electric charge chemical potentials that characterize freeze-out conditions in a heavy ion collision and confirm that in the temperature range 150 MeV < T < 170 MeV the hadron resonance gas model provides good approximations for these parameters that agree with QCD calculations on the (5-15)% level. We then show that a comparison of lattice QCD results for ratios of up to third order cumulants of electric charge fluctuations with experimental results allows to extract the freeze-out baryon chemical potential and the freeze-out temperature.

preprint2012arXiv

Recent Developments in Lattice Studies for Quarkonia

After discussing results of dilepton rates and electrical conductivity obtained from continuum extrapolated results of light quark correlation functions in quenched QCD I will give a review on recent developments in lattice QCD studies for quarkonia in the quark gluon plasma. Recent progress in the extraction of spectral properties from lattice QCD calculations of hadronic correlation functions will be discussed. Besides medium modifications of bound states and their dissociation in the plasma I will focus on transport coefficients, like heavy quark diffusion constants extracted from different correlation functions on the lattice. Present limitations and future perspectives for studies of quarkonia and related transport coefficients on the lattice will be discussed.

preprint2011arXiv

Heavy Quark diffusion from lattice QCD spectral functions

We analyze the low frequency part of charmonium spectral functions on large lattices close to the continuum limit in the temperature region $1.5\lesssim T/T_c\lesssim 3$ as well as for $T \simeq 0.75T_c$. We present evidence for the existence of a transport peak above $T_c$ and its absence below $T_c$. The heavy quark diffusion constant is then estimated using the Kubo formula. As part of the calculation we also determine the temperature dependence of the signature for the charmonium bound state in the spectral function and discuss the fate of charmonium states in the hot medium.

preprint2011arXiv

Meson screening masses from lattice QCD with two light and the strange quark

We present results for screening masses of mesons built from light and strange quarks in the temperature range of approximately between 140 MeV to 800 MeV. The lattice computations were performed with 2+1 dynamical light and strange flavors of improved (p4) staggered fermions along a line of constant physics defined by a pion mass of about 220 MeV and a kaon mass of 500 MeV. The lattices had temporal extents Nt = 4, 6 and 8 and aspect ratios of Ns / Nt \geq 4. At least up to a temperature of 140 MeV the pseudo-scalar screening mass remains almost equal to the corresponding zero temperature pseudo-scalar (pole) mass. At temperatures around 3Tc (Tc being the transition temperature) the continuum extrapolated pseudo-scalar screening mass approaches very close to the free continuum result of 2 πT from below. On the other hand, at high temperatures the vector screening mass turns out to be larger than the free continuum value of 2 πT. The pseudo-scalar and the vector screening masses do not become degenerate even for a temperature as high as 4Tc. Using these mesonic spatial correlation functions we have also investigated the restoration of chiral symmetry and the effective restoration of the axial symmetry. We have found that the vector and the axial-vector screening correlators become degenerate, indicating chiral symmetry restoration, at a temperature which is consistent with the QCD transition temperature obtained in previous studies. On the other hand, the pseudo-scalar and the scalar screening correlators become degenerate only at temperatures larger than 1.3Tc, indicating that the effective restoration of the axial symmetry takes place at a temperature larger than the QCD transition temperature.

preprint2011arXiv

Thermal dilepton rate and electrical conductivity: An analysis of vector current correlation functions in quenched lattice QCD

We calculate the vector current correlation function for light valence quarks in the deconfined phase of QCD. The calculations have been performed in quenched lattice QCD at T=1.45 Tc for four values of the lattice cut-off on lattices up to size 128^3x48. This allows to perform a continuum extrapolation of the correlation function in the Euclidean time interval tau*T -in [0.2, 0.5], which extends to the largest temporal separations possible at finite temperature, to better than 1% accuracy. In this interval, at the value of the temperature investigated, we find that the vector correlation function never deviates from the free correlator for massless quarks by more than 9%. We also determine the first two non-vanishing thermal moments of the vector meson spectral function. The second thermal moment deviates by less than 7% from the free value. With these constraints, we then proceed to extract information on the spectral representation of the vector correlator and discuss resulting consequences for the electrical conductivity and the thermal dilepton rate in the plasma phase.

preprint2011arXiv

Towards a non-perturbative measurement of the heavy quark momentum diffusion coefficient

We report on a lattice investigation of heavy quark diffusion within pure SU(3) plasma above the deconfinement transition, with the quarks treated to leading order in the heavy mass expansion. Using a multilevel algorithm, several volumes and lattice spacings, as well as tree-level improvement and perturbative renormalization, we measure the relevant "colour-electric" Euclidean correlator, finding that it clearly exceeds its perturbative counterpart. Even without analytic continuation, this suggests that at temperatures just above the critical one, non-perturbative interactions felt by the heavy quarks are stronger than within the weak-coupling expansion. After introducing rough modelling of the spectral shape, diffusion coefficients down to D ~ 0.5/T appear possible.

preprint2011arXiv

Width of the QCD transition in a Polyakov-loop DSE model

We consider the pseudocritical temperatures for the chiral and deconfinement transitions within a Polyakov-loop Dyson-Schwinger equation approach which employs a nonlocal rank-2 separable model for the effective gluon propagator. These pseudocritical temperatures differ by a factor of two when the quark and gluon sectors are considered separately, but get synchronized and become coincident when their coupling is switched on. The coupling of the Polyakov-loop to the chiral quark dynamics narrows the temperature region of the QCD transition in which chiral symmetry and deconfinement is established. We investigate the effect of rescaling the parameter T_0 in the Polyakov-loop potential on the QCD transition for both the logarithmic and polynomial forms of the potential. While the critical temperatures vary in a similar way, the width of the transition is stronger affected for the logarithmic potential. For this potential the character of the transition changes from crossover to a first order one when T_0 < 210 MeV, but it remains crossover in the whole range of relevant T_0 values for the polynomial form.

preprint2010arXiv

Charmonium correlation and spectral functions at finite temperature

We study the properties of charmonium states at finite temperature in quenched QCD on isotropic lattices. We measured charmonium correlators using non-perturbatively $\cO(a)$ improved clover fermions on fine ($a=0.01$ fm) lattices with a relatively large size of $128^{3}\times 96$, $128^3\times48$, $128^3\times32$ and $128^3\times24$ at $0.73~T_c$, $1.46~T_c$, $2.20~T_c$ and $2.93~T_c$, respectively. Our analysis suggests that $\Jpsi$ is melted already at $1.46~T_c$ and $η_c$ starts to dissolve at $1.46~T_c$ and does not exist at higher temperatures. We also identify the heavy quark transport contribution at the spectral function level for the first time.

preprint2010arXiv

Equation of State for physical quark masses

We calculate the QCD equation of state for temperatures corresponding to the transition region with physical mass values for two degenerate light quark flavors and a strange quark using an improved staggered fermion action (p4-action) on lattices with temporal extent N_tau=8. We compare our results with previous calculations performed at twice larger values of the light quark masses as well as with results obtained from a resonance gas model calculation. We also discuss the deconfining and chiral aspects of the QCD transition in terms of renormalized Polyakov loop, strangeness fluctuations and subtracted chiral condensate. We show that compared to the calculations performed at twice larger value of the light quark mass the transition region shifts by about 5 MeV toward smaller temperatures

preprint2010arXiv

The phase boundary for the chiral transition in (2+1)-flavor QCD at small values of the chemical potential

We determine the chiral phase transition line in (2+1)-flavor QCD for small values of the light quark chemical potential. We show that for small values of the chemical potential the curvature of the phase transition line can be deduced from an analysis of scaling properties of the chiral condensate and its susceptibilities. To do so we extend earlier studies of the magnetic equation of state in (2+1)-flavor QCD to finer lattice spacings, aT=1/8. We use these universal scaling properties of the chiral order parameter to extract the curvature of the transition line at two values of the cut-off, aT=1/4 and 1/8. We find that cut-off effects are small for the curvature parameter and determine the transition line in the chiral limit to leading order in the light quark chemical potential. We obtain Tc(μ_q)/Tc(0) = 1 - 0.059(2)(4) (μ_q/T)^2 +O(μ_q^4).

preprint2009arXiv

Baryon Number, Strangeness and Electric Charge Fluctuations in QCD at High Temperature

We analyze baryon number, strangeness and electric charge fluctuations as well as their correlations in QCD at high temperature. We present results obtained from lattice calculations performed with an improved staggered fermion action (p4-action) at two values of the lattice cut-off with almost physical up and down quark masses and a physical value for the strange quark mass. We compare these results, with an ideal quark gas at high temperature and a hadron resonance gas model at low temperature. We find that fluctuations and correlations are well described by the former already for temperatures about 1.5 times the transition temperature. At low temperature qualitative features of the lattice results are quite well described by a hadron resonance gas model. Higher order cumulants, which become increasingly sensitive to the light pions, however show deviations from a resonance gas in the vicinity of the transition temperature.

preprint2009arXiv

Charmonium correlators and spectral functions at finite temperature

We study charmonium correlators and spectral functions in quenched QCD, using Clover improved Wilson fermions on very fine (0.015 fm) isotropic lattices at 0.75 Tc and 1.5 Tc. We use a new approach to distinguish the zero mode contribution from the other contributions. Once this is removed, we find that the ratios of correlators to reconstructed correlators remain almost unity at all distances. The ground state peaks of spectral functions obtained at 0.75 Tc are reliable and robust. The present accuracy and limited number of points in the temporal direction at 1.5 Tc do not allow for a reliable conclusion about a possible melting of charmonium states in the QGP.

preprint2007arXiv

Color Screening and Quark-Quark Interactions in Finite Temperature QCD

We analyze the screening of static diquark sources in 2-flavor QCD and compare results with the screening of static quark-antiquark pairs. We show that a two quark system in a fixed color representations is screened at short distances like a single quark source in the same color representation whereas at large distances the two quarks are screened independently. At high temperatures we observe that the relative strength of the interaction in diquark and quark-antiquark systems, respectively, obeys Casimir scaling. We use this result to examine the possible existence of heavy quark-quark bound states in the high temperature phase of QCD. We find support for the existence of $bb$ states up to about $2T_c$ while $cc$ states are unlikely to be formed above $T_c$.

preprint2007arXiv

Screening at finite temperature and density

We present lattice QCD results on heavy quark free energies, extract from its temperature dependence entropy and internal energy contributions, and discuss the onset of medium effects that lead to screening of static quark-antiquark sources in a thermal medium. Most results are obtained in (2+1)-flavour QCD on a line of constant physics with almost realistic quark masses and compared to previous results from 2-flavor QCD as well as pure gauge theory. Furthermore, we discuss results on the density dependence of screening masses that have been obtained using a leading order Taylor expansion in the baryon chemical potential.

preprint2005arXiv

Heavy quark free energies and screening at finite temperature and density

We study the free energies of heavy quarks calculated from Polyakov loop correlation functions in full 2-flavour QCD using the p4-improved staggered fermion action. A small but finite Baryon number density is included via Taylor expansion of the fermion determinant in the Baryo-chemical potential mu. For temperatures above Tc we extract Debye screening masses from the large distance behaviour of the free energies and compare their mu-dependence to perturbative results.

preprint2005arXiv

Heavy quark potential and quarkonia dissociation rates

Quenched lattice data for the quark-antiquark interaction (in terms of heavy quark free energies) in the color singlet channel at finite temperatures are fitted and used within the nonrelativistic Schroedinger equation formalism to obtain binding energies and scattering phase shifts for the lowest eigenstates in the charmonium and bottomonium systems in a hot gluon plasma. The partial dissociation rate due to the Bhanot-Peskin process is calculated using different assumptions for the gluon distribution function, including free massless gluons, massive gluons, and massive damped gluons. It is demonstrated that a temperature dependent gluon mass has an essential influence on the heavy quarkonia dissociation, but that this process alone is insufficient to describe the heavy quarkonia dissociation rates.

preprint2005arXiv

Heavy two- and three-quark free energies at finite $T$

We investigate the behaviour of heavy $QQ$- and $QQQ$-free energies in SU(3) pure gauge theory at finite temperature. After stating some general properties of the free energies, we discuss the $QQ$-antitriplet and -sextet channels as well as free energies in different $QQQ$-channels. We perform a detailed finite volume scaling analysis of $QQ$-free energies below and above $T_c$. Moreover, we also discuss the relation of the $QQ$-antitriplet to the $Q\bar Q$-singlet free energies. A comparison of the effective running coupling indicates that the $r$-dependence of the free energies differs by a factor of two in agreement with the different Casimir factors. We observe that the baryonic free energies can be decomposed into a sum of two particle interactions.

preprint2005arXiv

The QCD equation of state for two flavours at non-zero chemical potential

We present results of a simulation of 2 flavour QCD on a $16^3\times4$ lattice using p4-improved staggered fermions with bare quark mass $m/T=0.4$. Derivatives of the thermodynamic grand canonical partition function $Z(V,T,μ_u,μ_d)$ with respect to chemical potentials $μ_{u,d}$ for different quark flavours are calculated up to sixth order, enabling estimates of the pressure and the quark number density as well as the chiral condensate and various susceptibilities as functions of $μ_{u,d}$ via Taylor series expansion. Results are compared to high temperature perturbation theory as well as a hadron resonance gas model. We also analyze baryon as well as isospin fluctuations and discuss the relation to the chiral critical point in the QCD phase diagram. We moreover discuss the dependence of the heavy quark free energy on the chemical potential.

preprint2004arXiv

Static quark anti-quark free energy and the running coupling at finite temperature

We analyze the free energy of a static quark anti-quark pair in quenched QCD at short and large distances. From this we deduce running couplings, g^2(r,T), and determine the length scale that separates at high temperature the short distance perturbative regime from the large distance non-perturbative regime in the QCD plasma phase. Ambiguities in the definition of a coupling beyond the perturbative regime are discussed in their relation to phenomenological considerations on heavy quark bound states in the quark gluon plasma. Our analysis suggests that it is more appropriate to characterize the non-perturbative properties of the QCD plasma phase close to T_c in terms remnants of the confinement part of the QCD force rather than a strong Coulombic force.

preprint2003arXiv

The Equation of State for Two Flavor QCD at Non-zero Chemical Potential

We present results of a simulation of QCD on a 4x16^3 lattice with 2 continuum flavors of p4-improved staggered fermion with mass m/T=0.4. Derivatives of the thermodynamic grand potential with respect to quark chemical potential mu_q up to fourth order are calculated, enabling estimates of the pressure, quark number density and associated susceptibilities as functions of mu_q via Taylor series expansion. Discretisation effects associated with various staggered fermion formulations are discussed in some detail. In addition it is possible to estimate the radius of convergence of the expansion as a function of temperature. We also discuss the calculation of energy and entropy densities which are defined via mixed derivatives of the thermodynamic grand potential with respect to the bare couplings and quark masses.

preprint2001arXiv

Short distance physics with heavy quark potentials

We present lattice studies of heavy quark potentials in the quenched approximation of QCD at finite temperatures. Both, the color singlet and color averaged potentials are calculated. While the potentials are well known at large distances, we give a detailed analysis of their short distance behavior (from 0.015 fm to 1 fm) near the critical temperature. At these distances we expect that the T-dependent potentials go over into the zero temperature potential. Indeed, we find evidences that the temperature influence gets suppressed and the potentials starts to become a unique function of the underlying distance scale. We use this feature to normalize the heavy quark potentials at short distances and extract the free energy of the quark system in a gluonic heat bath.

preprint1998arXiv

Thermodynamics of two-colour QCD

We discuss the thermodynamics of two-colour QCD with four flavours of staggered quarks on 8^3x4 and 16^3x4 lattices. In our simulations we use the Naik action for the fermions and a (1,2) tree-level improved gauge action. We analyze the deconfinement and chiral phase transitions for four different quark masses (m=0.1,0.05,0.025,0.015). Contrary to three-colour QCD the peak in the Polyakov loop susceptibility decreases with decreasing quark mass. This reflects an early breaking of the string in the heavy quark potential, which we verify explicitly by calculating the heavy quark potential at finite temperature using Polyakov loop correlations.