BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Stefan Schaefer

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

ResearcherUnclaimed source record
hep-lathep-phcond-mat.stat-mechcond-mat.dis-nnhep-th

Catalog footprint

What is connected

31works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

31 published item(s)

preprint2020arXiv

Frequency-splitting estimators for single-propagator traces

In these proceedings we address the computation of quark-line disconnected diagrams in lattice QCD. The evaluation of these diagrams is required for many phenomenologically interesting observables, but suffers from large statistical errors due to the vacuum and random-noise contributions to their variances. Motivated by a theoretical analysis of the variances, we introduce a new family of stochastic estimators of single-propagator traces built upon a frequency splitting combined with a hopping expansion of the quark propagator, and test their efficiency in two-flavour QCD with pions as light as 190 MeV. The use of these estimators reduces the cost of the computation by one to two orders of magnitude over standard estimators depending on the fermion bilinear. As a concrete application, we show the impact of these findings on the computation of the hadronic vacuum polarization contribution to the muon anomalous magnetic moment.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A multilevel algorithm for flow observables in gauge theories

We study the possibility of using multilevel algorithms for the computation of correlation functions of gradient flow observables. For each point in the correlation function an approximate flow is defined which depends only on links in a subset of the lattice. Together with a local action this allows for independent updates and consequently a convergence of the Monte Carlo process faster than the inverse square root of the number of measurements. We demonstrate the feasibility of this idea in the correlation functions of the topological charge and the energy density.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Domain decomposition, multi-level integration and exponential noise reduction in lattice QCD

We explore the possibility of computing fermionic correlators on the lattice by combining a domain decomposition with a multi-level integration scheme. The quark propagator is expanded in series of terms with a well defined hierarchical structure. The higher the order of a term, the (exponentially) smaller its magnitude, the less local is its dependence on the gauge field. Once inserted in a Wick contraction, the gauge-field dependence of the terms in the resulting series can be factorized so that it is suitable for multi-level Monte Carlo integration. We test the strategy in quenched QCD by computing the disconnected correlator of two flavor-diagonal pseudoscalar densities, and a nucleon two-point function. In either cases we observe a significant exponential increase of the signal-to-noise ratio.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

New extended interpolating operators for hadron correlation functions

New extended interpolating operators made of quenched three dimensional fermions are introduced in the context of lattice QCD. The mass of the 3D fermions can be tuned in a controlled way to find a better overlap of the extended operators with the states of interest. The extended operators have good renormalisation properties and are easy to control when taking the continuum limit. Moreover the short distance behaviour of the two point functions built from these operators is greatly improved. The operators have been numerically implemented and a comparison to point sources and Jacobi smeared sources has been performed on the new CLS configurations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The determination of $α_s$ by the ALPHA collaboration

We review the ALPHA collaboration strategy for obtaining the QCD coupling at high scale. In the three-flavor effective theory it avoids the use of perturbation theory at $α> 0.2$ and at the same time has the physical scales small compared to the cutoff $1/a$ in all stages of the computation. The result $Λ_\overline{MS}^{(3)}=332(14)$~MeV is translated to $α_\overline{MS}(m_Z)=0.1179(10)(2)$ by use of (high order) perturbative relations between the effective theory couplings at the charm and beauty quark "thresholds". The error of this perturbative step is discussed and estimated as $0.0002$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The large $N$ limit of the topological susceptibility of Yang-Mills gauge theory

We present a precise computation of the topological susceptibility $χ_{_\mathrm{YM}}$ of SU$(N)$ Yang-Mills theory in the large $N$ limit. The computation is done on the lattice, using high-statistics Monte Carlo simulations with $N=3, 4, 5, 6$ and three different lattice spacings. Two major improvements make it possible to go to finer lattice spacing and larger $N$ compared to previous works. First, the topological charge is implemented through the gradient flow definition; and second, open boundary conditions in the time direction are employed in order to avoid the freezing of the topological charge. The results allow us to extrapolate the dimensionless quantity $t_0^2χ_{_\mathrm{YM}}$ to the continuum and large $N$ limits with confidence. The accuracy of the final result represents a new quality in the verification of large $N$ scaling.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The topological susceptibility in the large-N limit of SU(N) Yang-Mills theory

We compute the topological susceptibility of the SU(N) Yang-Mills theory in the large-N limit with a percent level accuracy. This is achieved by measuring the gradient-flow definition of the susceptibility at three values of the lattice spacing for N=3,4,5,6. Thanks to this coverage of parameter space, we can extrapolate the results to the large-N and continuum limits with confidence. Open boundary conditions are instrumental to make simulations feasible on the finer lattices at the larger N.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Simulation of QCD with N_f=2+1 flavors of non-perturbatively improved Wilson fermions

We describe a new set of gauge configurations generated within the CLS effort. These ensembles have N_f=2+1 flavors of non-perturbatively improved Wilson fermions in the sea with the Luescher-Weisz action used for the gluons. Open boundary conditions in time are used to address the problem of topological freezing at small lattice spacings and twisted-mass reweighting for improved stability of the simulations. We give the bare parameters at which the ensembles have been generated and how these parameters have been chosen. Details of the algorithmic setup and its performance are presented as well as measurements of the pion and kaon masses alongside the scale parameter t_0.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the extraction of spectral quantities with open boundary conditions

We discuss methods to extract decay constants, meson masses and gluonic observables in the presence of open boundary conditions. The ensembles have been generated by the CLS effort and have 2+1 flavors of O(a)-improved Wilson fermions with a small twisted-mass term as proposed by Lüscher and Palombi. We analyse the effect of the associated reweighting factors on the computation of different observables.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Topological susceptibility and the sampling of field space in $N_f=2$ lattice QCD simulations

We present a measurement of the topological susceptibility in two flavor QCD. In this observable, large autocorrelations are present and also sizable cutoff effects have to be faced in the continuum extrapolation. Within the statistical accuracy of the computation, the result agrees with the expectation from leading order chiral perturbation theory.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Charm quark mass and D-meson decay constants from two-flavour lattice QCD

We present a computation of the charm quark's mass and the leptonic D-meson decay constants f_D and f_{D_s} in two-flavour lattice QCD with non-perturbatively O(a) improved Wilson quarks. Our analysis is based on the CLS configurations at two lattice spacings (a=0.065 and 0.048 fm, where the lattice scale is set by f_K) and pion masses ranging down to ~ 190 MeV at L*m_pi > 4, in order to perform controlled continuum and chiral extrapolations with small systematic uncertainties.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Improvement of N_f=3 lattice QCD with Wilson fermions and tree-level improved gauge action

We determine the parameter c_SW required for O(a)-improvement of the three flavor Wilson fermion action together with the tree-level Symanzik improved gauge action. The standard improvement condition is employed for a range of couplings. Additionally, we perform a check of the volume independence of c_SW and provide a preliminary estimate of the lattice spacing at our largest values of g_0^2.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Lattice QCD with open boundary conditions and twisted-mass reweighting

Lattice QCD simulations at small lattice spacings and quark masses close to their physical values are technically challenging. In particular, the simulations can get trapped in the topological charge sectors of field space or may run into instabilities triggered by accidental near-zero modes of the lattice Dirac operator. As already noted in ref. [1], the first problem is bypassed if open boundary conditions are imposed in the time direction, while the second can potentially be overcome through twisted-mass determinant reweighting [2]. In this paper, we show that twisted-mass reweighting works out as expected in QCD with open boundary conditions and 2+1 flavours of O(a) improved Wilson quarks. Further algorithmic improvements are tested as well and a few physical quantities are computed for illustration.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The strange quark mass and Lambda parameter of two flavor QCD

We complete the non-perturbative calculations of the strange quark mass and the Lambda parameter in two flavor QCD by the ALPHA collaboration. The missing lattice scale is determined via the kaon decay constant, for whose chiral extrapolation complementary strategies are compared. We also give a value for the scale r_0 in physical units as well as an improved determination of the renormalization constant Z_A.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Lattice QCD without topology barriers

As the continuum limit is approached, lattice QCD simulations tend to get trapped in the topological charge sectors of field space and may consequently give biased results in practice. We propose to bypass this problem by imposing open (Neumann) boundary conditions on the gauge field in the time direction. The topological charge can then flow in and out of the lattice, while many properties of the theory (the hadron spectrum, for example) are not affected. Extensive simulations of the SU(3) gauge theory, using the HMC and the closely related SMD algorithm, confirm the absence of topology barriers if these boundary conditions are chosen. Moreover, the calculated autocorrelation times are found to scale approximately like the square of the inverse lattice spacing, thus supporting the conjecture that the HMC algorithm is in the universality class of the Langevin equation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Non-renormalizability of the HMC algorithm

In lattice field theory, renormalizable simulation algorithms are attractive, because their scaling behaviour as a function of the lattice spacing is predictable. Algorithms implementing the Langevin equation, for example, are known to be renormalizable if the simulated theory is. In this paper we show that the situation is different in the case of the molecular-dynamics evolution on which the HMC algorithm is based. More precisely, studying the phi^4 theory, we find that the hyperbolic character of the molecular-dynamics equations leads to non-local (and thus non-removable) ultraviolet singularities already at one-loop order of perturbation theory.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Strange quark mass and Lambda parameter by the ALPHA collaboration

We determine f_K for lattice QCD in the two flavor approximation with non-perturbatively improved Wilson fermions. The result is used to set the scale for dimensionful quantities in CLS/ALPHA simulations. To control its dependence on the light quark mass, two different strategies for the chiral extrapolation are applied. Combining f_K and the bare strange quark mass with non-perturbative renormalization factors and step scaling functions computed in the Schroedinger Functional, we determine the RGI strange quark mass and the Lambda parameter in units of f_K.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Testing trivializing maps in the Hybrid Monte Carlo algorithm

We test a recent proposal to use approximate trivializing maps in a field theory to speed up Hybrid Monte Carlo simulations. Simulating the CP^{N-1} model, we find a small improvement with the leading order transformation, which is however compensated by the additional computational overhead. The scaling of the algorithm towards the continuum is not changed. In particular, the effect of the topological modes on the autocorrelation times is studied.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Algorithms for lattice QCD: progress and challenges

The development of improved algorithms for QCD on the lattice has enabled us to do calculations at small quark masses and get control over the chiral extrapolation. Also finer lattices have become possible, however, a severe slowing down associated with the topology of the gauge fields has been observed. This may prevent simulations of lattices fine enough for controlling the continuum extrapolation. This conference contribution introduces the basic concepts behind contemporary lattice algorithms, the current knowledge about their slowing down towards the continuum and its consequences for future lattice simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Autocorrelations in Hybrid Monte Carlo Simulations

Simulations of QCD suffer from severe critical slowing down towards the continuum limit. This problem is known to be prominent in the topological charge, however, all observables are affected to various degree by these slow modes in the Monte Carlo evolution. We investigate the slowing down in high statistics simulations and propose a new error analysis method, which gives a realistic estimate of the contribution of the slow modes to the errors.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Critical slowing down and error analysis in lattice QCD simulations

We study the critical slowing down towards the continuum limit of lattice QCD simulations with Hybrid Monte Carlo type algorithms. In particular for the squared topological charge we find it to be very severe with an effective dynamical critical exponent of about 5 in pure gauge theory. We also consider Wilson loops which we can demonstrate to decouple from the modes which slow down the topological charge. Quenched observables are studied and a comparison to simulations of full QCD is made. In order to deal with the slow modes in the simulation, we propose a method to incorporate the information from slow observables into the error analysis of physical observables and arrive at safer error estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Wilson fermions in the epsilon regime

We extend the epsilon-expansion of continuum chiral perturbation theory to nonzero lattice spacing in the framework of Wilson Chiral Perturbation Theory. We distinguish various regimes by defining the relative power counting of the quark mass m and the lattice spacing a. We observe that for m ~ a Lambda^2_QCD, the explicit breaking of chiral symmetry in Wilson fermions is still driven by the quark mass and lattice corrections are highly suppressed. The lattice spacing effects become more pronounced for smaller quark masses and may lead to non-trivial corrections of the continuum results at next-to-leading order. We compute these corrections for standard current and density correlation functions. A fit to lattice data shows that these corrections are small, as expected.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Epsilon regime calculations with reweighted clover fermions

We perform fully dynamical simulations at small quark masses by reweighting in the quark mass, calculating the weight factors stochastically. This approach avoids some of the technical difficulties associated with direct simulations. We find that the weight factors fluctuate only moderately on nHYP smeared dynamical Wilson-clover ensembles, and demonstrate that the overlap between the original and reweighted configurations is large both for short and long distance observables. We could successfully reweight 16^4, (1.85fm)^4 volume configurations from m_q 20 MeV to m_q~5 MeV, and 24^4, (2.77fm)^4 configurations from m_q~8MeV to m_q~MeV quark masses, reaching the epsilon-regime. Using the pseudoscalar and axialvector correlators we predict the low energy constants Sigma and F and study their volume and mass dependence.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Hypercubic Smeared Links for Dynamical Fermions

We investigate a variant of hypercubic gauge link smearing where the SU(3) projection is replaced with a normalization to the corresponding unitary group. This smearing is differentiable and thus suitable for use in dynamical fermion simulations using molecular dynamics type algorithms. We show that this smearing is as efficient as projected hypercubic smearing in removing ultraviolet noise from the gauge fields. We test the normalized hypercubic smearing in dynamical improved (clover) Wilson and valence overlap simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

New findings for topological excitations in SU(3) lattice gauge theory

We probe the SU(3) vacuum using eigenvectors of the Dirac operator with an arbitrary phase for the temporal boundary condition. We consider configurations with topological charge |Q| = 1 near the QCD phase transition and at low temperatures on a torus. For all our ensembles we show that the zero-mode of the Dirac operator changes its position as one changes the phase of the boundary condition. For ensembles near the QCD phase transition our results closely resemble the behavior of zero-modes for Kraan - van Baal solutions of the classical Yang-Mills equations where the individual lumps are interpreted as monopoles. Our findings near T_c and on the torus show that for both cases an excitation with topological charge |Q| = 1 is built from several separate lumps.

0 citations0 reviewsNo code linkedOpen access