Source author record

Stefan Krieg

Stefan Krieg appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

30works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

30 published item(s)

preprint2020arXiv

The Semimetal-Mott Insulator Quantum Phase Transition of the Hubbard Model on the Honeycomb Lattice

We take advantage of recent improvements in the grand canonical Hybrid Monte Carlo algorithm, to perform a precision study of the single-particle gap in the hexagonal Hubbard model, with on-site electron-electron interactions. After carefully controlled analyses of the Trotter error, the thermodynamic limit, and finite-size scaling with inverse temperature, we find a critical coupling of $U_c/κ=3.834(14)$ and the critical exponent $zν=1.185(43)$. Under the assumption that this corresponds to the expected anti-ferromagnetic Mott transition, we are also able to provide a preliminary estimate $β=1.095(37)$ for the critical exponent of the order parameter. We consider our findings in view of the $SU(2)$ Gross-Neveu, or chiral Heisenberg, universality class. We also discuss the computational scaling of the Hybrid Monte Carlo algorithm, and possible extensions of our work to carbon nanotubes, fullerenes, and topological insulators.

preprint2016arXiv

Computing the nucleon Dirac radius directly at $Q^2=0$

We describe a lattice approach for directly computing momentum derivatives of nucleon matrix elements using the Rome method, which we apply to obtain the isovector magnetic moment and Dirac radius. We present preliminary results calculated at the physical pion mass using a 2HEX-smeared Wilson-clover action. For removing the effects of excited-state contamination, the calculations were done at three source-sink separations and the summation method was used.

preprint2015arXiv

High-precision calculation of the strange nucleon electromagnetic form factors

We report a direct lattice QCD calculation of the strange nucleon electromagnetic form factors $G_E^s$ and $G_M^s$ in the kinematic range $0 \leq Q^2 \lesssim 1.2\: {\rm GeV}^2$. For the first time, both $G_E^s$ and $G_M^s$ are shown to be nonzero with high significance. This work uses closer-to-physical lattice parameters than previous calculations, and achieves an unprecedented statistical precision by implementing a recently proposed variance reduction technique called hierarchical probing. We perform model-independent fits of the form factor shapes using the $z$-expansion and determine the strange electric and magnetic radii and magnetic moment. We compare our results to parity-violating electron-proton scattering data and to other theoretical studies.

preprint2015arXiv

QCD thermodynamics with continuum extrapolated Wilson fermions II

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

preprint2014arXiv

Adaptive Aggregation Based Domain Decomposition Multigrid for the Lattice Wilson Dirac Operator

In lattice QCD computations a substantial amount of work is spent in solving discretized versions of the Dirac equation. Conventional Krylov solvers show critical slowing down for large system sizes and physically interesting parameter regions. We present a domain decomposition adaptive algebraic multigrid method used as a precondtioner to solve the "clover improved" Wilson discretization of the Dirac equation. This approach combines and improves two approaches, namely domain decomposition and adaptive algebraic multigrid, that have been used seperately in lattice QCD before. We show in extensive numerical test conducted with a parallel production code implementation that considerable speed-up over conventional Krylov subspace methods, domain decomposition methods and other hierarchical approaches for realistic system sizes can be achieved.

preprint2014arXiv

Charmonium spectral functions from 2+1 flavour lattice QCD

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

preprint2014arXiv

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

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

preprint2014arXiv

Recent results on the Equation of State of QCD

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

preprint2014arXiv

Spectral functions of charmonium with 2+1 flavours of dynamical quarks

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.

preprint2013arXiv

Determination of SU(2) ChPT LECs from 2+1 flavor staggered lattice simulations

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.

preprint2013arXiv

Lattice QCD at the physical point meets SU(2) chiral perturbation theory

We perform a detailed, fully-correlated study of the chiral behavior of the pion mass and decay constant, based on 2+1 flavor lattice QCD simulations. These calculations are implemented using tree-level, O(a)-improved Wilson fermions, at four values of the lattice spacing down to 0.054 fm and all the way down to below the physical value of the pion mass. They allow a sharp comparison with the predictions of SU(2) chiral perturbation theory (χPT) and a determination of some of its low energy constants. In particular, we systematically explore the range of applicability of NLO SU(2) χPT in two different expansions: the first in quark mass (x-expansion), and the second in pion mass (ξ-expansion). We find that these expansions begin showing signs of failure around M_π=300 MeV for the typical percent-level precision of our N_f=2+1 lattice results. We further determine the LO low energy constants (LECs), F=88.0 \pm 1.3\pm 0.3 and B^\msbar(2 GeV)=2.58 \pm 0.07 \pm 0.02 GeV, and the related quark condensate, Σ^\msbar(2 GeV)=(271\pm 4\pm 1 MeV)^3, as well as the NLO ones, l_3=2.5 \pm 0.5 \pm 0.4 and l_4=3.8 \pm 0.4 \pm 0.2, with fully controlled uncertainties. We also explore the NNLO expansions and the values of NNLO LECs. In addition, we show that the lattice results favor the presence of chiral logarithms. We further demonstrate how the absence of lattice results with pion masses below 200 MeV can lead to misleading results and conclusions. Our calculations allow a fully controlled, ab initio determination of the pion decay constant with a total 1% error, which is in excellent agreement with experiment.

preprint2013arXiv

Nucleon form factors with light Wilson quarks

We present nucleon observables - primarily isovector vector form factors - from calculations using 2+1 flavors of Wilson quarks. One ensemble is used for a dedicated high-precision study of excited-state effects using five source-sink separations between 0.7 and 1.6 fm. We also present results from a larger set of calculations that include an ensemble with pion mass 149 MeV and box size 5.6 fm, which nearly eliminates the uncertainty associated with extrapolation to the physical pion mass. The results show agreement with experiment for the vector form factors, which occurs only when excited-state contributions are reduced. Finally, we show results from a subset of ensembles that have pion mass 254 MeV with varying temporal and spatial box sizes, which we use for a controlled study of finite-volume effects and a test of the "$m_πL=4$" rule of thumb.

preprint2013arXiv

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

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.

preprint2012arXiv

Aggregation-based Multilevel Methods for Lattice QCD

In Lattice QCD computations a substantial amount of work is spent in solving the Dirac equation. In the recent past it has been observed that conventional Krylov solvers tend to critically slow down for large lattices and small quark masses. We present a Schwarz alternating procedure (SAP) multilevel method as a solver for the Clover improved Wilson discretization of the Dirac equation. This approach combines two components (SAP and algebraic multigrid) that have separately been used in lattice QCD before. In combination with a bootstrap setup procedure we show that considerable speed-up over conventional Krylov subspace methods for realistic configurations can be achieved.

preprint2012arXiv

Nucleon structure with pion mass down to 149 MeV

We present isovector nucleon observables: the axial, tensor, and scalar charges and the Dirac radius. Using the BMW clover-improved Wilson action and pion masses as low as 149 MeV, we achieve good control over chiral extrapolation to the physical point. Our analysis is done using three different source-sink separations in order to identify excited-state effects, and we make use of the summation method to reduce their size.

preprint2012arXiv

Systematic errors in partially-quenched QCD plus QED lattice simulations

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

preprint2011arXiv

Correlations and fluctuations from lattice QCD

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

preprint2011arXiv

Excited state contamination in nucleon structure calculations

Among the sources of systematic error in nucleon structure calculations is contamination from unwanted excited states. In order to measure this systematic error, we vary the operator insertion time and source-sink separation independently. We compute observables for three source-sink separations between 0.93 fm and 1.39 fm using clover-improved Wilson fermions and pion masses as low as 150 MeV. We explore the use of a two-state model fit to subtract off the contribution from excited states.

preprint2011arXiv

Fluctuations of conserved charges at finite temperature from lattice QCD

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

preprint2011arXiv

QCD thermodynamics with Wilson fermions

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

preprint2011arXiv

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

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

preprint2011arXiv

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

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

preprint2010arXiv

Is there still any Tc mystery in lattice QCD? Results with physical masses in the continuum limit III

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.

preprint2010arXiv

N_f=2+1 flavour equation of state

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

preprint2010arXiv

QCD transition temperature: full staggered result

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

preprint2010arXiv

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

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

preprint2010arXiv

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

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

preprint2010arXiv

The QCD equation of state with dynamical quarks

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