We report new results from our finite size scaling analysis of 4d compact pure U(1) gauge theory with Wilson action. Investigating several cumulants of the plaquette energy within the Borgs-Kotecky finite size scaling scheme we find strong evidence for a first-order phase transition and present a high precision value for the critical coupling in the thermodynamic limit.Comment: Lattice2002(Spin

We investigate optimal choices for the (outer) iteration method to use when solving linear systems with Neuberger's overlap operator in QCD. Different formulations for this operator give rise to different iterative solvers, which are optimal for the respective formulation. We compare these methods in theory and practice to find the overall optimal one. For the first time, we apply the so-called SUMR method of Jagels and Reichel to the shifted unitary version of Neuberger's operator, and show that this method is in a sense the optimal choice for propagator computations. When solving the "squared" equations in a dynamical simulation with two degenerate flavours, it turns out that the CG method should be used.

The computational costs of calculating the matrix sign function of the overlap operator together with fundamental numerical problems related to the discontinuity of the sign function in the kernel eigenvalues are the major obstacles towards simulations with dynamical overlap fermions using the Hybrid Monte Carlo algorithm. In a previous paper of the present series we introduced optimal numerical approximation of the sign function and have developed highly advanced preconditioning and relaxation techniques which speed up the the inversion of the overlap operator by nearly an order of magnitude.In this fourth paper of the series we construct an HMC algorithm for overlap fermions. We approximate the matrix sign function using the Zolotarev rational approximation, treating the smallest eigenvalues of the Wilson operator exactly within the fermionic force. Based on this we derive the fermionic force for the overlap operator. We explicitly solve the problem of the Dirac delta-function terms arising through zero crossings of eigenvalues of the Wilson operator. The main advantage of this scheme is that its energy violations scale better than O(∆τ 2 ) and thus are comparable with the violations of the standard leapfrog algorithm over the course of a trajectory. We explicitly prove that our algorithm satisfies reversibility and area conservation. We present test results from our algorithm on 4 4 , 6 4 , and 8 4 lattices.

We describe results of a high-statistics finite size scaling analysis of 4d compact U(1) lattice gauge theory with Wilson action at the phase transition point. Using a multicanonical hybrid Monte Carlo algorithm we generate data samples with more than 150 tunneling events between the metastable states of the system, on lattice sizes up to 18^4. We performed a first analysis within the Borgs-Kotecky finite size scaling scheme. As a result, we report evidence for a first-order phase transition with a plaquette energy gap, G=0.02667(20), at a transition coupling, beta_T=1.011128(11).Comment: Lattice 2000 (Topics in Gauge Theories),6 pages, 6 figures, LaTe

We demonstrate that substantial progress can be achieved in the study of the phase structure of 4-dimensional compact QED by a joint use of hybrid Monte Carlo and multicanonical algorithms, through an efficient parallel implementation. This is borne out by the observation of considerable speedup of tunnelling between the metastable states, close to the phase transition, on the Wilson line. We estimate that the creation of adequate samples (with order 100 flip-flops) becomes a matter of half a year's runtime at 2 Gflops sustained performance for lattices of size up to 244 .

No abstract

We demonstrate that substantial progress can be achieved in the study of the phase structure of 4-dimensional compact QED by a joint use of hybrid Monte Carlo and multicanonical algorithms, through an efficient parallel implementation. This is borne out by the observation of considerable speedup of tunnelling between the metastable states, close to the phase transition, on the Wilson line. Our approach leads to a general parallelization scheme for the efficient stochastic sampling of systems where (a part of) the Hamiltonian involves the total action or energy in each update step.Comment: LATTICE(Algorithm) 3 pages,1 figur

No abstract