Researcher profile

Sandor D. Katz

Sandor D. Katz contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
3topics
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

8 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.