title: Behind the success of the quark model

abstract: The ground-state three-quark (3Q) potential $V_{\rm 3Q}^{\rm g.s.}$ and the excited-state 3Q potential $V_{\rm 3Q}^{\rm e.s.}$ are studied using SU(3) lattice QCD at the quenched level. For more than 300 patterns of the 3Q systems, the ground-state potential $V_{\rm 3Q}^{\rm g.s.}$ is investigated in detail in lattice QCD with $12^3\times 24$ at $\beta=5.7$ and with $16^3\times 32$ at $\beta=5.8, 6.0$. As a result, the ground-state potential $V_{\rm 3Q}^{\rm g.s.}$ is found to be well described with Y-ansatz within the 1%-level deviation. From the comparison with the Q-$\rm\bar Q$ potential, we find the universality of the string tension as $\sigma_{\rm 3Q}\simeq\sigma_{\rm Q\bar Q}$ and the one-gluon-exchange result as $A_{\rm 3Q}\simeq\frac12 A_{\rm Q\bar Q}$. The excited-state potential $V_{\rm 3Q}^{\rm e.s.}$ is also studied in lattice QCD with $16^3\times 32$ at $\beta=5.8$ for 24 patterns of the 3Q systems.The energy gap between $V_{\rm 3Q}^{\rm g.s.}$ and $V_{\rm 3Q}^{\rm e.s.}$, which physically means the gluonic excitation energy, is found to be about 1GeV in the typical hadronic scale, which is relatively large compared with the excitation energy of the quark origin. This large gluonic excitation energy justifies the great success of the simple quark model.


title: Static Three-Quark Potential in SU(3) Lattice QCD Monte-Carlo Simulation

abstract: The static three-quark(3Q) potential is investigated in the SU(3) lattice QCD at the quenched level, using the standard Wilson action with $\beta =5.7 (a\simeq 0.2{\rm fm})$ and $12^3\times 24$. With the 3Q Wilson loop, the ground-state 3Q potential $V_{\rm 3Q}$ is calculated for 16 patterns of the 3Q configuration, using the smearing technique, which enhances the ground-state overlap of the 3Q operator. With the accuracy better than a few %, the lattice QCD data for $V_{\rm 3Q}$ are well described by a sum of the Coulomb term, a constant and the linear term which is proportional to the minimal length $L_{\rm min}$ of the flux-tube linking the 3 quarks. Our results seem to support the Y-type flux picture rather than $\Delta$-type one. From the comparison of $V_{\rm 3Q}$ with the ${\rm Q}$-${\rm \bar{Q}}$ potential, we find the universality of the string tension, $\sigma_{\rm 3Q}\simeq \sigma_{\rm Q\bar{Q}}$ as well as the OGE relation on the Coulomb coefficient, $A_{\rm 3Q}\simeq {1/2}A_{\rm Q\bar{Q}}$.


title: Thermal Width Broadening of the 0++ Glueball Spectrum near Tc

abstract: We study the 0++ glueball correlator constructed with SU(3) anisotropic quenched lattice QCD at various temperature taking into account the possible existence of the thermal width in the ground-state peak. For this purpose, we adopt the Breit-Wigner ansatz, analysing the lattice data obtained with 5,500-9,900 gauge configurations at each T. The results indicate the significant thermal width broadening as Gamma(Tc) \sim 300 MeV with a reduction in the peak center as Delta omega_0(Tc) \sim 100 MeV in the vicinity of the critical temperature Tc.


title: Dynamical Regge Calculus as Lattice Quantum Gravity

abstract: We propose a hybrid model of simplicial quantum gravity by performing at once dynamical triangulations and Regge calculus. A motive for the hybridization is to give a dynamical description of topology-changing processes of Euclidean spacetime. In addition, lattice diffeomorphisms as invariance of the simplicial geometry are generated by certain elementary moves in the model. We attempt also a lattice-theoretic derivation of the black hole entropy using the symmetry. Furthermore, numerical simulations of 3D pure gravity are carried out,exhibiting a large hysteresis between two phases. We also measure geometric properties of Euclidean `time slice' based on a geodesic distance, resulting in a fractal structure in the strong-coupling phase. Our hybrid model not only reproduces numerical results consistent with those of dynamical triangulations and Regge calculus, but also opens a possibility of studying quantum black hole physics on the lattice.


title: Vortex critical behavior at the de-confinement phase transition

abstract: The de-confinement phase transition in SU(2) Yang-Mills theory is revisited in the vortex picture. Defining the world sheets of the confining vortices by maximal center projection, the percolation properties of the vortex lines in the hypercube consisting of the time axis and two spatial axis are studied. Using the percolation cumulant, the temperature for the percolation transition is seen to be in good agreement with the critical temperature of the thermal transition. The finite size scaling function for the cumulant is obtained. The critical index of the finite size scaling function is consistent with the index of the 3D Ising model.


title: A perturbative determination of O(a) boundary improvement coefficients for the Schr\"odinger Functional coupling at 1-loop with improved gauge actions

abstract: We determine O($a$) boundary improvement coefficients up to 1-loop level for the Schr\"odinger Functional coupling with improved gauge actions including plaquette and rectangle loops. These coefficients are required to implement 1-loop O($a$) improvement in full QCD simulations for the coupling with the improved gauge actions. To this order, lattice artifacts of step scaling function for each improved gauge action are also investigated. In addition, passing through the SF scheme, we estimate the ratio of $\Lambda$-parameters between the improved gauge actions and the plaquette action more accurately.


title: Pion and Kaon masses in Staggered Chiral Perturbation Theory

abstract: We show how to compute chiral logarithms that take into account both the $\cO(a^2)$ taste-symmetry breaking of staggered fermions and the fourth-root trick that produces one taste per flavor. The calculation starts from the Lee-Sharpe Lagrangian generalized to multiple flavors. An error in a previous treatment by one of us is explained and corrected. The one loop chiral logarithm corrections to the pion and kaon masses in the full (unquenched), partially quenched, and quenched cases are computed as examples.


title: Lattice chiral symmetry with hopping interactions

abstract: We formulate Dirac fermions on a (1+1)-dimensional lattice based on a Hamiltonian formalism. The species doubling problem of the lattice fermion is resolved by introducing hopping interactions that mix left- and right-handed fermions around the momentum boundary. Approximate chiral symmetry is realized on the lattice. The deviation of the fermion propagator from the continuum one is small.


title: Advanced Algorithms for the Simulation of Gauge Theories with Dynamical Fermionic Degrees of Freedom

abstract: This thesis is focused on the implementation and the application of a novel kind of algorithm which is expected to overcome the limitations of older schemes. This new algorithm is named Multiboson Method. It allows to simulate an arbitrary number of dynamical fermion flavors, which can in principle have different masses. It will be shown that it exhibits better scaling properties for light fermions than other methods. Therefore, it has the potential to become the method of choice.


title: Heavy Quark Free Energies and Screening in SU(2) Gauge Theory

abstract: We study the properties of the free energy of infinitely heavy quark anti-quark pair in SU(2) gauge theory. By means of lattice Monte Carlo simulations we calculated the free energies in the singlet, triplet and color averaged channels, both in the confinement and in the deconfinement phase. The singlet and triplet free energies are defined in Coulomb gauge which is equivalent to their gauge invariant definitions recently introduced by Philipsen. We analyzed the short and the long distance behavior, making comparisons with the zero temperature case. The temperature dependence of the electric screening mass is carefully investigated. The order of the deconfining transition is manifest in the results near T_c and it allows a reliable test of a recently proposed method to renormalize the Polyakov loop.


title: Moments of Nucleon Generalized Parton Distributions in Lattice QCD

abstract: Calculation of moments of generalized parton distributions in lattice QCD requires more powerful techniques than those previously used to calculate moments of structure functions. Hence, we present a novel approach that exploits the full information content from a given lattice configuration by measuring an overdetermined set of lattice observables to provide maximal statistical constraints on the generalized form factors at a given virtuality, t. In an exploratory investigation using unquenched QCD configurations at intermediate sea quark masses, we demonstrate that our new technique is superior to conventional methods and leads to reliable numerical signals for the n=2 flavor singlet generalized form factors up to 3 GeV^2. The contribution from connected diagrams in the flavor singlet sector to the total quark angular momentum is measured to an accuracy of the order of one percent.


title: B --> pi and B --> K transitions in partially quenched chiral perturbation theory

abstract: We study the properties of the B-->pi and B-->K transition form factors in partially quenched QCD by using the approach of partially quenched chiral perturbation theory combined with the static heavy quark limit. We show that the form factors change almost linearly when varying the value of the sea quark mass, whereas the dependence on the valence quark mass contains both the standard and chirally divergent (quenched) logarithms. A simple strategy for the chiral extrapolations in the lattice studies with Nsea=2 is suggested. It consists of the linear extrapolations from the realistically accessible quark masses, first in the sea and then in the valence quark mass. From the present approach, we estimate the uncertainty induced by such extrapolations to be within 5%.


title: A lattice study of the two-dimensional Wess Zumino model

abstract: We present results from a numerical simulation of the two-dimensional Euclidean Wess-Zumino model. In the continuum the theory possesses N=1 supersymmetry. The lattice model we employ was analyzed by Golterman and Petcher in \cite{susy} where a perturbative proof was given that the continuum supersymmetric Ward identities are recovered without finite tuning in the limit of vanishing lattice spacing. Our simulations demonstrate the existence of important non-perturbative effects in finite volumes which modify these conclusions. It appears that in certain regions of parameter space the vacuum state can contain solitons corresponding to field configurations which interpolate between different classical vacua. In the background of these solitons supersymmetry is partially broken and a light fermion mode is observed. At fixed coupling the critical mass separating phases of broken and unbroken supersymmetry appears to be volume dependent. We discuss the implications of our results for continuum supersymmetry breaking.


title: The pseudo-Goldstone spectrum of 2-colour QCD at finite density

abstract: We examine the spectrum of 2-colour lattice QCD with 4 continuum flavours at a finite chemical potential ($\mu$) for quark-number, on a $12^3 \times 24$ lattice. First we present evidence that the system undergoes a transition to a state with a diquark condensate, which spontaneously breaks quark number at $\mu=m_\pi/2$, and that this transition is mean field in nature. We then examine the 3 states that would be Goldstone bosons at $\mu=0$ for zero Dirac and Majorana quark masses. The predictions of chiral effective Lagrangians give a good description of the behaviour of these masses for $\mu < m_\pi/2$. Except for the heaviest of these states, these predictions diverge from our measurements, once $\mu$ is significantly greater than $m_\pi/2$. However, the qualitative behaviour of these masses, indicates that the physics is very similar to that predicted by these effective Lagrangians, and there is some indication that at least part of these discrepancies is due to saturation, a lattice artifact.


title: Lattice Study of the Massive Schwinger Model with $\theta$ Term under L\"{u}scher's "Admissibility" Condition

abstract: L\"uscher's ``admissibility'' condition on the gauge field space plays an essential role in constructing lattice gauge theories which has exact chiral symmetries. We apply the gauge action proposed by L\"uscher with the domain-wall fermion action to the numerical simulation of the massive Schwinger model. We find this action can generate configurations in each topological sector separately without any topology changes. By developing a new method to sum over different topological sectors, we calculate meson masses in nonzero-$\theta$ vacuum.


title: Theta-vacuum systems via real action simulations

abstract: Inspired by the results of the Ising model within an imaginary external magnetic field, we introduce a transformation in quantum systems with a theta-vacuum term that amounts to a rescaling of z=cos(theta/2). Making use of this transformation we are able to determine the order parameter as a function of theta. The approach is successfully tested in models with both broken and unbroken CP symmetry at theta=pi.


title: A Ginsparg-Wilson Approach to Lattice Chern-Simons Theory

abstract: The concept of lattice modified symmetry formulations is adapted to the parity symmetry of gauge fields and applied to the pure Abelian Chern-Simons action in three dimensions. We derive an analogue of the Ginsparg-Wilson relation for the parity anti-symmetry, which is motivated from the perfect lattice action, and which we denote as the Chern-Simons-Ginsparg-Wilson relation (CSGWR). In addition to the overlap type solutions, we construct explicitly simple and local polynomial solutions to the CSGWR. We show that these actions are exactly invariant under a lattice modified parity transformation. That transformation is local as well, and it turns into the standard parity transformation in the continuum limit.


title: The Equation of State for Two Flavor QCD at Non-zero Chemical Potential

abstract: We present results of a simulation of QCD on a 4x16^3 lattice with 2 continuum flavors of p4-improved staggered fermion with mass m/T=0.4. Derivatives of the thermodynamic grand potential with respect to quark chemical potential mu_q up to fourth order are calculated, enabling estimates of the pressure, quark number density and associated susceptibilities as functions of mu_q via Taylor series expansion. Discretisation effects associated with various staggered fermion formulations are discussed in some detail. In addition it is possible to estimate the radius of convergence of the expansion as a function of temperature. We also discuss the calculation of energy and entropy densities which are defined via mixed derivatives of the thermodynamic grand potential with respect to the bare couplings and quark masses.


title: Hybrid configuration content of heavy S-wave mesons

abstract: We use the non-relativistic expansion of QCD (NRQCD) on the lattice to study the lowest hybrid configuration contribution to the ground state of heavy S-wave mesons. Using lowest-order lattice NRQCD to create the heavy-quark propagators, we form a basis of ``unperturbed'' S-wave and hybrid states. We then apply the lowest-order coupling of the quark spin and chromomagnetic field at an intermediate time slice to create ``mixed'' correlators between the S-wave and hybrid states. From the resulting amplitudes, we extract the off-diagonal element of our two-state Hamiltonian. Diagonalizing this Hamiltonian gives us the admixture of hybrid configuration within the meson ground state. The present effort represents a continuation of previous work: the analysis has been extended to include lattices of varying spacings, source operators having better overlap with the ground states, and the pseudoscalar (along with the vector) channel. Results are presented for bottomonium ($\Upsilon$, $\eta_b^{}$) using three different sets of quenched lattices. We also show results for charmonium ($J/\psi$, $\eta_c^{}$) from one lattice set, although we note that the non-relativistic approximation is not expected to be very good in this case.


title: Reducing Residual-Mass Effects for Domain-Wall Fermions

abstract: It has been suggested to project out a number of low-lying eigenvalues of the four-dimensional Wilson--Dirac operator that generates the transfer matrix of domain-wall fermions in order to improve simulations with domain-wall fermions. We investigate how this projection method reduces the residual chiral symmetry-breaking effects for a finite extent of the extra dimension. We use the standard Wilson as well as the renormalization--group--improved gauge action. In both cases we find a substantially reduced residual mass when the projection method is employed. In addition, the large fluctuations in this quantity disappear.


title: Infrared Features of the Landau Gauge QCD

abstract: The infrared features of Landau gauge QCD are studied by the lattice simulation of $\beta=6.0, 16^4, 24^4, 32^4$ and $\beta=6.4, 32^4, 48^4$. We adopt two definitions of the gauge field; 1) $U-$linear 2) $\log U$ and measured the gluon propagator and ghost propagator. Infrared singularity of the gluon propagator is less than that of tree level result but the gluon propagator at 0 momentum remains finite. The infrared singularity of ghost propagator is stronger than the tree level. The QCD running coupling measured by using the gluon propagator and the ghost propagator has a maximum $\alpha_s(p)\simeq 1$ at around $p=0.5GeV$ and decreases as $p$ approaches 0. The data are analyzed in use of formula of the principle of minimal sensitivity(PMS), the effective charge method and the contour-improved perturbation method, which suggest necessity of the resummation of perturbation series in the infrared region together with existence of the infrared fixed point. Kugo-Ojima parameter saturates at about -0.8 in contrast to the theoretically expected value -1.


title: Axial anomaly in the reduced model: Higher representations

abstract: The axial anomaly arising from the fermion sector of $\U(N)$ or $\SU(N)$ reduced model is studied under a certain restriction of gauge field configurations (the ``$\U(1)$ embedding'' with $N=L^d$). We use the overlap-Dirac operator and consider how the anomaly changes as a function of a gauge-group representation of the fermion. A simple argument shows that the anomaly vanishes for an irreducible representation expressed by a Young tableau whose number of boxes is a multiple of $L^2$ (such as the adjoint representation) and for a tensor-product of them. We also evaluate the anomaly for general gauge-group representations in the large $N$ limit. The large $N$ limit exhibits expected algebraic properties as the axial anomaly. Nevertheless, when the gauge group is $\SU(N)$, it does not have a structure such as the trace of a product of traceless gauge-group generators which is expected from the corresponding gauge field theory.


title: Thermal dilepton rates and in-medium hadron properties from lattice QCD

abstract: Recent progress of lattice investigations in thermal physics is summarized in this contribution. Hadronic spectral functions can be reconstructed from correlation functions in Euclidean time based on the Maximum Entropy Method without a priori assumptions on the spectral shape. The thermal modifications of hadron properties are investigated in the scalar, pseudo-scalar, vector and axial-vector channels near and above the deconfinement transition temperature for charmonium systems as well as for light quarks. Moreover, the reconstructed vector meson spectral function allows to extract the thermal cross section for the production of dilepton pairs at vanishing momentum.


title: Impossibility of spontaneously breaking local symmetries and the sign problem

abstract: Elitzur's theorem stating the impossibility of spontaneous breaking of local symmetries in a gauge theory is reexamined. The existing proofs of this theorem rely on gauge invariance as well as positivity of the weight in the Euclidean partition function. We examine the validity of Elitzur's theorem in gauge theories for which the Euclidean measure of the partition function is not positive definite. We find that Elitzur's theorem does not follow from gauge invariance alone. We formulate a general criterion under which spontaneous breaking of local symmetries in a gauge theory is excluded. Finally we illustrate the results in an exactly solvable two dimensional abelian gauge theory.


title: Non-perturbative renormalisation and improvement of the local vector current for quenched and unquenched Wilson fermions

abstract: By considering the local vector current between nucleon states and imposing charge conservation, we determine its renormalisation constant and quark mass improvement coefficient for Symanzik $O(a)$ improved Wilson fermions. The computation is first performed for quenched fermions (and for completeness also with unimproved fermions) and compared against known results. The two-flavour unquenched case is then considered.


title: Factorization Method for Simulating QCD at Finite Density

abstract: We propose a new method for simulating QCD at finite density. The method is based on a general factorization property of distribution functions of observables, and it is therefore applicable to any system with a complex action. The so-called overlap problem is completely eliminated by the use of constrained simulations. We test this method in a Random Matrix Theory for finite density QCD, where we are able to reproduce the exact results for the quark number density.


title: Light quark masses, chiral condensate and quark-gluon condensate in quenched lattice QCD with exact chiral symmetry

abstract: We determine several quantities in quenched lattice QCD with exact chiral symmetry. For 100 gauge configurations generated with Wilson gauge action at $ \beta = 6.0 $ on the $ 16^3 \times 32 $ lattice, we compute quenched quark propagators for 13 bare quark masses. The pion decay constant is extracted from the pion propagator, and from which the inverse lattice spacing is determined to be $ a^{-1} = 1.979(6) $ GeV. The parameters ($ C, \delta, B $) in the pseudoscalar meson mass formula in quenched chiral perturbation theory (q$\chi$PT) to one-loop order are determined. Further, we measure the index (topological) susceptibility of these 100 gauge configurations, $ \chi_t = (175 \pm 6 {MeV})^4 $, from which we obtain an estimate of the mass of $ \eta' $ in q$\chi$PT, and the coefficient of quenched chiral logarithm, both in good agreement with the values determined from the pion masses, as well as with the theoretical estimates. With our values of $ C, \delta, B $, the experimental inputs of pion and kaon masses, and the pion decay constant, we determine the light quark masses: $ m_{u,d} = 4.1 \pm 0.3 $ MeV, and $ m_s = 92 \pm 9 $ MeV, in the $ \bar{MS} $ scheme at scale $ \mu = 2 $ GeV. Also, we determine the quark condensate $<\bar q q > = -(250 \pm 3 {MeV})^3 $, and the quark-gluon condensate $ g < \bar q \sigma_{\mu\nu} F_{\mu\nu} q > = -(434 \pm 4 {MeV})^5 $, in the $ \bar{MS} $ scheme at scale 2 GeV.


title: QCD meson spectrum in the large $N_C$ limit

abstract: The low energy mass spectrum of QCD in the large $N_C$ limit is computed. The low-lying states in the meson spectrum are explicitly evaluated up to the fourth order in the strong coupling perturbative expansion. The 't Hooft limit is smooth and the meson masses are in very good agreement with the experimental values.


title: Heavy quark masses in the continuum limit of quenched Lattice QCD

abstract: We compute charm and bottom quark masses in the quenched approximation and in the continuum limit of lattice QCD. We make use of a step scaling method, previously introduced to deal with two scale problems, that allows to take the continuum limit of the lattice data. We determine the RGI quark masses and make the connection to the MSbar scheme. The continuum extrapolation gives us a value m_b^{RGI} = 6.73(16) GeV for the b-quark and m_c^{RGI} = 1.681(36) GeV for the c-quark, corresponding respectively to m_b^{MSbar}(m_b^{MSbar}) = 4.33(10) GeV and m_c^{MSbar}(m_c^{MSbar}) = 1.319(28) GeV. The latter result, in agreement with current estimates, is for us a check of the method. Using our results on the heavy quark masses we compute the mass of the Bc meson, M_{Bc} = 6.46(15) GeV.


title: Critical behaviour and scaling functions of the three-dimensional O(6) model

abstract: We numerically investigate the three-dimensional O(6) model on 12^3 to 120^3 lattices within the critical region at zero magnetic field, as well as at finite magnetic field on the critical isotherm and for several fixed couplings in the broken and the symmetric phase. We obtain from the Binder cumulant at vanishing magnetic field the critical coupling J_c=1.42865(3). The universal value of the Binder cumulant at this point is g_r(J_c)=-1.94456(10). At the critical coupling, the critical exponents \gamma=1.604(6), \beta=0.425(2) and \nu=0.818(5) are determined from a finite-size-scaling analysis. Furthermore, we verify predicted effects induced by massless Goldstone modes in the broken phase. The results are well described by the perturbative form of the model's equation of state. Our O(6)-result is compared to the corresponding Ising, O(2) and O(4) scaling functions. Finally, we study the finite-size-scaling behaviour of the magnetisation on the pseudocritical line.


title: The Supersymmetric Ward-Takahashi Identity in 1-Loop Lattice Perturbation Theory. I. General Procedure

abstract: The one-loop corrections to the lattice supersymmetric Ward-Takahashi identity (WTi) are investigated in the off-shell regime. In the Wilson formulation of the N=1 supersymmetric Yang-Mills (SYM) theory, supersymmetry (SUSY) is broken by the lattice, by the Wilson term and is softly broken by the presence of the gluino mass. However, the renormalization of the supercurrent can be realized in a scheme that restores the continuum supersymmetric WTi (once the on-shell condition is imposed). The general procedure used to calculate the renormalization constants and mixing coefficients for the local supercurrent is presented. The supercurrent not only mixes with the gauge invariant operator $T_\mu$. An extra mixing with other operators coming from the WTi appears. This extra mixing survives in the continuum limit in the off-shell regime and cancels out when the on-shell condition is imposed and the renormalized gluino mass is set to zero. Comparison with numerical results are also presented.


title: Monopole clusters at short and large distances

abstract: We present measurements of various geometrical characteristics of monopole clusters in SU(2) lattice gauge theory. The maximal Abelian projection is employed and both infinite, or percolating cluster and finite clusters are considered. In particular, we observe scaling for average length of segments of the percolating cluster between self-crossings, correlators of vacuum monopole currents, angular correlation between links along trajectories. Short clusters are random walks and their spectrum in length corresponds to free particles. At the hadronic scale, on the other hand, the monopole trajectories are no longer random walks. Moreover, we argue that the data on the density of finite clusters suggest that there are long-range correlations between finite clusters which can be understood as association of the clusters with two-dimensional surfaces, whose area scales.


title: Theta dependence of CP^9 model

abstract: We apply to the $CP^9$ model two recently proposed numerical techniques for simulation of systems with a theta term. The algorithms, successfully tested in the strong coupling limit, are applied to the weak coupling region. The results agree and errors have been evaluated and are at % level. The results scale well with the renormalization group equation and show that, for $CP^9$ in presence of a theta term, CP symmetry is spontaneously broken at $\theta=\pi$ in the continuum limit.


title: Beyond Quantum Field Theory: Chaotic Lattices?

abstract: We review the idea of chaotic quantization, based on the dynamics of classical lattice gauge systems as well as on non-abelian plasma physics in the infrared limit. The basic conjecture between Planck constant and properties of the five dimensional classical theory (h=aT) is demonstrated numerically for the U(1) lattice gauge group.


title: Heavy quark expansion parameters from lattice NRQCD

abstract: We present a lattice QCD calculation of the heavy quark expansion parameters $\mu_{\pi}^2$ and $\mu_G^2$ for heavy-light mesons and heavy-light-light baryons. The calculation is carried out on a 20$^3\times$48 lattice at $\beta$ = 6.0 in the quenched approximation, using the lattice NRQCD action for heavy quarks. We obtain the parameters $\mu_{\pi}^2$ and $\mu_G^2$ in two different methods: a direct calculation of the matrix elements and an indirect calculation through the mass spectrum, and confirm that the both methods give consistent results. We also discuss an application to the lifetime ratios.


title: Thermodynamics and in-medium hadron properties from lattice QCD

abstract: Non-perturbative studies of the thermodynamics of strongly interacting elementary particles within the context of lattice regularized QCD are being reviewed. After a short introduction into thermal QCD on the lattice we report on the present status of investigations of bulk properties. In particular, we discuss the present knowledge of the phase diagram including recent developments of QCD at non-zero baryon number density. We continue with the results obtained so far for the transition temperature as well as the temperature dependence of energy and pressure and comment on screening and the heavy quark free energies. A major section is devoted to the discussion of thermal modifications of hadron properties, taking special account of recent progress through the use of the maximum entropy method.


title: Gauge Theories with Nonunitary Parallel Transporters: A soluble Higgs model

abstract: It has been proposed to abandon the requirement that parallel transporters in gauge theories are unitary (or pseudoorthogonal). This leads to a geometric interpretation of Vierbein fields as parts of gauge fields, and nonunitary parallel transport in extra directions yields Higgs fields. In such theories, the holonomy group H is larger than the gauge group G. Here we study a 1-dimensional model with fermions which retains only the extra dimension, and which is soluble in the sense that its renormalization group flow may be exactly computed, with G=SU(2) and noncompact subgroup H of GL(2,C), or G=U(2), H=GL(2,C). In all cases the asymptotic behavior of the Higgs potential is computed, and with one possible exception for G=SU(2), H=GL(2,C), there is a flow of the action from an UV-fix point which describes a SU(2)-gauge theory with unitary parallel transporters, to a IR-fixpoint. We explain how a splitting of the masses of fermions of different flavor can arise through spontaneous symmetry breaking.


title: Numerical determination of monopole entropy in pure SU(2) QCD

abstract: We study numerically the length distributions of the infrared monopole clusters in pure SU(2) QCD. These distributions are Gaussian for all studied blocking steps of monopoles, lattice volumes and lattice coupling constant. We also investigate the monopole action for the infrared monopole clusters. The knowledge of both the length distribution and the monopole action allows us to determine the effective entropy of the monopole currents. The entropy is a descending function of blocking scale, indicating that the effective degrees of freedom of the extended monopoles are getting smaller as the blocking scale increases.


title: Localization in Lattice QCD

abstract: In this paper, we examine the phase diagram of quenched QCD with two flavors of Wilson fermions, proposing the following microscopic picture. The super-critical regions inside and outside the Aoki phase are characterized by the existence of a density of near-zero modes of the (hermitian) Wilson-Dirac operator, and thus by a non-vanishing pion condensate. Inside the Aoki phase, this density is built up from extended near-zero modes, while outside the Aoki phase, there is a non-vanishing density of exponentially localized near-zero modes, which occur in "exceptional" gauge-field configurations. Nevertheless, no Goldstone excitations appear outside the Aoki phase, and the existence of Goldstone excitations may therefore be used to define the Aoki phase in both the quenched and unquenched theories. We show that the density of localized near-zero modes gives rise to a divergent pion two-point function, thus providing an alternative mechanism for satisfying the relevant Ward identity in the presence of a non-zero order parameter. This divergence occurs when we take a "twisted" quark mass to zero, and we conclude that quenched QCD with Wilson fermions is well-defined only with a non-vanishing twisted mass. We show that this peculiar behavior of the near-zero-mode density is special to the quenched theory by demonstrating that this density vanishes in the unquenched theory outside the Aoki phase. We discuss the implications for domain-wall and overlap fermions constructed from a Wilson-Dirac kernel. We argue that both methods work outside the Aoki phase, but fail inside because of problems with locality and/or chiral symmetry, in both the quenched and unquenched theories.


title: Five-dimensional Lattice Gauge Theory as Multi-Layer World

abstract: A five-dimensional lattice space can be decomposed into a number of four-dimens ional lattices called as layers. The five-dimensional gauge theory on the lattice can be interpreted as four-dimensional gauge theories on the multi-layer with interactions between neighboring layers. In the theory, there exist two independent coupling constants; $\beta_4$ controls the dynamics inside a layer and $\beta_5$ does the strength of the inter-layer interaction.We propose the new possibility to realize the continuum limit of a five-dimensional theory using four-dimensional dynamics with large $\beta_4$ and small $\beta_5$. Our result is also related to the higher dimensional theory by deconstruction approach.


title: Low-lying meson spectrum of large $N_C$ strongly coupled lattice QCD

abstract: We compute the low energy mass spectrum of lattice QCD in the large $N_C$ limit. Expanding around a gauge-invariant ground state, which spontaneously breaks the discrete chiral symmetry, we derive an improved strong-coupling expansion and evaluate, for any value of $N_C$, the masses of the low-lying states in the unflavored meson spectrum. We then take the 't Hooft limit by rescaling $g^2 N_C\to g^2$; the 't Hooft limit is smooth and no arbitrary parameters are needed. We find, already at the fourth order of the strong coupling perturbation theory, a very good agreement between the results of our lattice computation and the known continuum values.


title: The static quark potential and scaling behavior of SU(3) l attice Yang-Mills theory

abstract: The potential between a static quark and antiquark in pure SU(3) Yang-Mills theory is evaluated non-perturbatively through computations on the lattice in the region from short to intermediate distances (0.05 fm < r < 0.8 fm). In the high energy regime the results are compared with the parameter-free prediction of perturbation theory obtained by solving the Renormalization Group equation at two and three loops. The choice of the renormalization scheme to define a running coupling turns out to be important for the accuracy of the perturbative prediction: by obtaining the running coupling through the force, perturbation theory is applicable up to alpha ~ 0.3, while from the static potential only up to alpha ~ 0.15. In the region where perturbation is supposed to be reliable, no large unexpected non-perturbative term is observed. In the second part of this work, universality and scaling behavior of different formulations of Yang-Mills theory on the lattice are discussed, in particular the Iwasaki and DBW2 actions. The deconfinement temperature T_c*r0 and the glueball masses (m_0++)*r0, (m_2++)*r0 are analyzed. Particular attention is dedicated to the violation of physical positivity which occur in these actions and the consequences in the extraction of physical quantities from Euclidean correlation functions.


title: A lattice calculation of vector meson couplings to the vector and tensor currents using chirally improved fermions

abstract: We present a quenched lattice calculation of $f_V^\perp/f_V$, the coupling of vector mesons to the tensor current normalized by the vector meson decay constant. The chirally improved lattice Dirac operator, which allows us to reach small quark masses, is used. We put emphasis on analyzing the quark mass dependence of $f_V^\perp/f_V$ and find only a rather weak dependence. Our results at the $\rho$ and $\phi$ masses agree well with QCD sum rule calculations and those from previous lattice studies.


title: Nucleon axial charge from quenched lattice QCD with domain wall fermions

abstract: We present a quenched lattice calculation of the nucleon isovector vector and axial-vector charges gV and gA. The chiral symmetry of domain wall fermions makes the calculation of the nucleon axial charge particularly easy since the Ward-Takahashi identity requires the vector and axial-vector currents to have the same renormalization, up to lattice spacing errors of order O(a^2). The DBW2 gauge action provides enhancement of the good chiral symmetry properties of domain wall fermions at larger lattice spacing than the conventional Wilson gauge action. Taking advantage of these methods and performing a high statistics simulation, we find a significant finite volume effect between the nucleon axial charges calculated on lattices with (1.2 fm)^3 and (2.4 fm)^3 volumes (with lattice spacing, a, of about 0.15 fm). On the large volume we find gA = 1.212 +/- 0.027(statistical error) +/- 0.024(normalization error). The quoted systematic error is the dominant (known) one, corresponding to current renormalization. We discuss other possible remaining sources of error. This theoretical first principles calculation, which does not yet include isospin breaking effects, yields a value of gA only a little bit below the experimental one, 1.2670 +/- 0.0030.


title: The Lattice $\Lambda$ Parameter in Domain Wall QCD

abstract: We evaluate the ratio of the scale parameter $\Lambda$ in domain wall QCD to the one in the continuum theory at one loop level incorporating the effect of massless quarks. We show that the Pauli-Villars regulator is required to subtract the unphysical massive fermion modes which emerge in the fermion loop contributions to the gluon self energy. Detailed results are presented as a function of the domain wall height $M$.


title: Vector-like Fermion Doubling and Pair-creation on a lattice

abstract: Based on the Pauli principle and Heisenberg uncertainty relationship, we show that the phenomenon of vectorlike fermion-doubling in quantum field theories on a lattice is related to the phenomenon of pair-creation of massless fermions, due to quantum-field fluctuations upon the zero-energy level of the ground state of regularized quantum field theories on the lattice.


title: Comparison of |Q|=1 and |Q|=2 gauge-field configurations on the lattice four-torus

abstract: It is known that exactly self-dual gauge-field configurations with topological charge |Q|=1 cannot exist on the untwisted continuum 4-torus. We explore the manifestation of this remarkable fact on the lattice 4-torus for SU(3) using advanced techniques for controlling lattice discretization errors, extending earlier work of De Forcrand et. al. for SU(2). We identify three distinct signals for the instability of |Q|=1 configurations, and show that these manifest themselves early in the cooling process, long before the would-be instanton has shrunk to a size comparable to the lattice discretization threshold. These signals do not appear for our |Q|=2 configurations. This indicates that these signals reflect the truly global nature of the instability, rather than local discretization effects. Monte-Carlo generated SU(3) gauge field configurations are cooled to the self-dual limit using an O(a^4)-improved gauge action chosen to have small but positive O(a^6) errors. This choice prevents lattice discretization errors from destroying instantons provided their size exceeds the dislocation threshold of the cooling algorithm. Lattice discretization errors are evaluated by comparing the O(a^4)-improved gauge-field action with an O(a^4)-improved action constructed from the square of an O(a^4)-improved lattice field-strength tensor, thus having different O(a^6) discretization errors. The number of action-density peaks, the instanton size and the topological charge of configurations is monitored. We observe a fluctuation in the total topological charge of |Q|=1 configurations, and demonstrate that the onset of this unusual behavior corresponds with the disappearance of multiple-peaks in the action density. At the same time discretization errors are minimal.


title: Observing string breaking with Wilson loops

abstract: An uncontroversial observation of adjoint string breaking is proposed, while measuring the static potential from Wilson loops only. The overlap of the Wilson loop with the broken-string state is small, but non-vanishing, so that the broken-string groundstate can be seen if the Wilson loop is long enough. We demonstrate this in the context of the (2+1)d SU(2) adjoint static potential, using an improved version of the Luscher-Weisz exponential variance reduction. To complete the picture we perform the more usual multichannel analysis with two basis states, the unbroken-string state and the broken-string state (two so-called gluelumps). As by-products, we obtain the temperature-dependent static potential measured from Polyakov loop correlations, and the fundamental SU(2) static potential with improved accuracy. Comparing the latter with the adjoint potential, we see clear deviations from Casimir scaling.


title: Lattice Gauge Fixing and the Violation of Spectral Positivity

abstract: Spectral positivity is known to be violated by some forms of lattice gauge fixing. The most notable example is lattice Landau gauge, where the effective gluon mass is observed to rise rather than fall with increasing distance. We trace this violation to the use of quenched auxiliary fields in the lattice gauge fixing process, and show that violation of spectral positivity is a general feature of quenching. We illustrate this with a simple quenched mass-mixing model in continuum field theory, and with a quenched form of the Ising model. For lattice gauge fixing associated with Abelian projection and lattice Landau gauge, we show that spectral positivity is violated by processes similar to those found in quenched QCD. For covariant gauges parametrized by a gauge-fixing parameter $\alpha $, the SU(2) gluon propagator is well described by a simple quenched mass-mixing formula. The gluon mass parameter appears to be independent of $\alpha$ for sufficiently large $\alpha$.


title: Lattice results with three quark flavours

abstract: There have been exciting advances recently in finite temperature calculations with three quark flavours and with non-zero chemical potential. The role of improved actions is explained and recent results from the Bielefeld group and MILC collaboration with improved Kogut-Susskind quarks are presented. Three new approaches to finite chemical potential are discussed.


title: Chirally improving Wilson fermions - I. O(a) improvement

abstract: We show that it is possible to improve the chiral behaviour and the approach to the continuum limit of correlation functions in lattice QCD with Wilson fermions by taking arithmetic averages of correlators computed in theories regularized with Wilson terms of opposite sign. Improved hadronic masses and matrix elements can be obtained by similarly averaging the corresponding physical quantities separately computed within the two regularizations. To deal with the problems related to the spectrum of the Wilson--Dirac operator, which are particularly worrisome when Wilson and mass terms are such as to give contributions of opposite sign to the real part of the eigenvalues, we propose to use twisted-mass lattice QCD for the actual computation of the quantities taking part to the averages. The choice $\pm \pi/2$ for the twisting angle is particularly interesting, as O($a$) improved estimates of physical quantities can be obtained even without averaging data from lattice formulations with opposite Wilson terms. In all cases little or no extra computing power is necessary, compared to simulations with standard Wilson fermions or twisted-mass lattice QCD.


title: Determination of the Improvement Coefficient c_SW up to One-Loop Order with the Conventional Perturbation Theory

abstract: We calculate the $O(a)$ improvement coefficient c_SW in the Sheikholeslami-Wohlert quark action for various improved gauge actions with six-link loops. We employ the conventional perturbation theory introducing the fictitious gluon mass to regularize the infrared divergence. Our results for some improved gauge actions are in agreement with those previously obtained with the Schr{\"o}dinger functional method.


title: Perturbative matching of staggered four-fermion operators with hypercubic fat links

abstract: We calculate the one-loop matching coefficients between continuum and lattice four-fermion operators for lattice operators constructed using staggered fermions and improved by the use of fattened links. In particular, we consider hypercubic fat links and SU(3) projected Fat-7 links, and their mean-field improved versions. We calculate only current-current diagrams, so that our results apply for operators whose flavor structure does not allow ``eye-diagrams''. We present general formulae, based on two independent approaches, and give numerical results for the cases in which the operators have the taste (staggered flavor) of the pseudo-Goldstone pion. We find that the one-loop corrections are reduced down to the 10-20% level, resolving the problem of large perturbative corrections for staggered fermion calculations of matrix elements.


title: Monte Carlo errors with less errors

abstract: We explain in detail how to estimate mean values and assess statistical errors for arbitrary functions of elementary observables in Monte Carlo simulations. The method is to estimate and sum the relevant autocorrelation functions, which is argued to produce more certain error estimates than binning techniques and hence to help toward a better exploitation of expensive simulations. An effective integrated autocorrelation time is computed which is suitable to benchmark efficiencies of simulation algorithms with regard to specific observables of interest. A Matlab code is offered for download that implements the method. It can also combine independent runs (replica) allowing to judge their consistency.


title: The apeNEXT project (Status report)

abstract: We present the current status of the apeNEXT project. Aim of this project is the development of the next generation of APE machines which will provide multi-teraflop computing power. Like previous machines, apeNEXT is based on a custom designed processor, which is specifically optimized for simulating QCD. We discuss the machine design, report on benchmarks, and give an overview on the status of the software development.


title: Glueball Regge Trajectories in (2+1) Dimensional Gauge Theories

abstract: We compute glueball masses for even spins ranging from 0 to 6, in the D=2+1 SU(2) lattice gauge theory. We do so over a wide range of lattice spacings, and this allows a well-controlled extrapolation to the continuum limit. When the resulting spectrum is presented in the form of a Chew-Frautschi plot we find that we can draw a straight Regge trajectory going through the lightest glueballs of spin 0, 2, 4 and 6. The slope of this trajectory is small and turns out to lie between the predictions of the adjoint-string and flux-tube glueball models. The intercept we find, \alpha_0 ~ -1, is much lower than is needed for this leading trajectory to play a `Pomeron-like' role of the kind it is often believed to play in D=3+1. We elaborate the Regge theory of high energy scattering in 2 space dimensions, and we conclude, from the observed low intercept, that high-energy glueball scattering is not dominated by the leading Regge pole exchange, but rather by a more complex singularity structure in the region 0 <= Re{\lambda} <= 1/2 of the complex angular momentum \lambda plane. We show that these conclusions do not change if we go to larger groups, SU(N>2), and indeed to SU(\infty), and we contrast all this with our very preliminary calculations in the D=3+1 SU(3) gauge theory.


title: Nonlinear Realization of Chiral Symmetry on the Lattice

abstract: We formulate lattice theories in which chiral symmetry is realized nonlinearly on the fermion fields. In this framework the fermion mass term does not break chiral symmetry. This property allows us to use the Wilson term to remove the doubler fermions while maintaining exact chiral symmetry on the lattice. Our lattice formulation enables us to address non-perturbative questions in effective field theories of baryons interacting with pions and in models involving constituent quarks interacting with pions and gluons. We show that a system containing a non-zero density of static baryons interacting with pions can be studied on the lattice without encountering complex action problems. In our formulation one can also decide non-perturbatively if the chiral quark model of Georgi and Manohar provides an appropriate low-energy description of QCD. If so, one could understand why the non-relativistic quark model works.


title: Chiral perturbation theory at O(a^2) for lattice QCD

abstract: We construct the chiral effective Lagrangian for two lattice theories: one with Wilson fermions and the other with Wilson sea fermions and Ginsparg-Wilson valence fermions. For each of these theories we construct the Symanzik action through order $a^2$. The chiral Lagrangian is then derived, including terms of order $a^2$, which have not been calculated before. We find that there are only few new terms at this order. Corrections to existing coefficients in the continuum chiral Lagrangian are proportional to $a^2$, and appear in the Lagrangian at order $a^2 p^2$ or higher. Similarly, O(4) symmetry breaking terms enter the Symanzik action at order $a^2$, but contribute to the chiral Lagrangian at order $a^2 p^4$ or higher. We calculate the light meson masses in chiral perturbation theory for both lattice theories. At next-to-leading order, we find that there are no order $a^2$ corrections to the valence-valence meson mass in the mixed theory due to the enhanced chiral symmetry of the valence sector.


title: Spectral Properties of the Overlap Dirac Operator in QCD

abstract: We discuss the eigenvalue distribution of the overlap Dirac operator in quenched QCD on lattices of size 8^{4}, 10^{4} and 12^{4} at \beta = 5.85 and \beta = 6. We distinguish the topological sectors and study the distributions of the leading non-zero eigenvalues, which are stereographically mapped onto the imaginary axis. Thus they can be compared to the predictions of random matrix theory applied to the \epsilon-expansion of chiral perturbation theory. We find a satisfactory agreement, if the physical volume exceeds about (1.2 fm)^{4}. For the unfolded level spacing distribution we find an accurate agreement with the random matrix conjecture on all volumes that we considered.


title: The QCDOC supercomputer: hardware, software, and performance

abstract: An overview is given of the QCDOC architecture, a massively parallel and highly scalable computer optimized for lattice QCD using system-on-a-chip technology. The heart of a single node is the PowerPC-based QCDOC ASIC, developed in collaboration with IBM Research, with a peak speed of 1 GFlop/s. The nodes communicate via high-speed serial links in a 6-dimensional mesh with nearest-neighbor connections. We find that highly optimized four-dimensional QCD code obtains over 50% efficiency in cycle accurate simulations of QCDOC, even for problems of fixed computational difficulty run on tens of thousands of nodes. We also provide an overview of the QCDOC operating system, which manages and runs QCDOC applications on partitions of variable dimensionality. Finally, the SciDAC activity for QCDOC and the message-passing interface QMP specified as a part of the SciDAC effort are discussed for QCDOC. We explain how to make optimal use of QMP routines on QCDOC in conjunction with existing C and C++ lattice QCD codes, including the publicly available MILC codes.


title: The Early Days of Lattice Gauge Theory

abstract: I discuss some of the historical circumstances that drove us to use the lattice as a non-perturbative regulator. This approach has had immense success, convincingly demonstrating quark confinement and obtaining crucial properties of the strong interactions from first principles. I wrap up with some challenges for the future.


title: A note on Neuberger's double pass algorithm

abstract: We analyze Neuberger's double pass algorithm for the matrix-vector multiplication R(H).Y (where R(H) is (n-1,n)-th degree rational polynomial of positive definite operator H), and show that the number of floating point operations is independent of the degree n, provided that the number of sites is much larger than the number of iterations in the conjugate gradient. This implies that the matrix-vector product $ (H)^{-1/2} Y \simeq R^{(n-1,n)}(H) \cdot Y $ can be approximated to very high precision with sufficiently large n, without noticeably extra costs. Further, we show that there exists a threshold $ n_T $ such that the double pass is faster than the single pass for $ n > n_T $, where $ n_T \simeq 12 - 25 $ for most platforms.


title: Pseudoscalar Decay Constants in Staggered Chiral Perturbation Theory

abstract: In a continuation of an ongoing program, we use staggered chiral perturbation theory to calculate the one-loop chiral logarithms and analytic terms in the pseudoscalar meson leptonic decay constants, $f_{\pi^+_5}$ and $f_{K^+_5}$. We consider the partially quenched, ``full QCD'' (with three dynamical flavors), and quenched cases.


title: Chiral perturbation theory with Wilson-type fermions including $a^2$ effects: $N_f=2$ degenerate case

abstract: We have derived the quark mass dependence of $m_{\pi}^2$, $m_{\rm AWI}$ and $f_{\pi}$, using the chiral perturbation theory which includes the $a^2$ effect associated with the explicit chiral symmetry breaking of the Wilson-type fermions, in the case of the $N_f=2$ degenerate quarks. Distinct features of the results are (1) the additive renormalization for the mass parameter $m_q$ in the Lagrangian, (2) $O(a)$ corrections to the chiral log ($m_q\log m_q$) term, (3) the existence of more singular term, $\log m_q$, generated by $a^2$ contributions, and (4) the existence of both $m_q\log m_q$ and $\log m_q$ terms in the quark mass from the axial Ward-Takahashi identity, $m_{\rm AWI}$. By fitting the mass dependence of $m_\pi^2$ and $m_{\rm AWI}$, obtained by the CP-PACS collaboration for $N_f=2$ full QCD simulations, we have found that the data are consistently described by the derived formulae. Resumming the most singular terms $\log m_q$, we have also derived the modified formulae, which show a better control over the next-to-leading order correction.


title: Analytic estimates for penguin operators in quenched QCD

abstract: Strong penguin operators are singlets under the right-handed flavor symmetry group SU(3)_R. However, they do not remain singlets when the operator is embedded in (partially) quenched QCD, but instead they become linear combinations of two operators with different transformation properties under the (partially) quenched symmetry group. This is an artifact of the quenched approximation. Each of these two operators is represented by a different set of low-energy constants in the chiral effective theory. In this paper, we give analytic estimates for the leading low-energy constants, in quenched and partially quenched QCD. We conclude that the effects of quenching on Q_6 are large.


title: Three-dimensional Gonihedric Potts model

abstract: We study, by the Mean Field and Monte Carlo methods, a generalized q-state Potts gonihedric model. The phase transition of the model becomes stronger with increasing $q.$ The value $k_c(q),$ at which the phase transition becomes second order, turns out to be an increasing function of $q.$


title: 10 D Euclidean dynamical triangulations

abstract: We investigate numerically 10 - dimensional Euclidean quantum gravity (with discretized Einstein - Hilbert action) in the framework of the dynamical triangulation approach. For the considered values of the gravitational coupling we observed two phases, the behavior of which is found to be similar to that of the crumpled and elongated phases of 3, 4 and 5 dimensional models. Surprisingly, (for the observed lattice sizes) the natural state of the 10 D system (when the Einstein - Hilbert action is turned off) is found to resemble branched polymer while in the low dimensional systems the natural state belongs to the crumpled phase.


title: Lattice QCD at Finite Density -- An introductory review

abstract: This is a pedagogical review of the lattice study of finite density QCD, which is intended to provide the minimum necessary contents, so that the paper may be used as the first reading for a newcomer to the field and also for those working in nonlattice communities. After a brief introduction to argue why finite density QCD can be a new attractive subject, we describe fundamental formulae which are necessary for the following sections. Then we survey lattice QCD simulations in small chemical potential regions, where several prominent works have been reported recently. Next, two-color QCD calculations are discussed, where we have a chance to glance at many new features of finite density QCD, and indeed recent simulations indicated quark pair condensation and the in-medium effect. Tables of SU(3) and SU(2) lattice simulations at finite baryon density are given. In the next section, we make a survey of several related works which may be a starting point of the new development in the future, although some works do not attract much attention now. Materials are described in a pedagogical manner. Starting from a simple 2-d model, we briefly discuss a lattice analysis of NJL model. We describe a non-perturbative analytical approach, i.e., strong coupling approximation method and some results. Canonical ensemble approach instead of usual grand canonical ensample may be another route to reach high density. We examine the density of state method and show this old idea includes recently proposed factorization method. An alternative method, complex Langevin equation and an interesting model, finite isospin model, are also discussed. In Appendix, we give several technical points which are useful in practical calculations.


title: Effective lattice theories for Polyakov loops

abstract: We derive effective actions for SU(2) Polyakov loops using inverse Monte Carlo techniques. In a first approach, we determine the effective couplings by requiring that the effective ensemble reproduces the single-site distribution of the Polyakov loops. The latter is flat below the critical temperature implying that the (untraced) Polyakov loop is distributed uniformly over its target space, the SU(2) group manifold. This allows for an analytic determination of the Binder cumulant and the distribution of the mean-field, which turns out to be approximately Gaussian. In a second approach, we employ novel lattice Schwinger-Dyson equations which reflect the SU(2) x SU(2) invariance of the functional Haar measure. Expanding the effective action in terms of SU(2) group characters makes the numerics sufficiently stable so that we are able to extract a total number of 14 couplings. The resulting action is short-ranged and reproduces the Yang-Mills correlators very well.


title: Fermions in 2D Lorentzian Quantum Gravity

abstract: We implement Wilson fermions on 2D Lorentzian triangulation and determine the spectrum of the Dirac-Wilson operator. We compare it to the spectrum of the corresponding operator in the Euclidean background. We use fermionic particle to probe the fractal properties of Lorentzian gravity coupled to $c=1/2$ and $c=4$ matter. We numerically determine the scaling exponent of the mass gap $M \sim N^{-1/d_H}$ to be $d_H=2.11(5)$, and $d_H=1.77(3)$ for $c=1/2$ and $c=4$, respectively


title: Kosterlitz Thouless Universality in Dimer Models

abstract: Using the monomer-dimer representation of strongly coupled U(N) lattice gauge theories with staggered fermions, we study finite temperature chiral phase transitions in (2+1) dimensions. A new cluster algorithm allows us to compute monomer-monomer and dimer-dimer correlations at zero monomer density (chiral limit) accurately on large lattices. This makes it possible to show convincingly, for the first time, that these models undergo a finite temperature phase transition which belongs to the Kosterlitz-Thouless universality class. We find that this universality class is unaffected even in the large N limit. This shows that the mean field analysis often used in this limit breaks down in the critical region.


title: Lattice extraction of $ K \to \pi \pi $ amplitudes to NLO in partially quenched and in full chiral perturbation theory

abstract: We show that it is possible to construct $\epsilon^\prime/\epsilon$ to NLO using partially quenched chiral perturbation theory (PQChPT) from amplitudes that are computable on the lattice. We demonstrate that none of the needed amplitudes require three-momentum on the lattice for either the full theory or the partially quenched theory; non-degenerate quark masses suffice. Furthermore, we find that the electro-weak penguin ($\Delta I=3/2$ and 1/2) contributions to $\epsilon^\prime/\epsilon$ in PQChPT can be determined to NLO using only degenerate ($m_K=m_\pi$) $K\to\pi$ computations without momentum insertion. Issues pertaining to power divergent contributions, originating from mixing with lower dimensional operators, are addressed. Direct calculations of $K\to\pi\pi$ at unphysical kinematics are plagued with enhanced finite volume effects in the (partially) quenched theory, but in simulations when the sea quark mass is equal to the up and down quark mass the enhanced finite volume effects vanish to NLO in PQChPT. In embedding the QCD penguin left-right operator onto PQChPT an ambiguity arises, as first emphasized by Golterman and Pallante. With one version (the "PQS") of the QCD penguin, the inputs needed from the lattice for constructing $K\to\pi\pi$ at NLO in PQChPT coincide with those needed for the full theory. Explicit expressions for the finite logarithms emerging from our NLO analysis to the above amplitudes are also given.


title: Nucleons Properties at Finite Lattice Spacing in Chiral Perturbation Theory

abstract: Properties of the proton and neutron are studied in partially-quenched chiral perturbation theory at finite lattice spacing. Masses, magnetic moments, the matrix elements of isovector twist-2 operators and axial-vector currents are examined at the one-loop level in a double expansion in the light-quark masses and the lattice spacing. This work will be useful in extrapolating the results of simulations using Wilson valence and sea quarks, as well as simulations using Wilson sea quarks and Ginsparg-Wilson valence quarks, to the continuum.


title: Scalar Mesons as qbar^2 q^2? Insight from the Lattice

abstract: I describe some insight obtained from a lattice calculation on the possibility that the light scalar mesons are \qbar^2 q^2 states rather than \qbar q. First I review some general features of \qbar^2 q^2 states in QCD inspired quark models. Then I describe a lattice QCD calculation of pseudoscalar meson scattering amplitudes, ignoring quark loops and quark annihilation, which finds indications that for sufficiently heavy quarks there is a stable four-quark bound state with JPC=0{++} and non-exotic flavor quantum numbers.


title: A Review Of Two Novel Numerical Methods in QFT

abstract: We outline two alternative schemes to perform numerical calculations in quantum field theory. In principle, both of these approaches are better suited to study phase structure than conventional Monte Carlo. The first method, Source Galerkin, is based on a numerical analysis of the Schwinger-Dyson equations using modern computer techniques. The nature of this approach makes dealing with fermions relatively straightforward, particularly since we can work on the continuum. Its ultimate success in non-trivial dimensions will depend on the power of a propagator expansion scheme which also greatly simplifies numerical calculation of traditional perturbation graphs. The second method extends Monte Carlo approaches by introducing a procedure to deal with rapidly oscillating integrals.


title: The spectrum of D_s mesons from lattice QCD

abstract: The spectrum of orbitally excited $D_s$ mesons is computed in the continuum limit of quenched lattice QCD. The results are consistent with the interpretation that the narrow resonance in the $D_s \pi^0$ channel discovered by the BABAR Collaboration is a $J^P=0^+$ $c\bar{s}$ meson. Furthermore, within statistical errors, the $1^+-1^-$ and the $0^+-0^-$ mass splittings are equal, in agreement with the chiral multiplet structure predicted by heavy hadron chiral effective theory. On our coarsest lattice we present results from the first study of orbitally excited $D_s$ mesons with two flavors of dynamical quarks, with mass slightly larger than the strange quark mass. These results are consistent with the quenched data.


title: Quark mass dependence of masses and decay constants of the pseudo-Goldstone bosons in QCD

abstract: The dependence of the pseudoscalar meson masses and decay constants on sea and valence quark masses is compared to next-to-leading order (NLO) Chiral Perturbation Theory (ChPT). The numerical simulations with two light dynamical quark flavors are performed with the Wilson-quark lattice action at gauge coupling beta=5.1 and hopping parameters kappa=0.176, 0.1765, 0.177 on a 16^4 lattice. O(a) lattice artifacts are taken into account by applying chiral perturbation theory for the Wilson lattice action. The values of the relevant combinations of Gasser-Leutwyler constants L_4, L_5, L_6 and L_8 are estimated.


title: Charge Radii of the Meson and Baryon Octets in Quenched and Partially Quenched Chiral Perturbation Theory

abstract: We calculate the electric charge radii of the SU(3) pseudoscalar mesons and the SU(3) octet baryons in quenched and partially quenched chiral perturbation theory. We work in the isospin limit, up to next-to-leading order in the chiral expansion, and to leading order in the heavy baryon expansion. The results are necessary for the extrapolation of future lattice calculations of meson and baryon charge radii. We also derive expressions for the nucleon and pion charge radii in SU(2) flavor away from the isospin limit.


title: Quenched charmonium spectrum

abstract: We study charmonium using the standard relativistic formalism in the quenched approximation, on a set of lattices with isotropic lattice spacings ranging from 0.1 to 0.04 fm. We concentrate on the calculation of the hyperfine splitting between eta_c and J/psi, aiming for a controlled continuum extrapolation of this quantity. The splitting extracted from the non-perturbatively improved clover Dirac operator shows very little dependence on the lattice spacing for $a \leq 0.1$ fm. The dependence is much stronger for Wilson and tree-level improved clover operators, but they still yield consistent extrapolations if sufficiently fine lattices, $a \leq 0.07$ fm ($a M(\eta_c) \leq 1$), are used. Our result for the hyperfine splitting is 77(2)(6) MeV (where Sommer's parameter, r_0, is used to fix the scale). This value remains about 30% below experiment. Dynamical fermions and OZI-forbidden diagrams both contribute to the remainder. Results for the eta_c and J/psi wave functions are also presented.


title: Heavy--light decay constants in the continuum limit of quenched lattice QCD

abstract: We compute the decay constants for the heavy--light pseudoscalar mesons in the quenched approximation and continuum limit of lattice QCD. Within the Schrodinger Functional framework, we make use of the step scaling method, which has been previously introduced in order to deal with the two scale problem represented by the coexistence of a light and a heavy quark. The continuum extrapolation gives us a value $f_{B_s} = 192(6)(4)$ MeV for the $B_s$ meson decay constant and $f_{D_s} = 240(5)(5)$ MeV for the $D_s$ meson.


title: Lattice Results on the Connected Neutron Charge Radius

abstract: We describe a calculation using quenched lattice QCD of the connected part of the neutron electric form factor for momentum transfers in the range $0.3 {\rm GeV^2} \stackrel{<}{\sim} -q^2 \stackrel{<}{\sim} 1.0 {\rm GeV^2}$. We extract the implied charge radius using a Galster parameterization and consider various ways of extrapolating the neutron charge radius value to the chiral limit. We find that the measured charge radii may be reconciled to experiment by standard phenomenology and lowest or next to lowest order contributions from chiral perturbation theory.


title: Analyticity and the phase diagram of QCD

abstract: Some consequences of the analyticity of the free energy (pressure) of QCD at finite chemical potential are deduced. These include a method for numerical exploration of the full phase diagram by a novel use of simulations at imaginary chemical potential to extract Yang-Lee zeroes of the grand-canonical partition function. We make use of, and comment on, CPT symmetries, positivity of non-linear susceptibilities and the finiteness of screening lengths. We also comment on the structure of zeroes expected for the usual picture of the phases of QCD, following a discussion of the physics of imaginary chemical potential.


title: Simulating the Schroedinger functional with two pseudo-fermions

abstract: We report on simulations with two flavors of O(a) improved degenerate Wilson fermions with Schroedinger functional boundary conditions. The algorithm which is used is Hybrid Monte Carlo with two pseudo-fermion fields as proposed by M. Hasenbusch. We investigate the numerical precision and sensitivity to reversibility violations of this algorithm. A gain of a factor two in CPU cost is reached compared with one pseudo-fermion field due to the larger possible step-size.


title: The Matter and the Pseudoscalar Densities in Lattice QCD

abstract: The matter and the pseudoscalar densities inside a hadron are calculated via gauge-invariant equal-time correlation functions. A comparison is made between the charge charge and the matter density distributions for the pion, the rho, the nucleon and the $\Delta^+$ within the quenched theory, and with two flavours of dynamical quarks.


title: The contribution of O(alpha) radiative corrections to the renormalised anisotropy and application to general tadpole improvement schemes: addendum to "One loop calculation of the renormalised anisotropy for improved anisotropic gluon actions on a lattice" [hep-lat/0208010]

abstract: General O(alpha) radiative corrections to lattice actions may be interpreted as counterterms that give additive contributions to the one-loop renormalisation of the anisotropy. The effect of changing the radiative coefficients is thus easily calculable. In particular, the results obtained in a previous paper for Landau mean link improved actions apply in any tadpole improvement scheme. We explain how this method can be exploited when tuning radiatively improved actions. Efficient methods for self-consistently tuning tadpole improvement factors are also discussed.


title: General chiral gauge theories on the lattice

abstract: We still extend the large class of Dirac operators decribing massless fermions on the lattice found recently, only requiring that such operators decompose into Weyl operators. After deriving general relations and constructions of operators, we study the basis representations of the chiral projections. We then investigate correlation functions of Weyl fermions for any value of the index, stressing the related conditions for basis transformations and their consequences, and getting the precise behaviors under gauge transformations and CP transformations. Various further developments include considerations of the explicit form of the effective action and of a representation of the general correlation functions in terms of alternating multilinear forms. For comparison we also consider gauge-field variations and their respective applications. Finally we compare with continuum perturbation theory.


title: Supersymmetry on a Euclidean Spacetime Lattice II: Theories with Eight Supercharges

abstract: We formulate Euclidean spacetime lattices whose continuum limits are supersymmetric Yang-Mills theories with eight supercharges in two and three dimensions. The lattice actions are themselves supersymmetric.


title: Quenched spectroscopy with fixed-point and chirally improved fermions

abstract: We present results from quenched spectroscopy calculations with the parametrized fixed-point and the chirally improved Dirac operators. Both these operators are approximate solutions of the Ginsparg-Wilson equation and have good chiral properties. This allows us to work at small quark masses and we explore pseudoscalar-mass to vector-mass ratios down to 0.28. We discuss meson and baryon masses, their scaling properties, finite volume effects and compare our results with recent large scale simulations. We find that the size of quenching artifacts of the masses is strongly correlated with their experimentally observed widths and that the gauge and hadronic scales are consistent.


title: Comment on: ``The Supersymmetric Ward-Takahashi Identity in 1-Loop Lattice Perturbation Theory. I. General Procedure'' by A. Feo; hep-lat/0305020

abstract: We comment on several issues concerning both the physics and authorship of article hep-lat/0305020.


title: On the scaling of computational particle physics codes on cluster computers

abstract: Many appplications in computational science are sufficiently compute-intensive that they depend on the power of parallel computing for viability. For all but the "embarrassingly parallel" problems, the performance depends upon the level of granularity that can be achieved on the computer platform. Our computational particle physics applications require machines that can support a wide range of granularities, but in general, compute-intensive state-of-the-art projects will require finely grained distributions. Of the different types of machines available for the task, we consider cluster computers. The use of clusters of commodity computers in high performance computing has many advantages including the raw price/performance ratio and the flexibility of machine configuration and upgrade. Here we focus on what is usually considered the weak point of cluster technology; the scaling behaviour when faced with a numerically intensive parallel computation. To this end we examine the scaling of our own applications from numerical quantum field theory on a cluster and infer conclusions about the more general case.


title: Analysis of the Axial Anomaly on the Lattice with O(a)-improved Wilson Action

abstract: Flavor singlet and non-singlet axial Ward identities are investigated using the Wilson formulation of lattice QCD with Clover O(a)-improvement, which breaks explicitly chiral symmetry. The matching at one-loop order of all the relevant renormalization constants with the continuum \bar{MS} scheme is presented. Our calculations include: 1) the contributions arising from the Clover term of the action; 2) the complete mixing of the gluon operator G \tilde{G} with the divergence of the singlet axial current; 3) the use of both local and extended definitions of the fermionic bilinear operators. A definition of the gluon operator G \tilde{G} on the lattice outside the chiral limit is proposed. Our definition takes into account the possible power-divergent mixing with the pseudoscalar density, generated by the breaking of chiral symmetry. A non-perturbative procedure for the evaluation of such mixing constant is worked out. Finally, the renormalization properties of the composite insertion of the topological charge operator relevant for the lattice calculation of the neutron electric dipole moment, induced by the strong CP-violating term of the QCD Lagrangian, are discussed.


title: The high temperature phase transition in SU(N) gauge theories

abstract: We calculate the continuum value of the deconfining temperature in units of the string tension for SU(4), SU(6) and SU(8) gauge theories, and we recalculate its value for SU(2) and SU(3). We find that the $N$-dependence for $2 \leq N \leq 8$ is well fitted by $T_c/\sqrt{sigma} = 0.596(4) + 0.453(30)/N^2$, showing a rapid convergence to the large-N limit. We confirm our earlier result that the phase transition is first order for $N \geq 3$ and that it becomes stronger with increasing $N$. We also confirm that as $N$ increases the finite volume corrections become rapidly smaller and the phase transition becomes visible on ever smaller volumes. We interpret the latter as being due to the fact that the tension of the domain wall that separates the confining and deconfining phases increases rapidly with $N$. We speculate on the connection to Eguchi-Kawai reduction and to the idea of a Master Field.


title: N to Delta electromagnetic transition form factors from Lattice QCD

abstract: The magnetic dipole, the electric quadrupole and the Coulomb quadrupole amplitudes for the transition $\gamma N\to \Delta$ are evaluated both in quenched lattice QCD at $\beta=6.0$ and using two dynamical Wilson fermions simulated at $\beta=5.6$. The dipole transition form factor is accurately determined at several values of momentum transfer. On the lattices studied in this work, the electric quadrupole amplitude is found to be non-zero yielding a negative value for the ratio, $ R_{EM}$, of electric quadrupole to magnetic dipole amplitudes at three values of momentum transfer.


title: Diquark condensation at strong coupling

abstract: The possibility of diquark condensation at sufficiently large baryon chemical potential and zero temperature is analyzed in QCD at strong coupling. In agreement with other strong coupling analysis, it is found that a first order phase transition separates a low density phase with chiral symmetry spontaneously broken from a high density phase where chiral symmetry is restored. In none of the phases diquark condensation takes place as an equilibrium state, but, for any value of the chemical potential, there is a metastable state characterized by a non-vanishing diquark condensate. The energy difference between this metastable state and the equilibrium state decreases with the chemical potential and is minimum in the high density phase. The results indicate that there is attraction in the quark-quark sector also at strong coupling, and that the attraction is more effective at high baryon density, but for infinite coupling it is not enough to produce diquark condensation. It is argued that the absence of diquark condensation is not a peculiarity of the strong coupling limit, but persists at sufficiently large finite couplings.


title: The QCD Phase Diagram for Three Degenerate Flavors and Small Baryon Density

abstract: We present results for the phase diagram of three flavor QCD for \mu_B ~ 500 MeV. Our simulations are performed with imaginary chemical potential \mu_I for which the fermion determinant is positive. Physical observables are then fitted by truncated Taylor series and continued to real chemical potential. We map out the location of the critical line T_c(\mu_B) with an accuracy up to terms of order (\mu_B/T)^6. We also give first results on a determination of the critical endpoint of the transition and its quark mass dependence. Our results for the endpoint differ significantly from those obtained by other methods, and we discuss possible reasons for this.


title: Lattice HQET with exponentially improved statistical precision

abstract: We introduce an alternative discretization for static quarks on the lattice retaining the O(a) improvement properties of the Eichten-Hill action. In this formulation, statistical fluctuations are reduced by a factor which grows exponentially with Euclidean time, x_0. For the first time, B-meson correlation functions are computed with good statistical precision in the static approximation for x_0>1 fm. At lattice spacings a \approx 0.1 fm, a \approx 0.08 fm and a \approx 0.07 fm the B_s-meson decay constant is determined in static and quenched approximations. A correction due to the finite mass of the b-quark is estimated by combining these static results with a recent determination of F_Ds.


title: Axial Anomaly and Index of the Overlap Hypercube Operator

abstract: The overlap hypercube fermion is constructed by inserting a lattice fermion with hypercubic couplings into the overlap formula. One obtains an exact Ginsparg-Wilson fermion, which is more complicated than the standard overlap fermion, but which has improved practical properties and is of current interest for use in numerical simulations. Here we deal with conceptual aspects of the overlap hypercube Dirac operator. Specifically, we evaluate the axial anomaly and the index, demonstrating that the correct classical continuum limit is recovered. Our derivation is non-perturbative and therefore valid in all topological sectors. At the non-perturbative level this result had previously only been shown for the standard overlap Dirac operator with Wilson kernel. The new techniques which we develop to accomplish this are of a general nature and have the potential to be extended to overlap Dirac operators with even more general kernels.


title: Chiral Lagrangian Parameters for Scalar and Pseudoscalar Mesons

abstract: The results of a high-statistics study of scalar and pseudoscalar meson propagators in quenched lattice QCD are presented. For two values of lattice spacing, $\beta=5.7$ ($a \approx .18$ fm) and 5.9 ($a \approx .12$ fm), we probe the light quark mass region using clover improved Wilson fermions with the MQA pole-shifting ansatz to treat the exceptional configuration problem. The quenched chiral loop parameters $m_0$ and $\alpha_{\Phi}$ are determined from a study of the pseudoscalar hairpin correlator. From a global fit to the meson correlators, estimates are obtained for the relevant chiral Lagrangian parameters, including the Leutwyler parameters $L_5$ and $L_8$. Using the parameters obtained from the singlet and nonsinglet pseudoscalar correlators, the quenched chiral loop effect in the nonsinglet scalar meson correlator is studied. By removing this QCL effect from the lattice correlator, we obtain the mass and decay constant of the ground state scalar, isovector meson $a_0$.


title: The fermion determinant in (4,4) 2d lattice super-Yang-Mills

abstract: We find that the fermion determinant is not generally positive in a class of lattice actions recently constructed by Cohen et al. [hep-lat/0307012]; these are actions that contain an exact lattice supersymmetry and have as their target (continuum) theory (4,4) 2-dimensional super-Yang-Mills. We discuss the implications of this finding for lattice simulations and give some preliminary results for the phase of the determinant in the phase-quenched ensemble.


title: Spectroscopy of Doubly-Charmed Baryons in Lattice QCD

abstract: We present results for masses of spin-1/2 and spin-3/2 double-charm baryons in quenched lattice QCD, from an exploratory study using a non-perturbatively improved clover action at beta=6.2. We have studied local operators and we observe, after appropriate projections, a good signal for the ground states. We also present results for single-charmed baryons and spin-splittings for both double- and single-charmed states.


title: Quark propagator and vertex: systematic corrections of hypercubic artifacts from lattice simulations

abstract: This is the first part of a study of the quark propagator and the vertex function of the vector current on the lattice in the Landau gauge and using both Wilson-clover and overlap actions. In order to be able to identify lattice artifacts and to reach large momenta we use a range of lattice spacings. The lattice artifacts turn out to be exceedingly large in this study. We present a new and very efficient method to eliminate the hypercubic (anisotropy) artifacts based on a systematic expansion on hypercubic invariants which are not SO(4) invariant. A simpler version of this method has been used in previous works. This method is shown to be significantly more efficient than the popular ``democratic'' methods. It can of course be applied to the lattice simulations of many other physical quantities. The analysis indicates a hierarchy in the size of hypercubic artifacts: overlap larger than clover and propagator larger than vertex function. This pleads for the combined study of propagators and vertex functions via Ward identities.


title: Meson and Baryon Spectroscopy on a Lattice

abstract: I review the results of hadron spectroscopy calculations from lattice QCD for an intended audience of low energy hadronic physicists. I briefly introduce the ideas of numerical lattice QCD. The various systematic errors, such as the lattice spacing and volume dependence, in lattice QCD calculations are discussed. In addition to the discussion of the properties of ground state hadrons, I also review the small amount of work done on the spectroscopy of excited hadrons and the effect of electromagnetic fields on hadron masses. I also discuss the attempts to understand the physical mechanisms behind hadron mass splittings.


title: Chaotic Quantization of Four-Dimensional U(1) Lattice Gauge Theory

abstract: We demonstrate that the quantized U(1) lattice gauge theory in four Euclidean dimensions can be obtained as the long time limit of the corresponding classical U(1) gauge theory in 4+1 dimensions. The Planck constant hbar is related to the excitation energy and the lattice constant of this classical template.


title: Hamiltonian Study of Improved $U(1$ Lattice Gauge Theory in Three Dimensions

abstract: A comprehensive analysis of the Symanzik improved anisotropic three-dimensional U(1) lattice gauge theory in the Hamiltonian limit is made. Monte Carlo techniques are used to obtain numerical results for the static potential, ratio of the renormalized and bare anisotropies, the string tension, lowest glueball masses and the mass ratio. Evidence that rotational symmetry is established more accurately for the Symanzik improved anisotropic action is presented. The discretization errors in the static potential and the renormalization of the bare anisotropy are found to be only a few percent compared to errors of about 20-25% for the unimproved gauge action. Evidence of scaling in the string tension, antisymmetric mass gap and the mass ratio is observed in the weak coupling region and the behaviour is tested against analytic and numerical results obtained in various other Hamiltonian studies of the theory. We find that more accurate determination of the scaling coefficients of the string tension and the antisymmetric mass gap has been achieved, and the agreement with various other Hamiltonian studies of the theory is excellent. The improved action is found to give faster convergence to the continuum limit. Very clear evidence is obtained that in the continuum limit the glueball ratio $M_{S}/M_{A}$ approaches exactly 2, as expected in a theory of free, massive bosons.


title: Vortex structures in pure SU(3) lattice gauge theory

abstract: The structures of confining vortices which underlie pure SU(3) Yang-Mills theory are studied by means of lattice gauge theory. Vortices and Z_3 monopoles are defined as dynamical degrees of freedom of the Z_3 gauge theory which emerges by center gauge fixing and by subsequent center projection. It is observed for the first time for the case of SU(3) that these degrees of freedom are sensible in the continuum limit: the planar vortex density and the monopole density properly scales with the lattice spacing. By contrast to earlier findings concerning the gauge group SU(2), the effective vortex theory only reproduces 62% of the full string tension. On the other hand, however, the removal of the vortices from the lattice configurations yields ensembles with vanishing string tension. SU(3) vortex matter which originates from Laplacian center gauge fixing is also discussed. Although these vortices recover the full string tension, they lack a direct interpretation as physical degrees of freedom in the continuum limit.


title: Wess-Zumino-Witten term on the lattice

abstract: We construct the Wess-Zumino-Witten (WZW) term in lattice gauge theory by using a Dirac operator which obeys the Ginsparg-Wilson relation. Topological properties of the WZW term known in the continuum are reproduced on the lattice as a consequence of a non-trivial topological structure of the space of admissible lattice gauge fields. In the course of this analysis, we observe that the gauge anomaly generally implies that there is no basis of a Weyl fermion which leads to a single-valued expectation value in the fermion sector. The lattice Witten term, which carries information of a gauge path along which the gauge anomaly is integrated, is separated from the WZW term and the multivaluedness of the Witten term is shown to be related to the homotopy group $\pi_{2n+1}(G)$. We also discuss the global $\SU(2)$ anomaly on the basis of the WZW term.


title: Correlation lengths and scaling functions in the three-dimensional O(4) model

abstract: We investigate numerically the transverse and longitudinal correlation lengths of the three-dimensional O(4) model as a function of the external field H. From our data we calculate the scaling function of the transverse correlation length, and that of the longitudinal correlation length for T>T_c. We show that the scaling functions do not only describe the critical behaviours of the correlation lengths but encompass as well the predicted Goldstone effects, in particular the H^{-1/2}-dependence of the transverse correlation length for T<T_c. In addition, we determine the critical exponent delta=4.824(9) and several critical amplitudes from which we derive the universal amplitude ratios R_{chi}=1.084(18), Q_c=0.431(9), Q_2^T=4.91(8), Q_2^L=1.265(24) and U_{xi}^c=1.99(1). The last result supports a relation between the longitudinal and transverse correlation functions, which was conjectured to hold below T_c but seems to be valid also at T_c.


title: Electromagnetic Properties of the Baryon Decuplet in Quenched and Partially Quenched Chiral Perturbation Theory

abstract: We calculate electromagnetic properties of the decuplet baryons in quenched and partially quenched chiral perturbation theory. We work at next-to-leading order in the chiral expansion, leading order in the heavy baryon expansion, and obtain expressions for the magnetic moments, charge radii, and electric quadrupole moments. The quenched calculation is shown to be pathological since only quenched chiral singularities are present at this order. We present the partially quenched results for both the SU(2) and SU(3) flavor groups and use the isospin limit in the latter. These results are necessary for proper extrapolation of lattice calculations of decuplet electromagnetic properties.


title: An Abelian effective action reproducing screening and confinement in quenched SU(2) QCD

abstract: In an Abelian projection SU(2) gluodynamics contains Abelian gauge fields (diagonal degrees of freedom) and Abelian matter fields (off-diagonal degrees). The matter fields are essential for the breaking of the adjoint string. We obtain numerically the effective action of the Abelian fields in quenched SU(2) QCD and show that the Abelian matter fields provide an essential contribution to the total action even in the infrared region.


title: QCD Dirac operator at nonzero chemical potential: lattice data and matrix model

abstract: Recently, a non-Hermitian chiral random matrix model was proposed to describe the eigenvalues of the QCD Dirac operator at nonzero chemical potential. This matrix model can be constructed from QCD by mapping it to an equivalent matrix model which has the same symmetries as QCD with chemical potential. Its microscopic spectral correlations are conjectured to be identical to those of the QCD Dirac operator. We investigate this conjecture by comparing large ensembles of Dirac eigenvalues in quenched SU(3) lattice QCD at nonzero chemical potential to the analytical predictions of the matrix model. Excellent agreement is found in the two regimes of weak and strong non-Hermiticity, for several different lattice volumes.


title: Entropy of monopoles from percolating cluster in quenched SU(2) QCD

abstract: The length distribution and the monopole action of the infrared monopole clusters are studied numerically in quenched SU(2) QCD. We determine the effective entropy of the monopole currents which turns out to be a descending function of the blocking scale, indicating that the effective degrees of freedom of the extended monopoles are getting smaller as the blocking scale increases.


title: Parallel implementation of a lattice-gauge-theory code: studying quark confinement on PC clusters

abstract: We consider the implementation of a parallel Monte Carlo code for high-performance simulations on PC clusters with MPI. We carry out tests of speedup and efficiency. The code is used for numerical simulations of pure SU(2) lattice gauge theory at very large lattice volumes, in order to study the infrared behavior of gluon and ghost propagators. This problem is directly related to the confinement of quarks and gluons in the physics of strong interactions.


title: On the continuum limit of gauge-fixed compact U(1) lattice gauge theory

abstract: We investigate the continuum limit of a compact formulation of the lattice U(1) gauge theory in 4 dimensions using a nonperturbative gauge-fixed regularization. We find clear evidence of a continuous phase transition in the pure gauge theory for all values of the gauge coupling (with gauge symmetry restored). When probed with quenched staggered fermions with U(1) charge, the theory clearly has a chiral transition for large gauge couplings. We identify the only possible region in the parameter space where a continuum limit with nonperturbative physics may appear.


title: Blocked lattice monopoles in quenched SU(2) QCD and dual superconductor model

abstract: We study the action of the lattice monopoles in quenched SU(2) QCD in the Maximal Abelian projection. We relate the lattice action of the monopole currents to the monopole degrees of freedom of the continuum dual superconductor model and obtain the value of the monopole condensate.


title: A detailed study of the Abelian-projected SU(2) flux tube and its dual Ginzburg-Landau analysis

abstract: The color-electric flux tube of Abelian-projected (AP) SU(2) lattice gauge theory in the maximally Abelian gauge (MAG) is revisited. It is shown that the lattice Gribov copy effect in the MAG is crucial for the monopole-related parts of the flux-tube profiles. Taking into account both the gauge fixing procedure and the effect of finite quark-antiquark distance properly, the scaling property of the flux-tube profile is confirmed. The quantitative relation between the measured AP flux tube and the flux-tube solution of the U(1) dual Abelian Higgs (DAH) model is also discussed. The fitting of the AP flux tube in terms of the DAH flux tube indicates that the vacuum can be classified as weakly type-I dual superconductor.


title: Mass dependence of the hairpin vertex in quenched QCD

abstract: The pseudoscalar ``hairpin'' vertex (i.e. quark-disconnected vertex) plays a key role in quenched chiral perturbation theory. Direct calculations using lattice simulations find that it has a significant dependence on quark mass. I show that this mass dependence can be used to determine the quenched Gasser-Leutwyler constant L5. This complements the calculation of L5 using the mass dependence of the axial decay constant of the pion. In an appendix, I discuss power counting for quenched chiral perturbation theory and describe the particular scheme used in this paper.


title: Extending the Utility of Partially Quenched QCD

abstract: It has been proposed that partially quenched chiral perturbation theory can be used together with partially quenched lattice QCD to determine the low-energy constants of real QCD. We point out a complication of this approach and show how it can be resolved. Partial quenching changes the chiral symmetry group and introduces a new independent four-derivative operator. We explore the effects of this new operator. We describe tree-level scattering processes to which it contributes, and show how they can be used, in principle, to deterimine its unknown coefficient. We also calculate its lowest-order (one-loop) contribution to charged meson masses and decay constants. In addition, we present the form of the analytic next-to-next-to-leading order meson mass and decay constant corrections, which are needed for lattice fits.


title: Phase structure of (2+1)d strongly coupled lattice gauge theories

abstract: We study the chiral phase transition in (2+1)d strongly coupled U(N) lattice gauge theories with staggered fermions. We show with high precision simulations performed directly in the chiral limit that these models undergo a Berezinski-Kosterlitz-Thouless (BKT) transition. We also show that this universality class is unaffected even in the large N limit.