Search

We investigate static quark anti-quark operators based on trial states formed from eigenvectors of the covariant three-dimensional lattice Laplace operator. We test the method by computing the static quark-anti-quark potential and comparing results to standard Wilson loop measurements. The new method is efficient not only for on-axis, but also for many off-axis quark-anti-quark separations when a fine spatial resolution is required. We further improve the ground-state overlap by using multiple eigenvector pairs, weighted with Gaussian profile functions of the eigenvalues, providing a variational basis. The method presented here can be applied to potential functions for all possible excitations of a gluonic string with fixed ends, hybrid or tetra-quark potentials, as well as static-light systems and allows visualization of the spatial distribution of the Laplace trial states.

We study the charmonium spectrum on an ensemble with two heavy dynamical quarks with a mass at half the physical charm quark mass. Operators for different quantum numbers are used in the framework of distillation with different smearing profiles to increase the overlap with ground and excited states. The use of exact distillation, large statistics and the absence of light quarks gives robust results for the charmonium spectrum. We also present preliminary results for the glueball spectrum in this theory., Comment: 10 pages, 8 figures, Proceedings of the 38th International Symposium on Lattice Field Theory, LATTICE2021 26th-30th July, 2021, Zoom/Gather@Massachusetts Institute of Technology

Eigenvalues of the Hermitian Wilson–Dirac operator are of special interest in several lattice QCD simulations, e.g., for noise reduction when evaluating all-to-all propagators. In this paper we present a Davidson-type eigensolver that utilizes the structural properties of the Hermitian Wilson–Dirac operator Q to compute eigenpairs of this operator corresponding to small eigenvalues. The main idea is to exploit a synergy between the (outer) eigensolver and its (inner) iterative scheme which solves shifted linear systems. This is achieved by adapting the multigrid DD- α AMG algorithm to a solver for shifted systems involving the Hermitian Wilson–Dirac operator. We demonstrate that updating the coarse grid operator using eigenvector information obtained in the generalized Davidson method is crucial to achieve good performance when calculating many eigenpairs, as our study of the local coherence shows. We compare our method with the commonly used software-packages PARPACK and PRIMME in numerical tests, where we are able to achieve significant improvements, with speed-ups of up to one order of magnitude and a near-linear scaling with respect to the number of eigenvalues. For illustration we compare the distribution of the small eigenvalues of Q on a 64 × 3 2 3 lattice with what is predicted by the Banks–Casher relation in the infinite volume limit.

In this review we summarize the ongoing effort to study extra-dimensional gauge theories with lattice simulations. In these models the Higgs field is identified with extra-dimensional components of the gauge field. The Higgs potential is generated by quantum corrections and is protected from divergencies by the higher dimensional gauge symmetry. Dimensional reduction to four dimensions can occur through compactification or localization. Gauge-Higgs unification models are often studied using perturbation theory. Numerical lattice simulations are used to go beyond these perturbative expectations and to include non-perturbative effects. We describe the known perturbative predictions and their fate in the strongly-coupled regime for various extra-dimensional models., Comment: Review published on the International Journal of Modern Physics A for the special issue "Lattice gauge theory beyond QCD". 22 pages and 7 figures

Most of today's lattice simulations are performed in the isospin symmetric limit of the light quark sector. Mass reweighting is a technique to include effects of isospin breaking in the sea quarks at moderate numerical cost. We will give a summary of our recent results on fine lattices with light quark masses and will show how light quark masses can be extracted by introducing suitable tuning conditions for the bare mass parameters. In general the reweighting factor introduces additional fluctuations and thus increases the statistical uncertainties. In the case of isospin reweighting this factor is a ratio of fermion determinants. The stochastic evaluation of the determinants potentially leads to stochastic noise in observables. We show the quark mass and the volume dependence of these fluctuations., Comment: 7 pages, 4 figures, 32st International Symposium on Lattice Field Theory - LATTICE 2014, Corrected (rescaled) x-axis of figures (1) and (2)

It is widely believed that the fermion determinant cannot be treated in global acceptance-rejection steps of gauge link configurations that differ in a large fraction of the links. However, for exact factorizations of the determinant that separate the ultraviolet from the infrared modes of the Dirac operator it is known that the latter show less variation under changes of the gauge field compared to the former. Using a factorization based on recursive domain decomposition allows for a hierarchical algorithm that starts with pure gauge updates of the links within the domains and ends after a number of filters with a global acceptance-rejection step. Ratios of determinants have to be treated stochastically and we construct techniques to reduce the noise. We find that the global acceptance rate is high on moderate lattice sizes and demonstrate the effectiveness of the hierarchical filter., Comment: 36 pages, 9 figures; improved version to be published in Comput.Phys.Commun., new results for the topological charge presented in Figure 7

In this paper, we develop a higher order symmetric partitioned Runge-Kutta method for a coupled system of differential equations on Lie groups. We start with a discussion on partitioned Runge-Kutta methods on Lie groups of arbitrary order. As symmetry is not met for higher orders, we generalize the method to a symmetric partitioned Runge-Kutta (SPRK) scheme. Furthermore, we derive a set of coefficients for convergence order 4. The SPRK integration method can be used, for example, in simulations of quantum field theories. Finally, we compare the new SPRK scheme numerically with the St\"ormer-Verlet scheme, one of the state-of-the-art schemes used in this subject., Comment: 18 pages, 2 figures, 2 tables

We present in detail a technique to extract the potential between a static quark and anti-quark pair from Wilson loops measured on dynamical configurations. This technique is based on HYP smearing and leads to an exponential improvement of the noise-to-signal ratio of Wilson loops. We explain why the correct continuum potential is obtained and show numerical evidence that the cut-off effects are small. We present precise results for the non-perturbative potential. As applications, we determine the scale r_0/a and study the shape of the static potential in the range of distances around r_0, where it can be compared with phenomenological potential models., Comment: 26 pages, 8 figures, 2 tables; new figure 2, update of figure 7, published version

We report on simulations with two flavors of O(a) improved degenerate Wilson fermions with Schroedinger functional boundary conditions. The algorithm which is used is Hybrid Monte Carlo with two pseudo-fermion fields as proposed by M. Hasenbusch. We investigate the numerical precision and sensitivity to reversibility violations of this algorithm. A gain of a factor two in CPU cost is reached compared with one pseudo-fermion field due to the larger possible step-size., Comment: 20 pages, 6 figures, changes in section 3, figures 4 and 6 replaced, addition of appendix A

A study of principle is conducted on the inclusion of the fermionic determinant as a Metropolis acceptance correction. It is carried out in the 2-D Schwinger model to prepare later applications to the Schr"odinger functional. A mixed stochastic/determistic acceptance step is found that allows to include some problematic modes in a way to avoid the collapse of the acceptance rate due to fluctuations., Comment: Lattice2003(machine), 3 pages, 3 Figures, talk by U.W

We present first results for the step scaling function sigma_P of the renormalization factor Z_P of the pseudoscalar density. The simulations are performed within the framework of the Schroedinger functional with two flavors of O(a) improved Wilson fermions. The knowledge of sigma_P is required to compute the renormalization group invariant quark masses. We also study the performance of a variant of the HMC algorithm using two pseudofermion fields., Comment: 3 pages, 2 figures, Lattice2002(spectrum)

Smearing the gauge links of dynamical configurations removes small scale unphysical vacuum fluctuations und thus improves the chiral properties of lattice fermions. We present a new algorithm for the simulation of dynamical fermions coupled via smeared links based on the standard pure gauge overrelaxation and heatbath updatings. Smeared links play a fundamental role in making this algorithm effective. At fixed lattice spacing the computational cost of the algorithm has an extra volume factor due to the finite volume of the lattice region which can be updated. As the continuum limit is approached the physical volume of the updated region remains constant. We simulated four flavors of staggered fermions coupled via hypercubic (HYP) smeared links. The simulation cost of the new algorithm on 10 fm^4 volumes is a factor 2-8 larger than with the standard Hybrid Monte Carlo but the improved properties of the HYP action allow to gain a factor 2 in the lattice spacing. The new algorithm could be applicable to simulations of more complicated chiral fermionic actions, like overlap or perfect actions., Comment: 9 pages, 1 figure

We present a method for finding lines of constant physics in the confinement phase of the SU(2) Higgs model on the lattice. The model is considered at finite values of the cut-off where it behaves like an effective field theory with three independent couplings. In particular, a renormalised quantity sensitive to a variation of the bare Higgs quartic self-coupling is constructed from generalised Binder cumulants. Numerical results for the non-perturbative matching of the bare parameters of the model between beta=2.2 and beta=2.4 are presented., Comment: 10 pages, 4 figures, small changes in the introduction