Search

We present a lattice QCD calculation of the hadronic light-by-light scattering contribution to the anomalous magnetic moment of the muon using $N_f=2+1+1$ flavors of staggered quarks with masses tuned to their physical values. Our final result, in the continuum limit, reads $a_{\mu}^{\mathrm{hlbl}} = 125.5(11.6)_{\mathrm{stat}}(0.4)_{\mathrm{syst}} \times 10^{-11}$ where the first error is statistical and the second is systematic. Light, strange and charm-quark contributions are considered. In addition to the connected and leading disconnected contributions, we also include an estimate of the sub-leading disconnected diagrams. Our result is compatible with previous lattice QCD and data-driven dispersive determinations., Comment: 28 pages, 18 figures

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

The precision with which hadronic vacuum polarization (HVP) is obtained determines how accurately important observables, such as the muon anomalous magnetic moment, a_\mu, or the low-energy running of the electromagnetic coupling, \alpha, are predicted. The two most precise approaches for determining HVP are: dispersive relations combined with e+e- to hadrons cross-section data, and lattice QCD. However, the results obtained in these two approaches display significant tensions, whose origins are not understood. Here we present a framework that sheds light on this issue and, if the two approaches can be reconciled, allows them to be combined. Via this framework, we test the hypothesis that the tensions can be explained by modifying the R-ratio in different intervals of center-of-mass energy sqrt(s). As ingredients, we consider observables that have been precisely determined in both approaches. These are the leading hadronic contributions to a_\mu, to the so-called intermediate window observable and to the running of \alpha between spacelike virtualities 1GeV^2 and 10GeV^2 (for which only a preliminary lattice result exists). Our tests take into account all uncertainties and correlations, as well as uncertainties on uncertainties in the lattice results. Among our findings, the most striking is that results obtained in the two approaches can be made to agree for all three observables by modifying the \rho peak in the experimental spectrum. In particular, we find that this requires a common ~5\% increase in the contributions of the peak to each of the three observables. This finding is robust against the presence or absence of one of the constraining observables. However, such an increase is much larger than the uncertainties on the measured R-ratio. We also discuss a variety of generalizations of the methods used here, as well as the limits in the information that can be extracted..., Comment: 38 pages, 8 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

We test here our recently introduced new lattice method for the $\beta$-function defined over infinite Euclidean space-time in the continuum from scale changes generated by infinitesimal or finite steps of the renormalized gauge coupling on the gradient flow. Harlander and Neumann calculated in this scheme the three-loop approximation to the continuum $\beta$-function. Our goal is the nonperturbative lattice implementation of the scheme which we tested originally in the chiral limit of the sextet model and in multi-flavor QCD with ten and twelve flavors of massless fermions. Results are reported here in the SU(3) Yang-Mills gauge sector without dynamical fermions and in ten-flavor QCD with massless femions. The three-loop gradient flow based $\beta$-function of Harlander and Neumann is used to connect the $\Lambda_{\overline{\rm MS}}$ scale of the SU(3) Yang-Mills gauge theory with the nonperturbative flow time scale $t_0$, or the equivalent Sommer scale $r_0$. Similarly, the $\Lambda_{\overline{\rm MS}}$ scale is connected with a selected nonperturbative scale in the ten-flavor theory, a pilot study of our new lattice based nonperturbative $\beta$-function for high precision determination of the strong coupling $\alpha_s$ at the Z-boson pole in QCD with three massless fermion flavors. This goal is an important alternative to results from the finite volume based step $\beta$-function of the Alpha collaboration. Work is ongoing on direct application of the method to QCD with three massless fermion flavors., Comment: To appear in the Proceedings of the 39th Annual International Symposium on Lattice Field Theory - LATTICE2022 8-13 Aug, 2022; 9 pages, 8 figures. arXiv admin note: text overlap with arXiv:2203.15847

QCD with infinite heavy quark masses exhibits a first-order thermal transition which is driven by the spontaneous breaking of the global $\mathcal{Z}_3$ center symmetry. We analyze the corresponding order parameter, namely the Polyakov loop and its moments, and show, with a rigorous finite size scaling, that in the continuum limit the transition is of first order. We show that the use of a parallel tempering algorithm can significantly reduce the large auto-correlation times which are mainly caused by the supercritical slowing down. As a result, we calculate the transition temperature $w_0 T_c$ with per-mill precision, and the latent heat, carrying out controlled continuum and infinite volume extrapolations., Comment: 8 pages, presentation at the 39th International Symposium on Lattice Field Theory, 8th-13th August 2022, University of Bonn, Germany

Since the release of the 2015 Long Range Plan in Nuclear Physics, major events have occurred that reshaped our understanding of quantum chromodynamics (QCD) and nuclear matter at large densities, in and out of equilibrium. The US nuclear community has an opportunity to capitalize on advances in astrophysical observations and nuclear experiments and engage in an interdisciplinary effort in the theory of dense baryonic matter that connects low- and high-energy nuclear physics, astrophysics, gravitational waves physics, and data science, Comment: 70 pages, 3 figures, White Paper for the Long Range Plan for Nuclear Science

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

Contribution from the USQCD Collaboration to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021)., Comment: 27 pp. main text, 4 pp. appendices, 29 pp. references, 1 p. index

New tests are applied to two $\beta$-functions of the much-discussed BSM model with ten massless fermion flavors in the fundamental representation of the SU(3) color gauge group. The renormalization scheme of the two $\beta$-functions is defined on the gauge field gradient flow in respective finite or infinite physical volumes at zero lattice spacing. Recently published results in the ten-flavor theory led to indicators of an infrared fixed point (IRFP) in the finite-volume step $\beta$-function in the strong coupling regime of the theory arXiv:2004.00754. We analyze our substantially extended set of ten-flavor lattice ensembles at strong renormalized gauge couplings and find no evidence or hint for IRFP in the finite-volume step $\beta$-function within controlled lattice reach. We also discuss new ten-flavor tests of the recently introduced lattice definition and algorithmic implementation of the $\beta$-function defined on the gradient flow of the gauge field over infinite Euclidean space-time in the continuum. Originally we introduced this new algorithm to match finite-volume step $\beta$-functions in massless near-conformal gauge theories with the infinite-volume $\beta$-function reached in the chiral limit from small fermion mass deformations of spontaneous chiral symmetry breaking. Results from the lattice analysis of the ten-flavor infinite-volume $\beta$-function are consistent with the absence of IRFP from our step $\beta$-function based analysis. We make important contact at weak coupling in infinite volume with gradient flow based three-loop perturbation theory, serving as a first pilot study toward the long-term goal of developing alternate approach to the determination of the strong coupling $\alpha_s$ at the Z-boson pole in QCD., Comment: To appear in the Proceedings of the 38th Annual International Symposium on Lattice Field Theory - LATTICE2021 26-30 July, 2021; 18 pages, 13 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

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 $\mu_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 $\mu_B$ from first principles., Comment: 10 pages, 3 figures; Contribution to the XXXIII International (ONLINE) Workshop on High Energy Physics "Hard Problems of Hadron Physics: Non-Perturbative QCD & Related Quests"; Based on 2108.09213 [hep-lat]. arXiv admin note: substantial text overlap with arXiv:2112.02134

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

Approaches to finite baryon density lattice QCD usually suffer from uncontrolled systematic uncertainties in addition to the well-known sign problem. We test a method - sign reweighting - that works directly at finite chemical potential and is yet free from any such uncontrolled systematics: with this approach the only problem is the sign problem itself. In practice the approach involves the generation of configurations with the positive fermionic weights given by the absolute value of the real part of the quark determinant, and a reweighting by a sign. There are only two sectors, +1 and -1 and as long as the average $\left\langle \pm \right\rangle \neq 0$ (with respect to the positive weight) this discrete reweighting has no overlap problem - unlike reweighting from $\mu=0$ - and the results are reliable. We also present results based on this algorithm on the phase diagram of lattice QCD with two different actions: as a first test, we apply the method to calculate the position of the critical endpoint with unimproved staggered fermions at $N_\tau=4$; as a second application, we study the phase diagram with 2stout improved staggered fermions at $N_\tau=6$. This second one is already a reasonably fine lattice - relevant for phenomenology. We demonstrate that the method penetrates the region of the phase diagram where the Taylor and imaginary chemical potential methods lose predictive power., Comment: 9 pages, 4 figures; Contribution to the Proceedings of The 38th International Symposium on Lattice Field Theory, LATTICE2021; Based on 2004.10800 and 2108.09213

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

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 $\mu_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., Comment: 11 pages, 6 figures

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

New results are reported from tests of a low-energy effective field theory (EFT) that includes a dilaton field to describe the emergent light scalar with ${ 0^{++} }$ quantum numbers in the strongly coupled near-conformal gauge theory with a massless fermion flavor doublet in the two-index symmetric (sextet) representation of the SU(3) color gauge group. In the parlor of walking --- based on the observed light scalar, the small $\beta$-function at strong coupling, and the large anomalous scale dimension of the chiral condensate --- the dilaton EFT hypothesis is introduced to test if it explains the slowly changing nearly scale invariant physics that connects the asymptotically free UV fixed point and the far-infrared scale of chiral symmetry breaking. The characteristic dilaton EFT signatures of scale symmetry breaking are probed in this report in the small Compton wavelength limit of Goldstone bosons relative to the size of the lattice volume (p-regime) and in the limit when the Goldstone wavelength exceeds the size of the volume ($\epsilon$-regime). Random matrix theory (RMT) analysis of the dilaton EFT is applied to the lowest part of the Dirac spectrum in the $\epsilon$-regime to directly test predictions for the fundamental EFT parameters. The predictions, sensitive to the choice of the dilaton potential, were limited before to the p-regime, using extrapolations from far above the chiral limit with untested uncertainties. The dilaton EFT analysis of the $\epsilon$-regime was first suggested in \cite{Fodor:2019vmw}, with some results presented at this conference and with our continued post-conference analysis added to stimulate discussions., Comment: To appear in the Proceedings of the 37th Annual International Symposium on Lattice Field Theory - LATTICE2019 16-22 June, 2019; 13 pages, 9 figures

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 updated results for the non-perturbative $\beta$-function of SU(3) gauge theories with $N_f = 12$ or 10 massless flavors in the fundamental rep or $N_f = 2$ in the sextet rep, measured with staggered fermions. New data at finer lattice spacing and our previously introduced method, the infinitesimal $\beta$-function, strengthen the case that the $N_f = 12$ model has no infrared fixed point up to $g^2 = 7.2$. We show how underestimated cutoff dependence in one domain wall study for $N_f = 10$ has been corrected, which is now consistent with staggered results showing a monotonically increasing $\beta$-function. A consistent theme is that too small volumes can lead to apparent fixed points which vanish towards the continuum limit. We also apply the infinitesimal $\beta$-function method to the $N_f = 10$ model, finding consistent behavior with the finite-step $\beta$-function. Ongoing simulations and analysis for the sextet model confirm our previous results from weak to strong coupling with a non-zero $\beta$-function throughout, in quantitative difference to Wilson fermion simulations~\cite{Hasenfratz:2015ssa}., Comment: 7 pages, 8 figures; Proceedings of the 37th International Symposium on Lattice Field Theory - Lattice 2019, Wuhan (China)

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

The dilaton low-energy effective field theory (EFT) of an emergent light scalar is probed in the paradigm of strongly coupled near-conformal gauge theories. These studies are motivated by models which exhibit small $\beta$-functions near the conformal window (CW), perhaps with slow scale-dependent walking and a light scalar with ${ 0^{++} }$ quantum numbers. We report our results from the hypothesis of a dilaton inspired EFT analysis with two massless fermions in the two-index symmetric (sextet) representation of the SU(3) color gauge group. With important caveats in our conclusions, conformal symmetry breaking entangled with chiral symmetry breaking would drive the near-conformal infrared behavior of the theory predicting characteristic dilaton signatures of the light scalar from broken scale invariance when probed on relevant scales of fermion mass deformations. From a recently reasoned choice of the dilaton potential in the EFT description~\cite{Golterman:2016lsd} we find an unexpectedly light dilaton mass in the chiral limit at $m_d/f_\pi = 1.56(28)$, set in units of the pion decay constant $f_\pi$. Subject to further statistical and systematic tests of continued post-conference analysis, this result is significantly lower than our earlier estimates from less controlled extrapolations of the light scalar (the $\sigma$-particle) to the massless fermion limit of chiral perturbation theory. We also discuss important distinctions between the dilaton EFT analysis and the linear $\sigma$-model without dilaton signatures. For comparative reasons, we comment on dilaton tests from recent work with fermions in the fundamental representation with $n_f=8$ flavors., Comment: 14 pages, 34 figures, Proceedings of the 36th International Symposium on Lattice Field Theory (Lattice 2018), July 22-28, 2018, East Lancing, USA; elimination of some fit redundancies with minor changes in related figures

SU(3) gauge theory with $N_f$ fermions in the fundamental representation serves as a theoretical testing ground for possible infrared conformal behavior, which could play a role in BSM composite Higgs models. We use lattice simulations to study the 10-flavor model, for which it has been claimed there is an infrared fixed point in the gauge coupling $\beta$-function. Our results suggest the opposite conclusion, namely we find no $\beta$-function fixed point in the explored range, with qualitative agreement with the 5-loop $\overline{MS}$ prediction. We comment on the inconsistency between our findings and other studies., Comment: 7 pages, 6 figures; Proceedings for the 36th Annual International Symposium on Lattice Field Theory, 22-28 July 2018, Michigan State University, East Lansing, Michigan, USA

We use lattice simulations to study SU(3) gauge theory with 13 massless fermions in the fundamental representation. We present evidence that the theory is conformal with a non-zero infrared fixed point in the gauge coupling. We use a newly-developed technique to calculate the mass anomalous dimension at the fixed point via step-scaling of the mode number, allowing us to take the continuum limit and compare to perturbative predictions. We comment on the relevance of these findings to the extended search for the conformal window in the fundamental representation and in particular 12 massless flavors., Comment: 7 pages, 7 figures; Proceedings for the 36th Annual International Symposium on Lattice Field Theory, 22-28 July 2018, Michigan State University, East Lansing, Michigan, USA

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

We compute the Euclidean correlators of the stress tensor in pure $SU(3)$ Yang-Mills theory at finite temperature at zero and finite spatial momenta with lattice simulations. We perform continuum extrapolations using $N_\tau=10,12,16,20$ lattices with renormalized anisotropy 2. We use these correlators to estimate the shear viscosity of the gluon plasma in the deconfined phase. For $T=1.5T_c$ we obtain $\eta/s=0.17(2)$.

We discuss near-conformal gauge theories beyond the standard model (BSM) where interesting results on the twelve-flavor $\beta$-function of massless fermions in the fundamental representation of the SU(3) color gauge group and dilaton tests of the light scalar with two massless fermions in the two-index symmetric tensor (sextet) representation can be viewed as parts of the same BSM paradigm under investigation. We report results from high precision analysis of the twelve-flavor $\beta$-function \cite{Fodor:2016zil} refuting its published IRFP \cite{Cheng:2014jba,Hasenfratz:2016dou}. We present our objections to recent claims \cite{Hasenfratz:2017mdh,Hasenfratz:2017qyr} for non-universal behavior of staggered fermions used in our analysis. We also report our first analysis of dilaton tests of the light $0^{++}$ scalar in the sextet model and comment on related post-conference developments. The dilaton test is the main thrust of this conference contribution including presentation #405 on the $n_f=12$ $\beta$-function and presentation #260 on dilaton tests of the sextet model. They are both selected from the near-conformal BSM paradigm., Comment: 14 pages, 10 figures; Proceedings of the 35th International Symposium on Lattice Field Theory, 18-24 June 2017, Granada, Spain

As part of our ongoing lattice study of SU(3) gauge theory with two flavors of fermions in the two-index symmetric representation (the sextet model), we present the current status of the pseudoscalar particle spectrum. We use a mixed action approach based on the gradient flow to control lattice artifacts, allowing a simultaneous extrapolation to the chiral and continuum limits. We find strong evidence that the pseudoscalar is a Goldstone boson state, with spontaneously broken chiral symmetry and a non-zero Goldstone decay constant in the chiral limit. In agreement with our study of the gauge coupling $\beta$ function, we find the sextet model appears to be a near-conformal gauge theory and serves as a prototype of the composite Higgs BSM template., Comment: 8 pages, 6 figures, Proceeding of the 35th International Symposium on Lattice Field Theory, 18-24 June 2017, Granada, Spain

We describe a new method to determine non-perturbatively the beta function of a gauge theory using lattice simulations in the p-regime of the theory. This complements alternative measurements of the beta function working directly at zero fermion mass and bridges the gap between the weak coupling perturbative regime and the strong coupling regime relevant to the mass spectrum of the theory. We apply this method to ${\mathrm {SU(3)} }$ gauge theory with two fermion flavors in the 2-index symmetric (sextet) representation. We find that the beta function is small but non-zero at the renormalized coupling value $g^2 = 6.7$, consistent with our previous independent investigation using simulations directly at zero fermion mass. The model continues to be a very interesting explicit realization of the near-conformal composite Higgs paradigm which could be relevant for Beyond Standard Model phenomenology., Comment: 8 pages, 7 figures; Proceedings of the 35th International Symposium on Lattice Field Theory, 18-24 June 2017, Granada, Spain

Results are reported for the beta-function of weakly coupled conformal gauge theories on the lattice, SU(3) with Nf=14 fundamental and Nf=3 sextet fermions. The models are chosen to be close to the upper end of the conformal window where perturbation theory is reliable hence a fixed point is expected. The study serves as a test of how well lattice methods perform in the weakly coupled conformal cases. We also comment on the 5-loop beta-function of two models close to the lower end of the conformal window, SU(3) with Nf=12 fundamental and Nf=2 sextet fermions., Comment: 8 pages, 20 figures, Proceedings of the 35th International Symposium on Lattice Field Theory, 18-24 June 2017, Granada, Spain

We report new results from high precision analysis of an important BSM gauge theory with twelve massless fermion flavors in the fundamental representation of the SU(3) color gauge group. The range of the renormalized gauge coupling is extended from our earlier work {Fodor:2016zil} to probe the existence of an infrared fixed point (IRFP) in the $\beta$-function reported at two different locations, originally in {Cheng:2014jba} and at a new location in {Hasenfratz:2016dou}. We find no evidence for the IRFP of the $\beta$-function in the extended range of the renormalized gauge coupling, in disagreement with {Cheng:2014jba,Hasenfratz:2016dou}. New arguments to guard the existence of the IRFP remain unconvincing {Hasenfratz:2017mdh}, including recent claims of an IRFP with ten massless fermion flavors {Chiu:2016uui,Chiu:2017kza} which we also rule out. Predictions of the recently completed 5-loop QCD $\beta$-function for general flavor number are discussed in this context., Comment: 7 pages, 9 figures

The first principle lattice QCD methods allow to calculate the thermodynamic observables at finite temperature and imaginary chemical potential. These can be compared to the predictions of various phenomenological models. We argue that Fourier coefficients with respect to imaginary baryochemical potential are sensitive to modeling of baryonic interactions. As a first application of this sensitivity, we consider the hadron resonance gas (HRG) model with repulsive baryonic interactions, which are modeled by means of the excluded volume correction. The Fourier coefficients of the imaginary part of the net-baryon density at imaginary baryochemical potential -- corresponding to the fugacity or virial expansion at real chemical potential -- are calculated within this model, and compared with the $N_t = 12$ lattice data. The lattice QCD behavior of the first four Fourier coefficients up to $T \simeq 185$ MeV is described fairly well by an interacting HRG with a single baryon-baryon eigenvolume interaction parameter $b \simeq 1$ fm$^3$, while the available lattice data on the difference $\chi_2^B - \chi_4^B$ of baryon number susceptibilities is reproduced up to $T \simeq 175$ MeV., Comment: 11 pages, 4 figures, new Fig. 2 depicting $\chi_2^B - \chi_4^B$, version accepted for publication in Physics Letters B

Fluctuations of conserved charges allow to study the chemical composition of hadronic matter. A comparison between lattice simulations and the Hadron Resonance Gas (HRG) model suggested the existence of missing strange resonances. To clarify this issue we calculate the partial pressures of mesons and baryons with different strangeness quantum numbers using lattice simulations in the confined phase of QCD. In order to make this calculation feasible, we perform simulations at imaginary strangeness chemical potentials. We systematically study the effect of different hadronic spectra on thermodynamic observables in the HRG model and compare to lattice QCD results. We show that, for each hadronic sector, the well established states are not enough in order to have agreement with the lattice results. Additional states, either listed in the Particle Data Group booklet (PDG) but not well established, or predicted by the Quark Model (QM), are necessary in order to reproduce the lattice data. For mesons, it appears that the PDG and the quark model do not list enough strange mesons, or that, in this sector, interactions beyond those included in the HRG model are needed to reproduce the lattice QCD results., Comment: 9 pages, 9 figures

We report new results on the conformal properties of an important strongly coupled gauge theory, a building block of composite Higgs models beyond the Standard Model. With twelve massless fermions in the fundamental representation of the SU(3) color gauge group, an infrared fixed point of the $\beta$-function was recently reported in the theory (Cheng:2014jba) with uncertainty in the location of the critical gauge coupling inside the narrow $[ 6.0

We analyze three sets of gauge ensembles in our extended physics program of a particularly important BSM gauge theory with a fermion doublet in the two-index symmetric (sextet) representation of the SU(3) BSM color gauge group. Our investigations include chiral symmetry breaking $\rm{(\chi SB)}$ in the p-regime and $\epsilon$-regime, the mass of the composite ${\rm 0^{++}}$ scalar, resonance spectroscopy, new physics from gauge anomaly constraints, and the role of stable sextet BSM baryons with Electroweak interactions in dark matter searches. Important new goals include studies of the ${\rm 0^{++}}$ scalar entangled with Goldstone dynamics in the p-regime and the $\epsilon$-regime, the resonance spectrum with particular attention to emerging LHC signals, like recent hints for diphoton excess at 750 GeV or diboson anomalies in the 2 TeV range. All results reported here are preliminary before journal publication including some post-conference material for the discussion., Comment: 20 pages, Proceedings of the 33rd International Symposium on Lattice Field Theory (LATTICE2015), Kobe, Japan, 14-18 July 2015. arXiv admin note: text overlap with arXiv:1502.00028; text overlap with arXiv:1211.1083 by other authors

We report tests and results from a new approach to the spectral density and the mode number distribution of the Dirac operator in lattice gauge theories. The algorithm generates the spectral density of the lattice Dirac operator as a continuous function over all scales of the complete eigenvalue spectrum. This is distinct from an earlier method where the integrated spectral density (mode number) was calculated efficiently for some preselected fixed range of the integration. The new algorithm allows global studies like the chiral condensate from the Dirac spectrum at any scale including the cutoff-dependent IR and UV range of the spectrum. Physics applications include the scale-dependent mass anomalous dimension, spectral representation of composite fermion operators, and the crossover transition from the $\epsilon$-regime of Random Matrix Theory to the p-regime in chiral perturbation theory. We present thorough tests of the algorithm in the 2-flavor sextet SU(3) gauge theory that we continue to pursue for its potential as a minimal realization of the composite Higgs scenario., Comment: 7 pages, Proceedings of the 33rd International Symposium on Lattice Field Theory (LATTICE2015), Kobe, Japan, 14-18 July 2015

We present a calculation of the leptonic decay-constant ratio $f_K/f_\pi$ in 2+1 flavor QCD. Our data set includes five lattice spacings and pion masses reaching down below the physical one. Special emphasis is placed on a careful study of all systematic uncertainties, especially the continuum extrapolation. Our result is perfectly compatible with the first-row unitarity constraint of the Standard Model., Comment: 19 pages, 6 figures, 3 tables; v2: added supplementary analysis, version published in Phys. Rev. D

The Electroweak interactions of a strongly coupled gauge theory are discussed with outlook beyond the Standard Model (BSM) under global and gauge anomaly constraints. The theory is built on a minimal massless fermion doublet of the SU(2) BSM flavor group (bsm-flavor) with a confining gauge force at the TeV scale in the two-index symmetric (sextet) representation of the BSM SU(3) color gauge group (bsm-color). The intriguing possibility of near-conformal sextet gauge dynamics could lead to the minimal realization of the composite Higgs mechanism with a light $0^{++}$ scalar, far separated from strongly coupled resonances of the confining gauge force in the 2-3 TeV range, distinct from Higgsless Technicolor. In previous publications we have presented results for the meson spectrum of the theory, including the light composite scalar, perhaps the emergent Higgs impostor. Here we discuss the critically important role of the baryon spectrum in the sextet model investigating its compatibility with what we know about thermal evolution of the early Universe including its galactic and terrestrial relics. For an important application, we report the first numerical results on the baryon spectrum of this theory from non-perturbative lattice simulations with baryon correlators in the staggered fermion implementation of the strongly coupled gauge sector. The quantum numbers of composite baryons and their spectroscopy from lattice simulations are required inputs for exploring dark matter contributions of the sextet BSM model, as outlined for future work., Comment: 23 pages, 3 figures