BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Ulli Wolff appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We show the equivalence of the 2D Ising model to standard free Euclidean lattice fermions of the Wilson Majorana type. The equality of the loop representations for the partition functions of both systems is established exactly for finite lattices with well-defined boundary conditions. The honeycomb lattice is particularly simple in this context and therefore discussed first and only then followed by the more familiar square lattice case.
We study the O(3) sigma model in $D=2$ on the lattice with a Boltzmann weight linearized in $β$ on each link. While the spin formulation now suffers from a sign-problem the equivalent loop model remains positive and becomes particularly simple. By studying the transfer matrix and by performing Monte Carlo simulations in the loop form we study the mass gap coupling in a step scaling analysis. The question addressed is, whether or not such a simplified action still has the right universal continuum limit. If the answer is affirmative this would be helpful in widening the applicability of worm algorithm methods.
We develop a strategy for the non-perturbative determination of the O(a)-improvement coefficient c_sw for Wilson fermions with massive sea quarks. The improvement condition is defined via the PCAC relation in the Schrödinger functional. It is imposed along a line of constant physics designed to be close to the correct mass of the charm quark. The numerical work uses the tree-level improved Lüscher-Weisz gauge action in N_f=3+1 Lattice QCD.
We review the long term project of the ALPHA collaboration to compute in QCD the running coupling constant and quark masses at high energy scales in terms of low energy hadronic quantities. The adapted techniques required to numerically carry out the required multiscale non-perturbative calculation with our special emphasis on the control of systematic errors are summarized. The complete results in the two dynamical flavor approximation are reviewed and an outlook is given on the ongoing three flavor extension of the programme with improved target precision.
We define a finite size renormalization scheme for $ϕ^4$ theory which in the thermodynamic limit reduces to the standard scheme used in the broken phase. We use it to re-investigate the question of triviality for the four dimensional infinite bare coupling (Ising) limit. The relevant observables all rely on two-point functions and are very suitable for a precise estimation with the worm algorithm. This contribution updates an earlier publication by analysing a much larger dataset.
We study the standard one-component $φ^4$-theory in four dimensions. A renormalized coupling is defined in a finite size renormalization scheme which becomes the standard scheme of the broken phase for large volumes. Numerical simulations are reported using the worm algorithm in the limit of infinite bare coupling. The cutoff dependence of the renormalized coupling closely follows the perturbative Callan Symanzik equation and the triviality scenario is hence further supported.
All-order strong coupling simulations have been used to derive precise energies of string states in the confined phase of three dimensional Z(2) lattice gauge theory. The behavior of the ground state energy is here compared with predictions of effective string theory. Our new data reported here are consistent with known universal terms of the long string length ($L_0$) expansion known from effective string models in the continuum limit. Our precision is however still not sufficient to disentangle non-univeral terms of order $L_0^{-7}
We present a Monte-Carlo algorithm for the simulation of the all-order strong coupling expansion of the Z2 gauge theory. This random surface ensemble is equivalent to the standard formulation, but allows to measure some quantities, like Polyakov loop correlators or excess free energies, with an accuracy that could not have been easily achieved with traditional simulation methods. One interesting application of the algorithm is the investigation of effective string theories.
We represent low dimensional quantum mechanical Hamiltonians by moderately sized finite matrices that reproduce the lowest O(10) boundstate energies and wave functions to machine precision. The method extends also to Hamiltonians that are neither Hermitian nor PT symmetric and thus allows to investigate whether or not the spectra in such cases are still real. Furthermore, the approach is especially useful for problems in which a position-dependent mass is adopted, for example in effective-mass models in solid-state physics or in the approximate treatment of coupled nuclear motion in molecular physics or quantum chemistry. The performance of the algorithm is demonstrated by considering the inversion motion of different isotopes of ammonia molecules within a position-dependent-mass model and some other examples of one- and two-dimensional Hamiltonians that allow for the comparison to analytical or numerical results in the literature.
We describe an explicit algorithm to factorize an even antisymmetric N^2 matrix into triangular and trivial factors. This allows for a straight forward computation of Pfaffians (including their signs) at the cost of N^3/3 flops.
Worm methods to simulate the Ising model in the Aizenman random current representation including a low noise estimator for the connected four point function are extended to allow for antiperiodic boundary conditions. In this setup several finite size renormalization schemes are formulated and studied with regard to the triviality of ϕ^4 theory in four dimensions. With antiperiodicity eliminating the zero momentum Fourier mode a closer agreement with perturbation theory is found compared to the periodic torus.
Worm algorithms have been very successful with the simulation of sigma models with fixed length spins which result from scalar field theories in the limit of infinite quartic coupling lambda. Here we investigate closer their algorithmic efficiency at finite and even vanishing lambda for the one component model in dimensions D = 2, 3, 4.
We present an algorithm in which the all-order strong coupling expansion of the Abelian U(1) gauge theory with Wilson plaquette action is sampled. In addition to the vacuum closed surface graphs of the partition function we propose to also allow for a class of defects (boundaries) related to Wilson loops in the ensemble. The efficiency of our scheme in estimating various observables is compared to a standard Metropolis algorithm.
We discuss the status of different determinations of alpha_s, motivating a precise and reliable computation from lattice QCD. In order to suppress perturbative errors, the non-perturbative computation has to reach high energy scales mu. Such results already exist in the SF-scheme for N_f=0,2 and N_f=3. We recently added the running with four massless flavours in a range of alpha from about 0.07 to 0.3 . It is based on our recent determination of the Sheikholeslami Wohlert coefficient in the four-flavour theory.
We reformulate the O(N) sigma model as a loop model whose configurations are the all-order strong coupling graphs of the original model. The loop configurations are represented by a pointer list in the computer and a Monte Carlo update scheme is proposed. Sample simulations are reported and the method turns out to be similarly efficient as the reflection cluster method, but it has greater potential for systematic generalization to other lattice field theories. A variant action suggested by the method is also simulated and leads to a rather extreme demonstration of the concept of universality of the scaling or continuum limit.
We exactly reformulate the lattice CP(N-1) spin model on a D dimensional torus as a loop model whose configurations correspond to the complete set of strong coupling graphs of the original system. A Monte Carlo algorithm is described and tested that samples the loop model with its configurations stored and manipulated as a linked list. Complete absence of critical slowing down and correspondingly small errors are found at D=2 for several observables including the mass gap. Using two different standard lattice actions universality is demonstrated in a finite size scaling study. The topological charge is identified in the loop model but not yet investigated numerically.
We give an overview on recently accomplished successful generalizations of `worm' or loop gas simulation methods to O(N) and CP(N-1) sigma models and to simple fermion models. Beside the advantage of (practically) eliminated critical slowing down we also explain additional opportunities to estimate some observables with extremely improved signal to noise levels.
We have calculated the step scaling function and the running coupling of QCD in the Schroedinger functional scheme with four flavors of O(a) improved Wilson quarks. Comparisons of our non-perturbative results with 2-loop and 3-loop perturbation theory as well as with non-perturbative data for only two flavors are made.
We study a renormalized coupling g and mass m in four dimensional phi^4 theory on tori with finite size z=mL. Precise numerical values close to the continuum limit are reported for z=1,2,4, based on Monte Carlo simulations performed in the equivalent all-order strong coupling reformulation. Ordinary renormalized perturbation theory is found to work marginally at z=2 and and to fail at z=1. By exactly integrating over the constant field mode we set up a renormalized expansion in z and compute three nontrivial orders. These results reasonably agree with the numerical data at small z. In the new expansion, the universal continuum limit exists as expected from multiplicative renormalizability. The triviality scenario is corroborated with significant precision.
We have determined the non-perturbative O(a)-improvement coefficient c_sw for four flavors of Wilson quarks with the plaquette gauge action in a range of beta >= 5. The data are fitted with several Pade approximation formulae to get an impression of the stability. A small extrapolation below beta=5 seems acceptable.
We present a non-perturbative computation of the running of the coupling alpha_s in QCD with two flavours of dynamical fermions in the Schroedinger functional scheme. We improve our previous results by a reliable continuum extrapolation. The Lambda-parameter characterizing the high-energy running is related to the value of the coupling at low energy in the continuum limit. An estimate of Lambda*r_0 is given using large-volume data with lattice spacings a from 0.07 fm to 0.1 fm. It translates into Lambda_{MSbar}^{(2)}=245(16)(16) MeV [assuming r_0=0.5 fm]. The last step still has to be improved to reduce the uncertainty.
Recent exact predictions for the massive scaling limit of the two dimensional XY-model are based on the equivalence with the sine-Gordon theory and include detailed results on the finite size behavior. The so-called step-scaling function of the mass gap is simulated with very high precision and found consistent with analytic results in the continuum limit. To come to this conclusion, an also predicted form of a logarithmic decay of lattice artifacts was essential to use for the extrapolation.
For the family of O(n) invariant nonlinear sigma-models we consider boundary conditions that are periodic up to an O(n) rotation. The helicity modulus is related to the change in free energy under variations of the corresponding angle. It defines a nonperturbative finite volume running coupling similar to the Schr"odinger functional for QCD. For the two-dimensional O(3)-model we investigate this quantity by analytical and numerical techniques. We establish its universal continuum relation to the finite volume massgap coupling at all scales and coupling strengths.
High precision measurements of the renormalized zero-momentum 4-point coupling g_R and of the Luscher-Weisz-Wolff running coupling gbar(L) = L*m(L) performed with two different lattice actions in the non-perturbative region confirm the earlier observations, that the cutoff effects look linear, in contrast to perturbative considerations. The use of different actions allows one to make a more reliable estimate on the continuum limit. The measurements were done for infinite volume correlation length up to 350.
The scaling of the Schr"odinger functional coupling is studied numerically and perturbatively for an SU(3) lattice gauge field coupled to an O(a) improved bosonic spinor field. This corresponds to QCD with minus two light flavours and is used as a numerically less costly test case for real QCD. A suitable algorithm is developed, and the influence of the matter fields on the continuum limit and the lattice artefacts are studied in detail.
Existing non-perturbative computations of the running coupling of quenched QCD in the Schroedinger functional scheme are extended to scales mu lying much deeper in the low-energy regime. We are able to reach 1/mu ~ 0.9 fm, where a significant deviation from its perturbative evolution is observed.
We report performance benchmarks for several algorithms that we have used to simulate the Schr"odinger functional with two flavors of dynamical quarks. They include hybrid and polynomial hybrid Monte Carlo with preconditioning. An appendix describes a method to deal with autocorrelations for nonlinear functions of primary observables as they are met here due to reweighting.
As a prerequisite to dynamical fermion simulations a detailed study of optimal parameters and scaling behavior is conducted for the quenched Schrödinger functional at fixed renormalized coupling. We compare standard hybrid overrelaxation techniques with local and global hybrid Monte Carlo. Our efficiency measure is designed to be directly relevant for the strong coupling constant as used by the ALPHA collaboration.
A non perturbative finite size scaling technique is used to study a running coupling in lattice Yang-Mills theory coupled to a bosonic Wilson spinor field in the Schrödinger functional scheme. This corresponds to two negative flavours. The scaling behaviour in this case is compared to quenched results and to QCD with two flavours. The continuum limit is confronted with renormalized perturbation theory.
The coefficients multiplying the counterterms required for O($a$) improvement of the action and the isovector axial current in lattice QCD are computed non-perturbatively, in the quenched approximation and for bare gauge couplings $g_0$ in the range $0 \leq g_0 \leq 1$. A finite-size method based on the Schrödinger functional is employed, which enables us to perform all calculations at zero or nearly zero quark mass. As a by-product the critical hopping parameter $κ_c$ is obtained at all couplings considered.
We compute the two loop coefficient in the relation between the lattice bare coupling and the running coupling defined through the Schroedinger functional for the case of pure SU(2) gauge theory. This result is needed as one computational component to relate the latter to the MSbar-coupling, and it allows us to implement O(a) improvement of the Schroedinger functional to two-loop order. In addition, the two-loop beta-function is verified in a perturbative computation on the lattice, and the behavior of an improved bare coupling is investigated beyond one loop.
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.
An overview is given over the recently developed and now widely used Monte Carlo algorithms with reduced or eliminated critical slowing down. The basic techniques are overrelaxation, cluster algorithms and multigrid methods. With these tools one is able to probe much closer than before the universal continuum behavior of field theories on the lattice.
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.