Source author record

Sandor D. Katz

Sandor D. Katz 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

28works
3topics
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

28 published item(s)

preprint2022arXiv

Lattice simulations of the QCD chiral transition at real $μ_B$

Most lattice studies of hot and dense QCD matter rely on extrapolation from zero or imaginary chemical potentials. The ill-posedness of numerical analytic continuation puts severe limitations on the reliability of such methods. We studied the QCD chiral transition at finite real baryon density with the more direct sign reweighting approach. We simulate up to a baryochemical potential-temperature ratio of $μ_B/T=2.7$, covering the RHIC Beam Energy Scan range, and penetrating the region where methods based on analytic continuation are unpredictive.This opens up a new window to study QCD matter at finite $μ_B$ from first principles.

preprint2022arXiv

Radius of convergence in lattice QCD at finite $μ_B$ with rooted staggered fermions

In typical statistical mechanical systems the grand canonical partition function at finite volume is proportional to a polynomial of the fugacity $e^{μ/T}$. The zero of this Lee-Yang polynomial closest to the origin determines the radius of convergence of the Taylor expansion of the pressure around $μ=0$. The computationally cheapest formulation of lattice QCD, rooted staggered fermions, with the usual definition of the rooted determinant, does not admit such a Lee-Yang polynomial. We argue that the radius of convergence is then bounded by the spectral gap of the reduced matrix of the unrooted staggered operator. This is a cutoff effect that potentially affects all estimates of the radius of convergence with the standard staggered rooting. We suggest a new definition of the rooted staggered determinant at finite chemical potential that allows for a definition of a Lee-Yang polynomial, and, therefore of the numerical study of Lee-Yang zeros. We also describe an algorithm to determine the Lee-Yang zeros and apply it to configurations generated with the 2-stout improved staggered action at $N_t = 4$. We perform a finite-volume scaling study of the leading Lee-Yang zeros and estimate the radius of convergence of the Taylor expansion extrapolated to an infinite volume. We show that the limiting singularity is not on the real line, thus giving a lower bound on the location of any possible phase transitions at this lattice spacing. In the vicinity of the crossover temperature at zero chemical potential, the radius of convergence turns out to be $μ_B/T \approx 2$ and roughly temperature independent. Our simulations are performed at strange quark chemical potential $μ_s=0$, but the method can be straightforwardly extended to strangeness chemical potential $μ_S=0$ or strangeness neutrality.

preprint2021arXiv

Lattice simulations of the QCD chiral transition at real baryon density

State-of-the-art lattice QCD studies of hot and dense strongly interacting matter currently rely on extrapolation from zero or imaginary chemical potentials. The ill-posedness of numerical analytic continuation puts severe limitations on the reliability of such methods. Here we use the more direct sign reweighting method to perform lattice QCD simulation of the QCD chiral transition at finite real baryon density on phenomenologically relevant lattices. This method does not require analytic continuation and avoids the overlap problem associated with generic reweighting schemes, so has only statistical but no uncontrolled systematic uncertainties for a fixed lattice setup. This opens up a new window to study hot and dense strongly interacting matter from first principles. We perform simulations up to a baryochemical potential-temperature ratio of $μ_B/T=2.5$ covering most of the RHIC Beam Energy Scan range in the chemical potential. We also clarify the connection of the approach to the more traditional phase reweighting method.

preprint2020arXiv

New approach to lattice QCD at finite density; results for the critical end point on coarse lattices

All approaches currently used to study finite baryon density lattice QCD suffer from uncontrolled systematic uncertainties in addition to the well-known sign problem. We formulate and test an algorithm, sign reweighting, that works directly at finite $μ= μ_B/3$ and is yet free from any such uncontrolled systematics. With this algorithm the {\em only} problem is the sign problem itself. This approach involves the generation of configurations with the positive fermionic weight $|{\rm Re\; det} D(μ)|$ where $D(μ)$ is the Dirac matrix and the signs ${\rm sign} \; ( {\rm Re\; det} D(μ) ) = \pm 1$ are handled by a discrete reweighting. Hence there are only two sectors, $+1$ and $-1$ and as long as the average $\langle\pm 1\rangle \neq 0$ (with respect to the positive weight) this discrete reweighting by the signs carries no overlap problem and the results are reliable. The approach is tested on $N_t = 4$ lattices with $2+1$ flavors and physical quark masses using the unimproved staggered discretization. By measuring the Fisher (sometimes also called Lee-Yang) zeros in the bare coupling on spatial lattices $L/a = 8, 10, 12$ we conclude that the cross-over present at $μ= 0$ becomes stronger at $μ> 0$ and is consistent with a true phase transition at around $μ_B/T \sim 2.4$.

preprint2020arXiv

The effect of stout smearing on the phase diagram from multiparameter reweigthing in lattice QCD

The phase diagram and the location of the critical endpoint (CEP) of lattice QCD with unimproved staggered fermions on a $N_t=4$ lattice was determined fifteen years ago with the multiparameter reweighting method by studying Fisher zeros. We first reproduce the old result with an exact algorithm (not known at the time) and with statistics larger by an order of magnitude. As an extension of the old analysis we introduce stout smearing in the fermion action in order to reduce the finite lattice spacing effects. First we show that increasing the smearing parameter $ρ$ the crossover at $μ= 0$ gets weaker, i.e., the leading Fisher zero gets farther away from the real axis. Furthermore as the chemical potential is increased the overlap problem gets severe sooner than in the unimproved case, therefore shrinking the range of applicability of the method. Nevertheless certain qualitative features remain, even after introducing the smearing. Namely, at small chemical potentials the Fisher zeros first get farther away from the real axis and later at around $aμ_q = 0.1 - 0.15$ they start to get closer, i.e., the crossover first gets weaker and later stronger as a function of $μ$. However, because of the more severe overlap problem the CEP is out of reach with the smeared action.

preprint2020arXiv

The QCD crossover at finite chemical potential from lattice simulations

We provide the most accurate results for the QCD transition line so far. We optimize the definition of the crossover temperature $T_c$, allowing for its very precise determination, and extrapolate from imaginary chemical potential up to real $μ_B \approx 300$ MeV. The definition of $T_c$ adopted in this work is based on the observation that the chiral susceptibility as a function of the condensate is an almost universal curve at zero and imaganiary $μ_B$. We obtain the parameters $κ_2=0.0153(18)$ and $κ_4=0.00032(67)$ as a continuum extrapolation based on $N_t=10,12$ and $16$ lattices with physical quark masses. We also extrapolate the peak value of the chiral susceptibility and the width of the chiral transition along the crossover line. In fact, both of these are consistent with a constant function of $μ_B$. We see no sign of criticality in the explored range.

preprint2020arXiv

Towards a reliable lower bound on the location of the critical endpoint

We perform the first direct determination of the position of the leading singularity of the pressure in the complex chemical potential $μ_B$ plane in lattice QCD using numerical simulations with 2-stout improved rooted staggered fermions. This provides a direct determination of the radius of convergence of the Taylor expansion of the pressure that does not rely on a finite-order truncation of the expansion. The analyticity issues in the complex $μ_B$ plane of the grand canonical partition function of QCD with rooted staggered fermions are solved with a careful redefinition of the fermion determinant that makes it a polynomial in the fugacity on any finite lattice, without changing the continuum limit of the observables. By performing a finite volume scaling study at a single coarse lattice spacing, we show that the limiting singularity is not on the real line in the thermodynamic limit, thus showing that the radius of convergence of the Taylor expansion gives a lower bound on the location of a possible phase transition. In the vicinity of the crossover temperature at zero chemical potential, the radius of convergence turns out to be $μ_B/T \approx 2$ and roughly temperature independent.

preprint2019arXiv

Magnetic catalysis and inverse catalysis for heavy pions

We investigate the QCD phase diagram for nonzero background magnetic fields using first-principles lattice simulations. At the physical point (in terms of quark masses), the thermodynamics of this system is controlled by two opposing effects: magnetic catalysis (enhancement of the quark condensate) at low temperature and inverse magnetic catalysis (reduction of the condensate) in the transition region. While the former is known to be robust and independent of the details of the interactions, inverse catalysis arises as a result of a delicate competition, effective only for light quarks. By performing simulations at different quark masses, we determine the pion mass above which inverse catalysis does not take place in the transition region anymore. Even for pions heavier than this limiting value - where the quark condensate undergoes magnetic catalysis - our results are consistent with the notion that the transition temperature is reduced by the magnetic field. These findings will be useful to guide low-energy models and effective theories of QCD.

preprint2016arXiv

Landau Levels in Lattice QCD

The spectrum of the two-dimensional continuum Dirac operator in the presence of a uniform background magnetic field consists of Landau levels, which are degenerate and separated by gaps. On the lattice the Landau levels are spread out by discretization artefacts, but a remnant of their structure is clearly visible (Hofstadter butterfly). If one switches on a non-Abelian interaction, the butterfly structure will be smeared out, but the lowest Landau level (LLL) will still be separated by a gap from the rest of the spectrum. In this talk we discuss how one can define the LLL in QCD and check how well certain physical quantities are approximated by taking into account only the LLL.

preprint2015arXiv

QCD thermodynamics with continuum extrapolated Wilson fermions II

We continue our investigation of 2+1 flavor QCD thermodynamics using dynamical Wilson fermions in the fixed scale approach. Two additional pion masses, approximately 440 MeV and 285 MeV, are added to our previous work at 545 MeV. The simulations were performed at 3 or 4 lattice spacings at each pion mass. The renormalized chiral condensate, strange quark number susceptibility and Polyakov loop is obtained as a function of the temperature and we observe a decrease in the light chiral pseudo-critical temperature as the pion mass is lowered while the pseudo-critical temperature associated with the strange quark number susceptibility or the Polyakov loop is only mildly sensitive to the pion mass. These findings are in agreement with previous continuum results obtained in the staggered formulation.

preprint2015arXiv

Static quark-antiquark pair free energy and screening masses: continuum results at the QCD physical point

We study the correlators of Polyakov loops, and the corresponding gauge invariant free energy of a static quark-antiquark pair in 2+1 flavor QCD at finite temperature. Our simulations were carried out on $N_t$ = 6, 8, 10, 12, 16 lattices using a Symanzik improved gauge action and a stout improved staggered action with physical quark masses. The free energies calculated from the Polyakov loop correlators are extrapolated to the continuum limit. For the free energies we use a two step renormalization procedure that only uses data at finite temperature. We also measure correlators with definite Euclidean time reversal and charge conjugation symmetry to extract two different screening masses, one in the magnetic, and one in the electric sector, to distinguish two different correlation lengths in the full Polyakov loop correlator. This conference contribution is based on the paper: JHEP 1504 (2015) 138

preprint2014arXiv

Charmonium spectral functions from 2+1 flavour lattice QCD

Finite temperature charmonium spectral functions in the pseudoscalar and vector channels are studied in lattice QCD with 2+1 flavours of dynamical Wilson quarks, on fine isotropic lattices (with a lattice spacing of 0.057 fm), with a non-physical pion mass of $m_π \approx$ 545 MeV. The highest temperature studied is approximately $1.4 T_c$. Up to this temperature no significant variation of the spectral function is seen in the pseudoscalar channel. The vector channel shows some temperature dependence, which seems to be consistent with a temperature dependent low frequency peak related to heavy quark transport, plus a temperature independent term at ω>0. These results are in accord with previous calculations using the quenched approximation.

preprint2014arXiv

Full result for the QCD equation of state with 2+1 flavors

We present a full result for the 2+1 flavor QCD equation of state. All the systematics are controlled, the quark masses are set to their physical values, and the continuum extrapolation is carried out. This extends our previous studies [JHEP 0601:089 (2006); 1011:077 (2010)] to even finer lattices and now includes ensembles with Nt = 6,8,10,12 up to Nt = 16. We use a Symanzik improved gauge and a stout-link improved staggered fermion action. Our findings confirm our earlier results. In order to facilitate the direct use of our equation of state we make our tabulated results available for download.

preprint2014arXiv

Recent results on the Equation of State of QCD

We report on a continuum extrapolated result (arXiv:1309.5258) for the equation of state (EoS) of QCD with $N_f=2+1$ dynamical quark flavors and discuss preliminary results obtained with an additional dynamical charm quark ($N_f=2+1+1$). For all our final results, the systematics are controlled, quark masses are set to their physical values, and the continuum limit is taken using at least three lattice spacings corresponding to temporal extents up to $N_t=16$.

preprint2014arXiv

The chiral transition as an Anderson transition

At low temperature the low-lying QCD Dirac spectrum obeys random matrix statistics. Recently we found that above $T_{c}$ the lowest part of the spectrum consists of localized modes that obey Poisson statistics. An interesting implication of this is that as the system crosses $T_{c}$ from above, the spectral statistics at $λ=0$ changes from Poisson to random matrix. Here we study this transition and its possible implications for the finite temperature transition of QCD-like theories.

preprint2012arXiv

Systematic errors in partially-quenched QCD plus QED lattice simulations

At the precision reached in current lattice QCD calculations, electromagnetic effects are becoming numerically relevant. Here, electromagnetic effects are included by superimposing $\mathrm{U}(1)$ degrees of freedom on $N_f = 2+1$ QCD configurations from the Budapest-Marseille-Wuppertal Collaboration. We present preliminary results for the electromagnetic corrections to light pseudoscalars mesons masses and discuss some of the associated systematic errors.

preprint2011arXiv

Correlations and fluctuations from lattice QCD

We present the new results of the Wuppertal-Budapest lattice QCD collaboration on flavor diagonal and non-diagonal quark number susceptibilities with 2+1 staggered quark flavors, in a temperature regime between 120 and 400 MeV. A Symanzik improved gauge and a stout-link improved staggered fermion action is utilized; the light and strange quark masses are set to their physical values. Lattices with $N_t=6, 8, 10, 12$ are used. We perform a continuum extrapolation of those observables for which the scaling regime is reached, and discretization errors are under control.

preprint2011arXiv

Fluctuations of conserved charges at finite temperature from lattice QCD

We present the full results of the Wuppertal-Budapest lattice QCD collaboration on flavor diagonal and non-diagonal quark number susceptibilities with 2+1 staggered quark flavors, in a temperature range between 125 and 400 MeV. The light and strange quark masses are set to their physical values. Lattices with Nt=6, 8, 10, 12, 16 are used. We perform a continuum extrapolation of all observables under study. A Symanzik improved gauge and a stout-link improved staggered fermion action is utilized. All results are compared to the Hadron Resonance Gas model predictions: good agreement is found in the temperature region below the transition.

preprint2011arXiv

QCD thermodynamics with Wilson fermions

QCD is investigated at finite temperature using Wilson fermions in the fixed scale approach. A 2+1 flavor stout and clover improved action is used at four lattice spacings allowing for control over discretization errors. The light quark masses in this first study are fixed to heavier than physical values. The renormalized chiral condensate, quark number susceptibility and the Polyakov loop is measured and the results are compared with the staggered formulation in the fixed N_t approach. The Wilson results at the finest lattice spacing agree with the staggered results at the highest N_t.

preprint2011arXiv

SU(2) chiral perturbation theory low-energy constants from staggered 2+1 flavor simulations

We measure the pion mass and decay constant on ensembles generated by the Wuppertal-Budapest Collaboration, and extract the NLO low-energy constants l_3 and l_4 of SU(2) chiral perturbation theory. The data are generated in 2+1 flavor simulations with Symanzik glue and 2-fold stout-smeared staggered fermions, with pion masses varying from 135 MeV to 400 MeV, lattice scales between 0.7 GeV and 2.0 GeV, and m_s kept at its physical value. Furthermore, by excluding the lightest mass points, we are able to test the reliability of SU(2) chPT as a tool to extrapolate towards the physical point from higher pion masses.

preprint2011arXiv

Transition temperature and the equation of state from lattice QCD, Wuppertal-Budapest results

The QCD transition is studied on lattices up to $N_t=16$. The chiral condensate is presented as a function of the temperature, and the corresponding transition temperature is extracted. The equation of state is determined on lattices with $N_t=6,8,10$ and at some temperature values with $N_t=12$. The pressure and the trace anomaly are presented as functions of the temperature in the range 100 ...1000 MeV . Using the same configurations we determine the continuum extrapolated phase diagram of QCD on the $μ-T$ plane for small to moderate chemical potentials. Two transition lines are defined with two quantities, the chiral condensate and the strange quark number susceptibility.

preprint2010arXiv

N_f=2+1 flavour equation of state

We conclude our investigation on the QCD equation of state (EoS) with 2+1 staggered flavors and one-link stout improvement. We extend our previous study [JHEP 0601:089 (2006)] by choosing even finer lattices. These new results [for details see arXiv:1007.2580] support our earlier findings. Lattices with N_t=6,8 and 10 are used, and the continuum limit is approached by checking the results at N_t=12. A Symanzik improved gauge and a stout-link improved staggered fermion action is taken; the light and strange quark masses are set to their physical values. Various observables are calculated in the temperature (T) interval of 100 to 1000~MeV. We compare our data to the equation of state obtained by the "hotQCD" collaboration.

preprint2010arXiv

QCD transition temperature: full staggered result

We conclude our investigations on the QCD cross-over transition temperatures with 2+1 staggered flavours and one-link stout improvement. We extend our previous two studies [Phys. Lett. B643 (2006) 46, JHEP 0906:088 (2009)] by choosing even finer lattices (N_t=16) and we work again with physical quark masses. These new results [for details see JHEP 1009:073,2010] support our earlier findings. We compare them with the published results of the hotQCD collaboration. All these results are confronted with the predictions of the Hadron Resonance Gas model and Chiral Perturbation Theory for temperatures below the transition region. Our results can be reproduced by using the physical spectrum. The findings of the hotQCD collaboration can be recovered only by using a distorted spectrum. This analysis provides a simple explanation for the observed discrepancy in the transition T between our and the hotQCD collaborations.

preprint2010arXiv

Recent results on QCD thermodynamics: lattice QCD versus Hadron Resonance Gas model

We present our most recent investigations on the QCD cross-over transition temperatures with 2+1 staggered flavours and one-link stout improvement [JHEP 1009:073, 2010]. We extend our previous two studies [Phys. Lett. B643 (2006) 46, JHEP 0906:088 (2009)] by choosing even finer lattices ($N_t$=16) and we work again with physical quark masses. All these results are confronted with the predictions of the Hadron Resonance Gas model and Chiral Perturbation Theory for temperatures below the transition region. Our results can be reproduced by using the physical spectrum in these analytic calculations. A comparison with the results of the hotQCD collaboration is also discussed.

preprint2010arXiv

Scaling study for 2 HEX smeared fermions: hadron and quark masses

The goal of this study is to investigate the scaling behaviour of our 2 HEX action. For this purpose, we compute the $N_f=3$ spectrum and compare the results to our 6 EXP action. We find a large scaling window up to $\sim 0.15\,\mathrm{fm}$ along with small scaling corrections at the 2%-level and full compatibility with our previous study. As a second important observable to be tested for scaling, we chose the non-perturbatively renormalized quenched strange quark mass. Here we find a fairly flat scaling with a broad scaling range up to $\simeq 0.15\,\mathrm{fm}$ and perfect agreement with the literature.

preprint2010arXiv

The QCD equation of state with dynamical quarks

The present paper concludes our investigation on the QCD equation of state with 2+1 staggered flavors and one-link stout improvement. We extend our previous study [JHEP 0601:089 (2006)] by choosing even finer lattices. Lattices with $N_t=6,8$ and 10 are used, and the continuum limit is approached by checking the results at $N_t=12$. A Symanzik improved gauge and a stout-link improved staggered fermion action is utilized. We use physical quark masses, that is, for the lightest staggered pions and kaons we fix the $m_π/f_K$ and $m_K/f_K$ ratios to their experimental values. The pressure, the interaction measure, the energy and entropy density and the speed of sound are presented as functions of the temperature in the range $100 ...1000 \textmd{MeV}$. We give estimates for the pion mass dependence and for the contribution of the charm quark. We compare our data to the equation of state obtained by the "hotQCD" collaboration.

preprint2007arXiv

The density of states method at non-zero chemical potential

We study the QCD phase diagram by first principle lattice calculations at so far unreached high densities. For this purpose we employ the density of states method. Unimproved staggered fermions, which describe four quark flavors in the continuum are used in this analysis. Though the method is quite expensive, small lattices show an indication for a triple-point connecting three different phases on the phase diagram.