Search

The QCD cross-over line in the temperature ($T$) -- baryo-chemical potential ($\mu_B$) plane has been computed by several lattice groups by calculating the chiral order parameter and its susceptibility at finite values of $\mu_B$. In this work we focus on the deconfinement aspect of the transition between hadronic and Quark Gluon Plasma (QGP) phases. We define the deconfinement temperature as the peak position of the static quark entropy ($S_Q(T,\mu_B)$) in $T$, which is based on the renormalized Polyakov loop. We extrapolate $S_Q(T,\mu_B)$ based on high statistics finite temperature ensembles on a $16^3\times 8$ lattice to finite density by means of a Taylor expansion to eighth order in $\mu_B$ (NNNLO) along the strangeness neutral line. For the simulations the 4HEX staggered action was used with 2+1 flavors at physical quark masses. In this setup the phase diagram can be drawn up to unprecedentedly high chemical potentials. Our results for the deconfinement temperature are in rough agreement with phenomenological estimates of the freeze-out curve in relativistic heavy ion collisions. In addition, we study the width of the deconfinement crossover. We show that up to $\mu_B \approx 400$ MeV, the deconfinement transition gets broader at higher densities, disfavoring the existence of a deconfinement critical endpoint in this range. Finally, we examine the transition line without the strangeness neutrality condition and observe a hint for the narrowing of the crossover towards large $\mu_B$., Comment: 14 pages, 8 figures

We present a lattice calculation of the $\eta$ and $\eta^{\prime}$ meson masses at the physical point and in the continuum limit, based on $N_f = 2+1+1$ flavors of rooted staggered quarks. Our analysis includes gauge ensembles at the physical pion and kaon masses spread over six lattice spacings in the range [0.064-0.1315]~fm. Our main results read $m_{\eta} = 543.5(5.6)~$MeV and $m_{\eta^{\prime}} = 986(38)~$MeV, consistent with the experimental values. This is an important numerical test that supports the validity of the fourth root procedure used in the staggered quark formalism. This calculation was the first step towards extracting the pseudoscalar transition form factors of the $\eta$ and $\eta^{\prime}$ mesons that play a crucial role in the hadronic light-by-light contribution to the muon $g-2$., Comment: 19 pages, 12 figures

The crossover from hadronic to quark matter is understood to be both a deconfinement as well as a chiral symmetry restoring transition. Here, we study observables related to both aspects using lattice simulations: the Polyakov loop and its derivatives and the chiral condensate and its derivatives. At zero baryochemical potential, and infinite volume, the chiral and deconfinement crossover temperatures almost agree. However, chiral and deconfinement related observables have a qualitatively different chemical potential and volume dependence. In general, deconfinement related observables have a milder volume dependence. Furthermore, while the deconfinement transition appears to get broader with increasing $\mu_B$, the width as well as the strength of the chiral transition is approximately constant. Our results are based on simulations at zero and imaginary chemical potentials using 4stout-improved staggered fermions with $N_\tau=12$ time-slices and physical quark masses., Comment: 10 pages, 8 figures

An efficient way to study the QCD phase diagram at small finite density is to extrapolate thermodynamical observables from imaginary chemical potential. The phase diagram features a crossover line starting from the transition temperature already determined at zero chemical potential. In this work we focus on the Taylor expansion of this line up to $\mu^4$ contributions. We present the continuum extrapolation of the crossover temperature based on different observables at several lattice spacings., Comment: Proceedings to Quark Matter Conference 2019

The effect of a finite volume presents itself both in heavy ion experiments as well as in recent model calculations. The magnitude is sensitive to the proximity of a nearby critical point. We calculate the finite volume effects at finite temperature in continuum QCD using lattice simulations and set the focus on the vicinity of the chiral crossover. We investigate the impact of finite volumes at zero and small chemical potentials on the QCD transition through the chiral observables.

Fluctuations play a key role in the study of QCD phases. Lattice QCD is a valuable tool to calculate them, but going to high orders is challenging. Up to the fourth order, continuum results are available since 2015. We present the first continuum results for sixth order baryon fluctuations for temperatures between $T=130 - 200$ MeV, and eighth order at $T=145$ MeV in a fixed volume. We show that for $T \leq 145$ MeV, relevant for criticality search, finite volume effects are under control. Our results are in sharp contrast with well known results in the literature obtained at finite lattice spacing., Comment: 5 pages, 2 figures (main text) + 5 pages, 7 figures (supplemental material)

Research on the QCD phase diagram with lattice field theory methods is dominated by the use of rooted staggered fermions, as they are the computationally cheapest discretization available. We show that rooted staggered fermions at a nonzero baryochemical potential $\mu_B$ predict a sharp rise in the baryon density at low temperatures and $\mu_B \gtrsim 3 m_\pi/2$, where $m_\pi$ is the Goldstone pion mass. We elucidate the nature of the non-analyticity behind this sharp rise in the density by a comparison of reweighting results with a Taylor expansion of high order. While at first sight this non-analytic behavior becomes apparent at the same position where the pion condensation transition takes place in the phase-quenched theory, the nature of the non-analyticity in the two theories appears to be quite different: While at nonzero isospin density the data are consistent with a genuine thermodynamic (branch-point) singularity, the results at nonzero baryon density point to an essential singularity at $\mu_B=0$. The effect is absent for four flavors of degenerate quarks, where rooting is not used. For the two-flavor case, we show numerical evidence that the magnitude of the effect diminishes on finer lattices. We discuss the implications of this technical complication on future studies of the QCD phase diagram., Comment: 12 pages, 6 figures

In this paper we present a first ab-initio calculation of the $\pi^0$, $\eta$ and $\eta^{\prime}$ transition form factors performed with physical light-quark masses. We provide a complete parametrization of the form factors that includes both single and double-virtual kinematics. Our results are compared with experimental measurements of the form factors in the space-like region and with the measured two-photon decay widths. In a second step, our parametrizations of the transition form factors are used to compute the dominant pseudoscalar-pole contributions to the hadronic light-by-light scattering in the muon $g-2$. Our final result reads $a_{\mu}^{\rm hlbl, ps-pole} = (85.1 \pm 5.2) \times 10^{-11}$. Although the pion-pole is dominant, we confirm that, together, the $\eta$ and $\eta^{\prime}$ provide roughly half of its contribution., Comment: 29 pages, 16 figures. Typos corrected. Add figure 6

Phase transitions in a non-perturbative regime can be studied by ab initio Lattice Field Theory methods. The status and future research directions for LFT investigations of Quantum Chromo-Dynamics under extreme conditions are reviewed, including properties of hadrons and of the hypothesized QCD axion as inferred from QCD topology in different phases. We discuss phase transitions in strong interactions in an extended parameter space, and the possibility of model building for Dark Matter and Electro-Weak Symmetry Breaking. Methodological challenges are addressed as well, including new developments in Artificial Intelligence geared towards the identification of different phases and transitions., Comment: 94 pages, 23 figures, GGI workshop "Phase Transitions in Particle Physics" review - v2: minor revisions and typo corrected, matches accepted version on PPNP

We report on our progress toward the computation of the $\pi^0$, $\eta$ and $\eta^{\prime}$ transition form factors using staggered quarks on $N_f=2+1+1$ gauge ensembles generated by the Budapest-Marseille-Wuppertal collaboration. These form factors are essential ingredients to evaluate the pseudoscalar-pole contributions to the hadronic light-by-light scattering in the muon $g-2$. Preliminary results for the pseudoscalar-pole contributions are presented, at finite lattice spacing, for all three light mesons., Comment: 10 pages, 7 figures, contribution to the 39th International Symposium on Lattice Field Theory

The equation of state of the quark gluon plasma is a key ingredient of heavy ion phenomenology. In addition to the traditional Taylor method, several novel approximation schemes have been proposed with the aim of calculating it at finite baryon density. In order to gain a pragmatic understanding of the limits of these schemes, we compare them to direct results at $\mu_B>0$, using reweighting techniques free from an overlap problem. We use 2stout improved staggered fermions with 8 time-slices and cover the entire RHIC BES range in the baryochemical potential, up to $\mu_B/T=3$., Comment: 7 pages, 3 figures

We calculate a resummed equation of state with lattice QCD simulations at imaginary chemical potentials. This work presents a generalization of the scheme introduced in 2102.06660 to the case of non-zero $\mu_S$, focusing on the line of strangeness neutrality. We present results up to $\mu_B/T \leq 3.5$ on the strangeness neutral line $\left\langle S \right\rangle = 0$ in the temperature range $130 \rm{MeV} \leq T \leq 280 \rm{MeV}$. We also extrapolate the finite baryon density equation of state to small non-zero values of the strangeness-to-baryon ratio $R=\left\langle S \right\rangle / \left\langle B \right\rangle$. We perform a continuum extrapolation using lattice simulations of the 4stout-improved staggered action with 8, 10, 12 and 16 timeslices., Comment: 15 pages, 12 figures; v2: contains ancillary files with tabulated data

In recent years there has been much progress on the investigation of the QCD phase diagram with lattice QCD. This talk will focus on the developments in the last few years. Especially the addition of external influences and extended ranges of $T$ and $\mu$ yield an increasing number of interesting results, a subset of which will be discussed. Many of these conditions are important for the understanding of both the QCD transition in the early universe and heavy ion collision experiments which are conducted for example at the LHC and RHIC. This offers many exciting opportunities for comparisons between theory and experiment., Comment: conference proceedings for Lattice 2021, this article draws from arXiv:2010.15503

We present preliminary results from our calculation of the pseudoscalar transition form factors of the $\eta$ and $\eta'$ mesons using staggered quarks on $N_f=2+1+1$ gauge ensembles generated by the Budapest-Marseille-Wuppertal collaboration. These transition form factors are an important input for the hadronic light-by-light contribution to the muon $(g-2)$. We first elaborate on the extraction of the masses of the $\eta$ and $\eta'$ mesons, that mix under the dynamics of QCD. Thereafter, we show our preliminary results for the pseudoscalar transition form factors, focusing on the $\eta$ meson in the absence of mixing., Comment: 9 pages, 7 figures, The 38th International Symposium on Lattice Field Theory, LATTICE2021 26th-30th July, 2021, Zoom/Gather@Massachusetts Institute of Technology

QCD with heavy dynamical quarks exhibits a first order thermal transition which is driven by the spontaneous breaking of the global $\mathcal{Z}_3$ center symmetry. Decreasing the quark masses weakens the transition until the corresponding latent heat vanishes at the critical mass. We explore the heavy mass region with three flavors of staggered quarks and analyze the Polyakov loop and its moments in a finite volume scaling study. We calculate the heavy critical mass in the three flavor theory in the infinite volume limit for $N_t=8$., Comment: 8 pages, presentation at the 38th International Symposium on Lattice Field Theory, 26th-30th July 2021, Massachusetts Institute of Technology, USA

The hadron resonance gas (HRG) model and its extensions are often used to describe the hadronic phase of strongly interacting matter. In our work we use lattice-QCD simulations with temporal extents of $N_\tau=8,10$ and $12$ to quantify corrections to the ideal HRG. Firstly, we determine a number of subleading fugacity expansion coefficients of the QCD free energy via a two-dimensional scan on the imaginary baryon number chemical potential ($\mu_B$) - strangeness chemical potential ($\mu_S$) plane. Using the aforementioned coefficients, we also extrapolate ratios of baryon number and strangeness fluctuations and correlations to finite chemical potentials via a truncated fugacity expansion. Our results extrapolated along the crossover line $T_\mathrm{c}(\mu_B)$ at strangeness neutrality are able to reproduce trends of experimental net-proton fluctuations measured by the STAR Collaboration., Comment: 10 pages, 4 figures, Contribution to the 38th International Symposium on Lattice Field Theory, LATTICE2021 26th-30th July, 2021

The equation of state of Quantum Chromodynamics (QCD) at finite density is currently known only in a limited range in the baryon chemical potential $\mu_B$. This is due to fundamental shortcomings of traditional methods such as Taylor expansion around $\mu_B=0$. In this contribution, we present an alternative scheme that displays substantially improved convergence over the Taylor expansion method. We calculate the alternative expansion coefficients in the continuum, and show our results for the thermodynamic observables up to $\mu_B/T\le3.5$., Comment: 8 pages, presentation at the 38th International Symposium on Lattice Field Theory, 26th-30th July 2021, Massachusetts Institute of Technology, USA

We present a new formalism for numerically treating the semiclassical gravitational collapse of a scalar quantum field in the radially symmetric case. Our formalism is time reversal invariant and the evolution of the scalar fields is unitary. We present some first results in the angular momentum $l=0$ approximation for an initially coherent state of a massless field and briefly discuss future prospects., Comment: 9 pages, 4 figures, talk presented at the 38th International Symposium on Lattice Field Theory - LATTICE2021, 26th-30th, July, 2021, Zoom/Gather@Massachusetts Institute of Technology

The hadron resonance gas (HRG) model is often believed to correctly describe the confined phase of QCD. This assumption is the basis of many phenomenological works on QCD thermodynamics and of the analysis of hadron yields in relativistic heavy ion collisions. We use first-principle lattice simulations to calculate corrections to the ideal HRG. Namely, we determine the sub-leading fugacity expansion coefficients of the grand canonical free energy, receiving contributions from processes like kaon-kaon or baryon-baryon scattering. We achieve this goal by performing a two dimensional scan on the imaginary baryon number chemical potential ($\mu_B$) - strangeness chemical potential ($\mu_S$) plane, where the fugacity expansion coefficients become Fourier coefficients. We carry out a continuum limit estimation of these coefficients by performing lattice simulations with temporal extents of $N_\tau=8,10,12$ using the 4stout-improved staggered action. We then use the truncated fugacity expansion to extrapolate ratios of baryon number and strangeness fluctuations and correlations to finite chemical potentials. Evaluating the fugacity expansion along the crossover line, we reproduce the trend seen in the experimental data on net-proton fluctuations by the STAR collaboration., Comment: 13 pages, 6 figures

In recent years there has been much progress on the investigation of the QCD phase diagram with lattice QCD simulations. In this review I focus on the developments in the last two years. Especially the addition of external influences or new parameter ranges yield an increasing number of interesting results. I discuss the progress for small, finite densities from both analytical continuation and Complex Langevin simulations, for heavy quark bound states (quarkonium), the dependence on the quark masses (Columbia plot) and the influence of a magnetic field. Many of these conditions are relevant for the understanding of both the QCD transition in the early universe and heavy ion collision experiments which are conducted for example at the LHC and RHIC., Comment: 59 pages, 31 figures

We present a method to study the semiclassical gravitational collapse of a radially symmetric scalar quantum field in a coherent initial state. The formalism utilizes a Fock space basis in the initial metric, is unitary and time reversal invariant up to numerical precision. It maintains exact compatibility of the metric with the expectation values of the energy momentum tensor in the scalar field coherent state throughout the entire time evolution. We find a simple criterion for the smallness of discretization effects, which is violated when a horizon forms. As a first example, we study the collapse of a specific state in the angular momentum $l=0$ approximation. Outside the simulated volume it produces a Schwarzschild metric with $r_s \sim 3.5 \ell_p$. We see behaviour that is compatible with the onset of horizon formation both in the semiclassical and corresponding classical cases in a regime where we see no evidence for large discretization artefacts. In our example setting, we see that quantum effects accelerate the possible horizon formation and move it radially outward. We find that this effect is robust against variations of the radial resolution, the time step, the volume, the initial position and shape of the inmoving state, the vacuum subtraction, the discretization of the time evolution operator and the integration scheme of the metric. We briefly discuss potential improvements of the method and the possibility of applying it to black hole evaporation. We also briefly touch on the extension of our formalism to higher angular momenta, but leave the details and numerics for a forthcoming publication., Comment: v2: 54 pages, 32 figures. Version to appear in PRD

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 $\mu_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 $\mu_B$. We obtain the parameters $\kappa_2=0.0153(18)$ and $\kappa_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 $\mu_B$. We see no sign of criticality in the explored range., Comment: 5 pages main text + 7 pages supplementary material

We present cross-correlators of QCD conserved charges at $\mu_B=0$ from lattice simulations and perform a Hadron Resonance Gas (HRG) model analysis to break down the hadronic contributions to these correlators. We construct a suitable hadronic proxy for the ratio $-\chi_{11}^{BS}/\chi_2^S$ and discuss the dependence on the chemical potential and experimental cuts. We then perform a comparison to preliminary STAR results and comment on a possible direct comparison of lattice and experiment., Comment: 6 pages, 3 figures, contribution to the proceedings from the 18th International Conference on Strangeness in Quark Matter (SQM 2019)

Like fluctuations, non-diagonal correlators of conserved charges provide a tool for the study of chemical freeze-out in heavy ion collisions. They can be calculated in thermal equilibrium using lattice simulations, and be connected to moments of event-by-event net-particle multiplicity distributions. We calculate them from continuum extrapolated lattice simulations at $\mu_B=0$, and present a finite-$\mu_B$ extrapolation, comparing two different methods. In order to relate the grand canonical observables to the experimentally available net-particle fluctuations and correlations, we perform a Hadron Resonance Gas (HRG) model analysis, which allows us to completely break down the contributions from different hadrons. We then construct suitable hadronic proxies for fluctuations ratios, and study their behavior at finite chemical potentials. We also study the effect of introducing acceptance cuts, and argue that the small dependence of certain ratios on the latter allows for a direct comparison with lattice QCD results, provided that the same cuts are applied to all hadronic species. Finally, we perform a comparison for the constructed quantities for experimentally available measurements from the STAR Collaboration. Thus, we estimate the chemical freeze-out temperature to 165 MeV using a strangeness-related proxy. This is a rather high temperature for the use of the Hadron Resonance Gas, thus, further lattice studies are necessary to provide first principle results at intermediate $\mu_B$., Comment: 20 pages, 17 figures, 1 table, version published in Phys.Rev.D

We discuss the reliability of available methods to constrain the location of the QCD critical endpoint with lattice simulations. In particular we calculate the baryon fluctuations up to $\chi^B_8$ using simulations at imaginary chemical potentials. We argue that they contain no hint of criticality., Comment: XXVIIth International Conference on Ultrarelativistic Nucleus-Nucleus Collisions (Quark Matter 2018)

We calculate several diagonal and non-diagonal fluctuations of conserved charges in a system of 2+1+1 quark flavors with physical masses, on a lattice with size $48^3\times12$. Higher order fluctuations at $\mu_B=0$ are obtained as derivatives of the lower order ones, simulated at imaginary chemical potential. From these correlations and fluctuations we construct ratios of net-baryon number cumulants as functions of temperature and chemical potential, which satisfy the experimental conditions of strangeness neutrality and proton/baryon ratio. Our results qualitatively explain the behavior of the measured cumulant ratios by the STAR collaboration., Comment: 22 pages, 12 figures