Source author record

Martin Lüscher

Martin Lüscher 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

24works
3topics
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

24 published item(s)

preprint2020arXiv

Master-field simulations of O(a)-improved lattice QCD: Algorithms, stability and exactness

In master-field simulations of lattice QCD, the expectation values of interest are obtained from a single or at most a few representative gauge-field configurations on very large lattices. If the light quarks are included, the generation of these fields using standard techniques is however challenging in view of various algorithmic instabilities and precision issues. Ways to overcome these problems are described here for the case of the O(a)-improved Wilson formulation of lattice QCD and the viability of the proposed measures is then checked in extensive simulations of the theory with 2+1 flavours of quarks.

preprint2016arXiv

The gradient flow coupling from numerical stochastic perturbation theory

Perturbative calculations of gradient flow observables are technically challenging. Current results are limited to a few quantities and, in general, to low perturbative orders. Numerical stochastic perturbation theory is a potentially powerful tool that may be applied in this context. Precise results using these techniques, however, require control over both statistical and systematic uncertainties. In this contribution, we discuss some recent algorithmic developments that lead to a substantial reduction of the cost of the computations. The matching of the ${\overline{\rm MS}}$ coupling with the gradient flow coupling in a finite box with Schrödinger functional boundary conditions is considered for illustration.

preprint2014arXiv

Properties and uses of the Wilson flow in lattice QCD

Theoretical and numerical studies of the Wilson flow in lattice QCD suggest that the gauge field obtained at flow time t>0 is a smooth renormalized field. The expectation values of local gauge-invariant expressions in this field are thus well-defined physical quantities that probe the theory at length scales on the order of sqrt(t). Moreover, by transforming the QCD functional integral to an integral over the gauge field at a specified flow time, the emergence of the topological (instanton) sectors in the continuum limit becomes transparent and is seen to be caused by a dynamical effect that rapidly separates the sectors when the lattice spacing is reduced from 0.1 fm to smaller values.

preprint2014arXiv

Step scaling and the Yang-Mills gradient flow

The use of the Yang-Mills gradient flow in step-scaling studies of lattice QCD is expected to lead to results of unprecedented precision. Step scaling is usually based on the Schrödinger functional, where time ranges over an interval [0,T] and all fields satisfy Dirichlet boundary conditions at time 0 and T. In these calculations, potentially important sources of systematic errors are boundary lattice effects and the infamous topology-freezing problem. The latter is here shown to be absent if Neumann instead of Dirichlet boundary conditions are imposed on the gauge field at time 0. Moreover, the expectation values of gauge-invariant local fields at positive flow time (and of other well localized observables) that reside in the center of the space-time volume are found to be largely insensitive to the boundary lattice effects.

preprint2014arXiv

Topology, the Wilson flow and the HMC algorithm

An old and apparently persistent problem in numerical lattice QCD is that the simulations tend to get trapped in a sector of fixed topological charge when the lattice spacing is taken to zero. The effect sets in very rapidly and may invalidate the simulation results in certain cases. In this talk, the issue is discussed using the Wilson flow as a tool. The flow has a simple scaling behaviour and allows one to understand how exactly the topological sectors emerge in the continuum limit. Further studies however suggest that the observed slowdown of the simulations at small lattice spacings is only partly caused by the emergence of the sectors.

preprint2013arXiv

Chiral symmetry and the Yang--Mills gradient flow

In the last few years, the Yang--Mills gradient flow was shown to be an attractive tool for non-perturbative studies of non-Abelian gauge theories. Here a simple extension of the flow to the quark fields in QCD is considered. As in the case of the pure-gauge gradient flow, the renormalizability of correlation functions involving local fields at positive flow times can be established using a representation through a local field theory in 4+1 dimensions. Applications of the extended flow in lattice QCD include non-perturbative renormalization and O(a) improvement as well as accurate calculations of the chiral condensate and of the pseudo-scalar decay constant in the chiral limit.

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.

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.

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.

preprint2011arXiv

Perturbative analysis of the gradient flow in non-abelian gauge theories

The gradient flow in non-abelian gauge theories on R^4 is defined by a local diffusion equation that evolves the gauge field as a function of the flow time in a gauge-covariant manner. Similarly to the case of the Langevin equation, the correlation functions of the time-dependent field can be expanded in perturbation theory, the Feynman rules being those of a renormalizable field theory on R^4 x [0,oo). For any matter multiplet and to all loop orders, we show that the correlation functions are finite, i.e. do not require additional renormalization, once the theory in four dimensions is renormalized in the usual way. The flow thus maps the gauge field to a one-parameter family of smooth renormalized fields.

preprint2010arXiv

Universality of the topological susceptibility in the SU(3) gauge theory

The definition and computation of the topological susceptibility in non-abelian gauge theories is complicated by the presence of non-integrable short-distance singularities. Recently, alternative representations of the susceptibility were discovered, which are singularity-free and do not require renormalization. Such an expression is here studied quantitatively, using the lattice formulation of the SU(3) gauge theory and numerical simulations. The results confirm the expected scaling of the susceptibility with respect to the lattice spacing and they also agree, within errors, with computations of the susceptibility based on the use of a chiral lattice Dirac operator.

preprint2009arXiv

Trivializing maps, the Wilson flow and the HMC algorithm

In lattice gauge theory, there exist field transformations that map the theory to the trivial one, where the basic field variables are completely decoupled from one another. Such maps can be constructed systematically by integrating certain flow equations in field space. The construction is worked out in some detail and it is proposed to combine the Wilson flow (which generates approximately trivializing maps for the Wilson gauge action) with the HMC simulation algorithm in order to improve the efficiency of lattice QCD simulations.

preprint2008arXiv

Fluctuations and reweighting of the quark determinant on large lattices

We propose to stabilise HMC simulations of lattice QCD with very light Wilson quarks by splitting the quark determinant into two factors and by treating the factor that includes the contribution of the low modes of the Dirac operator as a reweighting factor. In general, determinant reweighting becomes inefficient on large lattices, because the statistical fluctuations of quark determinants increase exponentially with the lattice volume. Random matrix theory and some numerical studies now suggest that the low-mode contribution to the determinant behaves differently, which allows factorisations to be devised that preserve the efficiency of the simulation on large lattices.

preprint2004arXiv

String excitation energies in SU(N) gauge theories beyond the free-string approximation

In the presence of a static quark--antiquark pair, the spectrum of the low-lying states in SU($N$) gauge theories is discrete and likely to be described, at large quark separations $r$, by an effective string theory. The expansion of the excitation energies in powers of $1/r$, which derives from the latter, involves an increasing number of unknown couplings that characterize the string self-interactions. Using open--closed string duality, we show that the possible values of the couplings are constrained by a set of algebraic relations. In particular, the corrections of order $1/r^2$ must vanish, while the $1/r^3$ terms (which we work out for the few lowest levels) depend on a single adjustable coupling only.

preprint2003arXiv

Lattice QCD and the Schwarz alternating procedure

A numerical simulation algorithm for lattice QCD is described, in which the short- and long-distance effects of the sea quarks are treated separately. The algorithm can be regarded, to some extent, as an implementation at the quantum level of the classical Schwarz alternating procedure for the solution of elliptic partial differential equations. No numerical tests are reported here, but theoretical arguments suggest that the algorithm should work well also at small quark masses.

preprint2002arXiv

Quark confinement and the bosonic string

Using a new type of simulation algorithm for the standard SU(3) lattice gauge theory that yields results with unprecedented precision, we confirm the presence of a $γ/r$ correction to the static quark potential at large distances $r$, with a coefficient $γ$ as predicted by the bosonic string theory. In both three and four dimensions, the transition from perturbative to string behaviour is evident from the data and takes place at surprisingly small distances.

preprint2001arXiv

Locality and exponential error reduction in numerical lattice gauge theory

In non-abelian gauge theories without matter fields, expectation values of large Wilson loops and loop correlation functions are difficult to compute through numerical simulation, because the signal-to-noise ratio is very rapidly decaying for increasing loop sizes. Using a multilevel scheme that exploits the locality of the theory, we show that the statistical errors in such calculations can be exponentially reduced. We explicitly demonstrate this in the SU(3) theory, for the case of the Polyakov loop correlation function, where the efficiency of the simulation is improved by many orders of magnitude when the area bounded by the loops exceeds 1 fm^2.

preprint1996arXiv

Non-perturbative determination of the axial current normalization constant in O(a) improved lattice QCD

A finite-size technique is employed to compute the normalization constant $Z_A$ of the isovector axial current in lattice QCD. The calculation is carried out in the quenched approximation for values of the bare gauge coupling $g_0$ ranging from 0 to 1. In the lattice action and the lattice expression for the axial current we include the counterterms required for O(a) improvement, with non-perturbatively determined coefficients. With little additional work the normalization constant $Z_V$ of the improved isospin current is also obtained.

preprint1992arXiv

Computation of the Running Coupling in the SU(2) Yang-Mills Theory

A finite-size scaling technique is applied to the SU(2) gauge theory (without matter fields) to compute a non-perturbatively defined running coupling alpha(q) for a range of momenta q given in units of the string tension K. We find that already at rather low q, the evolution of alpha(q) is well described by the 2-loop approximation to the Callan-Symanzik beta-function. At the highest momentum reached, q=20 sqrt(K), we obtain alpha_MSbar(q)=0.187 +/- 0.005 +/- 0.009 for the running coupling in the MSbar scheme of dimensional regularization.

preprint1992arXiv

The Schrödinger Functional - a Renormalizable Probe for Non-Abelian Gauge Theories

Following Symanzik we argue that the Schrödinger functional in lattice gauge theories without matter fields has a well-defined continuum limit. Due to gauge invariance no extra counter terms are required. The Schrödinger functional is, moreover, accessible to numerical simulations. It may hence be used to study the scaling properties of the theory and in particular the evolution of the renormalized gauge coupling from low to high energies. A concrete proposition along this line is made and the necessary perturbative analysis of the Schrödinger functional is carried through to 1-loop order.