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

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

In this study we employ staggered fermions to calculate the two-pion taste singlet states at rest. Leveraging the Clebsch-Gordan coefficients of the symmetry group associated with staggered fermions, we effectively compute the $\pi\pi$ contributions to the resting $\rho$-meson correlator. To discern the distinct energy states involved, we adopt a generalized eigenvalue problem-solving approach. This work will provide insight into the important role played by the two-pion contribution to the anomalous magnetic moment of the muon. In this paper we present our group theoretic considerations and preliminary results on the contribution of two-pion states to the rho meson correlation function., Comment: 8 pages, 2 tables, 5 figures

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)

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

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

In previous work, arXiv:1905.01909, we have calculated the $m_\varrho / f_\pi$ ratio in the chiral and continuum limit for $SU(3)$ gauge theory coupled to $N_f = 2,3,4,5,6$ fermions in the fundamental representation. The main result was that this ratio displays no statistically significant $N_f$-dependence. In the present work we continue the study of the $N_f$-dependence by extending the simulations to $N_f = 7, 8, 9, 10$. Along the way we also study in detail the $N_f$-dependence of finite volume effects on low energy observables and a particular translational symmetry breaking unphysical, lattice artefact phase specific to staggered fermions., Comment: 14 pages, 16 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

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

A common problem in lattice QCD simulations on the torus is the extremely long autocorrelation time of the topological charge, when one approaches the continuum limit. The reason is the suppressed tunneling between topological sectors. The problem can be circumvented by replacing the torus with a different manifold, so that the connectivity of the configuration space is changed. This can be achieved by using open boundary conditions on the fields, as proposed earlier. It has the side effect of breaking translational invariance strongly. Here we propose to use a non-orientable manifold, and show how to define and simulate lattice QCD on it. We demonstrate in quenched simulations that this leads to a drastic reduction of the autocorrelation time. A feature of the new proposal is, that translational invariance is preserved up to exponentially small corrections. A Dirac-fermion on a non-orientable manifold poses a challenge to numerical simulations: the fermion determinant becomes complex. We propose two approaches to circumvent this problem., Comment: 9 pages, 8 figures; v2: matches accepted version

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, Comment: 7 pages, 4 figures. Talk presented at the 33rd International Symposium on Lattice Field Theory (Lattice 2015), 14-18 July 2015, Kobe International Conference Center, Kobe, Japan

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., Comment: 19 pages, 13 figures, published version

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

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$., Comment: Conference proceedings: The 32nd International Symposium on Lattice Field Theory - Lattice 2014, June 23-28, 2014, Columbia University, New York, New York

Finite temperature charmonium spectral functions in the pseudoscalar(PS) and vector(V) channels are studied in lattice QCD with 2+1 flavours of dynamical Wilson quarks, on fine isotropic lattices (with a lattice spacing of 0.057fm), with a non-physical pion mass of 545MeV. The highest temperature studied is approximately 1.4Tc. Up to this temperature no significant variation of the spectral function is seen in the PS channel. The V 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 omega > 0. These results are in accord with previous calculations using the quenched approximation., Comment: Conference proceedings: The 32nd International Symposium on Lattice Field Theory - Lattice 2014 June 23-28, 2014 Columbia University, New York, New York This conference contribution draws heavily from the paper: arXiv:1401.5940 [hep-lat]

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_{\pi} \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 \omega>0. These results are in accord with previous calculations using the quenched approximation., Comment: 17 pages, 9 figures, 2 tables

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. In this study, all 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$. A Symanzik improved gauge and stout-link improved staggered fermion action is used. Our results are available online [ancillary file to arXiv:1309.5258]., Comment: Conference proceedings, 7 pages, 4 figures. Talk presented at 31st International Symposium on Lattice Field Theory (LATTICE 2013), July 29 - August 3, 2013, Mainz, Germany

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., Comment: 10 pages, 9 figures

By fitting pion masses and decay constants from 2+1 flavor staggered lattice simulations to the predictions of NLO and NNLO SU(2) chiral perturbation theory we determine the low-energy constants l_3 and l_4. The lattice ensembles were generated by the Wuppertal-Budapest collaboration and cover pion masses in the range of 135 to 435 MeV and lattice scales between 0.7 and 2.0 GeV. By choosing a suitable scaling trajectory, we were able to demonstrate that precise and stable results for the LECs can be obtained from continuum ChPT to NLO. The pion masses available in this work also allow us to study the applicability of using ChPT to extrapolate from higher masses to the physical pion mass., Comment: 8 pages, 8 figures, 1 table, talk presented at Xth Quark Confinement and the Hadron Spectrum, Munich, October 2012

We extract the next-to-leading-order low-energy constants \bar\ell_3 and \bar\ell_4 of SU(2) chiral perturbation theory, based on precise lattice data for the pion mass and decay constant on ensembles generated by the Wuppertal-Budapest Collaboration for QCD thermodynamics. These ensembles feature 2+1 flavors of two-fold stout-smeared dynamical staggered fermions combined with Symanzik glue, with pion masses varying from 135 to 435 MeV, lattice scales between 0.7 and 2.0 GeV, while m_s is kept fixed at its physical value. Moderate taste splittings and the scale being set through the pion decay constant allow us to restrict ourselves to the taste pseudoscalar state and to use formulas from continuum chiral perturbation theory. Finally, by dropping the data points near 135 MeV from the fits, we can explore the range of pion masses that is needed in SU(2) chiral perturbation theory to reliably extrapolate to the physical point., Comment: 40 pages, 22 figures, 3 tables; v2: expanded discussion, matches published version

We present an update on the QCD equation of state of the Wuppertal-Budapest Collaboration, extending our previous studies [JHEP 0601 (2006) 089, JHEP 1011 (2010) 077]. A Symanzik improved gauge and a stout-link improved staggered fermion action is utilized. We discuss partial quenching and present preliminary results for the fully dynamical charmed equation of state., Comment: Talk presented at the XXIX International Symposium on Lattice Field Theory, July 10-16, 2011, Lake Tahoe, California

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., Comment: 7 pages, 2 figures, The XXIX International Symposium on Lattice Field Theory, July 10-16, 2011, Squaw Valley, Lake Tahoe, California, USA

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., Comment: 7 pages, contribution to Lattice 2011, Squaw Valley, Lake Tahoe, California, July 2011

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., Comment: 7 pages, Talk presented at the XXIX International Symposium on Lattice Field Theory (Lattice 2011), July 10-16, 2011, Squaw Valley, Lake Tahoe, California, USA

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 $\mu-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., Comment: 4 pages, 2 figures, Proceedings for Quark Matter 2011

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., Comment: 4 pages, 3 figures, Proceedings for Quark Matter 2011

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., Comment: 4 pages, 3 figures, Presented at the 4th International Conference Hard Probes 2010

The dual condensate is a new QCD phase transition order parameter, which connnects confinement and chiral symmetry breaking as different mass limits. We discuss the relation between the fermion spectrum at general boundary conditions and the dual condensate and show numerical results for the latter from unquenched SU(3) lattice configurations., Comment: Talk given by Bo Zhang on 'Quark Confinement and the Hadron Spectrum IX', Madrid, Aug.2010

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., Comment: presented at the XXVIII. International Symposium on Lattice Field Theory, June 14-19,2010, Villasimius, Sardinia Italy

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., Comment: presented at the XXVIII. International Symposium on Lattice Field Theory, June 14-19,2010, Villasimius, Sardinia Italy

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., Comment: PoS for the XXVIII International Symposium on Lattice Field Theory, Lattice2010, 7 pages, 4 figures

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_\pi/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., Comment: Minor changes: Figure 1 added; Figure 15, Figure 17 and Table 5 changed. Accepted for publication in JHEP

The present paper concludes 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. The new results on this broad cross-over are in complete agreement with our earlier ones. We compare our findings 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 in these analytic calculations. The findings of the hotQCD collaboration can be recovered by using a distorted spectrum which takes into account lattice discretization artifacts and heavier than physical quark masses. This analysis provides a simple explanation for the observed discrepancy in the transition temperatures between our and the hotQCD collaborations., Comment: 25 pages, 10 figures and 3 tables

We determined the curvature of the phase transition line in the mu-T plane using a Taylor expansion in mu. The Polyakov loop and the strange quark number susceptibility were measured to locate the pseudocritical line. The analysis was carried out on Nt=4,6,8,10 lattices generated with a Symanzik improved gauge and stout-link improved (2+1) flavour staggered fermion action using physical quark masses., Comment: 7 pages, 5 figures; format adjusted

A quantitative comparison between the finite temperature behaviour of the staggered and Wilson fermion formulations are performed. The comparison is based on a physical quantity that is expected to be quite sensitive to the fermionic features of the action. For that purpose we use the height of the peak for $d\chi_s/dT$, where $\chi_s$ is the quark number susceptibility., Comment: 6 pages. Talk presented at Lattice 2007

The speed, bandwidth and cost characteristics of today's PC graphics cards make them an attractive target as general purpose computational platforms. High performance can be achieved also for lattice simulations but the actual implementation can be cumbersome. This paper outlines the architecture and programming model of modern graphics cards for the lattice practitioner with the goal of exploiting these chips for Monte Carlo simulations. Sample code is also given., Comment: 11 pages, 2 figures, sample code, new benchmarks added