title: Interpreting the lattice monopoles in the continuum terms

abstract: We review briefly current interpretation of the lattice monopoles, defined within the Maximal Abelian Projection, in terms of the continuum theory. We emphasize, in particular, that the lattice data, at the presently available lattices, indicate that the monopoles are associated with singular fields. This note is prompted by a recent analysis hep-lat/0211005 which is based on an implicit assumption that the fields are regular.


title: Non-Perturbative Determination of $c_{\rm SW}$ in Three-flavor Dynamical QCD

abstract: We present a fully non-perturbative determination of the $O(a)$ improvement coefficient $c_{\rm SW}$ in three-flavor dynamical QCD for the RG improved as well as the plaquette gauge actions, using the Schr\"odinger functional scheme. Results are compared with one-loop estimates at weak gauge coupling.


title: QCD and the eta prime Mass: Instantons or Confinement?

abstract: We argue that lattice calculations of the $\eta'$ mass in QCD with $N_c=2$ colors performed at non-zero baryon chemical potential can be used to study the mechanism responsible for the mass of the $\eta'$. QCD with two colors is an ideal laboratory because it exhibits confinement, chiral symmetry breaking and a would-be $U(1)_A$ Goldstone boson at all densities. Since the instanton density and the confinement scale vary with density in a very different way, instantons are clearly distinguishable from other possible mechanisms. There is an instanton prediction for the $\eta'$ mass at large density that can be compared to lattice results. The density dependence of the instanton contribution is a simple consequence of the integer topological charge carried by the instanton. We also argue that $N_c=3$ color QCD at finite isospin density can be used in order to study the origin of OZI-violation in the scalar sector.


title: Lattice Perturbation Theory

abstract: The consideration of quantum fields defined on a spacetime lattice provides computational techniques which are invaluable for studying gauge theories nonperturbatively from first principles. Perturbation theory is an essential aspect of computations on the lattice, especially for investigating the behavior of lattice theories near the continuum limit. Particularly important is its role in connecting the outcome of Monte Carlo simulations to continuum physical results. For these matchings the calculation of the renormalization factors of lattice matrix elements is required. In this review we explain the main methods and techniques of lattice perturbation theory, focusing on the cases of Wilson and Ginsparg-Wilson fermions. We will illustrate, among other topics, the peculiarities of perturbative techniques on the lattice, the use of computer codes for the analytic calculations and the computation of lattice integrals. Discussed are also methods for the computation of 1-loop integrals with very high precision. The review presents in a pedagogical fashion also some of the recent developments in this kind of calculations. The coordinate method of L\"uscher and Weisz is explained in detail. Also discussed are the novelties that Ginsparg-Wilson fermions have brought from the point of view of perturbation theory. Particular emphasis is given throughout the paper to the role of chiral symmetry on the lattice and to the mixing of lattice operators under renormalization. The construction of chiral gauge theories regularized on the lattice, made possible by the recent advances in the understanding of chiral symmetry, is also discussed. Finally, a few detailed examples of lattice perturbative calculations are presented.


title: The Second Moment of the Pion Light Cone Wave Function

abstract: We present a preliminary result for second moment of the light cone wave function of the pion. This parameter is the subject of a discrepancy between theoretical predictions (coming from lattice and sum rules) and a recent experimental result (that remarkably agrees with purely perturbative predictions). In this work we exploit lattice hypercubic symmetries to remove power divergences and, moreover, implement a full 1-loop matching for all the contributing operators.


title: The string spectrum from large Wilson loops

abstract: We look at energies of the low lying states of the hadronic string in three dimensional SU(2) lattice gauge theory by forming correlation matrices among different sources. We are able to go to previously inaccessible time separations. This is made possible by using a new algorithm proposed by L\"{u}scher and Weisz which lets us measure the exponentially small values of large Wilson loops with sufficient accuracy.


title: The Quark-Gluon Mixed Condensate g<\bar{q} sigma G q> in SU(3)c Quenched Lattice QCD

abstract: Using the SU(3)c lattice QCD with the Kogut-Susskind fermion at the quenched level, we study the quark-gluon mixed condensate g<\bar{q} sigma G q>, which is another chiral order parameter. For each current quark mass of m_q=21, 36, 52 MeV, we generate 100 gauge configurations in the 16^4 lattice with \beta = 6.0, and perform the measurement of the mixed condensate at 16 points in each gauge configuration. Using the 1600 data for each m_q, we find the ratio between the mixed condensate and the quark condensate, m_0^2 = g<\bar{q} sigma G q>/ <\bar{q}q> \simeq 2.5 GeV^2 at the lattice scale in the chiral limit. The large value of the mixed condensate suggests its importance in the operator product expansion in QCD.


title: Flavor Singlet Meson Mass in the Continuum Limit in Two-Flavor Lattice QCD

abstract: We present results for the mass of the eta-prime meson in the continuum limit for two-flavor lattice QCD, calculated on the CP-PACS computer, using a renormalization-group improved gauge action, and Sheikholeslami and Wohlert's fermion action with tadpole-improved csw. Correlation functions are measured at three values of the coupling constant beta corresponding to the lattice spacing a approx. 0.22, 0.16, 0.11 fm and for four values of the quark mass parameter kappa corresponding to mpi over mrho approx. 0.8, 0.75, 0.7 and 0.6. For each beta, kappa pair, 400-800 gauge configurations are used. The two-loop diagrams are evaluated using a noisy source method. We calculate eta-prime propagators using local sources, and find that excited state contributions are much reduced by smearing. A full analysis for the smeared propagators gives metaprime=0.960(87)+0.036-0.248 GeV, in the continuum limit, where the second error represents the systematic uncertainty coming from varying the functional form for chiral and continuum extrapolations.


title: Results for the Etaprime Mass from Two-Flavor Lattice QCD

abstract: We present results for the mass of the etaprime meson for two-flavor lattice QCD in the continuum limit, calculated on the CP-PACS computer, using an RG-improved gauge action and clover fermion action with tadpole-improved csw. Measurements are made at three couplings corresponding to a approx. 0.22, 0.16, 0.11 fm for four quark masses corresponding to mpi over mrho approx. 0.8, 0.75, 0.7, 0.6. Thw two-loop diagrams are evaluated using a noisy source method. Quark smearing for both one- and two- loop diagrams is successfully applied to obtain ground state signals in the etaprime channel. We obtain metaprime=0.960(87)+0.036-0.286GeV in the continuum limit, where the second error represents the systematic uncertainty coming from varying the functional form for chiral and continuum extrapolations.


title: Speeding up Lattice QCD simulations with clover-improved Wilson Fermions

abstract: We apply a recent proposal to speed up the Hybrid-Monte-Carlo simulation of systems with dynamical fermions to two flavour QCD with clover-improvement. The basic idea of our proposal is to split the fermion matrix into two factors with a reduced condition number each. In the effective action, for both factors a pseudo-fermion field is introduced. For our smallest quark masses we see a speed-up of more than a factor of two compared with the standard algorithm.


title: Finite-Volume Two-Pion Amplitudes in the I=0 Channel

abstract: We perform a calculation in one-loop chiral perturbation theory of the two-pion matrix elements and correlation functions of an I=0 scalar operator, in finite and infinite volumes for both full and quenched QCD. We show that major difficulties arise in the quenched theory due to the lack of unitarity. Similar problems are expected for quenched lattice calculations of $K \to \pi \pi$ amplitudes with $\Delta I=1/2$. Our results raise the important question of whether it is consistent to study $K\to\pi\pi$ amplitudes beyond leading order in chiral perturbation theory in quenched or partially quenched QCD.


title: Abelian representation for nonabelian Wilson loops and the Non - Abelian Stokes theorem on the lattice

abstract: We derive the Abelian - like expression for the lattice SU(N) Wilson loop in arbitrary irreducible representation. The continuum Abelian representation of the SU(N) Wilson loop (for the loop without selfintersections) that has been obtained by Diakonov and Petrov appears to be a continuum limit of this expression. We also obtain the lattice variant of a non - Abelian Stokes theorem and present the explicit expression for the matrix $\cal H$ used in the Diakonov - Petrov approach.


title: Mass Generation without Phase Coherence at Nonzero Temperature

abstract: We present results from numerical simulations of the 2+1d SU(2)xSU(2) Nambu - Jona-Lasinio model with N_f=4 fermion flavours at zero and nonzero temperature T. At zero temperature, critical exponents are extracted from the scaling of the order parameter and fermion mass and are found to be consistent with next-to-leading order predictions of the 1/N_f expansion. At nonzero temperature we observe fermion mass generation despite the lack of chiral symmetry breaking, which is forbidden by the Colemann-Mermin-Wagner theorem for all T>0. We study the effects of lattice discretisation and finite volume on the dynamically generated fermion mass. By study of the lattice dispersion relation we also show that in the hot phase there is no significant temperature induced modification to the speed of light. Studies of the equation of state are made by measuring the pressure as a function of temperature and comparison is made with large-N_f predictions.


title: Fine tuned vortices in lattice SU(2) gluodynamics

abstract: We report measurements of the action associated with center vortices in SU(2) pure lattice gauge theory. In the lattice units the excess of the action on the plaquettes belonging to the vortex is approximately a constant, independent on the lattice spacing 'a'. Therefore the action of the center vortex is of order 'A/a^2', where 'A' is its area. Since the area 'A' is known to scale in the physical units, the measurements imply that the suppression due to the surface action is balanced, or fine tuned to the entropy factor which is to be an exponential of 'A/a^2'.


title: Monopole creation operators as confinement--deconfinement order parameters

abstract: We study numerically two versions of the monopole creation operators proposed by Frohlich and Marchetti. The disadvantage of the old version of the monopole creation operator is due to visibility of the Dirac string entering the definition of the creation operator in the theories with coexisting electric and magnetic charges. This problem does not exist for the new creation operator which is rather complicated. Using the Abelian Higgs model with a compact gauge field we show that both definitions of the monopole creation operator can serve as order parameters for the confinement--deconfinement phase transition. The value of the monopole condensate for the old version depends on the length of Dirac string. However, as soon as the length is fixed the old operator certainly discriminates between the phases with condensed and non--condensed monopoles.


title: More on String Breaking in the 3D Abelian Higgs Model: the Photon Propagator

abstract: We study the Landau gauge photon propagator in the three--dimensional Abelian Higgs model with compact gauge field and fundamentally charged matter in the London limit. The total gauge field is split into singular and regular parts. On the confinement side of the string breaking crossover the momentum dependence of the total propagator is characterized by an anomalous dimension similarly to 3D compact QED. At the crossover and throughout the Higgs region the anomalous dimension disappears. This result perfectly agrees with recent observations that the monopole--antimonopole plasma leads to nonzero anomalous dimension and the presence of the matter fields causes monopole pairing into dipole bound states. The Yukawa mass characterizing the propagator part from regular gauge fields is non-vanishing at the Higgs side and coincides with the mass found for the total propagator. The regular gauge field without anomalous dimension becomes massless at the crossover and in the confinement region.


title: Quark-Gluon Mixed Condensate (anti-Q sigma G Q) from Lattice QCD

abstract: We study the quark-gluon mixed condensate g<\bar{q} sigma G q>, which is another chiral order parameter, in SU(3)c lattice QCD with the Kogut-Susskind fermion at the quenched level. Using 100 gauge configurations on the 16^4 lattice with \beta = 6.0, we measure the mixed condensate at 16 points in each gauge configuration for each current-quark mass of m_q=21, 36, 52 MeV. From the 1600 data for each m_q, we find the ratio between the mixed condensate and the quark condensate, m_0^2 = g<\bar{q} sigma G q> / <\bar{q}q> \simeq 2.5 GeV^2 at the lattice scale of 1/a \simeq 2 GeV in the chiral limit. The large value of the mixed condensate suggests its importance in the operator product expansions in QCD. We also show our preliminary results of the mixed condensate at finite temperature.


title: Charmonium near the deconfining transition on the lattice

abstract: We study the charmonium properties at finite temperature using quenched lattice QCD simulations. Although a simple potential model analysis indicates no bound state at $T>1.05T_c$, our analyses of the spatial correlation between quark and anti-quark and the spectral function indicate that a bound-state-like structure may survive even above $T_c$.


title: On the effects of (partial) quenching on penguin contributions to K-> pi pi

abstract: Recently, we pointed out that chiral transformation properties of strong penguin operators change in the transition from unquenched to (partially) quenched QCD. As a consequence, new penguin-like operators appear in the (partially) quenched theory, along with new low-energy constants, which should be interpreted as a quenching artifact. Here, we extend the analysis to the contribution of the new low-energy constants to the K^0 -> pi^+ pi^- amplitude, at leading order in chiral perturbation theory, and for arbitrary (momentum non-conserving) kinematics. Using these results, we provide a detailed discussion of the intrinsic systematic error due to this (partial) quenching artifact. We also give a simple recipe for the determination of the leading-order low-energy constant parameterizing the new operators in the case of strong $LR$ penguins.


title: Partially Quenched Chiral Perturbation Theory and a Massless Up Quark: A Lattice Calculation of the Light-Quark-Mass Ratio

abstract: The nontrivial topological structure of the QCD gauge vacuum generates a CP breaking term in the QCD Lagrangian. However, measurements of the neutron electric dipole moment have demonstrated that the term's coefficient is unnaturally small, a dilemma known as the strong CP problem. A massless up quark is one potential solution, as the term could then be absorbed through chiral rotations. Through the light-quark-mass ratio m_u / m_d, leading order Chiral Perturbation Theory appears to rule this scenario out. However, the Kaplan-Manohar ambiguity demonstrates that certain strong next-to-leading order corrections are indistinguishable from an up quark mass. Only a direct calculation of the Gasser-Leutwyler coefficient combination 2 L_8 - L_5 can resolve the issue. We carry out such a calculation, using partially quenched N_f = 3 staggered fermions and hypercubic blocking, and make a quantitative assessment of our systematic error. We find 2 L_8 - L_5 = (0.22 +/- 0.14) x 10^-3, which corresponds to a light-quark-mass ratio of m_u / m_d = 0.408 +/- 0.035. Thus, the massless-up-quark solution to the strong CP problem is ruled out. This is the first calculation of its type to use a physical number of light quarks, N_f = 3, and the first determination of 2 L_8 - L_5 to include a comprehensive study of statistical error.


title: The Thermal Glueball at Finite Temperature in SU(3) Anisotropic Lattice QCD

abstract: The thermal glueball is studied at finite temperature by using SU(3) anisotropic lattice QCD with beta_{lat}=6.25, the renormalized anisotropy xi = a_s/a_t = 4 over the lattice of the size 20^3xN_t with various N_t at the quenched level. While the narrow peak ansatz leads to the significant polemass reduction of about 300 MeV near the critical temperature T_c, Breit-Wigner ansatz which can take into account the possible appearance of the thermal width leads to the significant thermal width broadening of about 300 MeV with a modest reduction of the peak center of about 100 MeV.


title: The nonperturbative quark-gluon vertex

abstract: We show results for the quark-gluon vertex in the Landau gauge, using a mean-field improved Sheikholeslami-Wohlert fermion action. We compute all the three non-zero form factors of the vertex at zero gluon momentum, and compare them to the abelian vertex. The quark mass dependence of the vertex is also investigated and found to be negligible for the range of masses considered.


title: Numerical techniques for lattice QCD in the $\epsilon$--regime

abstract: In lattice QCD it is possible, in principle, to determine the parameters in the effective chiral lagrangian (including weak interaction couplings) by performing numerical simulations in the $\epsilon$--regime, i.e. at quark masses where the physical extent of the lattice is much smaller than the Compton wave length of the pion. The use of a formulation of the lattice theory that preserves chiral symmetry is attractive in this context, but the numerical implementation of any such approach requires special care in this kinematical situation due to the presence of some very low eigenvalues of the Dirac operator. We discuss a set of techniques (low-mode preconditioning and adapted-precision algorithms in particular) that make such computations numerically safe and more efficient by a large factor.


title: Uncovering Low-Dimensional Topological Structure in the QCD Vacuum

abstract: Recently, we have pointed out that sign-coherent 4-dimensional structures can not dominate topological charge fluctuations in QCD vacuum at all scales. Here we show that an enhanced lower-dimensional coherence is possible. In pure SU(3) lattice gauge theory we find that in a typical equilibrium configuration about 80% of space-time points are covered by two oppositely-charged connected structures built of elementary 3-dimensional coherent hypercubes. The hypercubes within the structure are connected through 2-dimensional common faces. We suggest that this coherence is a manifestation of a low-dimensional order present in the QCD vacuum. The use of a topological charge density associated with Ginsparg-Wilson fermions ("chiral smoothing") is crucial for observing this structure.


title: Correlators of left charges and weak operators in finite volume chiral perturbation theory

abstract: We compute the two-point correlator between left-handed flavour charges, and the three-point correlator between two left-handed charges and one strangeness violating \Delta I=3/2 weak operator, at next-to-leading order in finite volume SU(3)_L x SU(3)_R chiral perturbation theory, in the so-called epsilon-regime. Matching these results with the corresponding lattice measurements would in principle allow to extract the pion decay constant F, and the effective chiral theory parameter g_27, which determines the \Delta I = 3/2 amplitude of the weak decays K to \pi\pi as well as the kaon mixing parameter B_K in the chiral limit. We repeat the calculations in the replica formulation of quenched chiral perturbation theory, finding only mild modifications. In particular, a properly chosen ratio of the three-point and two-point functions is shown to be identical in the full and quenched theories at this order.


title: B_K from quenched overlap QCD

abstract: We present an exploratory calculation of the standard model Delta S=2 matrix element relevant for indirect CP violation in K -> pi pi decays. The computation is performed with overlap fermions in the quenched approximation at beta=6.0 on a 16^3x32 lattice. The resulting bare matrix element is renormalized non-perturbatively. Our preliminary result is B_K^{NDR}(2 GeV)=0.61(7), where the error does not yet include an estimate of systematic uncertainties.


title: Some remarks on cyclization and periodicity

abstract: Conditions when application of cyclization procedure does not lead to the periodic function are specified. Sufficient conditions under which obtained function is periodic are stated.


title: The measurement, fitting and interpretation of the radial distributions of Heavy-Light mesons calculated on a lattice with dynamical fermions

abstract: In our earlier work, the charge and matter radial distributions of heavy-light mesons were measured on a 16^3 times 24 lattice with a lattice spacing of a = 0.17 fm and a light quark mass about that of the strange quark. Several major improvements have now been made: 1) Dynamical fermions are used with a approximately 0.14 fm; 2) More gauge configurations are included (78 vs 20); 3) Off-axis, in addition to on-axis, insertions are made; 4) The data analysis is much more complete. In particular, distributions involving excited states are extracted.


title: Lattice calculation of the lowest order hadronic contribution to the muon anomalous magnetic moment

abstract: We present a quenched lattice calculation of the lowest order (alpha^2) hadronic contribution to the anomalous magnetic moment of the muon which arises from the hadronic vacuum polarization. A general method is presented for computing entirely in Euclidean space, obviating the need for the usual dispersive treatment which relies on experimental data for e^+e^- annihilation to hadrons. While the result is not yet of comparable accuracy to those state-of-the-art calculations, systematic improvement of the quenched lattice computation to this level of accuracy is straightforward and well within the reach of present computers. Including the effects of dynamical quarks is conceptually trivial, the computer resources required are not.


title: More on Electric and Magnetic Fluxes in SU(2)

abstract: The free energies of static charges and center monopoles are given by their fluxes. While electric fluxes show the universal behaviour of the deconfinement transition, the monopole free energies vanish in the thermodynamic limit at all temperatures and are thus irrelevant for the transition. Magnetic fluxes may, however, be used to measure the topological susceptibility without cooling.


title: Hadronic decay of a vector meson from the lattice

abstract: We explore the decay of a vector meson to two pseudoscalar mesons on the lattice with $N_f=2$ flavours of sea quark. Although we are working with quark masses that do not allow a physical decay, we show how the transition rate can be evaluated from the amplitude for $\rho \to \pi \pi$ and from the annihilation component of $\pi \pi \to \pi \pi$. We explore the decay amplitude for two different pion momenta and find consistent results. The coupling strength we find is in agreement with experiment. We also find evidence for a shift in the $\rho$ mass caused by mixing with two pion states.


title: Chiral fermion operators on the lattice

abstract: We only require generalized chiral symmetry and $\gamma_5$-hermiticity, which leads to a large class of Dirac operators describing massless fermions on the lattice, and use this framework to give an overview of developments in this field. Spectral representations turn out to be a powerful tool for obtaining detailed properties of the operators and a general construction of them. A basic unitary operator is seen to play a central r\^ole in this context. We discuss a number of special cases of the operators and elaborate on various aspects of index relations. We also show that our weaker conditions lead still properly to Weyl fermions and to chiral gauge theories.


title: Remarks on the reweighting method in the chemical potential direction

abstract: We comment on the reweighting method in the chemical potential $(\mu_{\rm q})$ direction. We study the fluctuation of the reweighting factor during Monte-Carlo steps. We find that it is the absolute value of the reweighting factor that mainly contributes to the shift of the phase transition line $(\beta_c)$ by the presence of $\mu_{\rm q}$. The phase fluctuation is a cause of the sign problem, but the effect on $\beta_c$ seems to be small. We also discuss $\beta_c$ for Iso-vector chemical potential and $\beta_c$ determined from simulations with imaginary chemical potential.


title: The profile of the broken string in the confined and deconfined phase in full QCD

abstract: We study the profile of the broken string (flux tube) in the maximally abelian gauge below and above the finite temperature phase transition in full QCD. In the deconfinement phase, the flux tube disappears and the electric field apperas to be Coulomb-like. In the confinement phase, but near $T_c$, at shorter distances a flux tube is formed like at zero temperature, while at larger distances the tube disappears similar to the deconfinement phase.


title: The flux distribution of the three quark system in SU(3)

abstract: We study the abelian color-flux distribution of the three quark system in the maximally abelian gauge on SU(3) lattices. The distribution of the color electric field suggests $Y Ansatz$, which might be interpreted through the dual superconductor picture as the result of the vacuum pressure in the confined phase. In order to clarify the flux structure, we investigate the color electric field in the three quark system also in the monopole part and in the photon part.


title: Lattice QCD study for the Quark-Gluon Mixed Condensate (anti-Q sigma G Q)

abstract: We study the quark-gluon mixed condensate g<\bar{q} sigma G q>, which is another chiral order parameter, using the SU(3)c lattice QCD with the Kogut-Susskind(KS) fermion at the quenched level. We generate 100 gauge configurations on the 16^4 lattice with \beta = 6.0, and measure the mixed condensate at 16 points in each gauge configuration for each current-quark mass of m_q=21, 36, 52 MeV. From the 1600 data for each m_q, we find the ratio between the mixed condensate and the quark condensate, m_0^2 = g<\bar{q} sigma G q> / <\bar{q}q> \simeq 2.5 GeV^2 at the lattice scale of 1/a \simeq 2 GeV in the chiral limit. The large value of the mixed condensate suggests its importance in the operator product expansions in QCD.


title: High Spin Glueballs from the Lattice

abstract: We discuss the principles underlying higher spin glueball calculations on the lattice. For that purpose, we develop numerical techniques to rotate Wilson loops by arbitrary angles in lattice gauge theories close to the continuum. As a first application, we compute the glueball spectrum of the SU(2) gauge theory in 2+1 dimensions for both parities and for spins ranging from 0 up to 4 inclusive. We measure glueball angular wave functions directly, decomposing them in Fourier modes and extrapolating the Fourier coefficients to the continuum. This allows a reliable labelling of the continuum states and gives insight into the way rotation symmetry is recovered. As one of our results, we demonstrate that the D=2+1 SU(2) glueball conventionally labelled as J^P = 0^- is in fact 4^- and that the lightest ``J=1'' state has, in fact, spin 3.


title: Vortices, Confinement and Higgs fields

abstract: We review lattice evidence for the vortex mechanism of quark confinement and study the influence of charged matter fields on the vortex distribution.


title: Some peculiarities of transition from discrete to continuum Fourier series in lattice theories

abstract: Transition from discrete to continuous Fourier series is studied for the functions becoming singular in the transition. Conditions are specified when summing replacement by integration is inadmissible.


title: New findings for topological excitations in SU(3) lattice gauge theory

abstract: We probe the SU(3) vacuum using eigenvectors of the Dirac operator with an arbitrary phase for the temporal boundary condition. We consider configurations with topological charge |Q| = 1 near the QCD phase transition and at low temperatures on a torus. For all our ensembles we show that the zero-mode of the Dirac operator changes its position as one changes the phase of the boundary condition. For ensembles near the QCD phase transition our results closely resemble the behavior of zero-modes for Kraan - van Baal solutions of the classical Yang-Mills equations where the individual lumps are interpreted as monopoles. Our findings near T_c and on the torus show that for both cases an excitation with topological charge |Q| = 1 is built from several separate lumps.


title: The BlueGene/L Supercomputer

abstract: The architecture of the BlueGene/L massively parallel supercomputer is described. Each computing node consists of a single compute ASIC plus 256 MB of external memory. The compute ASIC integrates two 700 MHz PowerPC 440 integer CPU cores, two 2.8 Gflops floating point units, 4 MB of embedded DRAM as cache, a memory controller for external memory, six 1.4 Gbit/s bi-directional ports for a 3-dimensional torus network connection, three 2.8 Gbit/s bi-directional ports for connecting to a global tree network and a Gigabit Ethernet for I/O. 65,536 of such nodes are connected into a 3-d torus with a geometry of 32x32x64. The total peak performance of the system is 360 Teraflops and the total amount of memory is 16 TeraBytes.


title: Convergence of chiral effective field theory

abstract: We formulate the expansion for the mass of the nucleon as a function of pion mass within chiral perturbation theory using a number of different ultra-violet regularisation schemes; including dimensional regularisation and various finite-ranged regulators. Leading and next-to-leading order non-analytic contributions are included through the standard one-loop Feynman graphs. In addition to the physical nucleon mass, the expansion is constrained by recent, extremely accurate, lattice QCD data obtained with two flavors of dynamical quarks. The extent to which different regulators can describe the chiral expansion is examined, while varying the range of quark mass over which the expansions are matched. Renormalised chiral expansion parameters are recovered from each regularisation prescription and compared. We find that the finite-range regulators produce consistent, model-independent results over a wide range of quark mass sufficient to solve the chiral extrapolation problem in lattice QCD.


title: Lattice effective theory and the phase transition at finite densities

abstract: The transition from the hadronic phase to the phase of color-superconductivity at large densities is addressed by an effective theory which incorporates the Yang-Mills dynamics in addition to the di-quark degree of freedom. A toy version of this theory is studied by lattice simulations. A first order phase transition separates the regime of broken color-electric flux tubes from the color superconducting phase. My findings suggest that the quark and gluon liberation occurs at the same critical chemical potential.


title: Massless QED_3 with explicit fermions

abstract: We study dynamical mass generation in QED in (2+1) dimensions using Hamiltonian lattice methods. We use staggered fermions, and perform simulations with explicit dynamical fermions in the chiral limit. We demonstrate that a recently developed method to reduce the fermion sign problem can successfully be applied to this problem. Our results are in agreement with both the strong coupling expansion and with Euclidean lattice simulations.


title: Direct simulations of small multi-fermion systems

abstract: I explore computer simulations of the dynamics of small multi-fermion lattice systems. The method is more general, but I concentrate on Hubbard type models where the fermions hop between a small number of connected sites. I use the natural mapping of fermion occupation numbers onto computer bits. Signs from fermion interchange are reduced to bit counting. The technique inherently requires computer resources growing exponentially with the system volume; so, it restricted to modestly small systems. Large volume results would require combining these techniques with further approximations, perhaps in a recursive renormalization group manner.


title: Non-perturbative Methods in Modal Field Theory

abstract: Several issues in the modal approach to quantum field theory are discussed. Within the formalism of spherical field theory, differential renormalization is presented and shown to result in a finite number of renormalization parameters. Computations of the massless Thirring model in 1+1 dimensions are presented using this approach. Diagonalization techniques in periodic field theory are demonstrated. Issues of very large Hilbert spaces are considered and several approaches are presented. The quasi sparse eigenvector (QSE) approach takes advantage of the relatively small number of basis states that typically contribute significantly to any particular eigenvector. Stochastic correction methods use Monte Carlo calculations to calculate higher order corrections to the quasi sparse result. The quasi sparse eigenvector method and stochastic error correction are applied to the Hubbard model. With U/t=4, the shift in the ground energy below the U=0 value is found to within 1% for the 8x8 Hubbard model with 25/64 filling.


title: Flux Tubes of Two- and Three-Quark System in Full QCD

abstract: We study the abelian color flux of two- and three-quark systems in the maximally abelian gauge in lattice QCD with dynamical fermions. We find that the abelian flux tube formed between quark and antiquark is very much the same as in quenched QCD up to quark separations of $R$$\sim$ 2fm. The profile of the color electric field in three-quark system suggests $Y$ ansatz, which might be interpreted as the result of the vacuum pressure in the confined phase. In order to clarify the flux structure, we investigate the color electric field of the three-quark system splittting the abelian gauge field into the monopole and photon parts.


title: Interactions between Lattice Hadrons

abstract: The effective residual interaction for a system of hadrons has a long tradition in theoretical physics. It has been mostly addressed in terms of boson exchange models. The aim of this review is to describe approaches based on lattice field theory and numerical simulation. At the present time this subject matter is in an exploratory stage. A large array of problems waits to be tackled, so that known features of hadron-hadron interactions will eventually be understood in a model-independent way. The lattice formulation, being capable of dealing with the nonperturbative regime, describes strong-interaction physics from first principles, i.e. quantum chromodynamics (QCD). Although the physics of hadron-hadron interactions may be intrinsically complicated, the methods used in lattice simulations are simple: For the most part they are based on standard mass calculations. This chapter addresses commonly used techniques, within QCD and also simpler lattice models, describes important results, and also gives some insight into numerical methods for multi-quark systems.


title: Lattice QCD, gauge fixing, and the transition to the perturbative regime

abstract: Perturbative QCD uses the Faddeev-Popov gauge-fixing procedure, which leads to ghosts and the local BRST invariance of the gauge-fixed perturbative QCD action. In the asymptotic regime, where perturbative QCD is relevant, Gribov copies can be neglected. In the nonperturbative regime, one must adopt either a nonlocal Gribov-copy free gauge (e.g., Laplacian gauge) or attempt to maintain local BRST invariance at the expense of admitting Gribov copies. These issues are explored and discussed. In addition, the relationship between recent Dyson-Schwinger based calculations of the infrared behavior of QCD Green's functions and the lattice calculation of these quantities is examined.


title: Light hadron spectroscopy with two flavors of $O(a)$-improved dynamical quarks

abstract: We present a high statistics study of the light hadron spectrum and quark masses in QCD with two flavors of dynamical quarks. Numerical simulations are carried out using the plaquette gauge action and the O(a)-improved Wilson quark action at \beta=5.2, where the lattice spacing is found to be a=0.0887(11)fm from \rho meson mass, on a 20^3\times 48 lattice. At each of five sea quark masses corresponding to m_{PS}/m_{V} \simeq 0.8-0.6, we generate 12000 trajectories using the symmetrically preconditioned Hybrid Monte Carlo algorithm. Finite spatial volume effects are investigated employing 12^3 \times 48, 16^3 \times 48 lattices. We also perform a set of simulations in quenched QCD with the same lattice actions at a similar lattice spacing to those for the full QCD runs. In the meson sector we find clear evidence of sea quark effects. The J parameter increases for lighter sea quark masses, and the full QCD meson masses are systematically closer to experiment than in quenched QCD. Careful finite-size studies are made to ascertain that these are not due to finite-size effects. Evidence of sea quark effects is less clear in the baryon sector due to larger finite-size effects. We also calculate light quark masses and find m_{ud}^{MS}(2GeV) =3.223(+0.046/-0.069)MeV and m_s^{MS}(2GeV)=84.5(+12.0/-1.7)MeV which are about 20% smaller than in quenched QCD.


title: The scaling region of the lattice O(N) sigma model at finite temperature

abstract: We present results from numerical studies of the finite temperature phase transition of the $(3+1)d$ O(N)-symmetric non-linear sigma model for $N=1,2$ and 3. We study the dependence of the width of the 3d critical region on $N$ and we show that the broken phase scaling region is much wider for N=2 and 3 than for N=1. We also compare the widths of the critical region in the low $T$ and high $T$ phases of the O(2) model and we show that the scaling region in the broken phase is much wider than in the symmetric phase. We also report results for the width of the scaling regions in the low $T$ phase$ (2+1)d$ Ising model and we show that the spatial correlation length has to be approximately twice the lattice temporal extent before the 2d scaling region is reached.


title: Opportunities, Challenges, and Fantasies in Lattice QCD

abstract: Some important problems in quantitative QCD will certainly yield to hard work and adequate investment of resources, others appear difficult but may be accessible, and still others will require essentially new ideas. Here I identify several examples in each class.


title: Interpolating the Free Energy Density Differences of Reweighting Methods

abstract: A discussion of the overlap problem of reweighting approaches to evaluating critical phenomenon in fermionic systems is motivated by highlighting the divergence of the joint probability density function of a general ratio. By identifying the bounds for which this integral can be expressed in closed form, we establish criteria for accurately mapping the joint ratio distribution of two disjoint ensembles through interpolation. The approach is applied to QCD with four staggered flavours to evaluate the critical line in the $\beta-\mu$ plane.


title: Heavy quark potential in lattice QCD at finite temperature

abstract: Results of the study of lattice QCD with two flavors of nonperturbatively improved Wilson fermions at finite temperature are presented. The transition temperature for m_{\pi}/m_{\rho} \sim 0.8 and lattice spacing a \sim 0.12 fm is determined. A two-exponent ansatz is successfully applied to describe the heavy quark potential in the confinement phase.


title: Finite temperature phase transition in lattice QCD with N_f=2 nonperturbatively improved Wilson fermions at N_t=8

abstract: The finite temperature lattice QCD with N_f=2 nonperturbatively improved Wilson fermions is studied on 16^3 8 lattice. Using abelian projection after fixing to MA gauge we determine the transition temperature for m_{\pi}/m_{\rho} \sim 0.8.


title: Dynamical generation of a gauge symmetry in the Double-Exchange model

abstract: It is shown that a bosonic formulation of the double-exchange model, one of the classical models for magnetism, generates dynamically a gauge-invariant phase in a finite region of the phase diagram. We use analytical methods, Monte Carlo simulations and Finite-Size Scaling analysis. We study the transition line between that region and the paramagnetic phase. The numerical results show that this transition line belongs to the Universality Class of the Antiferromagnetic RP(2) model. The fact that one can define a Universality Class for the Antiferromagnetic RP(2) model, different from the one of the O(N) models, is puzzling and somehow contradicts naive expectations about Universality.


title: Quark Mass Dependence of Nucleon Properties and Extrapolation from Lattice QCD

abstract: We summarize developments concerning the quark mass dependence of nucleon magnetic moments and the axial-vector coupling constant g_A. The aim is to explore the feasibility of chiral effective field theory methods for the extrapolation of lattice QCD results, from the relatively large quark masses that can be handled in such computations down to the physically relevant range.


title: The electrical conductivity and soft photon emissivity of the QCD plasma

abstract: The electrical conductivity in the hot phase of the QCD plasma is extracted from a quenched lattice measurement of the Euclidean time vector correlator for 1.5 < T/T_c < 3. The spectral density in the vicinity of the origin is examined using a method specially adapted to this region, and a peak at small energies is seen. The continuum limit of the electrical conductivity, and the closely related soft photon emissivity of the QCD plasma, are then extracted from a fit to the Fourier transform of the temporal vector correlator.


title: The Gribov problem and its solution from a toy model point of view

abstract: The standard Faddeev Popov gauge fixing procedure is put on solid grounds by taking into account a topological factor which corrects for the number of Gribov copies within the first Gribov horizon. Zwanziger's stochastic approach to gauge fixed Yang-Mills theory is briefly reviewed. A simple toy model is presented which illustrates both methods. Within the toy model, I show that a stochastic drift force can be constructed with which the gauged configurations are attracted by the fundamental modular region. The toy model shows that an action which gives rise to the drift force can be found. This makes a ``heat bath'' simulation possible, which is seen to be superior to the Langevin approach at the numerical level.


title: Evidence for Dyonic Structure of SU(2) Lattice Gauge Fields Below $T_c$

abstract: We report the observation of calorons with nontrivial holonomy and fractional topological charge objects in cooled lattice samples derived from SU(2) equilibrium ensembles at $T<T_c$.


title: In Medium Hadron Properties

abstract: We discuss modifications of hadron properties in the heat bath extracted from Euclidian correlation functions and spectral functions r econstructed with the Maximum Entropy Method. To investigate the cut-off dependencies we perform simulation on various lattice sizes and present the result of an analytic calculation of the mesonic spectral functions in the infinite temperature limit.


title: Confinement, deconfinement and the photon propagator in 3D cQED on the lattice

abstract: We report on a lattice study of the gauge boson propagator of 3D compact QED in Landau gauge at zero and finite temperature. Non-perturbative effects are reflected by the generation of a mass m and by an anomalous dimension alpha. These effects can be attributed to monopoles and are absent in the propagator of the regular part of the gauge field.


title: A hidden symmetry in the Standard Model

abstract: We found an additional symmetry hidden in the fermion and Higgs sectors of the Standard Model. It is connected to the centers of the SU(3) and SU(2) subgroups of the gauge group. A lattice regularization of the whole Standard Model is constructed that possesses this symmetry.


title: The Use of the Monte Carlo Hamiltonian

abstract: In order to solve quantum field theory in a non-perturbative way, Lagrangian lattice simulations have been very successful. Here we discuss a recently proposed alternative Hamiltonian lattice formulation - the Monte Carlo Hamiltonian. In order to show its working in the case of the scalar $\Phi^{4}_{1+1}$ model, we have computed thermodynamic functions like free energy, average energy, entropy and specific heat. We find good agreement between the results from the Monte Carlo Hamiltonian and standard Lagrangian lattice computations. However, the Monte Carlo Hamiltonian results show less fluctuations under variation of temperature. We address properties of the MC Hamiltonian, like a finite temperature window, and scaling properties. Also we discuss possible future applications - like quantum chaos in many-body systems, the non-perturbative computation of wave functions of elementary particles, as well as scattering amplitudes in high energy physics.


title: Numerical simulation with light Wilson-quarks

abstract: The computational cost of numerical simulations of QCD with light dynamical Wilson-quarks is estimated by determining the autocorrelation of various quantities. In test runs the expected qualitative behaviour of the pion mass and coupling at small quark masses is observed.


title: On the renormalization of the Polyakov loop

abstract: We discuss a non-perturbative renormalization of n-point Polyakov loop correlation functions by explicitly introducing a renormalization constant for the Polyakov loop operator on a lattice deduced from the short distance properties of 2-point correlators. We calculate this constant for the SU(3)gauge theory.


title: Renormalized quark-anti-quark free energy

abstract: We present results on the renormalized quark-anti-quark free energy in SU(3) gauge theory at finite temperatures. We discuss results for the singlet, octet and colour averaged free energies and comment on thermal relations which allow to extract separately the potential energy and entropy from the free energy.


title: Critical behaviour and Scaling functions for the three-dimensional O(6) spin model with external field

abstract: We numerically investigate the three-dimensional O(6) model on 12^3 to 120^3 lattices. From Binder's cumulant at vanishing magnetic field we obtain the critical coupling J_c=1.42865(5) and verify this value with the \chi^2-method. The universal value of Binder's cumulant at this point is g_r(J_c)=-1.94456(10). At the critical coupling we find the critical exponents \nu=0.818(5), \beta=0.425(2) and \gamma=1.604(6) from a finite size scaling analysis. We also determine the finite-size-scaling function on the critical line and the equation of state. Our O(6)-result for the equation of state is compared to the Ising, O(2) and O(4) results.


title: Matrix model correlation functions and lattice data for the QCD Dirac operator with chemical potential

abstract: We apply a complex chiral random matrix model as an effective model to QCD with a small chemical potential at zero temperature. In our model the correlation functions of complex eigenvalues can be determined analytically in two different limits, at weak and strong non-Hermiticity. We compare them to the distribution of the smallest Dirac operator eigenvalues from quenched QCD lattice data for small values of the chemical potential, appropriately rescaled with the volume. This confirms the existence of two different scaling regimes from lattice data.


title: Scaling behavior of the overlap quark propagator in Landau gauge

abstract: The properties of the momentum space quark propagator in Landau gauge are examined for the overlap quark action in quenched lattice QCD. Numerical calculations are done on three lattices with different lattice spacings and similar physical volumes to explore the approach of the quark propagator toward the continuum limit. We have calculated the nonperturbative momentum-dependent wave function renormalization function Z(p) and the nonperturbative mass function M(p) for a variety of bare quark masses and perform an extrapolation to the chiral limit. We find the behavior of Z(p) and M(p) are in reasonable agreement between the two finer lattices in the chiral limit, however the data suggest that an even finer lattice is desirable. The large momentum behavior is examined to determine the quark condensate.


title: Critical Slowing-Down in SU(2) Landau-Gauge-Fixing Algorithms at beta = infinity

abstract: We evaluate numerically and analytically the dynamic critical exponent $z$ for five gauge-fixing algorithms in SU(2) lattice Landau-gauge theory by considering the case $\beta = \infty$. Numerical data are obtained in two, three and four dimensions. Results are in agreement with those obtained previously at finite $\beta$ in two dimensions. The theoretical analysis, valid for any dimension $d$, helps us clarify the tuning of these algorithms. We also study generalizations of the overrelaxation algorithm and of the stochastic overrelaxation algorithm and verify that we cannot have a dynamic critical exponent $z$ smaller than 1 with these local algorithms. Finally, the analytic approach is applied to the so-called $\lambda$-gauges, again at $\beta = \infty$, and verified numerically for the two-dimensional case.


title: Coupling of the light vector meson to the vector and to the tensor current

abstract: We present results for the coupling of the light vector mesons to the tensor current, relative to the standard vector meson decay constants. From an O(a)-improved lattice study, performed at three values of the lattice spacing in the quenched approximation, our final values (in the continuum limit), in the MS-scheme at 2 GeV, are: (f^T/f)_rho =0.72(2)(+2/-0), (f^T/f)_K* =0.74(2), (f^T/f)_phi =0.76(1).


title: Chiral extrapolations of nucleon properties from lattice QCD

abstract: We report on recent work about the study of quark mass dependence of nucleon magnetic moments and axial-vector coupling constant. We examine the feasibility of chiral effective field theory methods for the extrapolation of lattice QCD data obtained at relative large pion masses down to the physical values.


title: A lattice study of $\Lambda_b$ semileptonic decay

abstract: We present results from a lattice study of the semileptonic decay $\Lambda_b \to \Lambda_c l \nu_l $. We use $O(a^2, \alpha_s a^2)$ improved quenched lattices of the MILC collaboration, with lattice spacing $\sim 0.13$ fm. For the valence quarks, the tadpole-improved clover action is used, with the Fermilab method employed for the heavy quarks. Form factors are extracted from the vector as well as the axial-vector part of the current.


title: The Confinement Problem in Lattice Gauge Theory

abstract: I review investigations of the quark confinement mechanism that have been carried out in the framework of SU(N) lattice gauge theory. The special role of Z(N) center symmetry is emphasized.


title: Lattice calculation of $1^{-+}$ hybrid mesons with improved Kogut-Susskind fermions

abstract: We report on a lattice determination of the mass of the exotic $1^{-+}$ hybrid meson using an improved Kogut-Susskind action. Results from both quenched and dynamical quark simulations are presented. We also compare with earlier results using Wilson quarks at heavier quark masses. The results on lattices with three flavors of dynamical quarks show effects of sea quarks on the hybrid propagators which probably result from coupling to two meson states. We extrapolate the quenched results to the physical light quark mass to allow comparison with experimental candidates for the $1^{-+}$ hybrid meson. The lattice result remains somewhat heavier than the experimental result, although it may be consistent with the $\pi_1(1600)$.


title: Supersymetry on the Noncommutative Lattice

abstract: Built upon the proposal of Kaplan et.al. [hep-lat/0206109], we construct noncommutative lattice gauge theory with manifest supersymmetry. We show that such theory is naturally implementable via orbifold conditions generalizing those used by Kaplan {\sl et.al.} We present the prescription in detail and illustrate it for noncommutative gauge theories latticized partially in two dimensions. We point out a deformation freedom in the defining theory by a complex-parameter, reminiscent of discrete torsion in string theory. We show that, in the continuum limit, the supersymmetry is enhanced only at a particular value of the deformation parameter, determined solely by the size of the noncommutativity.


title: Lattice chiral symmetry, CP-violation and Majorana fermions

abstract: A brief summary of lattice fermions defined by the general Ginsparg-Wilson algebra is first given. It is then shown that those general class of fermion operators have a conflict with CP invariance in chiral gauge theory and with the definition of Majorana fermions in the presence of chiral-symmetric Yukawa couplings. The same conclusion holds for the domain-wall fermion also.


title: Lattice QCD at non-vanishing density: phase diagram, equation of state

abstract: We propose a method to study lattice QCD at non-vanishing temperature (T) and chemical potential (\mu). We use n_f=2+1 dynamical staggered quarks with semi-realistic masses on L_t=4 lattices. The critical endpoint (E) of QCD on the Re(\mu)-T plane is located. We calculate the pressure (p), the energy density (\epsilon) and the baryon density (n_B) of QCD at non-vanishing T and \mu.


title: Lattice Supersymmetry and Topological Field Theory

abstract: It is known that certain theories with extended supersymmetry can be discretized in such a way as to preserve an exact fermionic symmetry. In the simplest model of this kind, we show that this residual supersymmetric invariance is actually a BRST symmetry associated with gauge fixing an underlying local shift symmetry. Furthermore, the starting lattice action is then seen to be entirely a gauge fixing term. The corresponding continuum theory is known to be a topological field theory. We look, in detail, at one example - supersymmetric quantum mechanics which possesses two such BRST symmetries. In this case, we show that the lattice theory can be obtained by blocking out of the continuum in a carefully chosen background metric. Such a procedure will not change the Ward identities corresponding to the BRST symmetries since they correspond to topological observables. Thus, at the quantum level, the continuum BRST symmetry is preserved in the lattice theory. Similar conclusions are reached for the two-dimensional complex Wess-Zumino model and imply that all the supersymmetric Ward identities are satisfied {\it exactly} on the lattice. Numerical results supporting these conclusions are presented.


title: Renormalization group improved action on anisotropic lattices

abstract: We study a block spin transformation in the SU(3) lattice gauge theory on anisotropic lattices to obtain Iwasaki's renormalization group improved action for anisotropic cases. For the class of actions with plaquette and $1\times2$ rectangular terms, we determine the improvement parameters as functions of the anisotropy $\xi= a_s/a_t$. We find that the program of improvement works well also on anisotropic lattices. From a study of an indicator which estimates the distance to the renormalized trajectory, we show that, for the range of the anisotropy $\xi \approx 1$--4, the coupling parameters previously determined for isotropic lattices improve the theory considerably.


title: Exotics

abstract: We review lattice QCD results for glueballs (including a discussion of mixing with scalar mesons), hybrid mesons and other exotic states (such as $B_s B_s$ molecules).


title: Abelian dominance and gluon propagators in the Maximally Abelian gauge of SU(2) lattice gauge theory

abstract: Propagators of the diagonal and the off-diagonal gluons are studied numerically in the Maximal Abelian gauge of SU(2) lattice gauge theory. It is found that in the infrared region the propagator of the diagonal gluon is strongly enhanced in comparison with the off--diagonal one. The enhancement factor is about 50 at our smallest momentum 325 MeV. We have also applied various fits to the propagator formfactors.


title: CP symmetry and lattice chiral gauge theories

abstract: The CP symmetry is a fundamental discrete symmetry in chiral gauge theory. Therefore this symmetry is expected to be kept also on the lattice. However, it has been pointed out by Hasenfratz that the chiral fermion action in Luscher's formulation of lattice chiral gauge theory is not invariant under CP transformation. In this thesis, we first review the method of constructing chiral gauge theory on the lattice. Then we generalize the analysis of Hasenfratz and show that CP symmetry is not manifestly implemented for the local and doubler-free Ginsparg-Wilson operator under rather general assumptions for chiral projection operators. We next calculate the fermion generating functional and precisely identify where the effects of this CP breaking appear in this formulation. The basics of lattice gauge theory are briefly summarized in Appendix.A.


title: The spectrum of screening masses near T_c: predictions from universality

abstract: We discuss the spectrum of screening masses in a pure gauge theory near the deconfinement temperature from the point of view of the dimensionally reduced model describing the spontaneous breaking of the center symmetry. Universality arguments can be used to predict the values of the mass ratios in the scaling region of the deconfined phase when the transition is of second order. One such prediction is that the scalar sector of the screening spectrum in SU(2) pure gauge theory contains a bound state of the fundamental excitation, corresponding through universality to the bound state found in the 3D Ising model and phi^4 theory in the broken symmetry phase. A Monte Carlo evaluation of the screening masses in the gauge theory confirms the validity of the prediction. We briefly discuss the possibility of using similar arguments for first order deconfinement transitions, and in particular for the physically relevant case of SU(3).


title: QCD at Fixed Topology

abstract: Since present Monte Carlo algorithms for lattice QCD may become trapped in a fixed topological charge sector, it is important to understand the effect of calculating at fixed topology. In this work, we show that although the restriction to a fixed topological sector becomes irrelevant in the infinite volume limit, it gives rise to characteristic finite size effects due to contributions from all $\theta$-vacua. We calculate these effects and show how to extract physical results from numerical data obtained at fixed topology.


title: Duality of gauge field singularities and the structure of the flux tube in Abelian-projected SU(2) gauge theory and the dual Abelian Higgs model

abstract: The structure of the flux-tube profile in Abelian-projected (AP) SU(2) gauge theory in the maximally Abelian gauge is studied. The connection between the AP flux tube and the classical flux-tube solution of the U(1) dual Abelian Higgs (DAH) model is clarified in terms of the path-integral duality transformation. This connection suggests that the electric photon and the magnetic monopole parts of the Abelian Wilson loop can act as separate sources creating the Coulombic and the solenoidal electric field inside a flux tube. The conjecture is confirmed by a lattice simulation which shows that the AP flux tube is composed of these two contributions.


title: Field theoretical approach to the study of theta dependence in Yang-Mills theories on the lattice

abstract: We discuss the extension of the field theoretical approach, already used in the lattice determination of the topological susceptibility, to the computation of further terms in the expansion of the ground state energy $F(\theta)$ around $\theta = 0$ in SU(N) Yang-Mills theories. In particular we determine the fourth order term in the expansion for SU(3) pure gauge theory and compare our results with previous cooling determinations. In the last part of the paper we make some considerations about the nature of the ultraviolet fluctuations responsible for the renormalization of the lattice topological charge correlation functions; in particular we propose and test an ansatz which leads to improved estimates of the fourth and higher order terms in the expansion of F(\theta).


title: Delta Baryon Magnetic Moments From Lattice QCD

abstract: Theoretical predictions for the magnetic moments of the physical Delta baryons are extracted from lattice QCD calculations. We utilize finite-range regulated effective field theory that is constructed to have the correct Dirac moment mass dependence in the region where the up and down quark masses are heavy. Of particular interest is the chiral nonanalytic behaviour encountered as the nucleon-pion decay channel opens. We find a Delta^++ magnetic moment (at the Delta pole) of 4.99 \pm 0.56 \mu_N. This result is within the Particle Data Group range of 3.7-7.5 \mu_N and compares well with the experimental result of Bosshard et al. of 4.52 \pm 0.51 \pm 0.45 \mu_N. The interplay between the different pion-loop contributions to the Delta^+ magnetic moment leads to the surprising result that the proton moment may exceed that of the Delta^+, contrary to conventional expectations.


title: Low-Dimensional Long-Range Topological Charge Structure in the QCD Vacuum

abstract: While sign-coherent 4-dimensional structures cannot dominate topological charge fluctuations in the QCD vacuum at all scales due to reflection positivity, it is possible that enhanced coherence exists over extended space-time regions of lower dimension. Using the overlap Dirac operator to calculate topological charge density, we present evidence for such structure in pure-glue SU(3) lattice gauge theory. It is found that a typical equilibrium configuration is dominated by two oppositely-charged sign-coherent connected structures (``sheets'') covering about 80% of space-time. Each sheet is built from elementary 3-d cubes connected through 2-d faces, and approximates a low-dimensional curved manifold (or possibly a fractal structure) embedded in the 4-d space. At the heart of the sheet is a ``skeleton'' formed by about 18% of the most intense space-time points organized into a global long-range structure, involving connected parts spreading over maximal possible distances. We find that the skeleton is locally 1-dimensional and propose that its geometrical properties might be relevant for understanding the possible role of topological charge fluctuations in the physics of chiral symmetry breaking.


title: Critical parameters and universal amplitude ratios of two-dimensional spin-S Ising models using high- and low-temperature expansions

abstract: For the study of Ising models of general spin S on the square lattice, we have combined our recently extended high-temperature expansions with the low-temperature expansions derived some time ago by Enting, Guttmann and Jensen. We have computed for the first time various critical parameters and improved the estimates of others. Moreover the properties of hyperscaling and of universality (spin S independence) of exponents and of various dimensionless amplitude combinations have been verified accurately. Assuming the validity of the lattice-lattice scaling, from our estimates of critical amplitudes for the square lattice we have also obtained estimates of the corresponding amplitudes for the spin S Ising model on the triangular, honeycomb, and kagome` lattices.


title: Partially quenched chiral perturbation theory and numerical simulations

abstract: The dependence of the pseudoscalar meson mass and decay constant is compared to one-loop Partially Quenched Chiral Perturbation Theory (PQChPT) in a numerical simulation with two light dynamical quarks. The characteristic behaviour with chiral logarithms is observed. The values of the fitted PQChPT-parameters are in a range close to the expectation in continuum in spite of the fact that the lattice spacing is still large, namely a=0.28 fm.


title: Anisotropic lattice actions for heavy quark

abstract: The anisotropic lattice fermion actions are investigated with the one-loop perturbative calculations aiming at constructing a formulation for heavy quark with controlled systematic uncertainties. For the heavy-light systems at rest the anisotropic lattice with small temporal lattice spacing $a_t$ suppresses the discretization error by a power of $a_t m_Q$ for a heavy quark of mass $m_Q$. We discuss the issue of large discretization errors, which scales as $a_s m_Q$ with $a_s$ the spatial lattice spacing. By performing one-loop calculations of the speed-of-light renormalization for several possible lattice actions in the limit of $a_t\to 0$, we show that one can eliminate the large systematic error on the anisotropic lattice.


title: Negative-parity baryon spectra in quenched anisotropic lattice QCD

abstract: We investigate the negative-parity baryon spectra in quenched lattice QCD. We employ the anisotropic lattice with standard Wilson gauge and O(a) improved Wilson quark actions at three values of lattice spacings with renormalized anisotropy xi=a_sigma/a_tau=4, where a_sigma and a_tau are spatial and temporal lattice spacings, respectively. The negative-parity baryons are measured with the parity projection. In particular, we pay much attention to the lowest SU(3) flavor-singlet negative-parity baryon, which is assigned as the Lambda(1405) in the quark model. For the flavor octet and decuplet negative-parity baryons, the calculated masses are close to the experimental values of corresponding lowest-lying negative-parity baryons. In contrast, the flavor-singlet baryon is found to be about 1.7 GeV, which is much heavier than the Lambda(1405). Therefore, it is difficult to identify the Lambda(1405) to be the flavor-singlet three-quark state, which seems to support an interesting picture of the penta-quark (uds qbar q) state or the N-Kbar molecule for the Lambda(1405).


title: The self-energy of improved staggered quarks

abstract: We calculate the fermion self-energy at O(alpha_s) for the Symanzik improved staggered fermion action used in recent lattice simulations by the MILC collaboration. We demonstrate that the algebraic approach to lattice perturbation theory, suggested by us recently, is a powerful tool also for improved lattice actions.


title: A model for string-breaking in QCD

abstract: We present a model for string breaking based on the existence of chromoelectric flux tubes. We predict the form of the long-range potential, and obtain an estimate of the string breaking length. A prediction is also obtained for the behaviour with temperature of the string breaking length near the deconfinement phase transition. We plan to use this model as a guide for a program of study of string breaking on the lattice.


title: A precise determination of the decay constant of the Ds-meson in quenched QCD

abstract: We present a lattice determination of the leptonic decay constant F_Ds of the Ds-meson using its mass and F_K as experimental input. Setting the scale with F_K=160 MeV, our final result is F_Ds = 252(9) MeV. The error contains all uncertainties apart from the quenched approximation. Setting the scale with the nucleon mass instead leads to an increase of 3(1) MeV of the leptonic decay constant which might indicate a small quenching error.


title: Supersymmetry on a Euclidean Spacetime Lattice I: A Target Theory with Four Supercharges

abstract: We formulate a Euclidean spacetime lattice whose continuum limit is (2,2) supersymmetric Yang-Mills theory in two dimensions, a theory which possesses four supercharges and an anomalous global chiral symmetry. The lattice action respects one exact supersymmetry, which allows the target theory to emerge in the continuum limit without fine-tuning. Our method exploits an orbifold construction described previously for spatial lattices in Minkowski space, and can be generalized to more complicated theories with additional supersymmetry and more spacetime dimensions.


title: Coulomb Energy, Vortices, and Confinement

abstract: We estimate the Coulomb energy of static quarks from a Monte Carlo calculation of the correlator of timelike link variables in Coulomb gauge. We find, in agreement with Cucchieri and Zwanziger, that this energy grows linearly with distance at large quark separations. The corresponding string tension, however, is several times greater than the accepted asymptotic string tension, indicating that a state containing only static sources, with no constituent gluons, is not the lowest energy flux tube state. The Coulomb energy is also measured on thermalized lattices with center vortices removed by the de Forcrand-D'Elia procedure. We find that when vortices are removed, the Coulomb string tension vanishes.


title: Non-perturbative renormalization of the static axial current in quenched QCD

abstract: We non-perturbatively calculate the scale dependence of the static axial current in the Schroedinger functional scheme by means of a recursive finite-size scaling technique, taking the continuum limit in each step. The bare current in the O(a) improved theory as well as in the original Wilson regularization is thus connected to the renormalization group invariant one. The latter may then be related to the current at the B-scale defined such that its matrix elements differ from the physical (QCD) ones by O(1/M). At present, a (probably small) perturbative uncertainty enters in this step. As an application, we renormalize existing unimproved data on F_B^{bare} and extrapolate to the continuum limit. We also study an interesting function h(d/L,u) derived from the Schroedinger functional amplitude describing the propagation of a static quark-antiquark pair.


title: Physical Nucleon Properties from Lattice QCD

abstract: We demonstrate that the extremely accurate lattice QCD data for the mass of the nucleon recently obtained by CP-PACS, combined with modern chiral extrapolation techniques, leads to a value for the mass of the physical nucleon which has a systematic error of less than one percent.


title: Fermi Surface Phenomena in the (2+1)d Four-Fermi Model

abstract: We study the Gross-Neveu model in 2+1 dimensions with a baryon chemical potential mu using both analytical and numerical methods. For mu greater than a critical value the model is chirally symmetric and has a Fermi surface with Fermi momentum ~mu. We have calculated the particle interaction in medium due to scalar meson exchange to leading order in 1/N, where N is the number of flavors, in the hard dense loop approach. The result has been used to calculate the relation between mu and the Fermi momentum and velocity in the resulting Fermi liquid to O(1/N). Simulation results from a 32^2x48 lattice for fermion and meson dispersion relations and meson wavefunctions are then presented, showing qualitative and in some cases quantitative agreement with analytic predictions. In particular, the simulations show clear evidence for the in-medium modification of the scalar propagator, oscillatory behaviour in the wavefunction consistent with a sharp Fermi surface, and tentative evidence for a massless pole in the vector meson channel resembling zero sound.


title: SU(2) Landau gluon propagator on a 140^3 lattice

abstract: We present a numerical study of the gluon propagator in lattice Landau gauge for three-dimensional pure-SU(2) lattice gauge theory at couplings beta = 4.2, 5.0, 6.0 and for lattice volumes V = 40^3, 80^3, 140^3. In the limit of large V we observe a decreasing gluon propagator for momenta smaller than p_{dec} = 350^{+ 100}_{- 50} MeV. Data are well fitted by Gribov-like formulae and seem to indicate an infra-red critical exponent kappa slightly above 0.6, in agreement with recent analytic results.


title: Exceptional Confinement in G(2) Gauge Theory

abstract: We study theories with the exceptional gauge group G(2). The 14 adjoint "gluons" of a G(2) gauge theory transform as {3}, {3bar} and {8} under the subgroup SU(3), and hence have the color quantum numbers of ordinary quarks, anti-quarks and gluons in QCD. Since G(2) has a trivial center, a "quark" in the {7} representation of G(2) can be screened by "gluons". As a result, in G(2) Yang-Mills theory the string between a pair of static "quarks" can break. In G(2) QCD there is a hybrid consisting of one "quark" and three "gluons". In supersymmetric G(2) Yang-Mills theory with a {14} Majorana "gluino" the chiral symmetry is Z(4)_\chi. Chiral symmetry breaking gives rise to distinct confined phases separated by confined-confined domain walls. A scalar Higgs field in the {7} representation breaks G(2) to SU(3) and allows us to interpolate between theories with exceptional and ordinary confinement. We also present strong coupling lattice calculations that reveal basic features of G(2) confinement. Just as in QCD, where dynamical quarks break the Z(3) symmetry explicitly, G(2) gauge theories confine even without a center. However, there is not necessarily a deconfinement phase transition at finite temperature.


title: Two flavors of dynamical quarks on anisotropic lattices

abstract: We report on our study of two-flavor full QCD on anisotropic lattices using $O(a)$-improved Wilson quarks coupled with an RG-improved glue. The bare gauge and quark anisotropies corresponding to the renormalized anisotropy $\xi=a_s/a_t = 2$ are determined as functions of $\beta$ and $\kappa$, which covers the region of spatial lattice spacings $a_s\approx 0.28$--0.16 fm and $m_{PS}/m_V\approx 0.6$--0.9. The calibrations of the bare anisotropies are performed with the Wilson loop and the meson dispersion relation at 4 lattice cutoffs and 5--6 quark masses. Using the calibration results we calculate the meson mass spectrum and the Sommer scale $r_0$. We confirm that the values of $r_0$ calculated for the calibration using pseudo scalar and vector meson energy momentum dispersion relation coincide in the continuum limit within errors. This work serves to lay ground toward studies of heavy quark systems and thermodynamics of QCD including the extraction of the equation of state in the continuum limit using Wilson-type quark actions.


title: A recurrence scheme for least-square optimized polynomials

abstract: A recurrence scheme is defined for the numerical determination of high degree polynomial approximations to functions as, for instance, inverse powers near zero. As an example, polynomials needed in the two-step multi-boson (TSMB) algorithm for fermion simulations are considered. For the polynomials needed in TSMB a code in C is provided which is easily applicable to polynomial degrees of several thousands.


title: Dynamical fermions as a global correction

abstract: In the simplified setting of the Schwinger model we present a systematic study on the simulation of dynamical fermions by global accept/reject steps that take into account the fermion determinant. A family of exact algorithms is developed, which combine stochastic estimates of the determinant ratio with the exploitation of some exact extremal eigenvalues of the generalized problem defined by the `old' and the `new' Dirac operator. In this way an acceptable acceptance rate is achieved with large proposed steps and over a wide range of couplings and masses.


title: Quark mass dependence of the nucleon axial-vector coupling constant

abstract: We study the quark mass expansion of the axial-vector coupling constant g_A of the nucleon. The aim is to explore the feasibility of chiral effective field theory methods for extrapolation of lattice QCD results - so far determined at relatively large quark masses corresponding to pion masses larger than 0.6 GeV - down to the physical value of the pion mass. We compare two versions of non-relativistic chiral effective field theory: One scheme restricted to pion and nucleon degrees of freedom only, and an alternative approach which incorporates explicit Delta(1230) resonance degrees of freedom. It turns out that, in order to approach the physical value of g_A in a leading-one-loop calculation, the inclusion of the explicit Delta(1230) degrees of freedom is crucial. With information on important higher order couplings constrained from analyses of inelastic pion production processes, a chiral extrapolation function for g_A is obtained, which works well from the chiral limit across the physical point into the region of present lattice data. The resulting enhancement of our extrapolation function near the physical pion mass is found to arise from an interplay between long- and short- distance physics.


title: Chiral Limit of Strongly Coupled Lattice Gauge Theories

abstract: We construct a new and efficient cluster algorithm for updating strongly coupled U(N) lattice gauge theories with staggered fermions in the chiral limit. The algorithm uses the constrained monomer-dimer representation of the theory and should also be of interest to researchers working on other models with similar constraints. Using the new algorithm we address questions related to the chiral limit of strongly coupled U(N) gauge theories beyond the mean field approximation. We show that the infinite volume chiral condensate is non-zero in three and four dimensions. However, on a square lattice of size $L$ we find $\sum_x < \bar\psi\psi(x) \bar\psi\psi(0) > \sim L^{2-\eta}$ for large $L$ where $\eta = 0.420(3)/N + 0.078(4)/N^2$. These results differ from an earlier conclusion obtained using a different algorithm. Here we argue that the earlier calculations were misleading due to uncontrolled autocorrelation times encountered by the previous algorithm.


title: Spectroscopy using the Anisotropic Clover Action

abstract: The calculation of the light-hadron spectrum in the quenched approximation to QCD using an anisotropic clover fermion action is presented. The tuning of the parameters of the action is discussed, using the pion and rho dispersion relation. The adoption of an anisotropic lattice provides clear advantages in the determination of the baryonic resonances, and in particular that of the so-called Roper resonance, the lightest radial excitation of the nucleon.


title: Gauge Theories on a 2+2 Anisotropic Lattice

abstract: The implementation of gauge theories on a four-dimensional anisotropic lattice with two distinct lattice spacings is discussed, with special attention to the case where two axes are finely and two axes are coarsely discretized. Feynman rules for the Wilson gauge action are derived and the renormalizability of the theory and the recovery of the continuum limit are analyzed. The calculation of the gluon propagator and the restoration of Lorentz invariance in on-shell states is presented to one-loop order in lattice perturbation theory for $SU(N_c)$ on both 2+2 and 3+1 lattices.


title: The pion form factor in improved lattice QCD

abstract: We calculate the electromagnetic form factor of the pion in lattice gauge theory. The non-perturbatively improved Sheikoleslami-Wohlert lattice action is used together with the O(a) improved current. The form factor is compared to results for other choices for the current and features of the structure of the pion deduced from the 'Bethe-Salpeter wave function' are discussed.


title: An Application of Feynman-Kleinert Approximants to the Massive Schwinger Model on a Lattice

abstract: A trial application of the method of Feynman-Kleinert approximants is made to perturbation series arising in connection with the lattice Schwinger model. In extrapolating the lattice strong-coupling series to the weak-coupling continuum limit, the approximants do not converge well. In interpolating between the continuum perturbation series at large fermion mass and small fermion mass, however, the approximants do give good results. In the course of the calculations, we picked up and rectified an error in an earlier derivation of the continuum series coefficients.


title: Aspects of Domain-Wall Fermion on the Lattice

abstract: Two transparent layers are introduced at the boundaries of the fifth dimension for the optimal domain-wall fermions. For the quark fields defined in terms of these two transparent layers, they obey the usual chiral projection rule in the continuum, independent of the gauge fields. Consequently, any observable constructed with the quark fields manifests the symmetries exactly as those of its counterpart in the continuum.


title: Lattice chiral fermions from continuum defects

abstract: We consider whether defects of co-dimension two could produce new lattice chiral fermions.


title: Status of B_K from Lattice QCD

abstract: A brief review of lattice calculations of the bag parameter B_K relevant for understanding indirect CP violation in the neutral kaon sector is given. A status report on current state-of-the-art calculations is presented as well as a discussion of the value of B_K exported to phenomenologists. This is a condensed and updated version of the review presented at the CKM Unitarity Triangle Workshop held at CERN during Feburary 13-16, 2002 (http://ckm-workshop.web.cern.ch/ckm-workshop/).


title: Renormalization of Anisotropy and Glueball Masses on Tadpole Improved Lattice Gauge Action

abstract: The Numerical calculations for tadpole-improved U(1) lattice gauge theory in three-dimensions on anisotropic lattices have been performed using standard path integral Monte Carlo techniques. Using average plaquette tadpole renormalization scheme, simulations were done with temporal lattice spacings much smaller than the spatial ones and results were obtained for the string tension, the renormalized anisotropy and scalar glueball masses. We find, by comparing the `regular' and `sideways' potentials, that tadpole improvement results in very little renormalization of the bare anisotropy and reduces the discretization errors in the static quark potential and in the glueball masses.


title: Continuous external momenta in non-perturbative lattice simulations: a computation of renormalization factors

abstract: We discuss the usage of continuous external momenta for computing renormalization factors as needed to renormalize operator matrix elements. These kind of external momenta are encoded in special boundary conditions for the fermion fields. The method allows to compute certain renormalization factors on the lattice that would have been very difficult, if not impossible, to compute with standard methods. As a result we give the renormalization group invariant step scaling function for a twist-2 operator corresponding to the average momentum of non-singlet quark densities.


title: Pressure and non-linear susceptibilities in QCD at finite chemical potentials

abstract: When the free energy density of QCD is expanded in a series in the chemical potential, mu, the Taylor coefficients are the non-linear quark number susceptibilities. We show that these depend on the prescription for putting chemical potential on the lattice, making all extrapolations in chemical potential prescription dependent at finite lattice spacing. To put bounds on the prescription dependence, we investigate the magnitude of the non-linear susceptibilities over a range of temperature, T, in QCD with two degenerate flavours of light dynamical quarks at lattice spacing 1/4T. The prescription dependence is removed in quenched QCD through a continuum extrapolation, and the dependence of the pressure, P, on mu is obtained.


title: The Determination of the Quark-Gluon Mixed Condensate (anti-Q sigma G Q) from Lattice QCD

abstract: We study the quark-gluon mixed condensate g<\bar{q} sigma G q>, using the SU(3)c lattice QCD with the Kogut-Susskind fermion at the quenched level. We generate 100 gauge configurations on the 16^4 lattice with \beta = 6.0, and perform the measurement of the mixed condensate at 16 points in each gauge configuration for each current quark mass of m_q=21, 36, 52 MeV. Using the 1600 data for each m_q, we find the ratio between the mixed condensate and the quark condensate, m_0^2 = g<\bar{q} sigma G q> / <\bar{q}q> \simeq 2.5 GeV^2 at the lattice scale of 1/a \simeq 2 GeV in the chiral limit. The large value of the mixed condensate suggests its importance in the operator product expansions in QCD. We study also chiral restoration at finite temperature in terms of the mixed condensate, which is another chiral order parameter. We present the lattice QCD results of the mixed condensate at finite temperature.


title: Parton distribution functions in the pion from lattice QCD

abstract: We analyze the moments of parton distribution functions in the pion calculated in lattice QCD, paying particular attention to their chiral extrapolation. Using the lowest three non-trivial moments calculated on the lattice, we assess the accuracy with which the x dependence of both the valence and sea quark distributions in the pion can be extracted. The resulting valence quark distributions at the physical pion mass are in fair agreement with existing Drell-Yan data, but the statistical errors are such that one cannot yet confirm (or rule out) the large-x behavior expected from hadron helicity conservation in perturbative QCD. However, one can expect that the next generation of calculations in lattice QCD will allow one to extract parton distributions with a level of accuracy comparable with current experiments.


title: SU(N) Glueball Masses in 2+1 Dimensions

abstract: We calculate the masses of the lowest lying eigenstates of improved SU(2), SU(3), SU(4) and SU(5) Hamiltonian lattice gauge theory (LGT) in 2+1 dimensions using an analytic variational approach. The ground state is approximated by a one plaquette trial state and mass gaps are calculated in the symmetric and antisymmetric sectors by minimising over a suitable basis of rectangular states. Analytic techniques are developed to handle the group integrals arising in the calculation.


title: Infinite temperature limit of meson spectral functions calculated on the lattice

abstract: We analyze the cut-off dependence of mesonic spectral functions calculated at finite temperature on Euclidean lattices with finite temporal extent. In the infinite temperature limit we present analytic results for lattice spectral functions calculated with standard Wilson fermions as well as a truncated perfect action. We explicitly determine the influence of `Wilson doublers' on the high momentum structure of the mesonic spectral functions and show that this cut-off effect is strongly suppressed when using an improved fermion action.


title: The Large N Glueball Mass Spectrum in 2+1 Dimensions

abstract: In this paper we explore the large N limit of the glueball mass spectrum for 2+1 dimensional pure gauge theory. We employ Hamiltonian lattice gauge theory (LGT) and analytic variational techniques to calculate glueball masses for finite values of N < 26. The results are then accurately extrapolated to infinite N. An empirical observation is discussed which links the infinite N glueball spectrum to a simple oscillator model.


title: Nucleon electromagnetic form factors on the lattice and in chiral effective field theory

abstract: We compute the electromagnetic form factors of the nucleon in quenched lattice QCD, using non-perturbatively improved Wilson fermions, and compare the results with phenomenology and chiral effective field theory.


title: Charmonium Spectrum from Quenched Lattice QCD with Tadpole Improvement Action

abstract: We report here our lattice simulation on the charmonium spectra in the quenched approximation. Because the full adjustment on the nonperturbative parameters such as $C_E$, $C_B$, $m_0a_s$ and $r_s$ needs many calculation time, we only adjust two of them, $m_0a_s$ and $r_s$($\xi_3$ plays this role in the paper) but with some rescale for mass splitting. After the rescale, we find that our results are in agreement with the experiment ones.


title: World-Line Path Integral Study of Supersymmetry Breaking in the Wess-Zumino Model

abstract: We study supersymmetry breaking in the lattice N=1 Wess-Zumino model by the world-line path integral algorithm. The ground state energy and supersymmetric Ward identities are exploited to support the expected symmetry breaking in finite volume. Non-Gaussian fluctuations of the topological charge are discussed and related to the infinite volume transition.


title: SU(N) Lattice Gauge Theory on a Single Cube

abstract: In this paper we study the viability of persuing analytic variational techniques for the calculation of glueball masses in 3+1 dimensional Hamiltonian lattice gauge theory (LGT) in the pure gauge sector. We discuss the major problems presented by a move from 2+1 to 3+1 dimensions and develop analytic techniques to approximate the integrals appearing in 3+1 dimensional variational glueball mass calculations. We calculate $0^{++}$ and $1^{+-}$ glueball masses on a lattice consisting of a single cube. Despite the use of a very simplistic model, promising signs of an approach to asymptotic scaling is displayed by the SU(N) $1^{+-}$ glueball mass as N is increased.


title: Large N reduction in the continuum three dimensional Yang-Mills theory

abstract: Numerical and theoretical evidence leads us to propose the following: Three dimensional Euclidean Yang-Mills theory in the planar limit undergoes a phase transition on a torus of side $l=l_c$. For $l>l_c$ the planar limit is $l$-independent, as expected of a non-interacting string theory. We expect the situation in four dimensions to be similar.


title: Color Confinement in Lattice Landau Gauge

abstract: The problem of color confinement in lattice Landau gauge QCD is briefly reviewed in the light of Kugo-Ojima criterion and Gribov-Zwanziger theory. The infrared properties of the Landau gauge QCD is studied by the measurement of the gluon propagator, ghost propagator, Kugo-Ojima parameter and the running coupling via $\beta=6.4,32^4,48^4$ lattice simulation. The data are analyzed in the principle of minimum sensitivity(PMS) in the $\widetilde{MOM}$ scheme. We observe the running coupling $\alpha_s(q)$ has a maximum about 1 and decreases near q=0. The Kugo-Ojima parameter, which is expected to be 1 was about 0.8.


title: Chiral Anomalies in the Reduced (or Matrix) Model

abstract: We show that, with an appropriate choice of a Dirac operator, there is a remnant of chiral anomalies in the reduced model in which there is no coordinate dependences of the gauge field. This result is obtained by exploring a topological nature of chiral anomalies associated to Ginsparg-Wilson-type lattice Dirac operators.


title: Accelerating the Hybrid Monte Carlo algorithm

abstract: An algorithm for separating the high- and low-frequency molecular dynamics modes in Hybrid Monte Carlo simulations of gauge theories with dynamical fermions is presented. The separation is based on splitting the pseudo-fermion action into two parts, as was initially proposed by Hasenbusch. We propose to introduce different evolution time-scales for each part. We test our proposal in realistic simulations of two-flavor O(a) improved Wilson fermions. A speed-up of more than a factor of three compared to the standard HMC algorithm is observed in a typical run.


title: N=(1,1) Super Yang-Mills on a (2+1) Dimensional Transverse Lattice with one Exact Supersymmetry

abstract: We present a formulation of N=(1,1), Super Yang-Mills theory in 2+1 dimensions using a transverse lattice methods that exactly preserves one supersymmetry. First, using a Lagrangian approach we obtain a standard transverse lattice formulation of the Hamiltonian. We then show that the Hamiltonian also can be written discretely as the square of a supercharge and that this produces a different result. Problems associated with the discrete realization of the full supercharge algebra are discussed. Numerically we solve for the bound states of the theory in the large N_c approximation and we find good convergence. We show that the number of fermion and boson massless bound states are closely related. Also we find that this theory admits winding states in the transverse direction and that their masses vary inversely with the winding number.


title: Spin-3/2 Nucleon and Delta Baryons in Lattice QCD

abstract: We present first results for masses of spin-3/2 N and Delta baryons in lattice QCD using Fat-Link Irrelevant Clover (FLIC) fermions. Spin-3/2 interpolating fields providing overlap with both spin-3/2 and spin-1/2 states are considered. In the isospin-1/2 sector, we observe, after appropriate spin and parity projection, a strong signal for the J^P=3/2^- state together with a weak but discernible signal for the 3/2^+ state with a mass splitting near that observed experimentally. We also find good agreement between the 1/2^+/- masses and earlier nucleon mass simulations with the standard spin-1/2 interpolating field. For the isospin-3/2 Delta states, clear mass splittings are observed between the various 1/2^+/- and 3/2^+/- channels, with the calculated level orderings in good agreement with those observed empirically.


title: Optimal lattice domain-wall fermions with finite Ns

abstract: I review the lattice formulations of vector-like gauge theories (e.g. QCD) with domain-wall fermions, and discuss how to optimize the chiral symmetry for any finite N_s (in the fifth dimension), as well as to eliminate its dependence on a_5.


title: Nonperturbative QCD, gauge-fixing, Gribov copies, and the lattice

abstract: Perturbative QCD uses the Faddeev-Popov gauge-fixing procedure, which leads to ghosts and the local BRST invariance of the gauge-fixed perturbative QCD action. In the asymptotic regime, where perturbative QCD is relevant, Gribov copies can be neglected. In the nonperturbative regime, one must adopt either a nonlocal Gribov-copy free gauge (e.g., Laplacian gauge) or attempt to maintain local BRST invariance at the expense of admitting Gribov copies. These issues are explored. In addition, we discuss the relationship between recent Dyson-Schwinger based model calculations of the infrared behavior of QCD Green's functions and the lattice calculation of these quantities.


title: High-Precision Lattice QCD Confronts Experiment

abstract: We argue that high-precision lattice QCD is now possible, for the first time, because of a new improved staggered quark discretization. We compare a wide variety of nonperturbative calculations in QCD with experiment, and find agreement to within statistical and systematic errors of 3% or less. We also present a new determination of alpha_msbar(Mz); we obtain 0.121(3). We discuss the implications of this breakthrough for phenomenology and, in particular, for heavy-quark physics.


title: Chiral Logs in Quenched QCD

abstract: The quenched chiral logs are examined on a $16^3 \times 28$ lattice with Iwasaki gauge action and overlap fermions. The pion decay constant $f_{\pi}$ is used to set the lattice spacing, $a = 0.200(3) {\rm fm}$. With pion mass as low as $\sim 180 {\rm MeV}$, we see the quenched chiral logs clearly in $m_{\pi}^2/m$ and $f_P$, the pseudoscalar decay constant. We analyze the data to determine how low the pion mass needs to be in order for the quenched one-loop chiral perturbation theory ($\chi$PT) to apply. With the constrained curve-fitting method, we are able to extract the quenched chiral log parameter $\delta$ together with other low-energy parameters. Only for $m_{\pi} \leq 300 {\rm MeV}$ do we obtain a consistent and stable fit with a constant $\delta$ which we determine to be 0.24(3)(4) (at the chiral scale $\Lambda_{\chi}=0.8 {\rm GeV}$). By comparing to the $12^3 \times 28$ lattice, we estimate the finite volume effect to be about 2.7% for the smallest pion mass. We also fitted the pion mass to the form for the re-summed cactus diagrams and found that its applicable region is extended farther than the range for the one-loop formula, perhaps up to $m_{\pi} \sim 500-600$ MeV. The scale independent $\delta$ is determined to be 0.20(3) in this case. We study the quenched non-analytic terms in the nucleon mass and find that the coefficient $C_{1/2}$ in the nucleon mass is consistent with the prediction of one-loop $\chi$PT\@. We also obtain the low energy constant $L_5$ from $f_{\pi}$. We conclude from this study that it is imperative to cover only the range of data with the pion mass less than $\sim 300 {\rm MeV}$ in order to examine the chiral behavior of the hadron masses and decay constants in quenched QCD and match them with quenched one-loop $\chi$PT\@.


title: Non-positive fermion determinants in lattice supersymmetry

abstract: We find that fermion determinants are not generally positive in a recent class of constructions with explicit lattice supersymmetry. These involve an orbifold of supersymmetric matrix models, and have as their target (continuum) theory (2,2) 2-dimensional super-Yang-Mills. The fermion determinant is shown to be identically zero for all boson configurations due to the existence of a zeromode fermion inherited from the "mother theory." Once this eigenvalue is factored out, the fermion determinant generically has arbitrary complex phase. We discuss the implications of this result for simulation of the models.


title: Lattice QCD and the Schwarz alternating procedure

abstract: A numerical simulation algorithm for lattice QCD is described, in which the short- and long-distance effects of the sea quarks are treated separately. The algorithm can be regarded, to some extent, as an implementation at the quantum level of the classical Schwarz alternating procedure for the solution of elliptic partial differential equations. No numerical tests are reported here, but theoretical arguments suggest that the algorithm should work well also at small quark masses.