title: $B$ and $B_c$ mesons with NRQCD and Clover actions

abstract: We present preliminary results from our study of the heavy-light spectrum and decay constants. For the heavy quark, we use NRQCD at various masses around and above the $b$ quark mass. For the first time, the heavy quark action and the heavy-light current consistently include corrections at second order in the non-relativistic expansion, as well as the leading finite $a$ corrections. The light quarks are simulated using a tadpole-improved Clover action at various masses in the strange and $c$ quark region.


title: Continuum Behaviour of Lattice QED, Discretized with One-Sided Lattice Differences, in One-Loop Order

abstract: A lattice action for QED is considered, where the derivatives in the Dirac operator are replaced by one-sided lattice differences. A systematic expansion in the lattice spacing of the one-loop contribution to the fermion self energy, vacuum polarization tensor and vertex function is carried out for an arbitrary choice of one-sided lattice differences. It is shown that only the vacuum polarization tensor possesses the correct continuum limit, while the fermion self energy and vertex function receive non-covariant contributions. A lattice action, discretized with a fixed choice of one-sided lattice differences therefore does not define a renormalizable field theory. The non-covariant contributions can however be eliminated by averaging the expressions over all possible choices of one-sided lattice differences.


title: Asymmetry Parameter Role in Description of Phase Structure of Lattice Gluodynamics at Finite Temperature

abstract: The role of lattice asymmetry parameter in description of the SU(2) gluodynamics phase structure at finite temperature is studied analytically. The fact that renormalization group relations which permit to remove the lattice asymmetry parameter from the thermodynamical quantities in the "naive" limit don't do the same in the approximation SU(N) = Z(N) keeping the "naive" limit results at the same time was point out. An additional condition which would fix asymmetry parameter is needed for this dependence removal.


title: Distribution of Instanton and Monopole Clustering

abstract: We study the relation between the instanton distribution and the monopole loop length in the SU(2) gauge theory with the abelian gauge fixing. We measure the monopole current from the multi-instanton ensemble on the $16^4$ lattice using the maximally abelian gauge. When the instanton density is dilute, there appear only small monopole loops. On the other hand, in the dense case, there appears one very long monopole loop, which is responsible for the confinement property, in each gauge configuration. We find a clear monopole clustering in the histogram of the monopole loop length from 240 gauge configurations.


title: Topological charge on the lattice: a field theoretical view of the geometrical approach

abstract: We construct sequences of ``field theoretical'' (analytical) lattice topological charge density operators which formally approach geometrical definitions in 2-d $CP^{N-1}$ models and 4-d $SU(N)$ Yang Mills theories. The analysis of these sequences of operators suggests a new way of looking at the geometrical method, showing that geometrical charges can be interpreted as limits of sequences of field theoretical (analytical) operators. In perturbation theory renormalization effects formally tend to vanish along such sequences. But, since the perturbative expansion is asymptotic, this does not necessarily lead to well behaved geometrical limits. It indeed leaves open the possibility that non-perturbative renormalizations survive.


title: Center Dominance and Z2 Vortices in SU(2) Lattice Gauge Theory

abstract: We find, in close analogy to abelian dominance in maximal abelian gauge, the phenomenon of center dominance in maximal center gauge for $SU(2)$ lattice gauge theory. Maximal center gauge is a gauge-fixing condition that preserves a residual $Z_2$ gauge symmetry; ``center projection'' is the projection of $SU(2)$ link variables onto $Z_2$ center elements, and ``center dominance'' is the fact that the center-projected link elements carry most of the information about the string tension of the full theory. We present numerical evidence that the thin $Z_2$ vortices of the projected configurations are associated with ``thick'' $Z_2$ vortices in the unprojected configurations. The evidence also suggests that the thick $Z_2$ vortices may play a significant role in the confinement process.


title: Simulations with Complex Measure

abstract: A method is proposed to handle the sign problem in the simulation of systems having indefinite or complex-valued measures. In general, this new approach, which is based on renormalisation blocking, is shown to yield statistical errors smaller than the crude Monte Carlo method using absolute values of the original measures. The improved method is applied to the 2D Ising model with temperature generalised to take on complex values. It is also adapted to implement Monte Carlo Renormalisation Group calculations of the magnetic and thermal critical exponents.


title: Effective Constraint Potential for Abelian Monopole in SU(2) Lattice Gauge Theory

abstract: We describe numerical calculation results for the probability distribution of the value of the monopole creation operator in the SU(2) lattice gluodynamics. We work in the maximal abelian projection. It occurs that at low temperatures, below the deconfinement phase transition, the maximum of the distribution is shifted from zero, which means that the effective constraint potential is of the Higgs type. Above the phase transition the minimum of the potential (the maximum of the monopole field distribution) is at the zero value of the monopole field. This fact confirms the existence of the abelian monopole condensate in the confinement phase of lattice gluodynamics, and agrees with the dual superconductor model of the confining vacuum.


title: A Wilson-Majorana Regularization for Lattice Chiral Gauge Theories

abstract: We discuss the regularization of chiral gauge theories on the lattice introducing only physical degrees of freedom. This is obtained by writing the Wilson term in a Majorana form, at the expense of the U(1) symmetry related to fermion number conservation. The idea of restoring chiral invariance in the continuum by introducing a properly chosen set of counterterms to be added to the tree level action is checked against one-loop perturbative calculations.


title: New results on topological susceptibility in SU(3) gauge theory

abstract: We survey recent lattice results on QCD topological properties. The behaviour of the topological susceptibility at the deconfining phase transition has been determined. This advance has been made possible by an i) an improvement of the topological charge operator and ii) a non-perturbative determination of renormalizations.


title: Heavy Light Weak Matrix Elements

abstract: I review the status of lattice calculations of heavy-light weak matrix elements, concentrating on semileptonic B decays to light mesons, B -> K* gamma, the B meson decay constant, f_B, and the mixing parameter B_B.


title: The Structure of Flux-tubes in SU(2)

abstract: The spatial distribution of the action and energy in the colour fields of flux-tubes is studied in lattice SU(2) field theory for static quarks at separations up to 1 fm. Special attention is paid to the structure of the colour fields associated with an excited flux-tube with E_u symmetry. We compare our results with hadronic string and flux-tube models. Sum rules are used to extract generalised beta-functions, to describe the expected colour field behaviour and to cross check the methods used.


title: Hubbard Model and Luscher fermions

abstract: We discuss numerical complexity of the L\"uscher algorithm applied to the Hubbard Model. In particular we present comparison to a certain algorithm based on direct computation of the fermionic determinant.


title: Curing O(a) Errors in 3-D Lattice SU(2) x U(1) Higgs Theory

abstract: We show how to make O(a) corrections in the bare parameters of 3-D lattice SU(2) times U(1) Higgs theory which remove O(a) errors in the match between the infrared behavior and the infrared behavior of the continuum theory. The corrections substantially improve the convergence of lattice data to a small a limit.


title: Lattice Monte Carlo Data versus Perturbation Theory

abstract: Differences between lattice Monte Carlo data and perturbation theory (for example the lack of asymptotic scaling) are usually associated with the `bad' behaviour of the bare lattice coupling g_0 due to the effects of large (and unknown) higher order terms in g_0. In this philosophy a new, renormalised coupling g' is defined with the aim of making the higher order coefficients of the perturbative series in g' as small as possible. In this paper an alternative scenario is discussed where lattice artifacts are proposed as the cause of the disagreement between Monte Carlo data and the g_0-perturbative series. We find that with the addition of a lattice artifact term, the usual asymptotic scaling expression in g_0 is in excellent agreement with Monte Carlo data. Lattice data studied includes the string tension, the hadronic scale r_0, the discrete beta function, M_rho, f_pi and the 1P-1S splitting in charmonium.


title: Comparing Different Improvement Programs for the N-Vector Model

abstract: We discuss the connection between various types of improved actions in the context of the two-dimensional sigma-model. We also discuss spectrum-improved actions showing that these actions do not have any improved behaviour. An O(a^2) on-shell improved action with all couplings defined on a plaquette and satisfying reflection positivity is also explicitly constructed.


abstract: The differences between lattice Monte Carlo data and perturbation theory are usually associated with the ``bad'' behaviour of the bare lattice coupling g_0 due to the effects of large (and unknown) higher order coefficients in the g_0 perturbative series. In this philosophy a new, renormalised coupling g' is defined with the aim of reducing the higher order coefficients of the perturbative series in g'. An improvement in the agreement between lattice data and this new perturbation series is generally observed. In this paper an alternative scenario is discussed where lattice artifacts are proposed as the cause of the disagreement between lattice data and the g_0-perturbative series. We find that this interpretation provides excellent agreement between lattice data and perturbation theory in g_0 corrected for lattice artifacts. We show that this viewpoint leads typically to an order of magnitude improvement in the agreement between lattice data and perturbation theory, compared to typical g' perturbation expansions. The success of this procedure leads to a determination of Lambda_MSbar^{N_f=0} of 220 +- 20 MeV. Lattice data studied includes quenched values of the string tension, the hadronic scale r_0, the discrete beta function Delta\beta, M_rho, f_pi and the 1P-1S splitting in charmonium. The new 3-loop term of the lattice beta- function has been incorporated in this study. A discussion of the implication of this result for lattice calculations is presented.


title: A new algorithm for numerical simulation of Langevin equations

abstract: Formulated is a new systematic method for obtaining higher order corrections in numerical simulation of stochastic differential equations (SDEs), i.e., Langevin equations. Random walk step algorithms within a given order of finite $\Delta t$, are obtained so as to reproduce within that order a corresponding transition density of the Fokker-Planck equations, in the weak Taylor approximation scheme. A great advantage of our method is its straightforwardness such that direct perturbative calculations produce the algorithm as an end result, so that the procedure is tractable by computer. Examples in general form for curved space cases as well as flat space cases are given in some order of approximations. Simulations are performed for specific examples of U(1) system and SU(2) systems, respectively.


title: The PZ method for estimating determinant ratios, with applications

abstract: We introduce a new method for estimating determinants or determinant ratios of large matrices, which combines the techniques of Pad\`{e} approximation with rational functions and $Z_{2}$ noise estimation of traces of large matrices. The method requires simultaneously solving several matrix equations, which can be conveniently accomplished using the MR method. We include some preliinary results, and indicate potential applications to non-Hermitian matrices, and Hybrid Monte Carlo without pseudofermions.


title: The Pseudo Specific Heat in SU(2) Gauge Theory : Finite Size Dependence and Finite Temperature Effects

abstract: We investigate the pseudo specific heat of SU(2) gauge theory near the crossover point on $4^4$ to $16^4$ lattices. Several different methods are used to determine the specific heat. The curious finite size dependence of the peak maximum is explained from the interplay of the crossover phenomenon with the deconfinement transition occurring due to the finite extension of the lattice. We find, that for lattices of size $8^4$ and larger the crossover peak is independent of lattice size at $\beta_{co}=2.23(2)$ and has a peak height of $C_{V,co}=1.685(10)$. We conclude therefore that the crossover peak is not the result of an ordinary phase transition. Further, the contributions to $C_V$ from different plaquette correlations are calculated. We find, that at the peak and far outside the peak the ratio of contributions from orthogonal and parallel plaquette correlations is different. To estimate the finite temperature influence on symmetric lattices far off the deconfinement transition point we calculate the modulus of the lattice average of the Polyakov loop on these lattices and compare it to predictions from a random walk model.


title: Hybrid Monte Carlo Without Pseudofermions

abstract: We introduce a dynamical fermion algorithm which is based on the hybrid Monte Carlo (HMC) algorithm, but without pseudofermions. The molecular dynamics steps in HMC are retained except the derivatives with respect to the gauge fields are calculated with the $Z_2$ noise. The determinant ratios are estimated with the Pa\`{d}e - $Z_2$ method. Finally, we use the Kennedy-Kuti linear accept/reject method for the Monte Carlo step which is shown to respect detailed balance. We comment on the comparison of this algorithm with the pseudofermion algorithm.


title: Multi-Grid Monte Carlo via $XY$ Embedding. II. Two-Dimensional $SU(3)$ Principal Chiral Model

abstract: We carry out a high-precision simulation of the two-dimensional $SU(3)$ principal chiral model at correlation lengths $\xi$ up to $\sim 4 \times 10^5$, using a multi-grid Monte Carlo (MGMC) algorithm and approximately one year of Cray C-90 CPU time. We extrapolate the finite-volume Monte Carlo data to infinite volume using finite-size-scaling theory, and we discuss carefully the systematic and statistical errors in this extrapolation. We then compare the extrapolated data to the renormalization-group predictions. The deviation from asymptotic scaling, which is $\approx 12%$ at $\xi \sim 25$, decreases to $\approx 2%$ at $\xi \sim 4 \times 10^5$. We also analyze the dynamic critical behavior of the MGMC algorithm using lattices up to $256 \times 256$, finding the dynamic critical exponent $z_{int,{\cal M}^2} \approx 0.45 \pm 0.02$ (subjective 68% confidence interval). Thus, for this asymptotically free model, critical slowing-down is greatly reduced compared to local algorithms, but not completely eliminated.


title: Stability of the Bulk Phase Diagram of the SU(2) Lattice Gauge Theory with Fundamental - Adjoint Action

abstract: Using improved mean field and strong coupling expansions we re-analyse the bulk phase diagram of the fundamental-adjoint action of the SU(2) Lattice Gauge Theory. We find that the qualitative features of the bulk phase diagram are robust and unchanged by the inclusion of higher order terms. On the other hand, some of the quantitative features, such as the location of the endpoint of the line of bulk phase transitions, seem to be strongly dependent on the higher terms of the strong coupling expansion.


title: Scaling in SU(3) Pure Gauge Theory with a Renormalization Group Improved Action

abstract: We study the scaling properties of the static quark potential and the ratio of the critical temperature $T_c$ to the square root of the string tension $\sigma$ in the SU(3) pure gauge theory using a renormalization group improved action. We first determine the critical coupling $\beta_c$ on lattices with temporal extension $N_t=3$, 4, and 6, and then calculate the static quark potential at the critical couplings on lattices at zero temperature. We note that the static quark potentials obtained are rotationally invariant with errors of at most 1 - 2 % in all the three cases, and that the potential $V(R)$ in physical units scales in the whole region of $R$ investigated. The values of $T_{c}/\sqrt{\sigma}$ for the three cases in the infinite volume limit are identical within errors. We estimate the value in the continuum limit to be $T_{c}/\sqrt{\sigma} = 0.656(4)$, which is slightly larger than the value in the continuum limit from the one-plaquette action, 0.629(3).


title: $B - \bar B$ Mixing in the HQET

abstract: We present a high statistics, quenched lattice calculation of the B-parameters $B_{B_d}$ and $B_{B_s}$, computed at lowest order in the HQET. The results were obtained using a sample of 600 quenched gauge field configurations, generated by Monte Carlo simulation at $\beta=6.0$ on a $24^{3}\times 40$ lattice. For the light quarks the SW-Clover action was used; the propagator of the lattice HQET was also tree-level improved. Our best estimate of the renormalization scale independent B-parameter is $\hat{B}_{B_d} = 1.03 \pm 0.06 \pm 0.18$. $\hat{B}_{B_d}$ has been obtained by using ``boosted'' perturbation theory to calculate the renormalization constants which relate the matrix elements of the lattice operators to the corresponding amplitudes in the continuum. Due to the large statistics, the errors in the extraction of the matrix elements of the relevant bare operators are rather small. The main systematic error, corresponding to $\pm 0.18$ in the above result, comes from the uncertainty in the evaluation of the renormalization constants, for which the one-loop corrections are rather large. The non-perturbative evaluation of these constants will help to reduce the final error. We also obtain $\hat{B}_{B_s}/\hat{B}_{B_d} = 1.01 \pm 0.01$ and $f^2_{B_s}\hat{B}_{B_s}/f^2_{B_d}\hat{B}_{B_d} = 1.38 \pm 0.07$.


title: Relating U(1) monopole configurations to SU(2) saddle-point configurations

abstract: We have studied field configurations of the 3-dimensional Georgi-Glashow model which interpolate between the $U(1)$ and the $SU(2)$ limits. In the intermediate region, these configurations contain 't-Hooft--Polyakov monopoles. We use cooling and extremization to find these configurations and investigate their evolution as we adiabatically move towards the $U(1)$ and the $SU(2)$ limits. We also evolve an $SU(2)$ saddle point solution towards the $U(1)$ limit to see the relation between the unstable solutions in the $SU(2)$ theory and the stable ones in the $U(1)$ theory.


title: A calculation of the $B_{B}$ parameter in the static limit

abstract: We calculate the $B_{B}$ parameter, relevant for $\overline{B}^0$ -- $B^0$ mixing, from a lattice gauge theory simulation at $\beta = 6.0$. The bottom quarks are simulated in the static theory, the light quarks with Wilson fermions. Improved smearing functions produced by a variational technique, MOST, are used to reduce statistical errors and minimize excited-state contamination of the ground-state signal. We obtain $B_B(4.33 GeV) = 0.98^{+4}_{-4}$ (statistical) $^{+3}_{-18}$ (systematic) which corresponds to $\widehat{B}_B = 1.40^{+6}_{-6}$ (statistical) $^{+4}_{-26}$ (systematic) for the one-loop renormalization-scheme-independent parameter. The systematic errors include the uncertainty due to alternative (less favored) treatments of the perturbatively-calculated mixing coefficients; this uncertainty is at least as large as residual differences between Wilson-static and clover-static results. Our result agrees with extrapolations of results from relativistic (Wilson) heavy quark simulations.


title: Light vector meson decay constants and the renormalization factor from a tadpole-improved action

abstract: The rho, K* and phi decay constants and the vector current renormalization factor are studied by using an O(a^2) classically-improved, tadpole-improved action. Tree-level calculations are used to show how the classical improvement of the action, involving next-nearest-neighbour timesteps, is transferred to the matrix elements. Simulations are performed on coarse lattices and compared to Wilson results from both coarse and fine lattices. The improved action data are found to resemble Wilson data obtained at 1/3 of the lattice spacing, which is the same degree of improvement that is seen by comparing the mass spectra.


title: Unstable Modes and Confinement in the Lattice Schr\"odinger Functional Approach

abstract: We analyze the problem of the Nielsen-Olesen unstable modes in the SU(2) lattice gauge theory by means of a recently introduced gauge-invariant effective action. We perform numerical simulations in the case of a constant Abelian chromomagnetic field. We find that for lattice sizes above a certain critical length the density of effective action shows a behaviour compatible with the presence of the unstable modes. We put out a possible relation between the dynamics of the unstable modes and the confinement.


title: The perfect action for non-degenerate staggered fermions

abstract: The perfect action of free staggered fermions is calculated by blocking from the continuum for degenerate and non-degenerate flavor masses. The symmetry structure, connecting flavor transformations and translations, is explained directly from the blocking scheme. It is convenient to use a modified Fourier transformation, respecting this connection, to treat the spin-flavor structure of the blockspins. The perfect action remains local in the non-degenerate case; it is explicitly calculated in two dimensions. I finally comment on the relation of the blocking scheme to the transition from Dirac-K\"ahler fermions to staggered fermions.


title: The Anti-Ferromagnetic Vacuum

abstract: Certain effective vertices may generate a non-homogeneous, periodic vacuum structure. The excitations above such a vacuum are studied in the framework of the $\phi^4$ and gauge models. The formation of the non-homogeneous vacuum is accompanied by the dynamical breakdown of the space-time inversion symmetry. Chiral transformation is introduced for bosons in close analogy with lattice fermions.


title: U(1) Lattice Gauge theory and its Dual

abstract: The three dimensional U(1) Lattice Gauge, in the weak coupling limit, is dual to a Discrete Gaussian model. We investigate this dual model and use it to calculate properties of the U(1) theory. We find that, because of the nature of the dual model, its advantages are outweighed by large autocorrelation times generated when the dual system becomes disordered.


title: Monte Carlo study of very weak first-order transitions in the three-dimensional Ashkin-Teller model

abstract: We propose numerical simulations of the Ashkin-Teller model as a foil for theoretical techniques for studying very weakly first-order phase transitions in three dimensions. The Ashkin-Teller model is a simple two-spin model whose parameters can be adjusted so that it has an arbitrarily weakly first-order phase transition. In this limit, there are quantities characterizing the first-order transition which are universal: we measure the relative discontinuity of the specific heat, the correlation length, and the susceptibility across the transition by Monte Carlo simulation.


title: Chiral zero modes on the domain-wall model in 4+1 dimensions

abstract: We investigate an original domain-wall model in 4+1 dimensions numerically in the presence of U(1) dynamical gauge field only in an extra dimension, corresponding to a weak coupling limit of 4-dimensional physical gauge coupling. Using a quenched approximation we carry out numerical simulation for this model at $\beta_{s} (= 1 / g^{2}_{s}) =$ 0.29 (``symmetric'' phase) and 0.5 (``broken'' phase), where $g_s$ is the gauge coupling constant of the extra dimension. In the broken phase, we found that there exists a critical value of a domain-wall mass $m_{0}^{c}$ which separates a region with a fermionic zero mode on the domain wall from the one without it in the same case of (2+1)-dimensional model. On the other hand, in the symmetric phase, our numerical data suggest that the chiral zero modes disappear in the infinite limit of 4-dimensional volume. From these results it seems difficult to construct the U(1) lattice chiral gauge theory via an original domain-wall formulation.


title: Domain walls, $Z(N)$ charge and $A_0$ condensate: a canonical ensemble study

abstract: The deconfinement phase transition is studied in the ensemble canonical with respect to triality. Since this ensemble implies a projection to the zero triality sector of the theory we introduce a quantity which is insensitive to $Z(N_c)$ symmetry but can reveal a critical behaviour in the theory with dynamical quarks. Further, we argue that in the canonical ensemble description of full QCD there exist domains of different $Z(N_c)$ phases which are degenerate and possess normal physical properties. This contradicts the predictions of the grand canonical ensemble. We propose a new order parameter to test the realization of the discrete $Z(N_c)$ symmetry at finite temperature and calculate it for the case of $Z(2)$ gauge fields coupled to fundamental fermions.


title: Another Look at Neural Multigrid

abstract: We present a new multigrid method called neural multigrid which is based on joining multigrid ideas with concepts from neural nets. The main idea is to use the Greenbaum criterion as a cost functional for the neural net. The algorithm is able to learn efficient interpolation operators in the case of the ordered Laplace equation with only a very small critical slowing down and with a surprisingly small amount of work comparable to that of a Conjugate Gradient solver. In the case of the two-dimensional Laplace equation with SU(2) gauge fields at beta=0 the learning exhibits critical slowing down with an exponent of about z = 0.4. The algorithm is able to find quite good interpolation operators in this case as well. Thereby it is proven that a practical true multigrid algorithm exists even for a gauge theory. An improved algorithm using dynamical blocks that will hopefully overcome the critical slowing down completely is sketched.


title: One-loop anisotropy for improved actions

abstract: We determine the one-loop correction to the anisotropy factor for the square Symanzik improved lattice action, extracted fromthe finite volume effective action for SU(N) gauge theories in the background of a zero-momentum gauge field. The result is smaller by approximately a factor 3 than the one-loop correction for the anisotropic Wilson action. We also comment on the Hamiltonian limit.


title: The Renormalization Group and Dynamical Triangulations

abstract: A block spin renormalization group approach is introduced which can be applied to dynamical triangulations in any dimension.


title: Continuum limit of finite temperature $\lambda\phi^4_3$ from lattice Monte Carlo

abstract: The $\phi^4_3$ model at finite temperature is simulated on the lattice. For fixed $N_t$ we compute the transition line for $N_s \to \infty$ by means of Finite Size Scaling techniques. The crossings of a Renormalization Group trajectory with the transition lines of increasing $N_t$ give a well defined limit for the critical temperature in the continuum. By considering different RG trajectories, we compute $T^c/g$ as a function of the renormalized parameters.


title: Finite Temperature Properties of the SO(3) Lattice Gauge Theory

abstract: We make a numerical study of the finite temperature properties of the SO(3) lattice gauge theory. As its symmetry properties are quite different from those of the SU(2) LGT, a different set of observables have to be considered in this model. We study several observables, such as, the plaquette square, the Z(2) monopole density, the fundamental and adjoint Wilson line, and the tiled Wilson line correlation function. Our simulations show that the Z(2) monopoles condense at strong coupling just as in the bulk system. This transition is seen at approximately the same location as in the bulk system. A surprising observation is the multiple valuedness of the adjoint Wilson line at high temperatures. At high temperatures, we observe long lived metastable states in which the adjoint Wilson line takes positive and negative values. The numerical values of other observables in these two states appear to be almost the same. We study these states using different methods and also make comparisons with the high temperature behaviour of the SU(2) LGT. Finally, we discuss various interpretations of our results and point out their relevance for the phase diagram of the SO(3) LGT at finite temperature.


title: The Abelianicity of Cooled SU(2) Lattice Configurations

abstract: We introduce a gauge-invariant measure of the local "abelianicity" of any given lattice configuration in non-abelian lattice gauge theory; it is essentially a comparison of the magnitude of field strength commutators to the magnitude of the field strength itself. This measure, in conjunction with the cooling technique, is used to probe the SU(2) lattice vacuum for a possible large-scale abelian background, underlying the local short-range field fluctuations. We do, in fact, find a substantial rise in abelianicity over 10 cooling steps or so, after which the abelianicity tends to drop again.


title: Is trivial the antiferromagnetic RP(2) model in four dimensions?

abstract: We study the antiferromagnetic RP(2) model in four dimensions. We find a second order transition with two order parameters, one ferromagnetic and the other antiferromagnetic. The antiferromagnetic sector has mean-field critical exponents and a renormalized coupling which goes to zero in the continuum limit. The exponents of the ferromagnetic channel are not the mean-field ones, but the difference can be interpreted as logarithmic corrections. We perform a detailed analysis of these corrections and conclude the triviality of the continuum limit of this model.


title: The critical points of lattice QCD with a non--zero quark density

abstract: We study the interplay of quark number density and chiral symmetry in lattice QCD. We suggest that both are controlled by the eigenvalue spectrum of the fermionic propagator matrix, which shapes the pattern of zeros of the partition function. The onset of the quark current would be triggered by the lowest lying eigenvalue, the chiral transition by the density of zeros, the two critical points being distinct in full QCD, and coincident in the quenched approximation. Our preliminary estimate for the critical point in full QCD in the infinite couling limit compares favourably with the predictions of the strong coupling expansions and of numerical simulations based on exact, alternative representations of the partition function. Several reasons of perplexity however remain, which are briefly discussed.


title: The exact equivalence of the one-flavour lattice Thirring model with Wilson fermions to a two-colour loop model

abstract: Within Euclidean lattice field theory an exact equivalence between the one-flavour 2D Thirring model with Wilson fermions and Wilson parameter $r = 1$ to a two-colour loop model on the square lattice is established. For non-interacting fermions this model reduces to an exactly solved loop model which is known to be a free fermion model. The two-colour loop model equivalent to the Thirring model can also be understood as a 4-state 49-vertex model.


title: A determination of the mass gap in the O(n) sigma model

abstract: We calculate the finite volume mass gap $M(L)$ at 3-loop level in the non-linear O($n$) $\sigma$-model in two dimensions in small volumes. By applying the Monte Carlo measurements of the running coupling $\bar g^2(L)=2nM(L)L/(n-1)$ by L\"uscher, Weisz and Wolff measured in units of the physical mass gap $m$, the result is used to determine $m$ in units of the $\Lambda$-parameter in the O(3) and O(4) models. Our determinations show good agreement with those by Hasenfratz, Maggiore and Niedermayer in both models. We note that this manuscript has been revised in our paper hep-lat/9810025 by using the corrected four-loop $\beta$-function on the lattice.


title: Weak Matrix Elements on the Lattice: Recent Developments in K-Physics

abstract: Some recent developments in the calculation of weak matrix elements on the lattice are reviewed. In particular, we concentrate on the application of MPSTV non-perturbative renormalization method to the operators relevant to kaon physics. The chiral behaviour of the $B_K$ parameter and the long-standing question of the $\Delta I=1/2$ rule are discussed.


title: Yang-Lee zeros of a random matrix model for QCD at finite density

abstract: We study the Yang-Lee zeros of a random matrix partition function with the global symmetries of the QCD partition function. We consider both zeros in the complex chemical potential plane and in the complex mass plane. In both cases we find that the zeros are located on a curve. In the thermodynamic limit, the zeros appear to merge to form a cut. The shape of this limiting curve can be obtained from a saddle-point analysis of the partition function. An explicit solution for the line of zeros in the complex chemical potential plane at zero mass is given in the form of a transcendental equation.


title: Comparison of Multi-quark Matrix Inversion Algorithms

abstract: We test iterative algorithms, MR, QMR$\gamma_5$ and BiCG$\gamma_5$, to compare their efficiency in matrix inversion with multi-quarks (shifted matrices) within one iteration process. Our results on the 8^3 x 12 and 16^3 x 24 show that MR admits multi-quark calculation with less memory requirement, whereas QMR is faster for the single quark calculation.


title: Improving Lattice Quark Actions

abstract: We explore the first stage of the Symanzik improvement program for lattice Dirac fermions, namely the construction of doubler-free, highly improved classical actions on isotropic as well as anisotropic lattices (where the temporal lattice spacing, a_t, is smaller than the spatial one). Using field transformations to eliminate doublers, we derive the previously presented isotropic D234 action with O(a^3) errors, as well as anisotropic D234 actions with O(a^4) or O(a_t^3, a^4) errors. Besides allowing the simulation of heavy quarks within a relativistic framework, anisotropic lattices alleviate potential problems due to unphysical branches of the quark dispersion relation (which are generic to improved actions), facilitate studies of lattice thermodynamics, and allow accurate mass determinations for particles with bad signal/noise properties, like glueballs and P-state mesons. We also show how field transformations can be used to completely eliminate unphysical branches of the dispersion relation. Finally, we briefly discuss future steps in the improvement program.


title: Hybrid mesons from quenched QCD

abstract: We use lattice methods to evaluate from first principles the spectrum of hybrid mesons produced by gluonic excitations in quenched QCD with quark masses near the strange quark mass. For the spin-exotic mesons with $J^{PC}=1^{-+},\ 0^{+-}$, and $2^{+-}$ which are not present in the quark model, we determine the lightest state to be $1^{-+}$ with mass of 2.0(2) GeV.


title: The U(1) Problem in Chiral Random Matrix Models

abstract: We show that conventional asymmetric chiral random matrix models (ChRMM), with a gaussian distribution in the asymmetry, provide for a screening of the topological charge and a resolution of the $U(1)$ problem in the unquenched approximation. Our exact results to order $1/N$ are in agreement with numerical estimates using large ensembles of asymmetric ChRMM with gaussian distributions.


title: Wilson--Yukawa Chiral Model on Lattice and Non-commutative Geometry

abstract: Based upon the mathematical formulas of Lattice gauge theory and non-commutative geometry differential calculus, we developed an approach of generalized gauge theory on a product of the spacetime lattice and the two discrete points(or a $Z_2$ discrete group). We introduce a differentiation for non-nearest-neighbour points and find that this differentiation may lead to the introduction of Wilson term in the free fermion Lagrangian on lattice. The Wilson-Yukawa chiral model on lattice is constructed by the generalized gauge theory and a toy model and Smit-Swift model are studied.


title: A Dedicated Computer for Ising-like Spin Glass Models

abstract: We present a parallel machine, based on programmable devices, dedicated to simulate spin glass models with Z2 variables and short range interaction. A working prototype is described for two lattices containing 312x312 spins each with an update time of 50 ns per spin. The final version of the three dimensional parallel machine is discussed with spin update time up to 312 ps.


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

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


title: Development in Lattice QCD

abstract: After a brief discussion of the promise and limitations of the lattice technique, I review lattice QCD results for several quantities of phenomenological interest. These are: matrix elements for heavy-to-light meson decays, leptonic decay constants f_B and f_D, the parameters B_B and B_K for neutral B and K meson mixing respectively, the strong coupling constant, light quark masses and the lightest scalar glueball mass.


title: Cellular automata and self-organized criticality

abstract: Cellular automata provide a fascinating class of dynamical systems capable of diverse complex behavior. These include simplified models for many phenomena seen in nature. Among other things, they provide insight into self-organized criticality, wherein dissipative systems naturally drive themselves to a critical state with important phenomena occurring over a wide range of length and time scales.


title: Self-organized criticality

abstract: I review the concept of self-organized criticality, wherein dissipative systems naturally drive themselves to a critical state with important phenomena occurring over a wide range of length and time scales. Several exact results are demonstrated for the Abelian sandpile.


title: Better Actions

abstract: We explain why compact U(1) confines and how to fix it. We show that plaquettes of negative trace carry most of the confinement signal in compact SU(2). We show how to perform noncompact gauge-invariant simulations without auxiliary fields. We suggest a way to simulate fermions without doublers.


title: Topology in 4D simplicial quantum gravity

abstract: We simulate 4d simplicial gravity for three topologis S4, S3xS1, (S1)^4 and show that the free energy for these three fixed topology ensembles is the same in the thermodynamic limit. We show, that the next-to-leading order corrections, at least away from the critical point, can be described by kinematic sources.


title: Light Quenched Hadron Spectrum and Decay Constants on different Lattices

abstract: In this paper we study O(2000) (quenched) lattice configurations from the APE collaboration, for different lattice volumes and for 6.0 \le beta \le 6.4 using both the Wilson and the SW-Clover fermion actions. We determine the light hadronic spectrum and meson decay constants and study the mesonic dispersion relation. We extract the hadronic variable J and the strange quark mass in the continuum at the next-to-leading order obtaining m_s^{MSbar}(mu=2 GeV) = 122 +/- 20 MeV. A study is made of their dependence on lattice spacing. We implement a newly developed technique to extract the inverse lattice spacing using data at the simulated values of the quark mass (i.e. at masses around the strange quark mass).


title: Reducing Discretization Errors in Lattice QCD Spectroscopy

abstract: The improved Wilson quark action - the clover action - is constructed to have smaller discretization errors than the normal Wilson quark action. We test this in a quenched spectroscopy computation on 6 lattice ensembles with spacings from 0.15 to 0.43 fm. To ensure that the dominant scaling violations come from the fermions we use an $O(a^2)$ improved 6-link $SU(3)$ pure gauge action. We find evidence that fermionic scaling violations are consistent with $O(a^2)$ for clover fermions and $O(a)$ with a nonnegligible $O(a^2)$ term for standard Wilson fermions. This latter mixed ansatz makes a reliable continuum extrapolation problematic for Wilson fermions. For clover fermions, on the other hand, we obtain accurate predictions for hadron masses in quenched continuum QCD. We find that the slopes of the scaling violations are roughly 200 MeV for both Wilson and clover fermions.


title: Matrix Elements of 4-fermion Operators with Quenched Wilson Fermions

abstract: We presents results for the matrix elements of a variety of four-fermion operators calculated using quenched Wilson fermions. Our simulations are done on 170 lattices of size $32^3 x 64$ at $\beta = 6.0$. We find $B_K=0.74 +- 0.04 +- 0.05$, $B_D= 0.78 +- 0.01$, $B_7^{3/2}= 0.58 +- 0.02 {+0.07 \atop -0.03}$, $B_8^{3/2}= 0.81 +- 0.03 {+.03 \atop -0.02}$, with all results being in the NDR scheme at $\mu=2 GeV$. We also calculate the B-parameter for the operator $\CQ_s$, which is needed in the study of the difference of B-meson lifetimes. Our best estimate is $B_s(NDR,\mu=1/a)=0.80 +- 0.01$. This is given at the lattice scale since the required 2-loop anomalous dimension matrix is not known. In all these estimates, the first error is statistical, while the second is due to the use of truncated perturbation theory to match continuum and lattice operators. Errors due to quenching and lattice discretization are not included. We also present new results for the perturbative matching coefficients, extending the calculation to all Lorentz scalar four-fermion operators, and using NDR as the continuum scheme.


title: Topology at the Deconfinement Transition Uncovered by Inverse Blocking in SU(2) Pure Gauge Theory with Fixed Point Action

abstract: Renormalization group transformations as discussed recently in deriving fixed point actions are used to analyse the vacuum structure near to the deconfinement temperature. Monte Carlo configurations are generated using the fixed point action. We compare equilibrium configurations with configurations obtained by inverse blocking from a coarser lattice. The absence of short range vacuum fluctuations in the latter does not influence the string tension. For the inversely blocked configurations we find the following: (i) the topological susceptibility chi_top is consistent with the phenomenological value in the confinement phase, (ii) chi_top drops across the deconfinement transition, (iii) density and size of instantons are estimated, (iv) the topological density is found to be correlated to Abelian monopole currents and (v) the density of spacelike monopole currents becomes a confinement order parameter.


title: Ab Initio Calculation of Relativistic Corrections to the Static Interquark potential I: SU(2) Gauge Theory

abstract: We test the capability of state-of-the-art lattice techniques for a precise determination of relativistic corrections to the static interquark potential, by use of SU(2) gauge theory. Emphasis is put on the short range structure of the spin dependent potentials, with lattice resolution a ranging from a approx 0.04 fm (at beta=2.74) down to a approx 0.02 fm (at beta=2.96) on volumes of 32^4 and 48^4 lattice sites. We find a new short range Coulomb-like contribution to the spin-orbit potential V_1'.


title: Quarkonium spin structure in lattice NRQCD

abstract: Numerical simulations of the quarkonium spin splittings are done in the framework of lattice nonrelativistic quantum chromodynamics (NRQCD). At leading order in the velocity expansion the spin splittings are of $O(M_Q v^4)$, where $M_Q$ is the renormalized quark mass and $v^2$ is the mean squared quark velocity. A systematic analysis is done of all next-to-leading order corrections. This includes the addition of $O(M_Q v^6)$ relativistic interactions, and the removal of $O(a^2 M_Q v^4)$ discretization errors in the leading-order interactions. Simulations are done for both S- and P-wave mesons, with a variety of heavy quark actions and over a wide range of lattice spacings. Two prescriptions for the tadpole improvement of the action are also studied in detail: one using the measured value of the average plaquette, the other using the mean link measured in Landau gauge. Next-to-leading order interactions result in a very large reduction in the charmonium splittings, down by about 60% from their values at leading order. There are further indications that the velocity expansion may be poorly convergent for charmonium. Prelimary results show a small correction to the hyperfine splitting in the Upsilon system.


title: Thermodynamics of the two--dimensional O(3) non linear sigma model with Fixed-Point lattice action

abstract: We calculate the free energy density of the two dimensional O(3) non linear sigma model over a large temperature region. At high temperatures the calculations could be done by perturbation theory whereas in the low temperature regime we used Monte-Carlo simulations. Using the Fixed--Point action the results do not show any remarkable cut-off effects and furthermore the behaviour of these effects is nearly the same over the whole temperature region. The results of the Fixed-Point action are then compared to those of the standard Wilson action.


title: Scaling Behavior in 4D Simplicial Quantum Gravity

abstract: Scaling relations in four-dimensional simplicial quantum gravity are proposed using the concept of the geodesic distance. Based on the analogy of a loop length distribution in the two-dimensional case, the scaling relations of the boundary volume distribution in four dimensions are discussed in three regions: the strong-coupling phase, the critical point and the weak-coupling phase. In each phase a different scaling behavior is found.


title: Abelian projection and studies of gauge-variant quantities in lattice QCD without gauge fixing

abstract: We suggest a new (dynamical) Abelian projection of the lattice QCD. It contains no gauge condition imposed on gauge fields so that Gribov copying is avoided. Configurations of gauge fields that turn into monopoles in the Abelian projection can be classified in a gauge invariant way. In the continuum limit, the theory respects the Lorentz invariance. A similar dynamical reduction of the gauge symmetry is proposed for studies of gauge-variant correlators (like a gluon propagator) in lattice QCD. Though the procedure is harder for numerical simulations, it is free of gauge fixing artifacts, like the Gribov horizon and copies.


title: QCD with domain wall quarks

abstract: We present lattice calculations in QCD using a variant of Kaplan fermions which retain the continuum SU(N)xSU(N) chiral symmetry on the lattice in the limit of an infinite extra dimension. In particular, we show that the pion mass and the four quark matrix element related to K_0-K_0-bar mixing have the expected behavior in the chiral limit, even on lattices with modest extent in the extra dimension, e.g. N_s=10.


title: Which Chiral Symmetry is Restored in High Temperature QCD?

abstract: Sigma models for the high temperature phase transition in quantum chromodynamics (QCD) suggest that at high temperature the SU(N_f) x SU(N_f) chiral symmetry becomes exact, but the anomalous axial U(1) symmetry need not be restored. In numerical lattice simulations, traditional methods for detecting symmetry restoration have sought multiplets in the screening mass spectrum. However, these methods were imprecise and the results, so far, incomplete. With improved statistics and methodology, we are now able to offer evidence for a restoration of the SU(2) x SU(2) chiral symmetry just above the crossover, but not of the axial U(1) chiral symmetry.


title: Quantum geometry of topological gravity

abstract: We study a c=-2 conformal field theory coupled to two-dimensional quantum gravity by means of dynamical triangulations. We define the geodesic distance r on the triangulated surface with N triangles, and show that dim[r^{d_H}]= dim[N], where the fractal dimension d_H = 3.58 +/- 0.04. This result lends support to the conjecture d_H = -2\alpha_1/\alpha_{-1}, where \alpha_{-n} is the gravitational dressing exponent of a spin-less primary field of conformal weight (n+1,n+1), and it disfavors the alternative prediction d_H = -2/\gamma_{str}. On the other hand, we find dim[l] = dim[r^2] with good accuracy, where l is the length of one of the boundaries of a circle with (geodesic) radius r, i.e. the length l has an anomalous dimension relative to the area of the surface. It is further shown that the spectral dimension d_s = 1.980 +/- 0.014 for the ensemble of (triangulated) manifolds used. The results are derived using finite size scaling and a very efficient recursive sampling technique known previously to work well for c=-2.


title: Critical properties and monopoles in U(1) lattice gauge theory

abstract: We present a detailed study of the properties of the phase transition in the four-dimensional compact U(1) lattice gauge theory supplemented by a monopole term, for values of the monopole coupling $\lambda$ such that the transition is of second order. By a finite size analysis we show that at $\lambda= 0.9$ the critical exponent is already characteristic of a second-order transition. Moreover, we find that this exponent is definitely different from the one of the Gaussian case. We further observe that the monopole density becomes approximately constant in the second-order region. Finally we reveal the unexpected phenomenon that the phase transition persists up to very large values of $\lambda$, where the transition moves to (large) negative $\beta$.


title: Questionable and unquestionable in the perturbation theory of non-Abelian models

abstract: We show, by explicit computation, that bare lattice perturbation theory in the two-dimensional O(n) nonlinear $\sigma$ models with superinstanton boundary conditions is divergent in the limit of an infinite number of points $|\Lambda|$. This is the analogue of David's statement that renormalized perturbation theory of these models is infrared divergent in the limit where the physical size of the box tends to infinity. We also give arguments which support the validity of the bare perturbative expansion of short-distance quantities obtained by taking the limit $|\Lambda|\to\infty$ term by term in the theory with more conventional boundary conditions such as Dirichlet, periodic, and free.


title: Gluinos on the lattice: quenched calculations

abstract: As a preparation for the numerical study of the SU(2) gauge theory with gluinos, the spectral properties of the fermion matrix and the masses of pseudoscalar and scalar states are investigated in the quenched approximation. The behaviour of the disconnected fermion diagrams on small lattices is also studied.


title: A Strong-Coupling Analysis of the Lattice $CP^{N-1}$ Models in the Presence of a $\theta$ Term

abstract: A $\theta$ term, which couples to topological charge, is added to the lattice $CP^{N-1}$ model. The strong-coupling character expansion is developed. The series for the free energy and mass gap are respectively computed to tenth order and fourth order. Several features of the strong-coupling analysis emerge. One is the loss of superconfinement. Another is that in the intermediate coupling constant region, there are indications of a transition to a deconfining phase when $\theta$ is sufficiently large. The transition is like the one which has been observed in Monte Carlo simulations of a similar lattice $CP^{N-1}$ action.


title: Improved Quark Actions for Lattice QCD

abstract: I present a brief summary of the status and prospects of improved Wilson-type quark actions for coarse lattice simulations. My conclusions are optimistic.


title: A Non-Perturbative Analysis of the Finite T Phase Transition in SU(2)xU(1) Electroweak Theory

abstract: The continuum 3d SU(2)$\times$U(1)+Higgs theory is an effective theory for a large class of 4d high-temperature gauge theories, including the minimal standard model and some of its supersymmetric extensions. We study the effects of the U(1) subgroup using lattice Monte Carlo techniques. When $g'^2/g^2$ is increased from the zero corresponding to pure SU(2)+Higgs theory, the phase transition gets stronger. However, the increase in the strength is close to what is expected perturbatively, and the qualitative features of the phase diagram remain the same as for $g'^2=0$. In particular, the first order transition still disappears for $m_H>m_{H,c}$. We measure the photon mass and mixing angle, and find that the mass vanishes in both phases within the statistical errors.


title: Perfect Lattice Actions for Staggered Fermions

abstract: We construct a perfect lattice action for staggered fermions by blocking from the continuum. The locality, spectrum and pressure of such perfect staggered fermions are discussed. We also derive a consistent fixed point action for free gauge fields and discuss its locality as well as the resulting static quark-antiquark potential. This provides a basis for the construction of (classically) perfect lattice actions for QCD using staggered fermions.


title: Three-Dimensional 3-State Potts Model Revisited With New Techniques

abstract: We report a fairly detailed finite-size scaling analysis of the first-order phase transition in the three-dimensional 3-state Potts model on cubic lattices with emphasis on recently introduced quantities whose infinite-volume extrapolations are governed `only' by exponentially small terms. In these quantities no asymptotic power series in the inverse volume are involved which complicate the finite-size scaling behaviour of standard observables related to the specific-heat maxima or Binder-parameter minima. Introduced initially for strong first-order phase transitions in q-state Potts models with ``large enough'' q, the new techniques prove to be surprisingly accurate for a q value as small as 3. On the basis of the high-precision Monte Carlo data of Alves `et al.' [Phys. Rev. B43 (1991) 5846], this leads to a refined estimate of $\beta_t = 0.550,565(10)$ for the infinite-volume transition point.


title: How to Measure the Fractal Geometry of the Relativistic Fermion Propagator

abstract: We study the geometry of propagation of relativistic fermions. We propose how to measure its quantum mechanical length. Numerical lattice results for the free propagator of Dirac-Wilson fermions yield Hausdorff dimension d_H=2 for the unit-matrix component and d_H=1 for any gamma-matrix component. A possible generalization when matter interacts with radiation is discussed.


title: Finite-Temperature Phase Structure of Lattice QCD with the Wilson Quark Action for Two and Four Flavors

abstract: We present further analyses of the finite-temperature phase structure of lattice QCD with the Wilson quark action based on spontaneous breakdown of parity-flavor symmetry. Results are reported on (i) an explicit demonstration of spontaneous breakdown of parity-flavor symmetry beyond the critical line, (ii) phase structure and order of chiral transition for the case of $N_f=4$ flavors, and (iii) approach toward the continuum limit.


title: Finite-Temperature QCD on the Lattice

abstract: Recent developments in finite-temperature studies of lattice QCD are reviewed. Topics include (i) tests of improved actions for the pure gauge system, (ii) scaling study of the two-flavor chiral transition and restoration of $U_A(1)$ symmetry with the Kogut-Susskind quark action, (iii) present understanding of the finite-temperature phase structure for the Wilson quark action. New results for finite-density QCD are briefly discussed.


title: Lattice Gauge Theory Simulations at Nonzero Chemical Potential in the Chiral Limit

abstract: We present a method of simulating lattice QCD at nonzero chemical potential in the chiral limit. By adding a weak four-fermi interaction to the standard staggered fermion SU(3) QCD action, we produce an algorithm in which the limit of massless fermions is well-behaved and physical. Using configurations at zero chemical potential, and an exact fugacity expansion of the fermion determinant, we can simulate QCD at nonzero chemical potential and evade the notorious problem of the complex action. Small lattice simulations give physical results: At strong gauge coupling the critical chemical potential \mu_c agrees with theoretical expectations and at weak gauge coupling \mu_c is nonzero in the low temperature confined phase of QCD and jumps to zero in the high temperature quark-gluon plasma phase. In all these simulations the quarks are exactly massless and there is a Goldstone pion.


title: Improvements of the local bosonic algorithm

abstract: We report on several improvements of the local bosonic algorithm proposed by M. Luescher. We find that preconditioning and over-relaxation works very well. A detailed comparison between the bosonic and the Kramers-algorithms shows comparable performance for the physical situation examined.


title: Krylov space solvers for shifted linear systems

abstract: We investigate the application of Krylov space methods to the solution of shifted linear systems of the form (A+\sigma) x - b = 0 for several values of \sigma simultaneously, using only as many matrix-vector operations as the solution of a single system requires. We find a suitable description of the problem, allowing us to understand known algorithms in a common framework and developing shifted methods basing on short recurrence methods, most notably the CG and the BiCGstab solvers. The convergence properties of these shifted solvers are well understood and the derivation of other shifted solvers is easily possible. The application of these methods to quark propagator calculations in quenched QCD using Wilson and Clover fermions is discussed and numerical examples in this framework are presented. With the shifted CG method an optimal algorithm for staggered fermions is available.


title: QCDF90: Lattice QCD with Fortran 90

abstract: We have used Fortran 90 to implement lattice QCD. We have designed a set of machine independent modules that define fields (gauge, fermions, scalars, etc...) and overloaded operators for all possible operations between fields, matrices and numbers. With these modules it is very simple to write high-level efficient programs for QCD simulations. To increase performances our modules also implements assignments that do not require temporaries, and a machine independent precision definition. We have also created a useful compression procedure for storing the lattice configurations, and a parallel implementation of the random generators. We have widely tested our program and modules on several parallel and single processor supercomputers obtaining excellent performances.


title: Infrared Renormalons and Finite Volume

abstract: We analyze the perturbative expansion of a condensate in the O(N) non-linear sigma model for large N on a two dimensional finite lattice. On an infinite volume this expansion is affected by an infrared renormalon. We extrapolate this analysis to the case of the gluon condensate of Yang-Mills theory and argue that infrared renormalons can be detected by performing perturbative studies even on relatively small lattices.


title: Lattice Formulation of the Standard Model

abstract: Combining the Kaplan surface mode approach for chiral fermions with added terms motivated by Eichten and Preskill suggests the possibility for a lattice regularization of the standard model which is finite, exactly gauge invariant, and only has physically desired states in its low energy spectrum. The conjectured scheme manifestly requires anomaly cancelation and explicitly contains baryon and lepton number violating terms.


title: Thermodynamics of Four-Flavour QCD with Improved Staggered Fermions

abstract: We have calculated the pressure and energy density in four-flavour QCD using improved fermion and gauge actions. We observe a strong reduction of finite cut-off effects in the high temperature regime, similar to what has been noted before for the SU(3) gauge theory. Calculations have been performed on $16^3\times 4$ and 16^4 lattices for two values of the quark mass, $ma = 0.05$ and 0.1. A calculation of the string tension at zero temperature yields a critical temperature $T_c/\sqrt{\sigma} = 0.407 \pm 0.010$ for the smaller quark mass value.


title: Structure of Critical Lines in Quenched Lattice QCD with the Wilson Quark Action

abstract: The structure of critical lines of vanishing pion mass for the Wilson quark action is examined in quenched lattice QCD. The numerical evidence is presented that critical lines spread into five branches beyond beta=5.6-5.7 at zero temperature. It is also shown that critical lines disappear in the deconfined phase for the case of finite temperatures.


title: Thermodynamics of lattice QCD with two light quark flavours on a 16^3 x 8 lattice II

abstract: We have extended our earlier simulations of the high temperature behaviour of lattice QCD with two light flavours of staggered quarks on a $16^3 \times 8$ lattice to lower quark mass (m_q=0.00625). The transition from hadronic matter to a quark-gluon plasma is observed at $6/g^2=5.49(2)$ corresponding to a temperature of $T_c \approx 140$MeV. We present measurements of observables which probe the nature of the quark-gluon plasma and serve to distinguish it from hadronic matter. Although the transition is quite abrupt, we have seen no indications that it is first order.


title: On the Phase Diagram of the SU(2) Adjoint Higgs Model in 2+1 Dimensions

abstract: The phase diagram is investigated for SU(2) lattice gauge theory in d=3, coupled to adjoint scalars. For small values of the quartic scalar coupling, lambda, the transition separating Higgs and confinement phases is found to be first-order, in agreement with earlier work by Nadkarni. The surface of second-order transitions conjectured by Nadkarni, however, is shown instead to correspond to crossover behaviour. This conclusion is based on a finite size analysis of the scalar mass and susceptibility. The nature of the phase transition at the termination of first-order behaviour is investigated and we find evidence for a critical point at which the scalar mass vanishes. The photon mass and confining string tension are measured and are found to be negligibly small in the Higgs phase. This is correlated with the very small density of magnetic monopoles in the Higgs phase. The string tension and photon mass rise rapidly as the crossover is traversed towards the symmetric phase.


title: Computer Aided Series Expansions for Critical Phenomena

abstract: Under quite general conditions critical phenomena can be described with high order linked cluster expansions. The coefficients of the series admit a graphical expansion that is generated with the aid of computers. Our generalization of linked cluster expansions from an infinite to a finite volume allows to perform a finite size scaling analysis. We also indicate a generalization to Dynamical Linked Cluster Expansions with possible applications to spin glasses and neural networks with coupled spin and interaction dynamics.


title: Four-dimensional Simulation of the Hot Electroweak Phase Transition with the SU(2) Gauge-Higgs Model

abstract: We study the finite-temperature electroweak phase transition of the minimal standard model within the four-dimensional SU(2) gauge-Higgs model. Monte Carlo simulations are performed for intermediate values of the Higgs boson mass in the range $50 \lesssim M_H \lesssim 100$GeV on a lattice with the temporal size $N_t=2$. The order of the transition is systematically examined using finite-size scaling methods. Behavior of the interface tension and the latent heat for an increasing Higgs boson mass is also investigated. Our results suggest that the first-order transition terminates around $M_H \sim 80$GeV.


title: Measures of critical exponents in the four dimensional site percolation

abstract: Using finite-size scaling methods we measure the thermal and magnetic exponents of the site percolation in four dimensions, obtaining a value for the anomalous dimension very different from the results found in the literature. We also obtain the leading corrections-to-scaling exponent and, with great accuracy, the critical density.


title: The equation of state for two flavor QCD at N_t=6

abstract: We calculate the two flavor equation of state for QCD on lattices with lattice spacing a=(6T)^{-1} and find that cutoff effects are substantially reduced compared to an earlier study using a=(4T)^{-1}. However, it is likely that significant cutoff effects remain. We fit the lattice data to expected forms of the free energy density for a second order phase transition at zero-quark-mass, which allows us to extrapolate the equation of state to m_q=0 and to extract the speed of sound. We find that the equation of state depends weakly on the quark mass for small quark mass.


title: Perturbation theory predictions and Monte Carlo simulations for the 2-d O(n) non-linear sigma-model

abstract: By using the results of a high-statistics (O(10^7) measurements) Monte Carlo simulation we test several predictions of perturbation theory on the O(n) non-linear sigma-model in 2 dimensions. We study the O(3) and O(8) models on large enough lattices to have a good control on finite-size effects. The magnetic susceptibility and three different definitions of the correlation length are measured. We check our results with large-n expansions as well as with standard formulae for asymptotic freedom up to 4 loops in the standard and effective schemes. For this purpose the weak coupling expansions of the energy up to 4 loops for the standard action and up to 3 loops for the Symanzik action are calculated. For the O(3) model we have used two different effective schemes and checked that they lead to compatible results. A great improvement in the results is obtained by using the effective scheme based on the energy at 3 and 4 loops. We find that the O(8) model follows very nicely (within few per mille) the perturbative predictions. For the O(3) model an acceptable agreement (within few per cent) is found.

categories: hep-lat comp-gas cond-mat hep-th nlin.CG nlin.SI solv-int


title: Computation of the one-loop Symanzik coefficients for the square action

abstract: We compute the one-loop coefficients for an alternative Symanzik improved pure gauge SU{N} lattice action (N=2 and N=3). For the standard Symanzik improved action we confirm previous results by L\"{u}scher and Weisz.


title: SU(N) gauge theories for all N in 3 and 4 dimensions

abstract: We compare the mass spectra and string tensions of SU(2), SU(3) and SU(4) gauge theories in 2+1 dimensions. We find that the ratios of masses are, to a first approximation, independent of N and that the remaining dependence can be accurately reproduced by a simple O(1/N^2) correction. This provides us with a prediction of these mass ratios for all SU(N) theories in 2+1 dimensions and demonstrates that these theories are `close' to N=\infty for N\geq 2. We also find that, when expressed in units of the dynamical length scale of the theory, the dimensionful coupling g^2 is proportional to 1/N at large N. We confirm that these theories are indeed confining in the limit N \to \infty. We describe preliminary calculations in 3+1 dimensions which indicate that the same will be true there.


title: SU(N) gauge theories for all N

abstract: We show that SU(N) gauge theories in 2+1 dimensions are close to N=\infty for N \geq 2. The dimensionful coupling, g^2, is proportional to 1/N, at large N, confirming the usual diagram-based expectation. Preliminary calculations in 3+1 dimensions indicate that the same is true there.


title: Gauge-ball spectrum of the four-dimensional pure U(1) gauge theory

abstract: We investigate the continuum limit of the gauge-ball spectrum in the four-dimensional pure U(1) lattice gauge theory. In the confinement phase we identify various states scaling with the correlation length exponent $\nu \simeq 0.35$. The square root of the string tension also scales with this exponent, which agrees with the non-Gaussian fixed point exponent recently found in the finite size studies of this theory. Possible scenarios for constructing a non-Gaussian continuum theory with the observed gauge-ball spectrum are discussed. The $0^{++}$ state, however, scales with a Gaussian value $\nu \simeq 0.5$. This suggests the existence of a second, Gaussian continuum limit in the confinement phase and also the presence of a light or possibly massless scalar in the non-Gaussian continuum theory. In the Coulomb phase we find evidence for a few gauge-balls, being resonances in multi-photon channels; they seem to approach the continuum limit with as yet unknown critical exponents. The maximal value of the renormalized coupling in this phase is determined and its universality confirmed.


title: Quantum geometry of 2d gravity coupled to unitary matter

abstract: We show that there exists a divergent correlation length in 2d quantum gravity for the matter fields close to the critical point provided one uses the invariant geodesic distance as the measure of distance. The corresponding reparameterization invariant two-point functions satisfy all scaling relations known from the ordinary theory of critical phenomena and the KPZ exponents are determined by the power-like fall off of these two-point functions. The only difference compared to flat space is the appearance of a dynamically generated fractal dimension d_h in the scaling relations. We analyze numerically the fractal properties of space-time for Ising and three-states Potts model coupled to 2d dimensional quantum gravity using finite size scaling as well as small distance scaling of invariant correlation functions. Our data are consistent with d_h=4, but we cannot rule out completely the conjecture d_H = -2\alpha_1/\alpha_{-1}, where \alpha_{-n} is the gravitational dressing exponent of a spin-less primary field of conformal weight (n+1,n+1). We compute the moments <L^n> and the loop-length distribution function and show that the fractal properties associated with these observables are identical, with good accuracy, to the pure gravity case.


title: Universal Amplitude Ratios in the 3D Ising Model

abstract: We present a high precision Monte Carlo study of various universal amplitude ratios of the three dimensional Ising spin model. Using state of the art simulation techniques we studied the model close to criticality in both phases. Great care was taken to control systematic errors due to finite size effects and correction to scaling terms. We obtain $C_+/C_-=4.75(3)$, $f_{+,2nd}/f_{-,2nd}=1.95(2)$ and $u^*=14.3(1)$. Our results are compatible with those obtained by field theoretic methods applied to the $\phi^4$ theory and high and low temperature series expansions of the Ising model. The mismatch with a previous Montecarlo study by Ruge et al. remains to be understood.


title: On Lattice QCD with Many Flavors

abstract: We discuss the confining and chiral-symmetry breaking properties of QCD with a large number of flavors $N_f$. In a Monte Carlo simulation of QCD with $N_f =16$ staggered fermions, we find clear evidence of a first order bulk phase transition which separates phases with broken and unbroken chiral symmetry. This is consistent with extrapolations of earlier studies with smaller $N_f$, and is also as expected from general arguments. Thus, even when the perturbative renormalization group flow has a new infrared stable fixed point near the origin, lattice artifacts induce chiral symmetry breaking, and presumably confinement, at sufficiently strong coupling.


title: Faster Fermions in the Tempered Hybrid Monte Carlo Algorithm

abstract: Tempering is used to change the quark mass while remaining in equilibrium between the trajectories of a standard hybrid Monte Carlo simulation of four flavours of staggered fermions. The algorithm is faster for small enough quark masses, and particularly so when more than one mass is required.


title: Random Walk Representation of the Lattice Fermionic Propagators and the Quark Model

abstract: A representation of the continuum fermionic propagator as a sum of directed random walks on a lattice is presented. Also a random walk representation for the lattice fermionic propagators is developed for the case of the naive, the Wilson, and the Kogut-Susskind fermions. For the naive fermions the phenomenon of fermion doubling appears as having $2^{D}$ distinct spin factors being associated with a single path in $D$-dimensions. In the case of the Wilson and the Kogut-Susskind fermions, in the naive continuum limit, the path integral representation coincides with the path integral representation for the continuum fermionic propagator. Using this representation the Green's functions of lattice QCD involving quark operators are written as a sum over the paths of valence quark, the gauge fields and the sea quarks being integrated out. Possible advantages of such a representation are illustrated by showing how one can use numerical simulations to obtain a heuristic insight into the relationship between QCD and the constituent quark model.


title: Monte Carlo Study of Pure-Phase Cumulants of 2D q-State Potts Models

abstract: We performed Monte Carlo simulations of the two-dimensional q-state Potts model with q=10, 15, and 20 to study the energy and magnetization cumulants in the ordered and disordered phase at the first-order transition point $\beta_t$. By using very large systems of size 300 x 300, 120 x 120, and 80 x 80 for q=10, 15, and 20, respectively, our numerical estimates provide practically (up to unavoidable, but very small statistical errors) exact results which can serve as a useful test of recent resummed large-q expansions for the energy cumulants by Bhattacharya `et al.' [J. Phys. I (France) 7 (1997) 81]. Up to the third order cumulant and down to q=10 we obtain very good agreement, and also the higher-order estimates are found to be compatible.


title: Topology of the SU(2) vacuum: a lattice study using improved cooling

abstract: We study the topological structure of the SU(2) vacuum at zero temperature: topological susceptibility, size, shape and distance distributions of the instantons. We use a cooling algorithm based on an improved action with scale invariant instanton solutions. This algorithm needs no monitoring or calibration, has an inherent cut off for dislocations and leaves unchanged instantons at physical scales. The physical relevance of our results is checked by studying the scaling and finite volume dependence. We obtain a susceptibility of (200(15) MeV)^4. The instanton size distribution is peaked around 0.43fm, and the distance distribution indicates a homogeneous, random spatial structure.


title: Four-Dimensional N=1 Supersymmetric Yang-Mills Theory on the Lattice without Fine-Tuning

abstract: We propose a method to formulate four-dimensional N=1 super Yang-Mills theory on the lattice without fine-tuning. We first show that four-dimensional Weyl fermion in a real representation, which is equivalent to Majorana fermion, can be formulated using the domain wall approach with an addition of a Majorana mass term only for the unwanted mirror fermion. This formalism has manifest gauge invariance. Fermion number conservation is violated only by the additional Majorana mass term for the mirror fermion and the violation is propagated to the physical fermion sector through anomalous currents. Due to this feature, the formalism, when applied to the gluino in the present case, ensures the restoration of supersymmetry in the continuum limit without fine-tuning, unlike the proposal by Curci and Veneziano.


title: The Svetitsky-Yaffe conjecture for the plaquette operator

abstract: According to the Svetitsky-Yaffe conjecture, a (d+1)-dimensional pure gauge theory undergoing a continuous deconfinement transition is in the same universality class as a d-dimensional statistical model with order parameter taking values in the center of the gauge group. We show that the plaquette operator of the gauge theory is mapped into the energy operator of the statistical model. For d=2, this identification allows us to use conformal field theory techniques to evaluate exactly the correlation functions of the plaquette operator at the critical point. In particular, we can evaluate exactly the plaquette expectation value in presence of static sources, which gives some new insight in the structure of the color flux tube in mesons and baryons.


title: Series Expansions for the Massive Schwinger Model in Hamiltonian lattice theory

abstract: It is shown that detailed and accurate information about the mass spectrum of the massive Schwinger model can be obtained using the technique of strong-coupling series expansions. Extended strong-coupling series for the energy eigenvalues are calculated, and extrapolated to the continuum limit by means of integrated differential approximants, which are matched onto a weak-coupling expansion. The numerical estimates are compared with exact results, and with finite-lattice results calculated for an equivalent lattice spin model with long-range interactions. Both the heavy fermion and the light fermion limits of the model are explored in some detail.


title: The Three Dimensional Thirring Model for Small N_f

abstract: We formulate the three dimensional Thirring model on a spacetime lattice and study it for various even numbers of fermion flavors N_f by Monte Carlo simulation. We find clear evidence for spontaneous chiral symmetry breaking at strong coupling, contradicting the predictions of the 1/N_f expansion. The critical point appears to correspond to an ultra-violet fixed point of the renormalisation group; a fit to a RG-inspired equation of state in the vicinity of the fixed point yields distinct critical exponents for N_f=2 and N_f=4, while no fit is found for N_f=6, suggesting there is a critical number N_fc<6 beyond which no chiral symmetry breaking occurs. The spectrum of the N_f=2 theory is studied; the states examined vary sharply but continuously across the transition.


title: Scaling of cluster fluctuations in two-dimensional q=5 and 7 state Potts models

abstract: The scaling behavior of fluctuations in cluster size is studied in q=5 and 7 state Potts models. This quantity exhibits scaling behavior on small lattices where the scaling of local operators like energy fluctuations and Binder cumulant can not be expected.


title: Study of Cluster Fluctuations in Two-dimensional q-State Potts Model

abstract: The two-dimensional Potts Model with 2 to 10 states is studied using a cluster algorithm to calculate fluctuations in cluster size as well as commonly used quantities like equilibrium averages and the histograms for energy and the order parameter. Results provide information about the variation of cluster sizes depending on the temperature and the number of states. They also give evidence for first-order transition when energy and the order parameter related measurables are inconclusive on small size lattices.


title: A toy model test of a new algorithm for bozonization of fermion determinants

abstract: Different aspects of bozonization algorithm proposed by Slavnov are tested by numerical simulations of a one dimensional toy model.


title: Monte Carlo Study of the 3D Thirring Model

abstract: I review three different non-perturbative approaches to the three dimensional Thirring model: the 1/N_f expansion, Schwinger-Dyson equations, and Monte Carlo simulation. Simulation results are presented to support the existence of a non-perturbative fixed point at a chiral symmetry breaking phase transition for N_f=2 and 4, but not for N_f=6. Spectrum calculations for $N_f=2$ reveal conventional level ordering near the transition.


title: Coexistence of monopoles and instantons for different topological charge definitions and lattice actions

abstract: We compute instanton sizes and study correlation functions between instantons and monopoles in maximum abelian projection within SU(2) lattice QCD at finite temperature. We compare several definitions of the topological charge, different lattice actions and methods of reducing quantum fluctuations. The average instanton size turns out to be $\sigma \approx 0.2$ fm. The correlation length between monopoles and instantons is $\zeta \approx 0.25$ fm and hardly affected by lattice artifacts as dislocations. We visualize several specific gauge field configurations and show directly that there is an enhanced probability for finding monopole loops in the vicinity of instantons. This feature is independent of the topological charge definition used.


title: A numerical test of the continuum index theorem on the lattice

abstract: The overlap formalism of chiral fermions provides a tool to measure the index, Q, of the chiral Dirac operator in a fixed gauge field background on the lattice. This enables a numerical measurement of the probability distribution, p(Q), in Yang-Mills theories. We have obtained an estimate for p(Q) in pure SU(2) gauge theory by measuring Q on 140 independent gauge field configurations generated on a 12^4 lattice using the standard single plaquette Wilson action at a coupling of beta=2.4. This distribution is in good agreement with a recent measurement [8] of the distribution of the topological charge on the same lattice using the same coupling and the same lattice gauge action. In particular we find <Q^2>=3.3(4) to be compared with <Q^2> = 3.9(5) found in [8]. The good agreement between the two distributions is an indication that the continuum index theorem can be carried over in a probabilistic sense on to the lattice.


title: Instantons and monopoles in lattice QCD

abstract: We analyze the interplay of topological objects in four-dimensional QCD on the lattice. The distributions of color magnetic monopoles in the maximum abelian gauge are computed around instantons in both pure and full QCD. We find an enhanced probability for monopoles inside the core of an instanton on gauge field average. This feature is independent of the topological charge definition used. For specific gauge field configurations we visualize the situation graphically. Moreover we investigate how monopole loops and instantons are correlated with the chiral condensate. Strong evidence is found that clusters of the quark condensate and topological objects coexist locally on individual configurations.


title: Possible evidence for the mass shift of $\eta$' meson at finite temperature

abstract: Topological charge distributions of SU(3) gauge theory at finite temperature are calculated on the lattice with high statistics in a manner of free from the uncertainty of $\beta$-function. Clear temperature dependence of the topological charge distribution is obtained. Combining the present result with the recent study about the relation between the masses of $\eta$ meson system and the gauge field topology, we obtain $m_{\eta'}(T)/m_{\eta'}(0) = 0.86 \pm 0.02$ and $0.69 \pm 0.01$ for $T=0.75 T_c$ and $T=0.93 T_c$, respectively. The result shows the first clear evidence of the mass shift of $\eta'$ meson at finite temperature.


title: Questionable Arguments for the Correctness of Perturbation Theory in Non-Abelian Models

abstract: We analyze the arguments put forward recently by Niedermayer et al in favor of the correctness of conventional perturbation theory in non-Abelian models and supposedly showing that our super-instanton counterexample was sick. We point out that within their own set of assumptions, the proof of Niedermayer et al regarding the correctness of perturbation theory is incorrect and provide a correct proof under more restrictive assumptions. We reply also to their claim that the S-matrix bootstrap approach of Balog et al supports the existence of asymptotic freedom in the O(3) model.


title: Fermionic Monte Carlo algorithms for lattice QCD

abstract: The increase with time of computer resources devoted to simulations of full QCD is spectacular. Yet the reduction of systematic errors is comparatively slow. This is due to the algorithmic complexity of the problem. I review, in elementary terms, the origin of this complexity, and estimate it for 3 exact fermion algorithms.


title: On the evaluation of universal non-perturbative constants in O(N) sigma models

abstract: We investigate the relation between on-shell and zero-momentum non-perturbative quantities entering the parametrization of the two-point Green's function of two-dimensional non-linear O(N) sigma models. We present accurate estimates of ratios of mass-scales and renormalization constants, obtained by an analysis of the strong-coupling expansion of the two-point Green's function. These ratios allow to connect exact known on-shell results with typical zero-momentum lattice evaluations. Our results are supported by the 1/N-expansion.


title: Z(3) Interfaces in Lattice Gauge Theory

abstract: A study is made of properties of the Z(3) interface which forms between the different ordered phases of pure SU(3) gauge theory above a critical temperature. The theory is simulated on a (2+1)-D lattice at various temperatures above this critical point. At high temperatures, the interface tension is shown to agree well with the prediction of perturbation theory. Near the critical temperature, the interface behaviour is characterised by various displacement moments, and modelled by an interacting scalar field theory. This thesis is provided for reference, as it gives full details of the computational and statistical methods outlined only briefly in preprints hep-lat/9605040 and hep-lat/9607005.


title: Monopoles and deconfinement transition in SU(2) lattice gauge theory

abstract: We investigate SU(2) lattice gauge theory in four dimensions in the maximally abelian projection. Studying the effects on different lattice sizes we show that the deconfinement transition of the fields and the percolation transition of the monopole currents in the three space dimensions are precisely related. To arrive properly at this result the uses of a mathematically sound characterization of the occurring networks of monopole currents and of an appropriate method of gauge fixing turn out to be crucial. In addition we investigate detailed features of the monopole structure in time direction.


title: On-shell improved lattice QCD with staggered fermions

abstract: By using Symanzik's improvement program, we study on-shell improved lattice QCD with staggered fermions. We find that there are as many as 15 independent lattice operators of dimension of six~(including both gauge and fermion operators) which must be added to the unimproved action to form an $O(a^2)$ improved action. Among them, the total number of dimension-6 gauge operators and fermion bilinears is 5. The other 10 terms are four-fermion operators. At the tree-level and tadpole-improved tree-level, all the 10 four fermion operators are absent.


title: Quark Confinement in C-periodic Cylinders at Temperatures above T_c

abstract: Due to the Gauss law, a single quark cannot exist in a periodic volume, while it can exist with C-periodic boundary conditions. In a C-periodic cylinder of cross section A = L_x L_y and length L_z >> L_x, L_y containing deconfined gluons, regions of different high temperature Z(3) phases are aligned along the z-direction, separated by deconfined- deconfined interfaces. In this geometry, the free energy of a single static quark diverges in proportion to L_z. Hence, paradoxically, the quark is confined, although the temperature T is larger than T_c. At T around T_c, the confined phase coexists with the three deconfined phases. The deconfined-deconfined interfaces can be completely or incompletely wet by the confined phase. The free energy of a quark behaves differently in these two cases. In contrast to claims in the literature, our results imply that deconfined-deconfined interfaces are not Euclidean artifacts, but have observable consequences in a system of hot gluons.


title: Chiral perturbation theory for K+ to pi+ pi0 decay in the continuum and on the lattice

abstract: In this paper we use one-loop chiral perturbation theory in order to compare lattice computations of the K+ to pi+ pi0 decay amplitude with the experimental value. This makes it possible to investigate three systematic effects that plague lattice computations: quenching, finite-volume effects, and the fact that lattice computations have been done at unphysical values of the quark masses and pion external momenta (only this latter effect shows up at tree level). We apply our results to the most recent lattice computation, and find that all three effects are substantial. We conclude that one-loop corrections in chiral perturbation theory help in explaining the discrepancy between lattice results and the real-world value. We also revisit B_K, which is closely related to the K+ to pi+ pi0 decay amplitude by chiral symmetry.


title: A Polynomial Hybrid Monte Carlo Algorithm

abstract: We present a simulation algorithm for dynamical fermions that combines the multiboson technique with the Hybrid Monte Carlo algorithm. We find that the algorithm gives a substantial gain over the standard methods in practical simulations. We point out the ability of the algorithm to treat fermion zeromodes in a clean and controllable manner.


title: Numerical study of plasmon properties in the SU(2)-Higgs model

abstract: Using the (effective) classical approximation, we compute numerically time-dependent correlation functions in the SU(2)-Higgs model around the electroweak phase transition, for $m_H \approx m_W$. The parameters of the classical model have been determined previously by the dimensional reduction relations for time-independent correlators. The $H$ and $W$ correlation functions correspond to gauge invariant fields. They show damped oscillatory behavior from which we extract frequencies $\om$ and damping rates $\gm$. In the Higgs phase the damping rates have roughly the values obtained in analytic calculations in the quantum theory. In the plasma phase (where analytic estimates for gauge invariant fields are not available), the damping rate associated with $H$ is an order of magnitude larger than in the Higgs phase, while the $W$ correlator appears to be overdamped, with a small rate. The frequency $\om_H$ shows a clear dip at the transition. The results are approximately independent of the lattice spacing, but this appears to be compatible with the lattice spacing dependence expected from perturbation theory.


title: Phase Structures of U(2) Gauge Theory with $\theta$-Term in 2 dimensions

abstract: U(2) lattice gauge theory with $\theta$-term in 2 space-time dimensions is investigated. It has non-Abelian real action and Abelian( U(1) type) imaginary action. The imaginary action is defined as the standard $\theta$-term. As the effect of renormalization group (RG) transformation, non-Abelian imaginary action is induced. After many steps of RG transformation, non-Abelian part will die away. After several steps of RG transformations, renormalized action approaches so called heat kernel action. Phase transition is found at $\theta =\pi$ only.


title: Pion loops in quenched Quantum Chromodynamics

abstract: We calculate the divergences of the generating functional of quenched Chiral Perturbation Theory to one loop for a generic number of flavours. The flavour number dependence of our result enlightens the mechanism of quark loop cancellation in the quenched effective theory for any Green function or S matrix element. We also apply our results to $\pi \pi$ scattering and evaluate the coefficient of the chiral log in the S-wave scattering lengths for the quenched case.


title: Measuring the interface tension when the electroweak phase transition becomes weak

abstract: We measure the interface tension near the phase transition endpoint of the 3d SU(2)-Higgs model. The tunnel correlation length method is used and compared to other approaches. A modified scaling behaviour for the mass gap as function of the transverse area is proposed.


title: The Higgs mechanism on the lattice

abstract: The lattice regularization of the Higgs sector of the standard model is summarized. The triviality bound and vacuum instability bound are described. The question of chiral gauge theories is discussed. Some aspects of the numerical simulations of the electroweak phase transition are considered.


title: A High Statistics Lattice Calculation of Heavy-Light Meson Decay Constants

abstract: We present a high statistics study of the D- and B-meson decay constants. The results were obtained by using the Clover and Wilson lattice actions at two different values of the lattice spacing $a$, corresponding to $\beta=6.0$ and 6.2. After a careful analysis of the systematic errors present in the extraction of the physical results, by assuming quite conservative discretization errors, we find $f_{D_s}=237 \pm 16$ MeV, $f_{D} = 221 \pm 17$ MeV ($f_{D_s}/f_D=1.07(4)$), $f_{B_s} = 205 \pm 35$ MeV, $f_{B} = 180 \pm 32$ MeV ($f_{B_s}/f_B=1.14(8)$), in good agreement with previous estimates.


title: QCD Thermodynamics with an Improved Lattice Action

abstract: We have investigated QCD with two flavors of degenerate fermions using a Symanzik-improved lattice action for both the gauge and fermion actions. Our study focuses on the deconfinement transition on an $N_t=4$ lattice. Having located the thermal transition, we performed zero temperature simulations nearby in order to compute hadronic masses and the static quark potential. We find that the present action reduces lattice artifacts present in thermodynamics with the standard Wilson (gauge and fermion) actions. However, it does not bring studies with Wilson-type quarks to the same level as those using the Kogut--Susskind formulation.


title: Three-dimensional lattice U(1) gauge-Higgs model at low $m_H$

abstract: We study the non-compact version of the U(1) gauge-Higgs model in three dimensions for $m_H = 30 GeV.$ We found that, using this formulation, rather modest lattices approach quite well the infinite volume behaviour.The phase transition is first order, as expected for this Higgs mass. The latent heat (in units of $T_{cr}^4$) is compatible with the predictions of the two-loop effective potential; it is an order of magnitude less than the corresponding SU(2) value. The transition temperature and $<\varphi^* \varphi>$ in units of the critical temperature are also compatible with the perturbative results.


title: Hubbard Model with L\"uscher fermions

abstract: We study the basic features of the two-dimensional quantum Hubbard Model at half-filling by means of the L\"uscher algorithm and the algorithm based on direct update of the determinant of the fermionic matrix. We implement the L\"uscher idea employing the transfer matrix formalism which allows to formulate the problem on the lattice in $(2+1)$ dimensions. We discuss the numerical complexity of the L\"uscher technique, systematic errors introduced by polynomial approximation and introduce some improvements which reduce long autocorrelations. In particular we show that preconditioning of the fermionic matrix speeds up the algorithm and extends the available range of parameters. We investigate the magnetic and the one-particle properties of the Hubbard Model at half-filling and show that they are in qualitative agreement with the existing Monte Carlo data and the mean-field predictions.


title: Fermion determinants in matrix models of QCD at nonzero chemical potential

abstract: The presence of a chemical potential completely changes the analytical structure of the QCD partition function. In particular, the eigenvalues of the Dirac operator are distributed over a finite area in the complex plane, whereas the zeros of the partition function in the complex mass plane remain on a curve. In this paper we study the effects of the fermion determinant at nonzero chemical potential on the Dirac spectrum by means of the resolvent, G(z), of the QCD Dirac operator. The resolvent is studied both in a one-dimensional U(1) model (Gibbs model) and in a random matrix model with the global symmetries of the QCD partition function. In both cases we find that, if the argument z of the resolvent is not equal to the mass m in the fermion determinant, the resolvent diverges in the thermodynamic limit. However, for z =m the resolvent in both models is well defined. In particular, the nature of the limit $z \rightarrow m$ is illuminated in the Gibbs model. The phase structure of the random matrix model in the complex m and \mu-planes is investigated both by a saddle point approximation and via the distribution of Yang-Lee zeros. Both methods are in complete agreement and lead to a well-defined chiral condensate and quark number density.