title: Localization in Lattice Gauge Theory and a New Multigrid Method

abstract: We show numerically that the lowest eigenmodes of the 2-dimensional Laplace-operator with SU(2) gauge couplings are strongly localized. A connection is drawn to the Anderson-Localization problem. A new Multigrid algorithm, capable to deal with these modes, shows no critical slowing down for this problem.


title: Spatial and Temporal Hadron Correlators below and above the Chiral Phase Transition

abstract: Hadronic correlation functions at finite temperature in QCD, with four flavours of dynamical quarks, have been analyzed both above and below the chiral symmetry restoration temperature. We have used both point and extended sources for spatial as well as temporal correlators. The effect of periodic temporal boundary conditions for the valence quarks on the spatial meson correlators has also been investigated. All our results are consistent with the existence of individual quarks at high temperatures. A measurement of the residual interaction between the quarks is presented.


title: The Block Spin Renormalization Group Approach and Two-Dimensional Quantum Gravity

abstract: A block spin renormalization group approach is proposed for the dynamical triangulation formulation of two-dimensional quantum gravity. The idea is to update link flips on the block lattice in response to link flips on the original lattice. Just as the connectivity of the original lattice is meant to be a lattice representation of the metric, the block links are determined in such a way that the connectivity of the block lattice represents a block metric. As an illustration, this approach is applied to the Ising model coupled to two-dimensional quantum gravity. The correct critical coupling is reproduced, but the critical exponent is obscured by unusually large finite size effects.


title: The random lattice as a regularization scheme

abstract: A semi-analytic method to compute the first coefficients of the renormalization group functions on a random lattice is introduced. It is used to show that the two-dimensional $O(N)$ non-linear $\sigma$-model regularized on a random lattice has the correct continuum limit. A degree $\kappa$ of ``randomness'' in the lattice is introduced and an estimate of the ratio $\Lambda_{random}/\Lambda_{regular}$ for two rather opposite values of $\kappa$ in the $\sigma$-model is also given. This ratio turns out to depend on $\kappa$.


title: Three-Dimensional Integrable Models and Associated Tangle Invariants

abstract: In this paper we show that the Boltzmann weights of the three-dimensional Baxter-Bazhanov model give representations of the braid group, if some suitable spectral limits are taken. In the trigonometric case we classify all possible spectral limits which produce braid group representations. Furthermore we prove that for some of them we get cyclotomic invariants of links and for others we obtain tangle invariants generalizing the cyclotomic ones.


title: On the exponential bound in four dimensional simplicial gravity

abstract: Simplicial quantum gravity has been proposed as a regularization for four dimensional quantum gravity. The partition function is constructed by performing a weighted sum over all triangulations of the 4-sphere. The model is well-defined only if the number of such triangulations consisting of $N$ simplexes is exponentially bounded. Numerical simulations seem so far to favor such a bound.


title: Fractal Wilson Loop - Area Law and Gauge Invariance in Next to Leading Order

abstract: We consider a fractal Wilson loop $<F_{P}>$ and present physical arguments why this should be a relevant observable in nature. We show for non-compact SU(2) lattice gauge theory in the next to leading order of strong coupling expansion that $<F_{P}>$ obeys an area law behavior and is gauge invariant.


title: Polyakov Loops and Finite-Size Effects of Hadron Masses in Lattice Full Q.C.D

abstract: The polarization of Polyakov type loops is responsible for the difference between quenched and unquenched finite size effects on the QCD mass spectrum. With a numerical simulation, using different sea quarks boundary conditions, we show that we can align the spatial Polyakov loops in a predefined direction. Starting from these results, we propose a procedure to partially remove the Polyakov type contributions in the meson propagators.


title: Volume dependence of the phase boundary in 4D dynamical triangulation

abstract: The number of configurations of the dynamical triangulation model of 4D euclidean quantum gravity appears to grow faster than exponentially with the volume, with the implication that the system would end up in the crumpled phase for any fixed $\kappa_2$ (inverse bare Newton constant). However, a scaling region is not excluded if we allow $\kappa_2$ to go to infinity together with the volume.


title: The Effect of Dynamical Gauge Field on the Chiral Fermion on a boundary

abstract: We study the effect of dynamical gauge field on the Kaplan's chiral fermion on a boundary in the strong gauge coupling limit in the extra dimension. To all orders of the hopping parameter expansion, we prove exact parity invariance of the fermion propagator on the boundary. This means that the chiral property of the boundary fermion, which seems to survive even in the presence of the gauge field from a perturbative point of view, is completely destroyed by the dynamics of the gauge field.


title: Antiferromagnetic Potts Models on the Square Lattice

abstract: We study the antiferromagnetic $q$-state Potts model on the square lattice for $q=3$ and $q=4$, using the Wang-Swendsen-Koteck\'y Monte Carlo algorithm and a new finite-size-scaling extrapolation method. For $q=3$ we obtain good control up to correlation length $\xi \sim 5000$; the data are consistent with $\xi(\beta) = A e^{2\beta} \beta (1 + a_1 e^{-\beta} + \ldots)$ as $\beta\to\infty$. For $q=4$ the model is disordered even at zero temperature.


title: Monte Carlo Methods for the Self-Avoiding Walk

abstract: This article is a pedagogical review of Monte Carlo methods for the self-avoiding walk, with emphasis on the extraordinarily efficient algorithms developed over the past decade.


title: Adiabatic Quantum Transport: Quantization and Fluctuations

abstract: Quillen's local index theorem is used to study the charge transport coefficients (adiabatic curvature) associated to the ground state of a Schr\"odinger operator for a charged (spinless) particle on a closed, multiply connected surface. The formula splits the adiabatic curvature into an explicit integral part and a fluctuating part which has a natural interpretation in terms of quantum chaos.


title: The SU(3) deconfining phase transition with Symanzik action

abstract: We report on the determination of the deconfining temperature in SU(3) pure gauge theory, using the Symanzik tree level improved action, on lattices of size 3 x 12^3, 4 x 16^3, 5 x 20^3, 6 x24^3. We find that the asymptotic scaling violation pattern is similar to the one observed using the Wilson action. We conclude that the irrelevant operators do not affect, in the range of couplings considered, the lattice beta function. An analysis based on an effective coupling formulation shows an apparent improvement.


title: Four Loop Result in $SU(3)$ Lattice Gauge Theory by a Stochastic Method: Lattice Correction to the Condensate

abstract: We describe a stochastic technique which allows one to compute numerically the coefficients of the weak coupling perturbative expansion of any observable in Lattice Gauge Theory. The idea is to insert the exponential representation of the link variables $U_\mu(x) \to \exp\{A_\mu(x)/\sqrt\beta\}$ into the Langevin algorithm and the observables and to perform the expansion in \beta^{-1/2}. The Langevin algorithm is converted into an infinite hierarchy of maps which can be exactly truncated at any order. We give the result for the simple plaquette of SU(3) up to fourth loop order (\beta^{-4}) which extends by one loop the previously known series.


title: Core -- a New Computational Technique for Lattice Systems

abstract: The COntractor REnormalization group (CORE) method, a new approach to solving Hamiltonian lattice systems, is introduced. The method combines contraction and variational techniques with the real-space renormalization group approach. It applies to lattice systems of infinite extent and is ideal for studying phase structure and critical phenomena. The CORE approximation is systematically improvable and can treat systems with dynamical fermions. The method is tested using the 1+1-dimensional Ising model.


title: Numerical Simulations and the Strength of the Electroweak Phase Transition

abstract: Numerical simulations are performed to study the finite temperature phase transition in the SU(2) Higgs model on the lattice. The strength of the first order phase transition is investigated by determining the latent heat and the interface tension on $L_t=2$ lattices. The values of the Higgs boson mass presently chosen are below 50 GeV. Our results are in qualitative agreement with two-loop resummed perturbation theory.


title: Scaling with a modified Wilson action which suppresses Z_2 artifacts in SU(2) lattice gauge theories

abstract: A modified Wilson action which suppresses plaquettes which take negative values is used to study the scaling behavior of the string tension. The use of the $\b_E$ scheme gives good agreement with asymptotic two loop results.


title: Finite-Size analysis of the 4-d abelian surface gauge model

abstract: We present the results of a finite-size analysis of the four dimensional abelian surface gauge model. This model is defined assigning abelian variables to the plaquettes of an hypercubical lattice, and is dual to the four dimensional Ising model. This last model is known to present a second order phase transition with mean field critical exponents. We have performed Monte Carlo simulations on several lattice sizes and high statistics. The analysis of the partition function zeroes and the specific heat scaling behaviour allowed us to estimate the critical coupling $\beta_c$ as well as the critical exponents $\nu$ and $\alpha$. Our results are consistent with the second order critical exponents $\nu = 1/2$ and $\alpha = 0$. The $\beta_c$ value is in perfect agreement with duality predictions from the 4-d Ising model. Nevertheless, the energy histograms show a seemingly non-vanishing double peak structure. The interface tension analysis suggests that this may be a finite size effect.


title: Exploratory Study of the 3-Gluon Vertex on the Lattice

abstract: We define and evaluate on the lattice the amputated 3-gluon vertex function in momentum space. We give numerical results for $16^3 \times 40$ and $24^3 \times 40$ quenched lattices at $\beta=6.0$. A good numerical signal is obtained, at the price of enforcing the gauge-fixing condition with high accuracy. By comparing results from two different lattice volumes, we try to investigate the crucial issue of finite volume effects. We also outline a method for the lattice evaluation of the QCD running coupling constant as defined from the 3-gluon vertex, while being aware that a realistic calculation will require larger $\beta$ values and very high statistics.


title: Detecting Dual Superconductivity in the Ground State of Gauge Theories - II

abstract: A monopole creation operator is constructed: its vacuum expectation value is an order parameter for dual superconductivity in that, if different from zero it signals spontaneous breaking of the U(1) symmetry corresponding to monopole charge conservation. The operator is tested on compact U(1) gauge theory on lattice. For SU(2) gauge theory it clearly demonstrates that confinement is produced by dual superconductivity.


title: Simulations of Dynamically Triangulated Gravity -- an Algorithm for Arbitrary Dimension

abstract: Recent models for discrete euclidean quantum gravity incorporate a sum over simplicial triangulations. We describe an algorithm for simulating such models in general dimensions. As illustration we show results from simulations in four dimensions


title: Gauge-Invariant Noncompact Lattice Simulations

abstract: Three techniques for performing gauge-invariant, noncompact lattice simulations of nonabelian gauge theories are discussed. In the first method, the action is not itself gauge invariant, but a kind of lattice gauge invariance is restored by random compact gauge transformations during the successive sweeps of the simulation. This method has been applied to pure $SU(2)$ gauge theory on a $12^4$ lattice, and Wilson loops have been measured at strong coupling, $\beta=0.5$. These Wilson loops display a confinement signal not seen in simulations performed earlier with the same action but without the random gauge transformations. In the second method, the action is gauge symmetrized by integrations over the group manifold. The third method is based upon a new, noncompact form of the action that is exactly invariant under lattice gauge transformations. The action is a natural discretization of the classical Yang-Mills action.


title: Radiative Corrections to the Kinetic Couplings in Nonrelativistic Lattice QCD

abstract: The heavy-quark mass and wave function renormalizations, energy shift, and radiative corrections to two important couplings, the so-called kinetic couplings, in nonrelativistic lattice QCD are determined to leading order in tadpole-improved perturbation theory. The scales at which to evaluate the running QCD coupling for these quantities, except the wave function renormalization, are obtained using the Lepage-Mackenzie prescription. When the bare quark mass is greater than the inverse lattice spacing, the kinetic coupling corrections are roughly 10\% of the tree-level coupling strengths; these corrections grow quickly as the bare quark mass becomes small. A need for computing the two-loop corrections to the energy shift and mass renormalization is demonstrated.


title: Fixed boundary conditions and phase transitions in pure gauge compact QED

abstract: We have simulated the pure gauge compact QED with fixed boundary conditions, on lattices from $6^4$ to $16^4$. We argue that a lattice with this fixed boundary imposition corresponds actually to a lattice with spherical topology. We have found the presence of a phase transition without any trace of discontinuity. Moreover, the specific heat and the Binder cumulant are qualitatively consistent with a second order phase transition. The implications of this observation on the nature of the compact lattice QED are discussed.


title: Monopole Percolation in pure gauge compact QED

abstract: The role of monopoles in quenched compact QED has been studied by measuring the cluster susceptibility and the order parameter $n_{max}/n_{tot}$ previously introduced by Hands and Wensley in the study of the percolation transition observed in non-compact QED. A correlation between these parameters and the energy (action) at the phase transition has been observed. We conclude that the order parameter $n_{max}/n_{tot}$ is a sensitive probe for studying the phase transition of pure gauge compact QED.


title: Nonassociative Algebras and Nonperturbative Field Theory for Hierarchical Models

abstract: Hierarchical renormalization group (RG) transformations are related to nonassociative algebras. These algebras serve as a new basic tool for a rigorous treatment of global RG flows and the search of nontrivial infrared fixed points. Convergent expansion methods are presented and analyzed in terms of algebra norms. It is shown that the infrared fixed points are solutions of a quadratic equation with an infinite number of unknowns. A continuous manifold of two dimensional periodic nontrivial fixed points is presented in terms of theta functions. Local Borel summability of the $\epsilon$- expansion for n-well fixed points is proved by algebraic methods.


title: Nonperturbative Studies of Quantum Gravity

abstract: One of several possibilities to construct a quantum theory of gravity is employing the Feynman path integral. This approach is plagued by some problems: the integration measure is not uniquely defined, the Einstein-Hilbert action unbounded, and perturbation theory nonrenormalizable. To make the path integral tractable one can approximate the continuous geometry of spacetime by a simplicial complex. The edge lengths of this lattice are considered as the dynamical degrees of freedom and Regge calculus is applied. In this work, numerical simulations using the Regge-Einstein action and a "compact" action show the occurence of a phase transition. The strength of this transition, separating a well-defined phase with finite expectation values from an ill-defined phase, is weaker for the compact action, which might be important for the continuum limit. To analyze the interaction mechanism of this formulation of quantum gravity, correlation functions between geometrical quantities like edge lengths, volume elements, and local curvatures have been computed. Our results for the two-point functions seem to prefer exchange particles with an effective mass. To ease treatment of quantum gravity a new approach is proposed consisting in a transformation of the path integral to the partition function of a spin system. This facilitates analytical and numerical calculations considerably. First results for the phase structure in two as in four dimensions are presented and indicate promising similarities to the original Regge theory.


title: $\pi NN$ and Pseudoscalar Form Factors from Lattice QCD

abstract: The $\pi NN$ form factor $g_{\pi NN}(q^2)$ is obtained from a quenched lattice QCD calculation of the pseudoscalar form factor $g_P(q^2)$ of the proton with pion pole dominance. We find that $g_{\pi NN}(q^2)$ fitted with the monopole form agrees well with the Goldberger-Treiman relation and is much preferred over the dipole form. The monopole mass is determined to be $0.75 \pm 0.14 {\rm GeV}$ which shows that $g_{\pi NN}(q^2)$ is rather soft. The extrapolated $\pi N$ coupling constant $g_{\pi NN} = 12.7 \pm 2.4$ is quite consistent with the phenomenological values. We also compare $g_{\pi NN}(q^2)$ with the axial form factor $g_A(q^2)$ to check the pion dominance in the induced pseudoscalar form factor $h_A(q^2)$ vis \`{a} vis chiral Ward identity.


title: Susceptibilities, the Specific Heat and a Cumulant in Two-Flavour QCD

abstract: We study the quark mass dependence of various response functions, which contribute to chiral susceptibilities and the specific heat in the staggered fermion formulation of two-flavour QCD. This yields information about the critical exponents $\alpha$, $\beta$ and $\delta$. In the case of the chiral susceptibility, obtained as derivative of the chiral order parameter with respect to the quark mass, we calculate all contributions. This allows to construct a cumulant of the order parameter, which is a scaling function and yields a direct determination of the critical exponent $\delta$. All our results are consistent with a second order phase transition.


title: Calculation and Interpretation of Hadron Correlation Functions in Lattice QCD

abstract: Several new developments in the calculation and interpretation of hadron density-density correlation functions are presented. The asymptotic behavior of correlation functions is determined from a tree diagram path integral. A method is developed to use this behavior to correct for leading image contributions on a finite periodic spatial lattice and to correct for the finite temporal extent of the lattice. Equal time correlation functions are shown to determine a sum of the ground state rms radius plus a polarization correction, and it is shown how to extract the hadron polarizability from unequal time correlation functions. Image-corrected correlation functions calculated in quenched lattice QCD are presented and used to determine the size of the pion and nucleon.


title: Monopole action from vacuum configurations in compact QED

abstract: It is possible to derive a monopole action from vacuum configurations obtained in Monte-Carlo simulations extending the method developed by Swendsen. We apply the method to compact QED both in the Villain and in the Wilson forms. The action of the natural monopoles in the Villain case is in fairly good agreement with that derived by the exact dual transformation. Comparing the monopole actions, we find (1) the DeGrand-Toussaint monopole definition may be useful for $\beta_V $ larger than about 0.5, (2) the Villain model well approximates the Wilson one for $\beta$ smaller than $\beta_c$ and (3) in the Wilson action the monopole condensation occurs in the confinement phase and $\beta_c$ may be explained by the energy-entropy balance of monopole loops like in the Villain case.


title: A Proposal of a Renormalization Group Transformation for Lattice Field Theories

abstract: We propose a new Real Space Renormalization Group transformation useful for Monte Carlo calculations in theories with global or local symmetries. From relaxation arguments we define the block-spin transformation with two tunable free parameters, adapted to the system's action. Varying them it is possible to place the fixed point very close to the simulation point. We show how the method works in a simple model with global symmetry: the three dimensional XY model.


title: Effective Scalar Field Theory for the Electroweak Phase Transition

abstract: We investigate an effective model for the finite temperature restoration phase transition of the electroweak theory. It is obtained by dimensional reduction of the $3+1$ dimensional full theory and by subsequent integration over all static gauge degrees of freedom. The resulting theory corresponds to a $3$-dimensional $O(4)$ ferromagnet containing cubic and quartic terms of the field in its potential function. Possible nonperturbative effects of a magnetic screening mass are parametrically included in the potential. We analyse the theory using mean field and numerical Monte Carlo (MC) simulation methods. At the value of the physical Higgs mass, $m_H=37~{\rm GeV}$, considered in the present investigation, we find a discontinuous symmetry restoring phase transition. We determine the critical temperature, order parameter jump, interface tension and latent heat characteristics of the transition. The Monte Carlo results indicate a somewhat weaker first order phase transition as compared to the mean field treatment, demonstrating that non-perturbative fluctuations of the Higgs field are relevant. This effect is especially important for the interface tension. Any observation of hard first order transition could result only from non-perturbative effects related to the gauge degrees of freedom.


title: Multigrid Monte Carlo Algorithms for SU(2) Lattice Gauge Theory: Two versus Four Dimensions

abstract: We study a multigrid method for nonabelian lattice gauge theory, the time slice blocking, in two and four dimensions. For SU(2) gauge fields in two dimensions, critical slowing down is almost completely eliminated by this method. This result is in accordance with theoretical arguments based on the analysis of the scale dependence of acceptance rates for nonlocal Metropolis updates. The generalization of the time slice blocking to SU(2) in four dimensions is investigated analytically and by numerical simulations. Compared to two dimensions, the local disorder in the four dimensional gauge field leads to kinematical problems.


title: Asymptotic scaling from strong coupling

abstract: Strong-coupling analysis of two-dimensional chiral models, extended to 15th order, allows for the identification of a scaling region where known continuum results are reproduced with great accuracy, and asymptotic scaling predictions are fulfilled. The properties of the large-$N$ second-order phase transition are quantitatively investigated.


title: Fermion production despite fermion number conservation

abstract: Lattice proposals for a nonperturbative formulation of the Standard Model easily lead to a global U(1) symmetry corresponding to exactly conserved fermion number. The absence of an anomaly in the fermion current would then appear to inhibit anomalous processes, such as electroweak baryogenesis in the early universe. One way to circumvent this problem is to formulate the theory such that this U(1) symmetry is explicitly broken. However we argue that in the framework of spectral flow, fermion creation and annihilation still in fact occurs, despite the exact fermion number conservation. The crucial observation is that fermions are excitations relative to the vacuum, at the surface of the Dirac sea. The exact global U(1) symmetry prohibits a state from changing its fermion number during time evolution, however nothing prevents the fermionic ground state from doing so. We illustrate our reasoning with a model in two dimensions which has axial-vector couplings, first using a sharp momentum cutoff, then using the lattice regulator with staggered fermions. The difference in fermion number between the time evolved state and the ground state is indeed in agreement with the anomaly. A study of the vacuum energy shows that the perturbative counterterm needed for restoration of gauge invariance is insufficient in a nonperturbative setting. For reference we also study a closely related model with vector couplings, the Schwinger model, and we examine the emergence of the $\theta$-vacuum structure of both theories.


title: Vortex Loops In The 3-d XY Model

abstract: We study the properties of vortex loops in the 3-d XY model. We find that the phase transition in this model is driven by the condensation of short loops into percolating loops. Paper comes as a uuencoded compressed tar-file, containing the complete PostScript file. Total file approx. 220 Kb.


title: Precision Upsilon Spectroscopy from Nonrelativistic Lattice QCD

abstract: The spectrum of the Upsilon system is investigated using the Nonrelativistic Lattice QCD approach to heavy quarks and ignoring light quark vacuum polarization. We find good agreement with experiment for the Upsilon(1S), Upsilon(2S), Upsilon(3S) and for the center of mass and fine structure of the chi_b states. The lattice calculations predict b-bbar D-states with center of mass at (10.20 +/- 0.07 +/- 0.03)GeV. Fitting procedures aimed at extracting both ground and excited state energies are developed. We calculate a nonperturbative dispersion mass for the Upsilon(1S) and compare with tadpole-improved lattice perturbation theory.


title: Failure of the Regge approach in two dimensional quantum gravity

abstract: Regge's method for regularizing euclidean quantum gravity is applied to two dimensional gravity. We use two different strategies to simulate the Regge path integral at a fixed value of the total area: A standard Metropolis simulation combined with a histogramming method and a direct simulation using a Hybrid Monte Carlo algorithm. Using topologies with genus zero and two and a scale invariant integration measure, we show that the Regge method does not reproduce the value of the string susceptibility of the continuum model. We show that the string susceptibility depends strongly on the choice of the measure in the path integral. We argue that the failure of the Regge approach is due to spurious contributions of reparametrization degrees of freedom to the path integral.


title: Canonical Demon Monte Carlo Renormalization Group

abstract: We describe a new method to compute renormalized coupling constants in a Monte Carlo renormalization group calculation. The method can be used for a general class of models, e.g., lattice spin or gauge models. The basic idea is to simulate a joint system of block spins and canonical demons. In contrast to the Microcanonical Renormalization Group invented by Creutz et al. our method does not suffer from systematical errors stemming from a simultaneous use of two different ensembles. We present numerical results for the $O(3)$ nonlinear $\sigma$-model.


title: The critical behaviour of Ising spins on 2D Regge lattices

abstract: We performed a high statistics simulation of Ising spins coupled to 2D quantum gravity on toroidal geometries. The tori were triangulated using the Regge calculus approach and contained up to $512^2$ vertices. We used a constant area ensemble with an added $R^2$ interaction term, employing the $dl/l$ measure. We find clear evidence that the critical exponents of the Ising phase transition are consistent with the static critical exponents and do not depend on the coupling strength of the $R^2$ interaction term. We definitively can exclude for this type of model a behaviour as predicted by Boulatov and Kazakov [Phys. Lett. {\bf B186}, 379 (1987)] for Ising spins coupled to dynamically triangulated surfaces.


title: Distribution of the color fields around static quarks: Flux tube profiles

abstract: We report detailed calculations of the profiles of energy and action densities in the quark-antiquark string in SU(2) lattice gauge theory.


title: Phase of the Wilson Line

abstract: This paper discusses the global $Z(N)$ symmetry of finite-temperature, $SU(N)$, pure Yang-Mills lattice gauge theory and the physics of the phase of the Wilson line expectation value. In the high $T$ phase, $\langle L \rangle$ takes one of $N$ distinct values proportional to the $Nth$ roots of unity in $Z(N)$, and the $Z(N)$ symmetry is broken. Only one of these is consistent with the usual interpretation $\langle L \rangle = e^{-F/T}$. This relation should be generalized to $\langle L \rangle = z e^{-F/T}$ with $z \in Z(N)$ so that it is consistent with the negative or complex values. In the Hamiltonian description, the {\em physical} variables are the group elements on the links of the spatial lattice. In a Lagrangian formulation, there are also group elements on links in the inverse-temperature direction from which the Wilson line is constructed. These are unphysical, auxiliary variables introduced to enforce the Gauss law constraints. The following results are obtained: The relation $\langle L \rangle=ze^{-F/T}$ is derived. The value of $z \in Z(N)$ is determined by the external field that is needed for the infinite-volume limit. There is a single physical, high-temperature phase, which is the same for all $z$. The global $Z(N)$ symmetry is not physical; it acts as the identity on all physical states. In the Hamiltonian formulation, the high-temperature phase is not distinguished by physical broken symmetry but rather by percolating flux.


title: Rough Interfaces Beyond the Gaussian Approximation

abstract: We compare predictions of the Capillary Wave Model beyond its Gaussian approximation with Monte Carlo results for the energy gap and the surface energy of the 3D Ising model in the scaling region. Our study reveals that the finite size effects of these quantities are well described by the Capillary Wave Model, expanded to two--loop order ( one order beyond the Gaussian approximation).


title: Nonuniformity of the $1/N$ Expansion for Two-Dimensional $O(N)$ Models

abstract: We point out that the $1/N$ expansion, which is widely invoked to infer properties of the $2D$ $O(N)$ models, is nonuniform in the temperature, i.e. with decreasing temperature the $1/N$ expansion truncated at a fixed order deviates more and more from the true answer. This fact precludes the use of the expansion to deduce low temperature properties such as asymptotic scaling for those models. By contrast, in the $1D$ $O(N)$ chains, there are no signs of such a nonuniformity.


title: Large Energy Cumulants in the 2D Potts Model and their Effects in Finite Size Analysis

abstract: We develop an ansatz for expressing the free energy of the two dimensional $q$-states Potts model for $q > 4$ near its first order phase transition point. We notice that for the moderate values of $ q \lesssim 15 $, the energy profile at the phase transition is not expressible as a sum of gaussians. We discuss how this affects the traditional finite size analysis of this phase transition. In particular, the dominant length scale governing the finite size corrections turns out to be much (${} \sim 6$ ~times) larger than the largest correlation length in the problem.


title: Compact U(1) Gauge Theory on Lattices with Trivial Homotopy Group

abstract: We study the pure gauge model on a lattice manifold with trivial fundamental homotopy group, homotopically equivalent to an $S_4$. Monopole loops may fluctuate freely on that lattice without restrictions due to the boundary conditions. For the original Wilson action on the hypertorus there is an established two-state signal in energy distribution functions which disappears for the new geometry. Our finite size scaling analysis suggests stringent upper bounds on possible discontinuities in the plaquette action. However, no consistent asymptotic finite size scaling behaviour is observed.


title: Multibondic Cluster Algorithm for Monte Carlo Simulations of First-Order Phase Transitions

abstract: Inspired by the multicanonical approach to simulations of first-order phase transitions we propose for $q$-state Potts models a combination of cluster updates with reweighting of the bond configurations in the Fortuin-Kastelein-Swendsen-Wang representation of this model. Numerical tests for the two-dimensional models with $q=7, 10$ and $20$ show that the autocorrelation times of this algorithm grow with the system size $V$ as $\tau \propto V^\alpha$, where the exponent takes the optimal random walk value of $\alpha \approx 1$.


title: Investigation of Laplacian Gauge Fixing for U(1) and SU(2) Gauge Fields

abstract: The Laplacian gauge on the lattice is investigated numerically using U(1) and SU(2) gauge fields. The problem of Gribov ambiguities is addressed and to asses the smoothness of the gauge fixed configurations, they are compared to configurations fixed to the Landau gauge. The results of these comparisons with the Landau gauge indicate that Laplacian gauge fixing works very well in practice and offers a viable alternative to Landau gauge fixing.


title: Comment on A. Patrascioiu and E. Seiler's paper ``Nonuniformity of the $1/N$ Expansion for Two-Dimensional $O(N)$ Models''

abstract: We remind that the non-uniformity in the temperature of the $1/N$ expansion for dimensionful quantities pointed out in Patrascioiu and Seiler's paper is not only compatible with but is predicted by asymptotic freedom, and it is present in the continuum as well as on the lattice. The convergence of the $1/N$ expansion for adimensional quantities is perfectly uniform.


title: Compact Lattice QED with Staggered Fermions and Chiral Symmetry Breaking

abstract: Different formulations of the $4d$ compact lattice QED with staggered fermions (standard Wilson and modified by suppression of lattice artifacts) are investigated by Monte Carlo simulations within the quenched approximation. We show that after suppressing lattice artifacts the system undergoes a phase transition from the Coulomb phase into a presumably weakly chirally broken phase only at (unphysical) negative $\beta$--values.


title: Simulation of Field Theories in Wavelet Representation

abstract: The field is expanded in a wavelet series and the wavelet coefficients are varied in a simulation of the 2D $\phi^4$ field theory. The drastically reduced autocorrelations result in a substantial decrease of computing requirements, compared to those in local Metropolis simulations. A large part of the improvement is shown to be the result of an additional freedom in the choice of the allowed range of change at the Metropolis update of wavelet components, namely the range can be optimized independently for all wavelet sizes.


title: Scaling in the Positive Plaquette Model and Universality in SU(2) Lattice Gauge Theory

abstract: We investigate universality, scaling, the beta-function and the topological charge in the positive plaquette model for SU(2) lattice gauge theory. Comparing physical quantities, like the critical temperature, the string tension, glueball masses, and their ratios, we explore the effect of a complete suppression of a certain lattice artifact, namely the negative plaquettes, for SU(2) lattice gauge theory. Our result is that this modification does not change the continuum limit, i.e., the universality class. The positive plaquette model and the standard Wilson formulation describe the same physical situation. The approach to the continuum limit given by the beta-function in terms of the bare lattice coupling, however, is rather different: the beta-function of the positive plaquette model does not show a dip like the model with standard Wilson action.


title: An Exploratory Study of Nucleon-Nucleon Scattering Lengths in Lattice QCD

abstract: An exploratory study is made of the nucleon-nucleon $s$-wave scattering lengths in quenched lattice QCD with the Wilson quark action. The $\pi$-$N$ and $\pi$-$\pi$ scattering lengths are also calculated for comparison. The calculations are made with heavy quarks corresponding to $m_\pi/m_\rho\approx 0.73-0.95$. The results show that the $N$-$N$ system has an attractive force in both spin-singlet and triplet channels, with their scattering lengths significantly larger than those for the $\pi$-$N$ and $\pi$-$\pi$ cases, a trend which is qualitatively consistent with the experiment. Problems toward a more realistic calculation for light quarks are discussed.


title: Probing with Penguins: A lattice calculation of the branching ratio for some of the exclusive modes of $b \to s \gamma$

abstract: We calculate the leading-order matrix element for exclusive decays of $b \to s\gamma$ in the quenched approximation of lattice QCD on a $24^3\times48$ lattice at $\beta=6.2$, using an O(a)-improved fermion action. The matrix element is used to extract the on-shell form factor $T_1(q^2=0)$ for $B \to K^*\gamma$ and $B_s \to \phi\gamma$, using two different assumptions about the $q^2$ dependence of the form factors for these decays. For $B \to K^*\gamma$, $T_1(q^2=0)$ is determined to be $0.159^{+34}_{-33}\pm 0.067$ or $0.124^{+20}_{-18} \pm 0.022$ in the two cases. We find the results to be consistent (in the Standard Model) with the CLEO experimental branching ratio of $BR(B\to K^*\gamma) = (4.5 \pm 1.5 \pm 0.9) \times 10^{-5}$.


title: Curvature and scaling in 4D dynamical triangulation

abstract: We study the average number of simplices $N'(r)$ at geodesic distance $r$ in the dynamical triangulation model of euclidean quantum gravity in four dimensions. We use $N'(r)$ to explore definitions of curvature and of effective global dimension. An effective curvature $R_V$ goes from negative values for low $\kappa_2$ (the inverse bare Newton constant) to slightly positive values around the transition $\kappa_2^c$. Far above the transition $R_V$ is hard to compute. This $R_V$ depends on the distance scale involved and we therefore investigate a similar explicitly $r$ dependent `running' curvature $R_{\rm eff}(r)$. This increases from values of order $R_V$ at intermediate distances to very high values at short distances. A global dimension $d$ goes from high values in the region with low $\kappa_2$ to $d=2$ at high $\kappa_2$. At the transition $d$ is consistent with 4. We present evidence for scaling of $N'(r)$ and introduce a scaling dimension $d_s$ which turns out to be approximately 4 in both weak and strong coupling regions. We discuss possible implications of the results, the emergence of classical euclidean spacetime and a possible `triviality' of the theory.


title: Manifestation of Sea Quark Effects in the Strong Coupling Constant in Lattice QCD

abstract: We demonstrate that sea quark effects of a magnitude expected from renormalization group considerations are clearly visible in the strong coupling constant measured in current full QCD simulations. Building on this result an estimate of $\alpha_{\overline{MS}}^{(5)}(M_Z)$ is made employing the charmonium $1S-1P$ mass splitting calculated on full QCD configurations generated with two flavors of dynamical Kogut-Susskind quarks to fix the scale.


title: Fractal Structure in Two-Dimensional Quantum Regge Calculus

abstract: We study the fractal structure of the surface in two-dimensional quantum Regge calculus by performing Monte Carlo simulation with up to 200,000 triangles. The result can be compared with the universal scaling function obtained analytically in the continuum limit of dynamical triangulation, which provides us with a definite criterion whether Regge calculus serves as a proper regularization of quantum gravity. When the scale-invariant measure is taken as the measure of the link-length integration, we observe the correct scaling behavior in the data for the type of loop attached to a baby universe. The data seem to converge to the universal scaling function as the number of triangles is increased. The data for the type of loop attached to the mother universe, on the other hand, shows no scaling behavior up to the present size.


title: The Phases and Triviality of Scalar Quantum Electrodynamics

abstract: The phase diagram and critical behavior of scalar quantum electrodynamics are investigated using lattice gauge theory techniques. The lattice action fixes the length of the scalar (``Higgs'') field and treats the gauge field as non-compact. The phase diagram is two dimensional. No fine tuning or extrapolations are needed to study the theory's critical behovior. Two lines of second order phase transitions are discovered and the scaling laws for each are studied by finite size scaling methods on lattices ranging from $6^4$ through $24^4$. One line corresponds to monopole percolation and the other to a transition between a ``Higgs'' and a ``Coulomb'' phase, labelled by divergent specific heats. The lines of transitions cross in the interior of the phase diagram and appear to be unrelated. The monopole percolation transition has critical indices which are compatible with ordinary four dimensional percolation uneffected by interactions. Finite size scaling and histogram methods reveal that the specific heats on the ``Higgs-Coulomb'' transition line are well-fit by the hypothesis that scalar quantum electrodynamics is logarithmically trivial. The logarithms are measured in both finite size scaling of the specific heat peaks as a function of volume as well as in the coupling constant dependence of the specific heats measured on fixed but large lattices. The theory is seen to be qualitatively similar to $\lambda\phi^{4}$. The standard CRAY random number generator RANF proved to be inadequate


title: Random Geometries and Real Space Renormalization Group

abstract: A method of ``blocking'' triangulations that rests on the self-similarity feature of dynamically triangulated random manifolds is proposed. The method is used to define the renormalization group for random geometries. As an illustration, the idea is applied to pure euclidean quantum gravity in two dimensions. Generalisation to more complicated systems and to higher dimensionalities of space-time appears straightforward.


title: Deconfinement Transition in Large N Lattice Gauge Theory

abstract: We study analytically the phase diagram of the pure $SU(N)$ lattice gauge theory at finite temperature, and we attempt to estimate the critical deconfinement temperature. We apply large $N$ techniques to the Wilson and to the Heat Kernel action, and we study the resulting models both in the strong coupling and in the weak coupling limits. Using the Heat Kernel action, we establish an interesting connection between the Douglas-Kazakov phase transition of two-dimensional QCD and the deconfining transition in $d$ dimensions. The analytic results obtained for the critical temperature compare well with Montecarlo simulations of the full theory in $(2+1)$ and in $(3+1)$ dimensions.


title: Regge gravity on general triangulations

abstract: We investigate quantum gravity in four dimensions using the Regge approach on triangulations of the four-torus with general, non-regular incidence matrices. We find that the simplicial lattice tends to develop spikes for vertices with low coordination numbers even for vanishing gravitational coupling. Different to the regular, hypercubic lattices almost exclusively used in previous studies, we find now that the observables depend on the measure. Computations with nonvanishing gravitational coupling still reveal the existence of a region with well-defined expectation values. However, the phase structure depends on the triangulation. Even with additional higher- order terms in the action the critical behavior of the system changes with varying (local) coordination numbers.


title: Can Sigma Models Describe Finite Temperature Chiral Transitions?

abstract: Large-N expansions and computer simulations indicate that the universality class of the finite temperature chiral symmetry restoration transition in the 3D Gross-Neveu model is mean field theory. This is a counterexample to the standard 'sigma model' scenario which predicts the 2D Ising model universality class. We trace the breakdown of the standard scenario (dimensional reduction and universality) to the absence of canonical scalar fields in the model. We point out that our results could be generic for theories with dynamical symmetry breaking, such as Quantum Chromodynamics.


title: String tension and monopoles in $T \neq 0$ SU(2) QCD

abstract: Monopole and photon contributions to abelian Wilson loops are calculated using Monte-Carlo simulations of finite-temperature $SU(2)$ QCD in the maximally abelian gauge. Long monopole loops alone are responsible for the behavior of the string tension in the confinement phase up to the critical $\beta_c$. Short monopole loops and photons do not contribute to the string tension. The abelian and the monopole spacial string tensions (both of which agree with the normal ones for $\beta < \beta_c$) show a $g^{4}(T) T^2$ scaling behavior in the deconfinement phase. The abelian spacial string tension is in agreement with the full one even in the deconfinement phase.


title: Low-Temperature Series for the Correlation Length in $d=3$ Ising Model

abstract: We extend low-temperature series for the second moment of the correlation function in $d=3$ simple-cubic Ising model from $u^{15}$ to $u^{26}$ using finite-lattice method, and combining with the series for the susceptibility we obtain the low-temperature series for the second-moment correlation length to $u^{23}$. An analysis of the obtained series by inhomogeneous differential approximants gives critical exponents $ 2\nu^{\prime} + \gamma^{\prime} \approx 2.55 $ and $ 2\nu^{\prime} \approx 1.27 $.


title: Ising Spins on Thin Graphs

abstract: The Ising model on ``thin'' graphs (standard Feynman diagrams) displays several interesting properties. For ferromagnetic couplings there is a mean field phase transition at the corresponding Bethe lattice transition point. For antiferromagnetic couplings the replica trick gives some evidence for a spin glass phase. In this paper we investigate both the ferromagnetic and antiferromagnetic models with the aid of simulations. We confirm the Bethe lattice values of the critical points for the ferromagnetic model on $\phi^3$ and $\phi^4$ graphs and examine the putative spin glass phase in the antiferromagnetic model by looking at the overlap between replicas in a quenched ensemble of graphs. We also compare the Ising results with those for higher state Potts models and Ising models on ``fat'' graphs, such as those used in 2D gravity simulations.


title: Properties of B-Mesons in Lattice QCD

abstract: The results of an extensive study of B-meson properties in quenched lattice QCD are presented. The studies are carried out in the static quark limit where the b-quark is taken to be infinitely massive. Our computations rely on a multistate smearing method introduced previously, with smearing functions generated from a relativistic lattice quark model. Systematic errors arising from excited state contamination, finite volume effects, and the chiral extrapolation for the light quarks are estimated. We obtain continuum results for the mass splitting M_{B_s}- M_{B_u} = 86 (+/-)12(stat) {+7/-9}(syst) MeV, the ratio of decay constants f_{B_s}/f_{B_u} = 1.22 (+/-)0.04(stat) (+/-)0.02 (syst). For the B-meson decay constant we separately exhibit the sizable uncertainties in the extrapolation to the continuum limit a -> 0 and higher order perturbative matching. We obtain f_{B} = 188 (+/-)23(stat) (+/-)15(syst) {+26/-0}(extrap) (+/-)14 (pert) MeV. ----- [Postscript version of paper available by anonymous ftp at fncrd6.fnal.gov. The file is fb.ps in subdirectory theory.]


title: Heating and small-size instantons in the $O(3) \:\sigma$ model on the lattice

abstract: We study the role of small-size instantons in the determination of the topological susceptibility of the 2-d $O(3) \: \sigma $ model on the lattice. In particular, we analyze how they affect the non-perturbative determination, by Monte Carlo techniques, of the renormalizations on the lattice. As a result, we obtain a high-precision non-perturbative determination of the mixing with the unity operator, finding good agreement with perturbative computations. We also obtain the size distribution of instantons in the physical vacuum up to very small values of the size in physical units, without observing any ultraviolet cut-off. Moreover, we show by analytical calculation that the mixing of the topological susceptibility with the action density is a negligible part of the whole non-perturbative signal.


title: Beautiful Baryons from Lattice QCD

abstract: We perform a lattice study of heavy baryons, containing one ($\Lambda_b$) or two $b$-quarks ($\Xi_b$). Using the quenched approximation we obtain for the mass of $\Lambda_b$ $$ M_{\Lambda_b}= 5.728 \pm 0.144 \pm 0.018 {\rm GeV}.$$ The mass splitting between the $\Lambda_b$ and the B-meson is found to increase by about 20\% if the light quark mass is varied from the chiral limit to the strange quark mass.


title: Non perturbative determination of the running coupling constant in quenched SU(2)

abstract: Through a finite size renormalization group technique we calculate the running coupling constant for quenched SU(2) with a few percent error over a range of energy varying by a factor thirty. The definition is based on ratio of correlations of Polyakov loops with twisted boundary conditions. The extrapolation to the continuum limit is governed by corrections due to lattice artifacts which are proportional to the square of the lattice spacing and appears rather smooth.


title: Supersymmetric corrections to $\epsilon^\prime /\epsilon at the leading order in QCD and QED

abstract: We study the corrections to $\epsilon^\prime /\epsilon$ in the minimal supersymmetric model at the leading order in QCD and QED. Supersymmetry can increase the standard model prediction for $\epsilon^\prime /\epsilon$ by at most 40\% for $m_t=174$ GeV, an enhancement which is indistinguishable from the present theoretical uncertainties. The most conspicuous effect of supersymmetry is a strong depletion of $\epsilon^\prime /\epsilon$. For certain choices of supersymmetric parameters, vanishing and even small negative values of $\epsilon^\prime /\epsilon$ can be obtained for the top quark in the CDF range.


title: Reply to M.Campostrini's and P.Rossi's Comment on our paper `Nonuniformity of the $1/N$ Expansion for Two-Dimensional $O(N)$ Models'

abstract: This reply tries to rectify some misunderstandings that are in our opinion contained in the Comment by Campostrini and Rossi, <hep-lat 99407008> on our paper <hep-lat 9407003>.


title: Chiral perturbation theory, finite size effects and the three-dimensional $XY$ model

abstract: We study finite size effects of the d=3 $XY$ model in terms of the chiral perturbation theory. We calculate by Monte Carlo simulations physical quantities which are, to order of $(1/L)^2$, uniquely determined only by two low energy constants. They are the magnetization and the helicity modulus (or the Goldstone boson decay constant) in infinite volume. We also pay a special attention to the region of the validity of the two possible expansions in the theory.


title: Pion-Nucleon Sigma Term in Lattice QCD

abstract: We calculate both connected and disconnected contribution to the $\pi$-$N$ $\sigma$-term in quenched lattice QCD with Wilson quark action on a $12^3\times 20$ lattice at $\beta=5.7$. The latter is evaluated with the aid of the variant wall source method, which was previously applied successfully for extraction of $\pi$-$\pi$ scattering lengths and $\eta^\prime$ meson mass. We found that the disconnected contribution is about twice larger than the connected one. The value for the full $\pi$-$N$ $\sigma$ term $\sigma=40-60$ MeV is consistent with the experimental estimates. The nucleon matrix element of the strange quark density $\bar s s$ is fairly large in our result.


title: Polyakov loops and monopoles in QCD

abstract: Monte-Carlo simulations of abelian projection of $T \neq 0$ pure lattice QCD show that 1) Polyakov loops written in terms of abelian link fields alone play a role of an order parameter of deconfinement transition, 2) the abelian Polyakov loops are decomposed into contributions from Dirac strings of monopoles and from photons, 3) vanishing of the abelian Polyakov loops in the confinement phase is due to the Dirac strings alone and the photons give a finite contribution in both phases. Moreover, these results appear to hold good with any abelian projection as seen from the studies in the maximally abelian gauge and in various unitary gauges.


title: Monopole action and condensation in SU(2) QCD

abstract: An effective monopole action for various extended monopoles is derived from vacuum configurations after abelian projection in the maximally abelian gauge in $SU(2)$ QCD. The action appears to be independent of the lattice volume. Moreover it seems to depend only on the physical lattice spacing of the renormalized lattice, not on $\beta$. Entropy dominance over energy of monopole loops is seen on the renormalized lattice with the spacing $b>b_c\simeq 5.2\times10^{-3} \Lambda_L^{-1}$. This suggests that monopole condensation always (for all $\beta$) occurs in the infinite-volume limit of lattice QCD.


title: The Sine Gordon Model: Perturbation Theory and Cluster Monte Carlo

abstract: We study the expansion of the surface thickness in the 2-dimensional lattice Sine Gordon model in powers of the fugacity z. Using the expansion to order z**2, we derive lines of constant physics in the rough phase. We describe and test a VMR cluster algorithm for the Monte Carlo simulation of the model. The algorithm shows nearly no critical slowing down. We apply the algorithm in a comparison of our perturbative results with Monte Carlo data.


title: A Lattice Calculation of the Heavy Quark Universal Form Factor

abstract: A preliminary computation of the Isgur-Wise universal form factor using a lattice formulation of the Heavy Quark Effective Theory (HQET) is described and compared with the recent data from ARGUS and CLEO on B -> D*lv decay.


title: Nucleon Structure from Lattice QCD

abstract: We report results on the nucleon structure obtained from the lattice quantum chromodynamics calculation. These include the axial, electromagnetic, $\pi NN$, and scalar form factors. The calculation is carried out at $\beta = 6$ on a $16^3 \times 24$ lattice with 24 quenched gauge configurations. The chiral limit results are extrapolated from several light quark cases. For the disconnected insertion (sea-quark contribution), we used the stochastic estimation with the $Z_2$ noise to calculate the diagonal and off-diagonal traces of the inverse matrices with a size of $10^6 \times 10^6$. It is found that the $Z_2$ noise is the optimal choice and its comparison with the Gaussian noise for our quark matrix is given. For the sea-quark contribution, we report results on the strange condensate in the nucleon and the $\pi N \sigma$ term.


title: Non-Perturbative Renormalization Group Flows in Two-Dimensional Quantum Gravity

abstract: Recently a block spin renormalization group approach was proposed for the dynamical triangulation formulation of two-dimensional quantum gravity. We use this approach to examine non-perturbatively a particular class of higher derivative actions for pure gravity.


title: Scaling Properties of the Energy Density in SU(2) Lattice Gauge Theory

abstract: The lattice data for the energy density of $SU(2)$ gauge theory are calculated with \nop~derivatives of the coupling constants. These derivatives are obtained from two sources : i) a parametrization of the \nop~beta function in accord with the measured critical temperature and $\Delta\beta-$values and ii) a \nop~calculation of the presssure. We then perform a detailed finite size scaling analysis of the energy density near $T_c$. It is shown that at the critical temperature the energy density is scaling as a function of $VT^3$ with the corresponding $3d$ Ising model critical exponents. The value of $\epsilon(T_c)/T^4_c$ in the continuum limit is estimated to be 0.256(23). In the high temperature regime the energy density is approaching its weak coupling limit from below, at $T/T_c \approx 2$ it has reached only about $70\%$ of the limit.


title: "Confinement Mechanism in Various Abelian Projections of $SU(2)$ Lattice Gluodynamics"

abstract: We show that the monopole confinement mechanism in lattice gluodynamics is a particular feature of the maximal abelian projection. We give an explicit example of the $SU(2) \rightarrow U(1)$ projection (the minimal abelian projection), in which the confinement is due to topological objects other than monopoles. We perform analytical and numerical study of the loop expansion of the Faddeev--Popov determinant for the maximal and the minimal abelian projections, and discuss the fundamental modular region for these projections.


title: Tricritical Phenomena in a Z(3) Lattice Gauge Theory

abstract: The Z(3) gauge model with double plaquette representation of the action on a generalized Bethe lattice of plaquettes is constructed. It is reduced to the spin-1 Blume-Emery-Griffiths (BEG) model. An Ising-type critical line of a second-order phase transition ending in the tricritical point is found.


title: Interrelation of a Z(3) Gauge Theory on the Flat Lattices and a Spin-1 BEG Model

abstract: The Z(3) gauge model with double plaquette representation of the action on the flat triangular and square lattices is constructed. It is reduced to the spin-1 Blume-Emery-Griffiths (BEG) model. An Ising-type critical line of a second-order phase transition is found.


title: Spectrum of the Dirac Operator and Multigrid Algorithm with Dynamical Staggered Fermions

abstract: Complete spectra of the staggered Dirac operator $\Dirac$ are determined in quenched four-dimensional $SU(2)$ gauge fields, and also in the presence of dynamical fermions. Periodic as well as antiperiodic boundary conditions are used. An attempt is made to relate the performance of multigrid (MG) and conjugate gradient (CG) algorithms for propagators with the distribution of the eigenvalues of~$\Dirac$. The convergence of the CG algorithm is determined only by the condition number~$\kappa$ and by the lattice size. Since~$\kappa$'s do not vary significantly when quarks become dynamic, CG convergence in unquenched fields can be predicted from quenched simulations. On the other hand, MG convergence is not affected by~$\kappa$ but depends on the spectrum in a more subtle way.


title: Compact lattice QED with Wilson fermions

abstract: We study the phase structure and the chiral limit of $4d$ compact lattice QED with Wilson fermions (both dynamical and quenched). We use the standard Wilson gauge action and also a modified one suppressing lattice artifacts. Different techniques and observables to locate the chiral limit are discussed.


title: Entropy and the Approach to the Thermodynamic Limit in Three-Dimensional Simplicial Gravity

abstract: We present numerical results supporting the existence of an exponential bound in the dynamical triangulation model of three-dimensional quantum gravity.Both the critical coupling and various other quantities show a slow power law approach to the infinite volume limit.


title: The (LATTICE) QCD Potential and Running Coupling: How to Accurately Interpolate between Multi-Loop QCD and the String Picture

abstract: We present a simple parameterization of a running coupling constant, defined via the static potential, that interpolates between 2-loop QCD in the UV and the string prediction in the IR. Besides the usual $\Lam$-parameter and the string tension, the coupling depends on one dimensionless parameter, determining how fast the crossover from UV to IR behavior occurs (in principle we know how to take into account any number of loops by adding more parameters). Using a new Ansatz for the LATTICE potential in terms of the continuum coupling, we can fit quenched and unquenched Monte Carlo results for the potential down to ONE lattice spacing, and at the same time extract the running coupling to high precision. We compare our Ansatz with 1-loop results for the lattice potential, and use the coupling from our fits to quantitatively check the accuracy of 2-loop evolution, compare with the Lepage-Mackenzie estimate of the coupling extracted from the plaquette, and determine Sommer's scale $r_0$ much more accurately than previously possible. For pure SU(3) we find that the coupling scales on the percent level for $\beta\geq 6$.


title: Topological Zero Modes in Monte Carlo Simulations

abstract: We present an improvement of global Metropolis updating steps, the instanton hits, used in a hybrid Monte Carlo simulation of the two-flavor Schwinger model with staggered fermions. These hits are designed to change the topological sector of the gauge field. In order to match these hits to an unquenched simulation with pseudofermions, the approximate zero mode structure of the lattice Dirac operator has to be considered explicitly.


title: Lee-Yang zeroes in the one flavour massive lattice Schwinger model

abstract: We study the partition function of the model formulated with Wilson fermions with only one species, both analytically and numerically. At strong coupling we construct the solution for lattice size up to $8\times 8$, a polynomial in the hopping parameter up to $O(\ka^{128})$. At $\be>0$ we evaluate the expectation value of the fermion determinant for complex values of $\ka$. From the Lee-Yang zeroes we find support for the existence of a line of phase transitions from $(\be=0, \ka\simeq 0.38)$ up to $(\be=\infty, \ka=1/4)$.


title: An Improved Estimator for the Correlation Function of 2D Nonlinear Sigma Models

abstract: I present a new improved estimator for the correlation function of 2D nonlinear sigma models. Numerical tests for the 2D XY model and the 2D O(3)-invariant vector model were performed. For small physical volume, i.e. a lattice size small compared to the to the bulk correlation length, a reduction of the statistical error of the finite system correlation length by a factor of up to 30 compared to the cluster-improved estimator was observed. This improvement allows for a very accurate determination of the running coupling proposed by M. L"uscher et al. for 2D O(N)-invariant vector models.


title: Complex-Temperature Singularities of the Susceptibility in the $d=2$ Ising Model. I. Square Lattice

abstract: We investigate the complex-temperature singularities of the susceptibility of the 2D Ising model on a square lattice. From an analysis of low-temperature series expansions, we find evidence that as one approaches the point $u=u_s=-1$ (where $u=e^{-4K}$) from within the complex extensions of the FM or AFM phases, the susceptibility has a divergent singularity of the form $\chi \sim A_s'(1+u)^{-\gamma_s'}$ with exponent $\gamma_s'=3/2$. The critical amplitude $A_s'$ is calculated. Other critical exponents are found to be $\alpha_s'=\alpha_s=0$ and $\beta_s=1/4$, so that the scaling relation $\alpha_s'+2\beta_s+\gamma_s'=2$ is satisfied. However, using exact results for $\beta_s$ on the square, triangular, and honeycomb lattices, we show that universality is violated at this singularity: $\beta_s$ is lattice-dependent. Finally, from an analysis of spin-spin correlation functions, we demonstrate that the correlation length and hence susceptibility are finite as one approaches the point $u=-1$ from within the symmetric phase. This is confirmed by an explicit study of high-temperature series expansions.


title: The High-Temperature Expansion of the Hierarchical Ising Model: From Poincar\'e Symmetry to an Algebraic Algorithm

abstract: We show that the hierarchical model at finite volume has a symmetry group which can be decomposed into rotations and translations as the familiar Poincar\'e groups. Using these symmetries, we show that the intricate sums appearing in the calculation of the high-temperature expansion of the magnetic susceptibility can be performed, at least up to the fourth order, using elementary algebraic manipulations which can be implemented with a computer. These symmetries appear more clearly if we use the 2-adic fractions to label the sites. We then apply the new algebraic methods to the calculation of quantities having a random walk interpretation. In particular, we show that the probability of returning at the starting point after $m$ steps has poles at $D=-2,-4,....-2m$ , where $D$ is a free parameter playing a role similar to the dimensionality in nearest neighbor models.


title: Noncompact Gauge-Invariant Lattice Simulations

abstract: We have applied a new gauge-invariant, noncompact, Monte Carlo method to simulate $U(1)$, $SU(2)$, and $SU(3)$ gauge theories on $8^4$ and $12^4$ lattices. The Creutz ratios of the Wilson loops agree with the exact results for $U(1)$ for $\beta \ge 1$ apart from a renormalization of the charge. The $SU(2)$ and $SU(3)$ Creutz ratios robustly display quark confinement at $\beta = 0.5$ and $\beta = 1.5$, respectively. At much weaker coupling, the $SU(2)$ and $SU(3)$ Creutz ratios agree with perturbation theory after a renormalization of the coupling constant.


title: Critical exponents of a three dimensional O(4) spin model

abstract: By Monte Carlo simulation we study the critical exponents governing the transition of the three-dimensional classical O(4) Heisenberg model, which is considered to be in the same universality class as the finite-temperature QCD with massless two flavors. We use the single cluster algorithm and the histogram reweighting technique to obtain observables at the critical temperature. After estimating an accurate value of the inverse critical temperature $\Kc=0.9360(1)$, we make non-perturbative estimates for various critical exponents by finite-size scaling analysis. They are in excellent agreement with those obtained with the $4-\epsilon$ expansion method with errors reduced to about halves of them.


title: Subtleties and Fancies in Gauge Theory Non Trivial Vacuum

abstract: The one loop effective potential for a non-Abelian gauge configuration is analyzed using the background field method. The Savvidy result and the non-Abelian ansatz, the other alternative possible background that generates a constant color magnetic field configuration, are compared. This second possibility is very interesting because it avoids the possible coordinate singularity, ${\rm Det}B_i^a=0$, and it is easy to implement in lattice simulations. We emphasize the interesting dependence of the potential by the gauge fixing parameter $\alpha$, when the loop expansion is performed around a non trivial background configuration. Finally, we point out some crucial differences in analyzing the vacuum structure between non-Abelian gauge theories and the cases of scalar and Abelian gauge theories.


title: Critical Exponents, Hyperscaling and Universal Amplitude Ratios for Two- and Three-Dimensional Self-Avoiding Walks

abstract: We make a high-precision Monte Carlo study of two- and three-dimensional self-avoiding walks (SAWs) of length up to 80000 steps, using the pivot algorithm and the Karp-Luby algorithm. We study the critical exponents $\nu$ and $2\Delta_4 -\gamma$ as well as several universal amplitude ratios; in particular, we make an extremely sensitive test of the hyperscaling relation $d\nu = 2\Delta_4 -\gamma$. In two dimensions, we confirm the predicted exponent $\nu = 3/4$ and the hyperscaling relation; we estimate the universal ratios $\<R_g^2\> / \<R_e^2\> = 0.14026 \pm 0.00007$, $\<R_m^2\> / \<R_e^2\> = 0.43961 \pm 0.00034$ and $\Psi^* = 0.66296 \pm 0.00043$ (68\% confidence limits). In three dimensions, we estimate $\nu = 0.5877 \pm 0.0006$ with a correction-to-scaling exponent $\Delta_1 = 0.56 \pm 0.03$ (subjective 68\% confidence limits). This value for $\nu$ agrees excellently with the field-theoretic renormalization-group prediction, but there is some discrepancy for $\Delta_1$. Earlier Monte Carlo estimates of $\nu$, which were $\approx\! 0.592$, are now seen to be biased by corrections to scaling. We estimate the universal ratios $\<R_g^2\> / \<R_e^2\> = 0.1599 \pm 0.0002$ and $\Psi^* = 0.2471 \pm 0.0003$; since $\Psi^* > 0$, hyperscaling holds. The approach to $\Psi^*$ is from above, contrary to the prediction of the two-parameter renormalization-group theory. We critically reexamine this theory, and explain where the error lies.


title: Finite-Size Scaling at $\xi/L \gg 1$

abstract: We present a simple and powerful method for extrapolating finite-volume Monte Carlo data to infinite volume, based on finite-size-scaling theory. We discuss carefully its systematic and statistical errors, and we illustrate it using three examples: the two-dimensional three-state Potts antiferromagnet on the square lattice, and the two-dimensional $O(3)$ and $O(\infty)$ $\sigma$-models. In favorable cases it is possible to obtain reliable extrapolations (errors of a few percent) even when the correlation length is 1000 times larger than the lattice.


title: Observing Long Colour Flux Tubes in SU(2) Lattice Gauge Theory

abstract: We present results of a high statistics study of the chromo field distribution between static quarks in SU(2) gauge theory on lattices of volumes 16^4, 32^4, and 48^3*64, with physical extent ranging from 1.3 fm up to 2.7 fm at beta=2.5, beta=2.635, and beta=2.74. We establish string formation over physical distances as large as 2 fm. The results are tested against Michael's sum rules. A detailed investigation of the transverse action and energy flux tube profiles is provided. As a by-product, we obtain the static lattice potential in unpreceded accuracy.


title: An Analysis of Four-quark Energies in SU(2) Lattice Monte Carlo using the Flux-tube Symmetry:

abstract: Energies of four-quark systems calculated by the static quenched SU(2) lattice Monte Carlo method are analyzed in $2\times 2$ bases for square, rectangle, tilted rectangle, linear and quadrilateral geometry configurations and in $3\times 3$ bases for a non-planar geometry configuration. For small interquark distances, a lattice effect is taken into account by considering perimeter dependent terms which are characterized by the cubic symmetry. It is then found that a parameter $f$ - that can be identified as a gluon field overlap factor - is rather well described by the form $exp(-[b_sE{\cal A}+\sqrt{b_s}F{\cal P}])$, where ${\cal A}$ and ${\cal P}$ are the area and perimeter mainly defined by the positions of the four quarks, $b_s$ is the string constant in the 2-quark potentials and $E,F$ are constants.


title: Core -- a New Method for Solving Hamiltonian Lattice Systems

abstract: The COntractor REnormalization group (CORE) approximation, a new method for solving Hamiltonian lattice systems, is introduced. The approach combines variational and contraction techniques with the real-space renormalization group approach and is systematically improvable. Since it applies to lattice systems of infinite extent, the method is suitable for studying critical phenomena and phase structure; systems with dynamical fermions can also be treated. The method is tested using the 1+1-dimensional Ising model.


title: Multigrid Methods for Propagators in Lattice Gauge Theories

abstract: Multigrid methods were invented for the solution of discretized partial differential equations in ordered systems. The slowness of traditional algorithms is overcome by updates on various length scales. In this article we discuss generalizations of multigrid methods for disordered systems, in particular for propagators in lattice gauge theories. A discretized nonabelian gauge theory can be formulated as a system of statistical mechanics where the gauge field degrees of freedom are $SU(N)$ matrices on the links of the lattice. These $SU(N)$ matrices appear as random coefficients in Dirac equations. We aim at finding an efficient method by which one can solve Dirac equations without critical slowing down. If this could be achieved, Monte Carlo simulations of Quantum Chromodynamics (the theory of the strong interaction) would be accelerated considerably. In principle, however, the methods discussed can be used in arbitrary space-time dimension and for arbitrary gauge group. Moreover, there are applications in multigrid Monte Carlo simulations, and for the definition of block spins and blocked gauge fields in Monte Carlo renormalization group studies. As a central results it was found that {\em the geometric multigrid method works in principle in arbitrarily disordered gauge fields.} Finally, an overview is given of other approaches to the propagator problem in lattice gauge theories.


title: Renormalization Group And Scaling Within The Microcanonical Fermionic Average Approach

abstract: The MFA approach for simulations with dynamical fermions in lattice gauge theories allows in principle to explore the parameters space of the theory (e.g. the $\beta, m$ plane for the study of chiral condensate in QED) without the need of computing the fermionic determinant at each point. We exploit this possibility for extracting both the renormalization group trajectories ("constant physics lines") and the scaling function, and we test it in the Schwinger Model. We discuss the applicability of this method to realistic theories.


title: Exploring Lattice Quantum Chromodynamics by Cooling

abstract: The effect of cooling on a number of observables is calculated in SU(2) lattice gauge theory. The static quark-antiquark potential and spin-dependent interactions are studied, and the topological charge is monitored. The chiral symmetry breaking order parameter $\langle \overline{\chi}\chi \rangle$ and meson correlators are calculated using staggered fermions. Interactions on the distance scale of a few lattice spacings are found to be essentially eliminated by cooling, including the spin-dependent potentials. $\langle \overline{\chi}\chi \rangle$ and meson correlators up to time separations of several lattice spacings relax very quickly to their free-field values. At larger times, there is evidence of a difference between the pseudoscalar and vector channels. A fit to the pseudoscalar correlation function yields ``mass'' values about $2/3$ (in lattice units) of the uncooled masses. These results raise the question of how to reconcile the large-time behavior of the hadron correlators with the fact that the spin-dependent potentials and $\langle \overline{\chi}\chi \rangle$ essentially disappear (in lattice units) after only a small amount of cooling.


title: Boundary condition for Staggered Fermion in Lattice Schr\"odinger Functional of QCD

abstract: The fermionic part of the Schr\"odinger functional of QCD is formulated in the lattice regularization with the staggered fermion. The boundary condition imposed on the staggered fermion field are examined in terms of the four-component Dirac spinor. The boundary terms are different from those of the Symanzik's theory in the flavor structure due to the species doubling. It is argued that, in the case of the homogeneous Dirichlet boundary condition, surface divergence does not occur if the link variables of gauge field are introduced on the original lattice, not on the blocked one. Its application to the numerical calculation of the running coupling constant in QCD is discussed.


title: Gribov Noise on the Lattice Axial Current Renormalisation Constant

abstract: We study the influence of Gribov copies, in the Landau gauge, on the lattice renormalisation constant of the axial current, obtained from a Ward identity on quark state correlation functions, with the Clover action, in quenched SU(3) gauge theory. A comparison between the gauge invariant determination of the renorma= lization constant and the gauge dependent one is discussed. We find that the residual gauge freedom associated to Gribov copies induces observable effects, which, at the level of numerical accuracy of our simulation, are included in the statistical uncertainty inherent in a Monte Carlo simulation. Doubling the statistics suggests that the fluctuation due to the lattice Gribov ambiguity scales down at least as fast as a pure statistical error.


title: Domain wall fermions in a waveguide: the phase diagram at large Yukawa coupling

abstract: In this paper we return to a model with domain wall fermions in a waveguide. This model contains a Yukawa coupling $y$ which is needed for gauge invariance. A previous paper left the analysis for large values of this coupling incomplete. We fill the gap by developing a systematic expansion suitable for large $y$, and using this, we gain an analytic understanding of the phase diagram and fermion spectrum. We find that in a sense all the species doublers come back for large $y$. The conclusion that no lattice chiral gauge theory can be obtained from this approach therefore remains valid.


title: Generalized Hybrid Monte-Carlo

abstract: We propose a modification of the Hybrid Monte-Carlo method to sample equilibrium distributions of continuous field models. The method allows an efficient implementation of Fourier acceleration and is shown to reduce completely critical slowing down for the Gaussian model, i. e., $z=0$.


title: Low--Temperature Series for Renormalized Operators: the Ferromagnetic Square--Lattice Ising Model.

abstract: A method for computing low--temperature series for renormalized operators in the two--dimensional Ising model is proposed. These series are applied to the study of the properties of the truncated renormalized Hamiltonians when we start at very low temperature and zero field. The truncated Hamiltonians for majority rule, Kadanoff transformation and decimation for $2 \times 2$ blocks depend on the how we approach the first--order phase--transition line. These Renormalization Group transformations are multi--valued and discontinuous at this first--order transition line when restricted to some finite--dimensional interaction space.


title: Pathological Behavior of Renormalization-Group Maps at High Fields and Above the Transition Temperature

abstract: We show that decimation transformations applied to high-$q$ Potts models result in non-Gibbsian measures even for temperatures higher than the transition temperature. We also show that majority transformations applied to the Ising model in a very strong field at low temperatures produce non-Gibbsian measures. This shows that pathological behavior of renormalization-group transformations is even more widespread than previous examples already suggested.


title: Simulating the Electroweak Phase Transition in the SU(2) Higgs Model

abstract: Numerical simulations are performed to study the finite temperature phase transition in the SU(2) Higgs model on the lattice. In the presently investigated range of the Higgs boson mass, below 50 GeV, the phase transition turns out to be of first order and its strength is rapidly decreasing with increasing Higgs boson mass. In order to control the systematic errors, we also perform studies of scaling violations and of finite volume effects.


title: Perfect Lattice Actions for the Gross-Neveu Model at large N

abstract: Fixed point actions for free and interacting staggered lattice fermions are constructed by iterating renormalization group transformations. At large N the fixed point action for the Gross-Neveu model is a perfect action in the sense of Hasenfratz and Niedermayer, i.e. cut-off effects are completely eliminated. In particular, the fermionic 1-particle energy spectrum of the lattice theory is identical with the one of the continuum even for arbitrarily small correlation lengths. The cut-off effects of the chiral condensate are eliminated using a perfect operator. (The paper is stored as a ps-file containing both the text and 5 figures.)


title: A Solution to the Strong CP Problem

abstract: One may argue that QCD solves the strong CP problem by itself, without having to introduce new symmetries and particles. To test this idea, a lattice simulation is performed. The problem is investigated in the CP$^3$ model first. It is found that the model has a first order phase transition in $\theta$ from a confining phase at small $\theta$ to a deconfining phase at large $\theta$, and that the critical value of $\theta$ decreases towards zero as $\beta$ is taken to infinity. This suggests that $\theta$ is tuned to zero in the continuum limit. Preliminary studies of the SU(2) Yang-Mills theory in four dimensions show a phase transition in $\theta$ as well, so that it is quite likely that the strong CP problem in QCD is solved along the same line.


title: Screening of the topological charge in a correlated instanton vacuum

abstract: Screening of the topological charge due to he fermion-induced interactions is an important phenomenon, closely related with the resolution of the strong CP and U(1) problems. We study the mechanism of such screening in a 'correlated instanton vacuum', as opposed to the 'random' one. Both scalar and pseudoscalar gluonic correlators are analyzed by means of an observable that minimizes finite size effects. Screening of the topological charge is established. This allows us to calculate the $\eta'$ mass without having to invert the Dirac operator. We suggest that this method might be used in lattice QCD calculations as well. Our results for the screening of the topological charge are in agreement with the chiral Ward identities, and the scalar gluonic correlator satisfies a low energy theorem first derived by Novikov et al. \cite{Novikov-etal}. We also propose to evaluate the topological susceptibility in the Witten-Veneziano formula not in an infinite box in an world $without$ fermions but in an infinitesimal box in a world $with$ fermions.


title: The Elusive Asymptotic Behavior of the High-Temperature Expansion of the Hierarchical Ising Model

abstract: We present a differential formulation of the recursion formula of the hierarchical model which provides a recursive method of calculation for the high-temperature expansion. We calculate the first 30 coefficients of the high temperature expansion of the magnetic susceptibility of the Ising hierarchical model with 12 significant digits. We study the departure from the approximation which consists of identifying the coefficients with the values they would take if a $[0,1]$ Pad\'e approximant were exact. We show that, when the order in the high-temperature expansion increases, the departure from this approximation grows more slowly than for nearest neighbor models. As a consequence, the value of the critical exponent $\gamma $ estimated using Pad\'e approximants converges very slowly and the estimations using 30 coefficients have errors larger than 0.05. A (presumably much) larger number of coefficients is necessary to obtain the critical exponents with a precision comparable to the precision obtained for nearest neighbor models with less coefficients. We also discuss the possibility of constructing models where a $[0,1]$ Pad\'e approximant would be exact.


title: Lattice $(\Phi^4)_4$ Effective Potential Giving Spontaneous Symmetry Breaking and the Role of the Higgs Mass

abstract: We present a critical reappraisal of the available results on the broken phase of $\lambda(\Phi^4)_4$ theory, as obtained from rigorous formal analyses and from lattice calculations. All the existing evidence is compatible with Spontaneous Symmetry Breaking but dictates a trivially free shifted field that becomes controlled by a quadratic hamiltonian in the continuum limit. As recently pointed out, this implies that the simple one-loop effective potential should become effectively exact. Moreover, the usual naive assumption that the Higgs mass-squared $m^2_h$ is proportional to its ``renormalized'' self-coupling $\lambda_R$ is not valid outside perturbation theory: the appropriate continuum limit has $m_h$ finite and vanishing $\lambda_R$. A Monte Carlo lattice computation of the $\lambda(\Phi^4)_4$ effective potential, both in the single-component and in the O(2)-symmetric cases, is shown to agree very well with the one-loop prediction. Moreover, its perturbative leading-log improvement (based on the concept of $\lambda_R$) fails to reproduce the Monte Carlo data. These results, while supporting in a new fashion the peculiar ``triviality'' of the $\lambda(\Phi^4)_4$ theory, also imply that, outside perturbation theory, the magnitude of the Higgs mass does not give a measure of the observable interactions in the scalar sector of the standard model.


title: Confining Forces and String Formation from the Lattice

abstract: We show the running coupling as derived from the SU(3) QQ potential and discuss preliminary results on spin dependent heavy quark potentials from high statistics lattice simulations of SU(2) gauge theory. The precision suffices to study scaling properties and lattice artifacts (at short distances). We identify the Coulomb like short range interaction as a mixed vector-scalar exchange. We measure flux tube formation between a static QQ pair over physical distances up to 2 fm, with a spatial resolution as small as .05 fm. Consistency with the string picture is found for separation larger than about 1 fm, with a half width of the profile of approximately .7 fm.


title: Real Space Renormalization Group Analysis of U(1)-Gauge Theory with theta-term in 2 Dimensions

abstract: U(1) lattice gauge theory with $\theta$-term is investigated by real space renormalization group approach. Flows of renormalized coupling constants are analyzed. For each $\theta$, renormalization flows converge to a single trajectory irrespective of bare coupling constants of real action. For $\theta \not= \pi$ the system is in the confinement phase controlled by an infrared fixed point. We found a phase transition at $\theta = \pi$. Imaginary part of the action given by $\theta$-term stays fixed under renormalization group transformations leading to deconfinement.


title: Universality and the Deconfinement Phase Transition in SU(2) Lattice Gauge Theory

abstract: We study the three dimensional fundamental-adjoint $SU(2)$ lattice gauge theory at non-zero temperatures by Monte Carlo simulations. On an $8^3 \times 2$ lattice, at $\beta_A = 1.1$, where $\beta_A$ is the adjoint coupling, we find no evidence of any transition at the location of a previously known bulk phase transition around $\beta = 1.33$. Moreover, the deconfinement transition at $ \beta_A = 1.1$ occurs at $\beta=1.20$ and is of first order for $\beta_A \ge 1.1$, thus implying a change of universality class from that of the Wilson action at $\beta_A=0$. Computations of the plaquette susceptibility and the temporal and spatial Polyakov loops on $8^3 \times 4$ and $16^3 \times 8$ lattices at $\beta_A = 1.1$ further support these conclusions and suggest that the previously claimed bulk transition around $\beta = 1.33$ is, in fact, the first order deconfinement transition. Simulations at larger $\beta_A$ and the measurements of the mass gaps from the correlation functions of temporal and spatial Polyakov loops also confirm the temperature dependent nature of the above transition. The consequences of our results on universality are discussed.


title: Low-Temperature Series Expansions for the Spin-1 Ising Model

abstract: The finite lattice method of series expansion has been used to extend low-temperature series for the partition function, order parameter and susceptibility of the spin-1 Ising model on the square lattice. A new formalism is described that uses two distinct transfer matrix approaches in order to significantly reduce computer memory requirements and which permits the derivation of the series to 79th order. Subsequent analysis of the series clearly confirms that the spin-1 model has the same dominant critical exponents as the spin-$\frac{1}{2}$ Ising model. Accurate estimates for both the critical temperature and non-physical singularities are obtained. In addition, evidence for a non-analytic confluent correction with exponent $\Delta_{1} = 1.1 \pm 0.1$ is found.


title: Spectrum of the Dirac Operator and Inversion Algorithms with Dynamical Staggered Fermions

abstract: Complete spectra of the staggered Dirac operator $\Dirac$ are determined in four-dimensional $SU(2)$ gauge fields with and without dynamical fermions. An attempt is made to relate the performance of multigrid and conjugate gradient algorithms for propagators with the distribution of the eigenvalues of~$\Dirac$.


title: Lattice Results for Heavy Light Matrix Elements

abstract: Lattice results for heavy light matrix elements are reviewed and some of their implications are very briefly discussed. Despite the fact that in most cases the lattice results for weak matrix elements at the moment have only a modest accuracy of about 20--30\% they already have important phenomenological repercussions; e.g.\ for $V_{td}/V_{ts}$, $x_s/x_d$ and $B\to K^\ast\gamma$. Presented at the XXVII International Conference on High Energy Physics, Glasgow, July 1994.


title: Chiral symmetry and lattice gauge theory

abstract: I review the problem of formulating chiral symmetry in lattice gauge theory. I discuss recent approaches involving an infinite tower of additional heavy states to absorb Fermion doublers. For hadronic physics this provides a natural scheme for taking quark masses to zero without requiring a precise tuning of parameters. A mirror Fermion variation provides a possible way of extending the picture to chirally coupled light Fermions. Talk presented at "Quark Confinement and the Hadron Spectrum," Como, Italy, 20-24 June 1994.


title: Tricritical Behavior of Two-Dimensional Scalar Field Theories

abstract: We compute by Monte Carlo numerical simulations the critical exponents of two-dimensional scalar field theories at the $\lambda\phi^6$ tricritical point. The results are in agreement with the Zamolodchikov conjecture based on conformal invariance.


title: Monte Carlo calculation for systems consisting of several coordinate patches

abstract: I investigate the time step dependence of Monte Carlo simulations for coordinate-spaces consisting of several patches. It is shown that a naive kinetic term does not necessarily converge to the same spectrum as a Hamiltonian calculation. Then an improved kinetic term is presented which allows one to connect the Monte Carlo and Rayleigh-Ritz results of intermediate volume SU(2) gauge theory.


title: The Gluon Propagator on a Large Volume, at $\beta=6.0$

abstract: We present the results of a high statistics lattice study of the gluon propagator, in the Landau gauge, at $\beta=6.0$. As suggested by previous studies, we find that, in momentum space, the propagator is well described by the expression $G(k^2)= \Big[ M^2 + Z\cdot k^2(k^2/\Lambda^2)^\eta\Big]^{-1} $. By comparing $G(k^2)$ on different volumes, we obtain a precise determination of the exponent $\eta=0.532(12)$, and verify that $M^2$ does not vanish in the infinite volume limit. The behaviour of $\eta$ and $M^2$ in the continuum limit is not known, and can only be studied by increasing the value of $\beta$.


title: An ``Improved" Lattice Study of Semi-leptonic Decays of D-Mesons

abstract: We present results of a lattice computation of the matrix elements of the vector and axial-vector currents which are relevant for the semi-leptonic decays $D \rightarrow K$ and $D \rightarrow K^*$. The computations are performed in the quenched approximation to lattice QCD on a $24^3 \times 48$ lattice at $\beta=6.2$, using an $O(a)$-improved fermionic action. In the limit of zero lepton masses the semi-leptonic decays $D \rightarrow K$ and $D \rightarrow K^*$ are described by four form factors: $f^{+}_K,V,A_1$ and $A_2$, which are functions of $q^2$, where $q^{\mu}$ is the four-momentum transferred in the process. Our results for these form factors at $q^2=0$ are: $f^+_K(0)=0.67 \er{7}{8}$ , $V(0)=1.01 \err{30}{13}$ , $A_1(0)=0.70 \err{7}{10}$ , $A_2(0)=0.66 \err{10}{15}$ , which are consistent with the most recent experimental world average values. We have also determined the $q^2$ dependence of the form factors, which we find to be reasonably well described by a simple pole-dominance model. Results for other form factors, including those relevant to the decays \dpi and \drho, are also given.


title: Geometrical Volume Effects in the Computation of the Slope of the Isgur-Wise Function

abstract: We use a method recently suggested for evaluating the slope of the Isgur-Wise function, at the zero-recoil point, on the lattice. The computations are performed in the quenched approximation to lattice QCD, on a $24^3 \times 48$ lattice at $\beta=6.2$, using an $O(a)$-improved action for the fermions. We have found unexpectedly large finite-volume effects in such a calculation. These volume corrections turned out to be purely geometrical and independent of the dynamics of the system. After the study of these effects on a smaller volume and for different quark masses, we give approximate expressions that account for them. Using these approximations we find $\xi^\prime(1)=-1.7 \pm 0.2$ and $\xi^\prime(1)=-1.4 +0.2-0.1$ for the slope of the Isgur-Wise function, for two mesons composed of a heavy quark slightly heavier and lighter, respectively, than the charm quark, and in both cases, a light antiquark whose mass is about that of the strange quark.


title: The beta function and equation of state for QCD with two flavors of quarks

abstract: We measure the pressure and energy density of two flavor QCD in a wide range of quark masses and temperatures. The pressure is obtained from an integral over the average plaquette or psi-bar-psi. We measure the QCD beta function, including the anomalous dimension of the quark mass, in new Monte Carlo simulations and from results in the literature. We use it to find the interaction measure, E-3p, yielding non-perturbative values for both the energy density E and the pressure p. uuencoded compressed PostScript file Revised version should work on more PostScript printers.


title: Embedding the N=2 Supersymmetric Yang-Mills Theory in the Adjoint Higgs-Yukawa Model on the Lattice

abstract: The N=2 supersymmetric Yang-Mills theory is formulated on the lattice. The feasibility of numerical simulations is discussed.


title: The continuum limit of $f_B$ from the lattice in the static approximation

abstract: We present an analysis of the continuum extrapolation of $f_B$ in the static approximation from lattice data. The method described here aims to uncover the systematic effects which enter in this extrapolation and has not been described before. Our conclusions are that we see statistical evidence for scaling of $f_B^{stat}$ for inverse lattice spacings $\gtap 2$ GeV but not for $\ltap 2$ GeV. We observe a lack of {\em asymptotic} scaling for a variety of quantities, including $f_B^{stat}$, at all energy scales considered. This can be associated with finite lattice spacing systematics. Once these effects are taken into account, we obtain a value of 230(35) MeV for $f_B^{stat}$ in the continuum where the error represents uncertainties due to both the statistics and the continuum extrapolation. In this method there is no error due to uncertainties in the renormalization constant connecting the lattice and continuum effective theories.


title: QCD Strings as Constrained Grassmannian Sigma Model:

abstract: We present calculations for the effective action of string world sheet in R3 and R4 utilizing its correspondence with the constrained Grassmannian sigma model. Minimal surfaces describe the dynamics of open strings while harmonic surfaces describe that of closed strings. The one-loop effective action for these are calculated with instanton and anti-instanton background, reprsenting N-string interactions at the tree level. The effective action is found to be the partition function of a classical modified Coulomb gas in the confining phase, with a dynamically generated mass gap.


title: Domain Wall Fermions and Chiral Gauge Theories

abstract: We review the status of the domain wall fermion approach to construct chiral gauge theories on the lattice. In this model an extra, fifth dimension is added and our 4-dimensional world lives on a domainwall induced by a soliton shaped mass defect that depends on the extra dimension only. We demonstrate that the domain wall model gives the correct anomaly structure when external gauge fields are used. We discuss two ways of adding dynamical gauge fields aiming at a lattice regularized chiral gauge theory. An approach is presented to keep the lattice infinite by regarding the fifth direction as the time-axis of a 4-dimensional Hamiltonian. Finally, a prospect to use domain wall fermions for simulating QCD is given.


title: Correlation Function at $\beta_t$ in the Disordered Phase of 2D Potts Models

abstract: We use Monte Carlo simulations to measure the spin-spin correlation function in the disordered phase of two-dimensional $q$-state Potts models with $q=10,15$, and $20$ at the first-order transition point $\beta_t$. To extract the correlation length $\xi_d$ from the exponential decay of the correlation function over several decades with the desired accuracy we make extensively use of cluster-update techniques and improved estimators. Our results for $\xi_d$ are compatible with an analytic formula. As a byproduct we also measure the energy moments in the disordered phase and find very good agreement with a recent large $q$ expansion at $\beta_t$.