title: Improvement and Analytic Techniques in Hamiltonian Lattice Gauge Theory

abstract: This thesis is concerned with two topics in Hamiltonian lattice gauge theory: improvement and the application of analytic techniques. On the topic of improvement, we develop a direct method for improving lattice Hamiltonians for gluons, in which linear combinations of gauge invariant lattice operators are chosen to cancel the lowest order discretisation errors. On the topic of analytic methods, we extend the techniques that have been used in 2+1 dimensional SU(2) variational calculations for many years, to the general case of SU(N). For this purpose a number of group integrals are calculated in terms of Toeplitz determinants. As generating functions these group integrals allow the calculation of all matrix elements appearing in variational glueball mass calculations in 2+1 dimensions. Making use of these analytic techniques, glueball masses in various symmetry sectors are calculated with N as large as 25 allowing extrapolations to infinite N to be performed. We finish with a feasibility study of applying analytic Hamiltonian techniques in 3+1 dimensional calculations.


title: Continuum limit of proton decay matrix elements in quenched lattice QCD

abstract: We present a lattice QCD calculation of the parameters \alpha and \beta which are necessary in the theoretical estimation of the proton lifetime in grand unified theories (GUTs) using chiral lagrangian approach. The simulation is carried out using the Wilson quark action at three gauge coupling constants in the quenched approximation. We obtain |\alpha(2GeV)|=0.0091(08)(^{+10}_{-19})GeV^3 and |\beta(2GeV)|=0.0098(08)(^{+10}_{-20})GeV^3 in the continuum limit where the first error is statistical and the second one is due to scale setting.


title: Lattice Extraction of $\epsilon^\prime/\epsilon$ using (PQ) Chiral Perturbation Theory

abstract: Results at NLO in (PQ)ChPT for extracting $\epsilon^\prime/\epsilon$ from lattice data are discussed.


title: Study of finite volume effects in the non-perturbative determination of $\csw$ with the SF method in full three-flavor lattice QCD

abstract: The non-perturbative $\csw$ determined by the Schr\"{o}dinger functional (SF) method with the RG-improved gauge action in dynamical $N_f=3$ QCD shows a finite volume effect when the numerical simulations are carried out at a constant lattice size $L/a$. We remove the unwanted finite volume effect by keeping physical lattice extent $L$ at a constant. The details of the method and the result obtained for non-perturbative $\csw$ with a constant $L$ are reported.


title: Gluon and ghost propagators from the viewpoint of general principles of quantized gauge field theories

abstract: We discuss the possible form of gluon and ghost propagators in the infrared region of Yang-Mills theory in the covariant gauge from the viewpoint of general principles of quantized gauge field theories.


title: Application of Maximum Entropy Method to Lattice Field Theory with a Topological Term

abstract: In Monte Carlo simulation, lattice field theory with a $\theta$ term suffers from the sign problem. This problem can be circumvented by Fourier-transforming the topological charge distribution $P(Q)$. Although this strategy works well for small lattice volume, effect of errors of $P(Q)$ becomes serious with increasing volume and prevents one from studying the phase structure. This is called flattening. As an alternative approach, we apply the maximum entropy method (MEM) to the Gaussian $P(Q)$. It is found that the flattening could be much improved by use of the MEM.


title: Finite temperature phase transition in full QCD with $N_f=2$ flavors of clover fermions at $N_t=8$ and 10

abstract: We present results for QCD with $N_f=2$ flavors of dynamical quarks using nonperturbatively improved Wilson fermions at finite temperature on $16^3 \times 8$ and $24^3 \times 10$ lattices. We determine the transition temperature in the range of quark masses $0.6 \lesssim m_\pi/m_\rho \leq 0.8$. After fixing the Maximal Abelian gauge we investigate the contribution of Abelian monopoles to the Polyakov loop, Polyakov loop susceptibility and confirm Abelian and monopole dominance in full QCD.


title: Topological susceptibility from the overlap

abstract: The chiral symmetry at finite lattice spacing of Ginsparg-Wilson fermionic actions constrains the renormalization of the lattice operators; in particular, the topological susceptibility does not require any renormalization, when using a fermionic estimator to define the topological charge. Therefore, the overlap formalism appears as an appealing candidate to study the continuum limit of the topological susceptibility while keeping the systematic errors under theoretical control. We present results for the SU(3) pure gauge theory using the index of the overlap Dirac operator to study the topology of the gauge configurations. The topological charge is obtained from the zero modes of the overlap and using a new algorithm for the spectral flow analysis. A detailed comparison with cooling techniques is presented. Particular care is taken in assessing the systematic errors. Relatively high statistics (500 to 1000 independent configurations) yield an extrapolated continuum limit with errors that are comparable with other methods. Our current value from the overlap is $\chi^{1/4} = 188 \pm 12 \pm 5 \MeV$.


title: QCD at zero baryon density

abstract: While the grand canonical partition function Z_{GC}(mu) with chemical potential mu explicitly breaks the Z_3 symmetry with the Dirac determinant, the canonical partition function at fixed baryon number Z_C(B) is manifestly Z_3-symmetric. We compare Z_{GC}(mu=0) and Z_C(B=0) formally and by numerical simulations, in particular with respect to properties of the deconfinement transition. Differences between the two ensembles, for physical observables characterising the phase transition, vanish with increasing lattice size. We show numerically that the free energy density is the same for both ensembles in the thermodynamic limit.


title: Effective string picture for confinement at finite temperature: theoretical predictions and high precision numerical results

abstract: The effective string picture of confinement is used to derive theoretical predictions for the interquark potential at finite temperature. At short distances, the leading string correction to the linear confining potential between a heavy quark-antiquark pair is the "L\"uscher term''. We assume a Nambu--Goto effective string action, and work out subleading contributions in an analytical way. We discuss the contribution given by a possible ``boundary term'' in the effective action, comparing these predictions with results from simulations of lattice $Z_2$ gauge theory in three dimensions, obtained with an algorithm which exploits the duality of the $Z_2$ gauge model with the Ising spin model.


title: Improving the locality of the overlap Dirac operator via approximate solutions of the Ginsparg-Wilson relation

abstract: We determine the free field hypercubic Dirac operator which is optimally close to satisfying the Ginsparg-Wilson relation. Inserting this operator into the overlap formula, we show that the analytic locality bound on the resulting overlap Dirac operator is substantially stronger than in the standard case. This improvement generally persists in gauge backgrounds when the plaquette variables are all close to unity.


title: Benchmarking computer platforms for lattice QCD applications

abstract: We define a benchmark suite for lattice QCD and report on benchmark results from several computer platforms. The platforms considered are apeNEXT, CRAY T3E, Hitachi SR8000, IBM p690, PC-Clusters, and QCDOC.


title: Non perturbative renormalization in coordinate space

abstract: We present an exploratory study of a gauge-invariant non-perturbative renormalization technique. The renormalization conditions are imposed on correlation functions of composite operators in coordinate space on the lattice. Numerical results for bilinears obtained with overlap and O(a)-improved Wilson fermions are presented. The measurement of the quark condensate is also discussed.


title: Hadronic matrix elements of proton decay on the lattice

abstract: We report on our on-going project to calculate proton decay matrix elements using domain-wall fermions on the lattice. By summarizing the history of the proton decay calculation on the lattice, we reveal the systematic errors of those calculations. Then we discuss our approach to tackle those uncertainties and show our preliminary results on the matrix elements.


title: Nucleon Decay with Domain-Wall Fermions

abstract: We report on our on-going project to calculate the nucleon decay matrix elements with domain-wall fermions. Operator mixing is discussed employing a non-perturbative renormalization. Bare matrix elements of all the possible decay modes induced by the dimension-six operators are calculated with the direct method, which are compared with the indirect calculation using chiral perturbation theory.


title: Deconfinement transition in 2+1-dimensional SU(4) lattice gauge theory

abstract: A missing piece is added to the Svetitsky-Yaffe conjecture. The spin model in the same universality class as the (2+1)d SU(4) theory, the 2d Ashkin-Teller model, has a line of continuously varying critical exponents. The exponents measured in the gauge theory correspond best to the Potts point on the Ashkin-Teller line.


title: Dynamical FLIC Fermions

abstract: The use of APE smearing or other blocking techniques in fermion actions can provide many advantages. There are many variants of these fat link actions in lattice QCD currently, such as FLIC fermions. The FLIC fermion formalism makes use of the APE blocking technique in combination with a projection of the blocked links back into the special unitary group. This reunitarisation is often performed using an iterative maximisation of a gauge invariant measure. This technique is not differentiable with respect to the gauge field and thus prevents the use of standard Hybrid Monte Carlo simulation algorithms. The use of an alternative projection technique circumvents this difficulty and allows the simulation of dynamical fat link fermions with standard HMC and its variants.


title: I=2 pion-pion scattering phase shift in the continuum limit calculated with two-flavor full QCD

abstract: We present a calculation of the scattering phase shift for the I=2 S-wave pion-pion system in the continuum limit with two-flavor full QCD. Calculations are made at three lattice spacings, using the finite volume method of L\"uscher in the center of mass frame, and its extension to the laboratory frame.


title: Maximum Entropy Method Approach to $\theta$ Term

abstract: In Monte Carlo simulations of lattice field theory with a $\theta$ term, one confronts the complex weight problem, or the sign problem. This is circumvented by performing the Fourier transform of the topological charge distribution $P(Q)$. This procedure, however, causes flattening phenomenon of the free energy $f(\theta)$, which makes study of the phase structure unfeasible. In order to treat this problem, we apply the maximum entropy method (MEM) to a Gaussian form of $P(Q)$, which serves as a good example to test whether the MEM can be applied effectively to the $\theta$ term. We study the case with flattening as well as that without flattening. In the latter case, the results of the MEM agree with those obtained from the direct application of the Fourier transform. For the former, the MEM gives a smoother $f(\theta)$ than that of the Fourier transform. Among various default models investigated, the images which yield the least error do not show flattening, although some others cannot be excluded given the uncertainty related to statistical error.


title: Chirally improving Wilson fermions

abstract: It is possible to improve the chiral behaviour and the approach to the continuum limit of correlation functions in lattice QCD with standard and twisted Wilson fermions by taking arithmetic averages of correlators computed in theories regularized with Wilson terms of opposite sign. To avoid the problem of the spurious zero-modes of the Wilson--Dirac operator, twisted-mass lattice QCD should be used for the actual computation of the correlators taking part in the averages. A particularly useful choice for the twisting angle is $\pm\pi/2$ where many physical quantities (e.g. masses and zero-momentum matrix elements) are automatically improved with no need of averaging.


title: Higher orders of the high-temperature expansion for the Ising model in three dimensions

abstract: The new algorithm of the finite lattice method is applied to generate the high-temperature expansion series of the simple cubic Ising model to $\beta^{50}$ for the free energy, to $\beta^{32}$ for the magnetic susceptibility and to $\beta^{29}$ for the second moment correlation length. The series are analyzed to give the precise value of the critical point and the critical exponents of the model.


title: Phase structure of lattice QCD for general number of flavors

abstract: We investigate the phase structure of lattice QCD for the general number of flavors in the parameter space of gauge coupling constant and quark mass, employing the one-plaquette gauge action and the standard Wilson quark action. Performing a series of simulations for the number of flavors $N_F=6$--360 with degenerate-mass quarks, we find that when $N_F \ge 7$ there is a line of a bulk first order phase transition between the confined phase and a deconfined phase at a finite current quark mass in the strong coupling region and the intermediate coupling region. The massless quark line exists only in the deconfined phase. Based on these numerical results in the strong coupling limit and in the intermediate coupling region, we propose the following phase structure, depending on the number of flavors whose masses are less than $\Lambda_d$ which is the physical scale characterizing the phase transition in the weak coupling region: When $N_F \ge 17$, there is only a trivial IR fixed point and therefore the theory in the continuum limit is free. On the other hand, when $16 \ge N_F \ge 7$, there is a non-trivial IR fixed point and therefore the theory is non-trivial with anomalous dimensions, however, without quark confinement. Theories which satisfy both quark confinement and spontaneous chiral symmetry breaking in the continuum limit exist only for $N_F \le 6$.


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 perturbatively determine O($a$) boundary improvement coefficients at 1-loop for the Schr\"odinger Functional coupling with improved gauge actions. These coefficients are required to implement the 1-loop O($a$) improvement in full QCD simulations for the coupling with the improved gauge actions. To this order, lattice artifacts of the step scaling function (SSF) for each improved gauge action are also investigated. Furthermore, we investigate the effect of the 1-loop O($a$) improvement to the SSF in numerical simulations of the pure SU(3) gauge theory.


title: A Perturbative Determination of Mass Dependent $O(a)$ Improvement Coefficients in a Relativistic Heavy Quark Action

abstract: We present the results for a perturbative determination of mass dependent improvement coefficients $\nu$, $r_s$, $c_E$ and $c_B$ in a relativistic heavy quark action, which we have designed to control $m_Qa$ errors by extending the on-shell $O(a)$ improvement program to the case of $m_Q \gg \Lambda_{\rm QCD}$, where $m_Q$ is the heavy quark mass. The parameters $\nu$ and $r_s$ are determined from the quark propagator and $c_E$ and $c_B$ are from the on-shell quark-quark scattering amplitude. We show that all the parameters, together with the quark wave function and the mass renormalization factors, are determined free from infrared divergences once their tree level values are properly tuned. Results of these parameters are shown as a function of $m_Q a$ for various improved gauge actions.


title: Heavy hadrons on an anisotropic lattice

abstract: Results from simulations of quarkonia and heavy-light mesons on an anisotropic lattice are presented. The improved quark action and action-parameter tuning used in this study are discussed.


title: Non-perturbative propagators and dimension 2 condensate in Yang-Mills theory

abstract: We have found ultraviolet asymptotic slutions of the Schwinger-Dyson equation for the gluon and ghost propagators which have simultaneously the perturbative logarithmic correction and the non-perturbative $1/p^2$ power correction. By including the perturbative corrections, the power correction reproduces exactly the leading OPE result suggesting the existence of dimension two condensate.


title: Solving the Schwinger-Dyson Equations for Gluodynamics in the Maximal Abelian Gauge

abstract: We derive the Schwinger-Dyson equations for the SU(2) Yang-Mills theory in the maximal Abelian gauge and solve them in the infrared asymptotic region. We find that the infrared asymptotic solutions for the gluon and ghost propagators are consistent with the hypothesis of Abelian dominance.


title: Infrared Features of the Lattice Landau Gauge QCD

abstract: The infrared properties of lattice Landau gauge QCD of SU(3) are studied by measuring gluon propagator, ghost propagator, QCD running coupling and Kugo-Ojima parameter of $\beta=6.0, 16^4,24^4,32^4$ and $\beta=6.4, 32^4,48^4$ lattices. We study Gribov copy problem by parallel tempering (PT) gauge fixing of $\beta=2.2, 16^4$ SU(2) configurations. Rather unexpected features of the unique gauge(fundamental modular gauge, FMG) are observed in relation of physical values vs the optimizing function.


title: Infrared Features of Lattice Landau Gauge QCD and the Gribov Copy Problem

abstract: Infrared features of gluon propagator, ghost propagator, QCD running coupling and the Kugo-Ojima parameter in lattice Landau gauge QCD are presented. The framework of PMS analysis suggests that there appear infrared, intermediate and ultraviolet regions specified by $\Lambda_{\bar{MS}}$, $\beta_0$ and $\beta_1$. The propagators and the running coupling of $q>1GeV$ are fitted by the $\widetilde{MOM}$ scheme using the factorization scale as the scale parameter. The running coupling data of $q<14GeV$ are fitted by the contour improved perturbation method. The Gribov copy problem is studied in SU(2), $\beta= 2.2$ samples by comparing the data gauge fixed by the parallel tempering method and the data gauge fixed by the straightforward gauge fixing. The Gribov noise effect turned out to be about 2%.


title: Chiral Dirac fermions on the lattice using Geometric Discretisation

abstract: A new approach to the problem of doubling is presented with the Dirac-Kahler (DK) theory as a starting point and using Geometric Discretisation providing us with a new way of extracting the Dirac field in the discrete setting of a hyper-cubic complex.


title: I=2 Pion Scattering Length from Two-Pion Wave Function

abstract: We present a preliminary report on a calculation of scattering length for I=2 $S$-wave two-pion system directly from the two-pion wave function. Results are compared with those calculated from the time dependence of two-pion four-point functions. Calculations are made with an RG-improved action for gluons and improved Wilson action for quarks at $a^{-1}=1.207(12)$ GeV on $20^3 \times 48$ and $24^3 \times 48$ lattices.


title: Exploring two non-perturbative definitions of $c_A$

abstract: We present two determinations of the coefficient $c_A$ in quenched QCD, needed to build the $O(a)$ improved axial current. The first condition used is the requirement that the PCAC quark mass, as a function of $x_0$, stays flat for a non-trivial spatial phase for the fermions in the Schr\"odinger functional. The second condition is that the PCAC relation for the ground-state and the first excited state at finite L give the same quark mass. Our results confirm previous findings that in the quenched theory the intrinsic $O(a)$ ambiguity of $c_A$ gets relevant around $\beta\simeq6.0$.


title: Features of SU(N) Gauge Theories

abstract: We review recent lattice results for the large $N$ limit of SU(N) gauge theories. In particular, we focus on glueball masses, topology and its relation to chiral symmetry breaking (relevant for phenomenology), on the tension of strings connecting sources in higher representations of the gauge group (relevant for models of confinement and as a comparative ground for theories beyond the Standard Model) and on the finite temperature deconfinement phase transition (relevant for RHIC-like experiments). In the final part we present open challenges for the future.


title: Large cutoff effects of dynamical Wilson fermions

abstract: We present and discuss results for cutoff effects in the PCAC masses and the mass dependence of r_0 for full QCD and various fermion actions. Our discussion of how one computes mass dependences - here of r_0 - is also relevant for comparisons with chiral perturbation theory.


title: Adventures in Coulomb Gauge

abstract: We study the phase structure of SU(2) gauge theories at zero and high temperature, with and without scalar matter fields, in terms of the symmetric/broken realization of the remnant gauge symmetry which exists after fixing to Coulomb gauge. The symmetric realization is associated with a linearly rising color Coulomb potential (which we compute numerically), and is a necessary but not sufficient condition for confinement.


title: Pion electromagnetic form-factor with domain wall fermions

abstract: Motivated by recent measurements at J-Lab, the pion electromagnetic form-factor is investigated with quenched domain wall fermions and a renormalization group improved gauge action called DBW2. We see that quark mass dependence of the form-factor with finite momentum transfers is rather small.


title: Calculation of CP Violation in Non-leptonic Kaon Decay on the Lattice

abstract: We give a progress report of our lattice calculation of direct and indirect CP violation in kaon decays, parametrized as $\epsilon'/\epsilon$ and $B_K$, which require non-perturbative calculation of the matrix elements of the Standard Model effective Hamiltonian.


title: Calcutation of kaon matrix elements in quenched domain-wall QCD with DBW2 gauge action

abstract: We give a progress report of our new $a^{-1}\approx 3$ GeV quenched calculation of kaon matrix elements with domain-wall fermion and DBW2 gauge action. Our smaller lattice spacing allows us to address the effect of charmed quark on the lattice. We show preliminary results of $B_K$ renormalized non-perturbatively and $K\to\pi$ matrix elements.


title: Structure of the baryonic flux tube in N_{f}=2 lattice QCD at finite temperature

abstract: We study the flux tube profile in the baryonic system in full QCD at finite temperature on $N_{t}=8$ lattice. We fix the maximally Abelian gauge and measure the monopole and the photon parts of the Abelian action density, the color electric field and the monopole current on both sides of the finite temperature transition. We demonstrate the disappearance of the flux tube in the high temperature phase.


title: Excited nucleon spectrum from lattice QCD with maximum entropy method

abstract: We study excited states of the nucleon in quenched lattice QCD with the spectral analysis using the maximum entropy method. Our simulations are performed on three lattice sizes $16^3\times 32$, $24^3\times 32$ and $32^3\times 32$, at $\beta=6.0$ to address the finite volume issue. We find a significant finite volume effect on the mass of the Roper resonance for light quark masses. After removing this systematic error, its mass becomes considerably reduced toward the direction to solve the level order puzzle between the Roper resonance $N'(1440)$ and the negative-parity nucleon $N^*(1535)$.


title: Light meson correlation functions near the deconfining transition on anisotropic lattices

abstract: We study the light hadron correlators near the deconfining transition by extracting the spectral function on quenched anisotropic lattices. We adopt the method successfully applied to the charmonium systems: the use of the smeared operators and the analysis method which is composed of the maximum entropy method and the $\chi^2$ fitting assuming several forms. The numerical simulations are performed on lattices of $a_s^{-1} \simeq 2$ GeV and $a_s/a_t=4$, with the clover quark around the strange quark mass. Towards applications to the plasma phase, we check the reliability of our procedure by examining how the results for the correlators below $T_c$ ($N_t=32$) changes under variation of input parameters such as the smearing function.


title: Large Mass and Chemical Potential Model: A Laboratory for QCD?

abstract: We use a model based on the hopping parameter expansion to study QCD at large $\mu$. We find interesting behavior in the region expected to show flavor-color locking.


title: Excitations of torelon

abstract: The excitations of gluonic flux tube in a periodic lattice are examined. Monte Carlo simulations from an anisotropic lattice are presented and the comparison with effective string models is discussed.


title: Non-perturbative renormalization of moments of parton distribution functions

abstract: We compute non-perturbatively the evolution of the twist-2 operators corresponding to the average momentum of non-singlet quark densities. The calculation is based on a finite-size technique, using the Schr\"odinger Functional, in quenched QCD. We find that a careful choice of the boundary conditions, is essential, for such operators, to render possible the computation. As a by-product we apply the non-perturbatively computed renormalization constants to available data of bare matrix elements between nucleon states.


title: D-Theory: Field Quantization by Dimensional Reduction of Discrete Variables

abstract: D-theory is an alternative non-perturbative approach to quantum field theory formulated in terms of discrete quantized variables instead of classical fields. Classical scalar fields are replaced by generalized quantum spins and classical gauge fields are replaced by quantum links. The classical fields of a d-dimensional quantum field theory reappear as low-energy effective degrees of freedom of the discrete variables, provided the (d+1)-dimensional D-theory is massless. When the extent of the extra Euclidean dimension becomes small in units of the correlation length, an ordinary d-dimensional quantum field theory emerges by dimensional reduction. The D-theory formulation of scalar field theories with various global symmetries and of gauge theories with various gauge groups is constructed explicitly and the mechanism of dimensional reduction is investigated.


title: The pion form factor from first principles

abstract: We calculate the electromagnetic form factor of the pion in quenched lattice QCD. The non-perturbatively improved Sheikoleslami-Wohlert lattice action is used together with the ${\mathcal O}(a)$ improved current. We calculate form factor for pion masses down to $m_\pi = 360 MeV$. We compare the mean square radius for the pion extracted from our form factors to the value obtained from the 'Bethe Salpeter amplitude'. Using (quenched) chiral perturbation theory, we extrapolate our results towards the physical pion mass.


title: A global optimization method for Landau gauge fixing in Lattice QCD

abstract: An algorithm for gauge fixing to the minimal Landau gauge in lattice QCD is described. The method, a combination of an evolutionary algorithm with a steepest descent method, is able to solve the problem of the nonperturbative gauge fixing. The performance of the combined algorithm is investigated on $8^4$, $\beta = 5.7$, and $16^4$, $\beta = 6.0$, lattice SU(3) gauge configurations.


title: Non-perturbative renormalization of vector and axial vector currents in quenched QCD for a renormalization group improved gauge action

abstract: Renormalization constants of vector ($Z_V$) and axial-vector ($Z_A$) currents are determined non-perturbatively in quenched QCD for an RG-improved gauge action and a tadpole-improved clover quark action using the Schr\"odinger functional method. Meson decay constants $f_\rho$ and $f_\pi$ show much better scaling when $Z_V$ and $Z_A$ estimated for infinite physical volume are used instead of $Z$-factors from tadpole-improved one-loop perturbation theory.


title: Electromagnetic Form Factors with FLIC Fermions

abstract: The Fat-Link Irrelevant Clover (FLIC) fermion action provides a new form of nonperturbative O(a) improvement and allows efficient access to the light quark-mass regime. FLIC fermions enable the construction of the nonperturbatively O(a)-improved conserved vector current without the difficulties associated with the fine tuning of the improvement coefficients. The simulations are performed with an O(a^2) mean-field improved plaquette-plus-rectangle gluon action on a 20^3 x 40 lattice with a lattice spacing of 0.128 fm, enabling the first simulation of baryon form factors at light quark masses on a large volume lattice. Magnetic moments, electric charge radii and magnetic radii are extracted from these form factors, and show interesting chiral nonanalytic behavior in the light quark mass regime.


title: Chiral Structure in Baryon Magnetic Moments

abstract: We highlight recent advances in finite-range regularised chiral effective field theory. Application to baryon magnetic moments show an unambiguous signal of quenched chiral physics in FLIC fermion simulations.


title: Numerical study of the equation of state for two flavor QCD at finite density

abstract: We discuss the equation of state for 2 flavor QCD at non-zero temperature and density. Derivatives of $\ln Z$ 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 a Taylor series expansion. It is found that the fluctuations in the quark number density increase in the vicinity of the phase transition temperature and the susceptibilities start to develop a pronounced peak as $\mu_q$ is increased. This suggests the presence of a critical endpoint in the $(T, \mu_q)$ plane.


title: Lattice QCD in the epsilon-regime and random matrix theory

abstract: In the epsilon-regime of QCD the main features of the spectrum of the low-lying eigenvalues of the (euclidean) Dirac operator are expected to be described by a certain universality class of random matrix models. In particular, the latter predict the joint statistical distribution of the individual eigenvalues in any topological sector of the theory. We compare some of these predictions with high-precision numerical data obtained from lattice QCD for a range of lattice spacings and volumes. While no complete matching is observed, the results agree with theoretical expectations at volumes larger than about 5 fm^4.


title: Dimensional reduction and a Z(3) symmetric model

abstract: We present first results from a numerical investigation of a Z(3) symmetric model based on dimensional reduction.


title: Free Meson Spectral Functions on the Lattice

abstract: We present results from an analytic calculation of thermal meson spectral functions in the infinite temperature (free field) limit. We compare spectral functions for various lattice fermion formulations used at present in studies of in-medium properties of hadrons based on the maximum entropy method (MEM). In particular, we will present a new calculation of spectral functions performed with extended quark sources.


title: One-loop calculation of mass dependent ${\cal O}(a)$ improvement coefficients for the relativistic heavy quarks on the lattice

abstract: We carry out the one-loop calculation of mass dependent ${\cal O}(a)$ improvement coefficients in the relativistic heavy quark action recently proposed, employing the ordinary perturbation theory with the fictitious gluon mass as an infrared regulator. We also determine renormalization factors and improvement coefficients for the axial-vector current at the one-loop level. It is shown that the improvement coefficients are infrared finite at the one-loop level if and only if the improvement coefficients in the action are properly tuned at the tree level.


title: Energy spectrum of a Heavy-Light meson

abstract: The energies of different angular momentum states of a heavy-light meson were measured on a lattice by Michael and Peisa in Phys. Rev. D 58, 34506 (1998). We have now repeated this study using a larger, quenched lattice to check that there the finite size effects are under control. Also the effect of quenching is studied by using dynamical fermion configurations. In the latter case the lattice spacing a is smaller, which allows one to make an estimate of a -> 0 corrections. The new measurements confirm the earlier calculated energy spectrum. From some theoretical considerations it is expected that, for high angular momentum states, the multiplets should be inverted compared with the Coulomb spectrum, i.e. L- should lie higher than L+ (see H. J. Schnitzer Phys. Lett. B 226, 171 (1989)). This happens when a pure scalar confining potential is present. Here the notation L+(-) means that the light quark spin couples to angular momentum L giving the total j=L+/- 1/2. Experimentally this inversion is not seen for P-wave mesons, and lattice measurements now suggest that there is no inversion in the D-wave either. Indeed we find that the spin-orbit splitting is very small for both values L=1 and 2. This almost exact cancellation need not be a coincidence, and could probably be understood by using straightforward symmetry arguments (P.R. Page, T. Goldman and J.N. Ginocchio, Phys. Rev. Lett. 86, 204 (2001)).


title: Extraction of parton distributions from lattice QCD

abstract: We review the calculation of moments of both the polarized and unpolarized parton distribution functions of the nucleon in lattice QCD, and in particular their extrapolation to the physical region. We also discuss the reconstruction of the x dependence of the valence quark distributions in the nucleon from a finite number of lattice moments.


title: The Factorization Method for Simulating Systems With a Complex Action

abstract: We propose a method for Monte Carlo simulations of systems with a complex action. The method has the advantages of being in principle applicable to any such system and provides a solution to the overlap problem. We apply it in random matrix theory of finite density QCD where we compare with analytic results. In this model we find non--commutativity of the limits $\mu\to 0$ and $N\to\infty$ which could be of relevance in QCD at finite density.


title: String effects in SU(2) lattice gauge theory

abstract: We discuss the effective string picture for the confining regime of lattice gauge theories at zero and finite temperature. We present results of extensive Monte Carlo simulations - performed with the Luscher and Weisz algorithm - for SU(2) Yang-Mills theory in 2+1 dimensions. We also address the issue of "string universality" by comparing our results with those obtained in other lattice gauge theories.


title: Preliminary results from a simulation of quenched QCD with overlap fermions on a large lattice

abstract: We simulate quenched QCD with the overlap Dirac operator. We work with the Wilson gauge action at beta=6 on an 18^3x64 lattice. We calculate quark propagators for a single source point and quark mass ranging from am_q=0.03 to 0.75. We present here preliminary results based on the propagators for 60 gauge field configurations.


title: A study of colour field distributions in the baryon

abstract: The distributions of chromo-electric and chromo-magnetic field associated with flux tubes in the baryon are studied in SU(3) lattice QCD. Maximal Abelian projection is used to reduce the statistical fluctuations. For a fixed source geometry, many different string configurations are possible. We investigated whether the string configuration, that is the choice of operator, biases the observed flux distribution.


title: Heavy quark masses from finite volume effects

abstract: I discuss the results of a new calculation of the charm and bottom quark masses in the quenched approximation and in the continuum limit of lattice QCD. The work has been done by the APE group at the ``Tor Vergata'' University making use of the step scaling method, previously introduced to deal with two scale problems, that allows to take the continuum limit of the lattice data. We have computed the RGI quark masses and then we have connected the results to the MSbar scheme. The continuum numbers are 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.


title: Heavy-light decay constants from the step scaling method

abstract: We discuss results for the heavy-light decay constants in the continuum limit of quenched lattice QCD from finite size scaling techniques. We disentangle the simultaneous presence of the different energy scales characterizing heavy-light physics by first performing simulations at the unphysical volume L_0=0.4 fm, and then evolving the results towards the infinite volume. We find f_{Bs}= 192(6)(4) MeV and f_{Ds}=240(5)(5) MeV. The approach has been developed by the APE group at the University of Rome ``Tor Vergata''.


title: Hadronic transitions from the lattice

abstract: I discuss strategies to determine hadronic decay couplings from lattice studies. As a check of the methods, I explore the decay of a vector meson to two pseudoscalar mesons with $N_f=2$ flavours of sea quark. Although we are working with quark masses that do not allow a physical decay, I show how the transition rate can be evaluated from the amplitude for $\rho \to \pi \pi$ and from the annihilation component of $\pi \pi \to \pi \pi$. I explore the decay amplitude for two different pion momenta and find consistent results. The coupling strength found is in agreement with experiment. I also find evidence for a shift in the $\rho$ mass caused by mixing with two pion states. I also present results for the decay of a hybrid meson, for the case of heavy valence quarks.


title: Dynamics of Monopoles and Flux Tubes in Two-Flavor Dynamical QCD

abstract: We investigate the confining properties of the QCD vacuum with $N_f=2$ flavors of dynamical quarks, and compare the results with the properties of the quenched theory. We use non-perturbatively $\mathcal{O}(a)$ improved Wilson fermions to keep cut-off effects small. We focus on color magnetic monopoles. Among the quantities we study are the monopole density and the monopole screening length, the static potential and the profile of the color electric flux tube. We furthermore derive the low-energy effective monopole action. Marked differences between the quenched and dynamical vacuum are found.


title: Unphysical Operators in Partially Quenched QCD

abstract: We point out that the chiral Lagrangian describing pseudo-Goldstone bosons in partially quenched QCD has one more four-derivative operator than that for unquenched QCD with three flavors. The new operator can be chosen to vanish in the unquenched sector of the partially quenched theory. Its contributions begin at next-to-leading order in the chiral expansion. At this order it contributes only to unphysical scattering processes, and we work out some examples. Its contributions to pseudo-Goldstone properties begin at next-to-next-to-leading order, and we determine their form. We also determine all the zero and two derivative operators in the $O(p^6)$ partially quenched chiral Lagrangian, finding three more than in unquenched QCD, and use these to give the general form of the analytic next-to-next-to-leading order contributions to the pseudo-Goldstone mass and decay constant. We discuss the general implications of such additional operators for the utility of partially quenched simulations


title: Glueball Matrix Elements on Anisotropic Lattices

abstract: The glueball-to-vacuum matrix elements of local gluonic operators in scalar, tensor, and pseudoscalar channels are investigated numerically on several anisotropic lattices with the spatial lattice spacing in the range 0.1fm -- 0.2fm. These matrix elements are needed to predict the glueball branching ratios in $J/\psi$ radiative decays which will help to identify the glueball states in experiments. Two types of improved local gluonic operators are constructed for a self-consistent check, and the finite volume effects are also studied. The lattice spacing dependence of our results is very small and the continuum limits are reliably extrapolated.


title: Lattice study of exotic S=+1 baryon

abstract: We propose S=+1 baryon interpolating operators, which are based on an exotic description of the antidecuplet baryon like diquark-diqaurk-antiquark. By using one of the new operators, the mass spectrum of the spin-1/2 pentaquark states is calculated in quenched lattice QCD at $\beta=6/g^2=6.2$ on a $32^3\times48$ lattice. It is found that the $J^P$ assignment of the lowest $\Theta(uudd \bar s)$ state is most likely $(1/2)^-$. We also calculate the mass of the charm analog of the $\Theta$ and find that the $\Theta_c(uudd \bar c)$ state lies much higher than the $DN$ threshold, in contrast to several model predictions.


title: Lattice QCD on Earth Simulator

abstract: We report on coding and performance of our polynomial hybrid Monte Carlo program on the Earth Simulator. At present the entire program achieves 25--40% efficiency. An analysis of overheads shows that a tuning of inter-node communications is required for further improvement.


title: A More Improved Lattice Action for Heavy Quarks

abstract: We extend the Fermilab formalism for heavy quarks to develop a more improved action. We give results of matching calculations of the improvement couplings at tree level. Finally, we estimate the discretization errors associated with the new action.


title: Consistency of lattice definitions of U(1) flux in Abelian projected SU(2) gauge theory

abstract: We reexamine the dual Abrikosov vortex under the requirement that the lattice averages of the fields satisfy exact Maxwell equations [ME]. The electric ME accounts for the total flux and the magnetic ME determines the shape of the confining string. This leads to unique and consistent definitions of flux and electric and magnetic currents at finite lattice spacing. The resulting modification of the standard DeGrand-Toussaint construction gives a magnetic current comprised of smeared monopoles.


title: A New Method for Simulating QCD at Finite Density

abstract: We propose a new method for simulating QCD at finite density, where interesting phases such as the color superconductivity phase is conjectured to appear. 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. The achieved system size is large enough to extract the thermodynamic limit. Our results provide a clear understanding of how the expected first order phase transition is induced by the imaginary part of the action. We also discuss the noncommutativity of the zero chemical potential limit and the thermodynamic limit, which is relevant to recent Monte Carlo studies at small chemical potential.


title: Lattice Superstring and Noncommutative Geometry

abstract: Recent developments in superstring theory and noncommutative geometry are deeply related to the idea of Eguchi-Kawai reduction in large N gauge theories which dates back to early 80s. After a general review on this subject including revived interests in solving planar QCD, we present some results in the superstring matrix model suggesting the dynamical generation of 4d space-time due to the collapse of the eigenvalue distribution. We then discuss interesting dynamical properties of field theories in noncommutative geometry, which have been revealed by Monte Carlo simulations of twisted reduced models. We conclude with a comment on the recent proposal for a lattice construction of supersymmetric gauge theories based on reduced models.


title: An Approach to Higher Dimensional Theories Based on Lattice Gauge Theory

abstract: A higher dimensional lattice space can be decomposed into a number of four-dimensional lattices called as layers. The higher dimensional gauge theory on the lattice can be interpreted as four-dimensional gauge theories on the multi-layer with interactions between neighboring layers. We propose the new possibility to realize the continuum limit of a five-dimensional theory based on the property of the phase diagram.


title: Confinement of Color: Recent Progress

abstract: Recent progress done in Pisa on the subject is presented. It is shown that dual superconductivity of the vacuum or absence of it is an intrinsic property of QCD vacuum, independent of the choice of the abelian projection. The order of the deconfining phase transition in $N_f=2$ QCD is studied as a key to understand the mechanism of confinement.


title: What is the order of the deconfining phase transition?

abstract: The techniques are discussed by which the order of the deconfining phase transition is investigated on the lattice. QCD with $N_f=2$ is a special case, which can provide information on the mechanism of confinement.


title: Confinement of Color: A Review

abstract: The status of our understanding of the mechanisms of color confinement is reviewed, in particular the results of numerical simulations on the lattice.


title: Classical simulations and particle production in heavy-ion collisions

abstract: The classical approximation may be applied to a number of problems in non-equilibrium field theory. The principles and limits of classical real-time lattice simulations are presented, with particular emphasis on the definition of particle numbers and energies and on applications to the earliest stages of heavy-ion collisions.


title: Consistent Definitions of Flux and Electric and Magnetic Current in Abelian Projected SU(2) Lattice Gauge Theory

abstract: Through the use of a lattice U(1) Ward-Takahashi identity, one can find a precise definition of flux and electric four-current that does not rely on the continuum limit. The magnetic four-current defined for example by the DeGrand-Toussaint construction introduces order a^2 errors in the field distributions. We advocate using a single definition of flux in order to be consistent with both the electric and magnetic Maxwell's equations at any lattice spacing. In a U(1) theory the monopoles are slightly smeared by this choice, i.e. are no longer associated with a single lattice cube. In Abelian projected SU(2) the consistent definition suggests further modifications. For simulations in the scaling window, we do not foresee large changes in the standard analysis of the dual Abrikosov vortex in the maximal Abelian gauge because the order a^2 corrections have small fluctuations and tend to cancel out. However in other gauges, the consequences of our definitions could lead to large effects which may help in understanding the choice of gauge. We also examine the effect of truncating all monopoles except for the dominant cluster on the profile of the dual Abrikosov vortex.


title: Large scale simulations of lattice gauge theory: key NSF project of China and recent developments

abstract: In a recent funded Key Project of National Science Foundation of China, we have planned to do large scale simulations of lattice Quantum Chromodynamics, using the parallel supercomputing facilities in mainland China. Here I briefly review our plan and progress in recent years.


title: Vacuum polarisation and the muon g-2 in lattice QCD

abstract: We measure the hadronic contribution to the vacuum polarisation tensor, and use it to estimate the hadronic contribution to (g-2)_mu, the muon anomalous magnetic moment.


title: Quark spectra and light hadron phenomenology from overlap fermions with improved gauge field action

abstract: We present first results from a simulation of quenched overlap fermions with improved gauge field action. Among the quantities we study are the spectral properties of the overlap operator, the chiral condensate and topological charge, quark and hadron masses, and selected nucleon matrix elements. To make contact with continuum physics, we compute the renormalization constants of quark bilinear operators in perturbation theory and beyond.


title: Full QCD Algorithms towards the Chiral Limit

abstract: I discuss the behaviour of algorithms for dynamical fermions as the sea-quark mass decreases. I focus on the Hybrid-Monte-Carlo (HMC) algorithm applied to two degenerate flavours of Wilson fermions. First, I briefly review the performance obtained in large scale HMC simulations. Then I discuss a modified pseudo-fermion action for the HMC simulation that has been introduced three years ago. I summarize recent results obtained with this pseudo-fermion action by the QCDSF and the ALPHA collaborations. I comment on alternatives to the HMC, like the Multiboson algorithm and variants of it.


title: On Dirac sheet configurations of SU(2) lattice fields

abstract: Finite temperature Euclidean SU(2) lattice gauge fields close to the deconfinement phase transition are subjected to cooling. We find relatively stable or absolutely stable configurations with an action below the one-instanton action $S_{inst}=2\pi^2$ both in the deconfinement and the confinement phases. In this paper we attempt an interpretation of these lowest action configurations. Their action is purely magnetic and amounts to $S/S_{inst} \approx N_t/N_s$, where $N_t$ ($N_s$) is the timelike (spacelike) lattice size, while the topological charge vanishes. In the confined phase part of the corresponding lattice configurations turns out to be absolutely stable with respect to the cooling process in which case Abelian projection reveals a homogeneous, purely Abelian magnetic field closed over the "boundary" in one of the spatial directions. Referring to the dyonic structure established for the confinement phase near $T_c$ and based on the observation made for this phase that such events below the instanton action $S_{inst}$ emerge from dyon-antidyon annihilation, the question of stability (metastability) is discussed for both phases. The hypothetically different dyonic structure of the deconfinement phase, inaccessible by cooling, could explain the metastability.


title: A Novel Ansatz for the Energy-Momentum Tensor on the Lattice

abstract: The comparison of structural analogies between the energy-momentum tensors in general relativity and in a gauge theory of Yang-Mills type is tentatively extended to lattice physics. These considerations are guiding to a new lattice model for the symmetric energy-momentum tensor $\Theta_{\mu\nu}$ of the pure Yang-Mills gauge sector, basing on half powers of the plaquette variable. The concept of non-trivial principal square roots of unitary matrices in lattice gauge theories can be epitomized to reconcile the pretension to a uniform construction principle for the components of $\Theta_{\mu\nu}$ with general qualitative thermodynamic demands concerning arguments in favour of a Wilson form for $<\Theta^\mu_\mu>$ and $\Theta_{44}$. SU(2) Monte Carlo results for the Euclidean expectation values on a 10**4 lattice are compared with that of competing hitherto existing lattice models for $\Theta_{\mu\nu}$.


title: Spontaneous symmetry breaking in strong-coupling lattice QCD at high density

abstract: We determine the patterns of spontaneous symmetry breaking in strong-coupling lattice QCD in a fixed background baryon density. We employ a next-nearest-neighbor fermion formulation that possesses the SU(N_f)xSU(N_f) chiral symmetry of the continuum theory. We find that the global symmetry of the ground state varies with N_f and with the background baryon density. In all cases the condensate breaks the discrete rotational symmetry of the lattice as well as part of the chiral symmetry group.


title: On the gauge dependence of the singlet and adjoint potentials

abstract: We study gauge dependence of the recently suggested definition of the singlet and adjoint potentials in SU(2) lattice gauge theory. We find that in the (time local) maximal tree axial gauge the singlet potential obtained from the gauge dependent correlator ${Tr} L(x)L^\dagger(y)$ differs from that computed in the Coulomb gauge. In the generalized Coulomb gauge we find the range of the parameter values in which the singlet potential differs from that in the Coulomb gauge.


title: The Study of the Continuum Limit of the Supersymmetric Ward-Takahashi Identity for N=1 Super Yang-Mills Theory

abstract: The one-loop corrections to the supersymmetric Ward-Takahashi identity (WTi) are investigated in the off-shell regime in the Wilson formulation of the discretized N=1 Super Yang-Mills (SYM) theory. The study of the continuum limit as well as the renormalization procedure for the supercurrent are presented.


title: Non-perturbative Heavy Quark Effective Theory

abstract: We explain how to perform non-perturbative computations in HQET on the lattice. In particular the problem of the subtraction of power-law divergences is solved by a non-perturbative matching of HQET and QCD. As examples, we present a full calculation of the mass of the b-quark in the combined static and quenched approximation and outline an alternative way to obtain the B-meson decay constant at lowest order. Since no excessively large lattices are required, our strategy can also be applied including dynamical fermions.


title: Non perturbative study of QCD on a 2+2 anisotropic lattice

abstract: We present preliminary results for a non perturbative determination of the parameters in the action required to restore Lorentz invariance in long distance correlators, using the static interquark potential. Comparison with analytical results is made and further applications are discussed.


title: Quenched and Unquenched Chiral Perturbation Theory in the \epsilon-Regime

abstract: The chiral limit of finite-volume QCD is the $\epsilon$-regime of the theory. We discuss how this regime can be used for determining low-energy observables of QCD by means of comparisons between lattice simulations and quenched and unquenched chiral perturbation theory. The quenched theory suffers in the $\epsilon$-regime from ``quenched finite volume logs'', the finite-volume analogs of quenched chiral logs.


title: A BCS Gapped Superfluid on the Lattice

abstract: We discuss results of numerical analyses of the 3+1 dimensional Nambu--Jona-Lasinio model at non-zero baryon density. By studying relevant order parameters, it is shown that the model goes through a transition between a vacuum with broken chiral symmetry, and a chirally restored phase at high baryon chemical potential (mu) with a diquark condensate roughly proportional to the area of the Fermi surface and a broken U(1) baryon number symmetry. The condensation of diquark pairs from the fermionic spectrum is shown to lead to an energy gap (Delta) about the Fermi energy (E_F), which is found to be approximately 15% of the vacuum fermion mass and roughly independent of mu in the chirally restored phase. The ground-state is believed, therefore, to be that of a BCS superfluid. We also discuss the effect of simulating on a finite volume, and any caveats this places on the conclusions made above.


title: The Casimir Energy Paradox of the QCD String

abstract: It is widely thought that the early onset of the asymptotic Casimir energy with unit conformal charge signals bosonic string formation of the confining flux connecting a static quark-antiquark pair in QCD. This is observed on a scale where most of the string eigenmodes do not exist and the few stable modes above the ground state are displaced. Hints for the resolution of this paradox are suggested.


title: Three flavor dynamical QCD project by CP-PACS/JLQCD

abstract: The CP-PACS and JLQCD Collaborations have been made systematic studies of lattice QCD carrying out both chiral and continuum extrapolations. Importance of dynamical quark effects has been clarified by a comparison of quenched QCD and two flavor full QCD simulations. In two flavor simulations, the dynamical effects of $u$ and $d$ quarks are taken into account, but the third quark $s$ is still treated in a quenched approximation. As the final step towards a fully realistic lattice simulation of QCD, two collaborations started a joint project of three flavor QCD, concentrating all big computers available. Based on a series of preparative studies of exact simulation algorithm and non-perturbative improvement coefficient, a large scale simulation of three flavor QCD has been started. I present the results of light hadron spectrum and light quark masses from the first production runs.


title: Results on improved KS dynamical configurations: spectrum, decay constants, etc

abstract: The MILC Collaboration has been producing ensembles of lattice configurations with three dynamical flavors for the past few years. There are now results for three lattice spacings for a variety of light and strange quark masses, ranging down to $m_l=0.1 m_s$, where $m_s$ is the dynamical strange quark mass and $m_l$ is the common mass of the $u$ and $d$ quarks. Recently, the Fermilab, HPQCD, MILC and UKQCD collaborations have presented a summary of results obtained using these lattices. Compared with quenched results, these new calculations show great improvement in agreement with experiment. This talk addresses the technical improvements that make these calculations possible and provides additional details of calculations not presented in the initial summary. We demonstrate that a wide range of hadronic observables can now be calculated to 2--3% accuracy.


title: Properties of charmonium in lattice QCD with 2+1 flavors of improved staggered sea quarks

abstract: We use the dynamical gluon configurations provided by the MILC collaboration in a study of the charmonium spectrum and psi leptonic width. We examine sea quark effects on mass splitting and on the leptonic decay matrix element for light masses as low as m_s/5, while keeping the strange quark mass fixed and the lattice spacing nearly constant.


title: Recent Developments of Domain-wall/Overlap Fermions for Lattice QCD

abstract: I review the lattice formulations of vector-like gauge theories (e.g. QCD) with domain-wall/overlap fermions, and discuss how to optimize the chiral symmetry for any finite $ N_s $ (sites in the fifth dimension). In this formulation, quark propagators in gauge background can be computed efficiently through the effective 4D lattice Dirac operator.


title: Higher-order perturbation theory for highly-improved actions

abstract: I review techniques and applications of higher-order perturbation theory for highly-improved lattice actions.


title: D_s spectrum and leptonic decays with Fermilab heavy quarks and improved staggered light quarks

abstract: We present preliminary results for the D_s meson spectrum and decay constants in unquenched lattice QCD. Simulations are carried out with 2+1 dynamical quarks using gauge configurations generated by the MILC collaboration. We use the ``asqtad'' a^2 improved staggered action for the light quarks, and the clover heavy quark action with the Fermilab interpretation. We compare our spectrum results with the newly discovered 0+ and 1+ states in the D_s system.


title: Phase Fluctuation of Fermion Determinant in Lattice QCD at Finite Density

abstract: Once the quark chemical potential $\mu$ is introduced in finite density QCD, the fermion determinant appeared in the path integral measure becomes complex. In order to investigate the phase effect of SU(3) lattice QCD (2-flavors), we calculated the fluctuation of the phase of $\det\Delta(\mu)$ on a $8^3\times4$ lattice at $\mu = 0.1$ and 0.2. Then we calculated the chiral condensate and Polyakov line in the no phase and reweighted cases. There is little difference between these two cases at $\mu = 0.1$ and 0.2. We consider a possible reason for this result below $\mu \le 0.27$ in terms of $Z_3$ symmetry.


title: Recent progress in calculating weak matrix elements using staggered fermions

abstract: We present a chronological review of the progress in calculating weak matrix elements using staggered fermions. We review the perturbative calculation of one-loop matching formula including both current-current diagrams and penguin diagrams using improved staggered fermions. We also present preliminary results of weak matrix elements relevant to CP violation calculated using the improved (HYP (II)) staggered fermions. Since the complete set of matching coefficients at the one-loop level became available recently, we have constructed lattice operators with all the $g^2$ corrections included. The main results include both $\Delta I = 3/2$ and $\Delta I = 1/2$ contributions.


title: Solution of the Dirac equation in lattice QCD using a domain decomposition method

abstract: Efficient algorithms for the solution of partial differential equations on parallel computers are often based on domain decomposition methods. Schwarz preconditioners combined with standard Krylov space solvers are widely used in this context, and such a combination is shown here to perform very well in the case of the Wilson--Dirac equation in lattice QCD. In particular, with respect to even-odd preconditioned solvers, the communication overhead is significantly reduced, which allows the computational work to be distributed over a large number of processors with only small parallelization losses.


title: On the phase diagram of quenched QCD with Wilson fermions

abstract: In this talk, I reported on recent work on the Aoki phase diagram for quenched QCD with two flavors of Wilson fermions. Part of this work was done in collaboration with Yigal Shamir, and a shorter account of this part appeared recently. In this write-up, I will therefore limit myself to reporting on work done with Steve Sharpe and Robert Singleton, Jr, which has not been published yet. We discuss the symmetries of quenched QCD, paying careful attention to non-perturbative issues. This allows us to derive an effective lagrangian which agrees with standard quenched chiral perturbation theory, but which can also be used to address questions of a non-perturbative nature.


title: Lattice measurement of the couplings g_infinity and g_Bstar_B_pi

abstract: We present the results of a quenched lattice QCD study of the coupling g_hat, in the static heavy quark limit. After combining this with our previous results obtained by using propagating heavy quarks with a mass around the physical charm quark, we are able to interpolate to the b-quark sector. Our results are g_hat_infinity = 0.48 +/- 0.03 +/- 0.11, g_hat_b = 0.58 +/- 0.06 +/- 0.10, and g_Bstar_ B_p = 47 +/- 5 +/- 8.


title: Lattice QCD at finite T and \mu

abstract: Recent results of lattice QCD at finite temperature and density are reviewed. At vanishing density the transition temperature, the equation of state and hadron properties are discussed both for the pure gauge theory and for dynamical staggered, Wilson and overlap fermions. The second part deals with finite density. There are recent results for full QCD at finite temperature and moderate density, while at larger densities QCD-like models are studied.


title: Simulation of Wilson fermion at finite isospin density

abstract: For QCD with Wilson fermions at isospin chemical potential we study the finite phase transition on an $8^3\times4$ lattice at $\kappa=0.15$. We use two gauge actions: Wilson action and DBW2 action. Both actions give the same results. The phase diagram is qualitatively similar to the one obtained for QCD at small baryon chemical potential. We also calculate the number density for various isospin chemical potentials.


title: The pion electromagnetic form factor

abstract: A ratio of lattice correlation functions is identified from which the pion form factor can be obtained directly. Preliminary results from quenched Wilson simulations are presented.


title: Lattice regularized chiral perturbation theory

abstract: Chiral perturbation theory can be defined and regularized on a spacetime lattice. A few motivations are discussed here, and an explicit lattice Lagrangian is reviewed. A particular aspect of the connection between lattice chiral perturbation theory and lattice QCD is explored through a study of the Wess-Zumino-Witten term.


title: Dynamical Domain Wall Fermions

abstract: We report on an exploratory study of N_f=2 dynamical domain wall fermions and the DBW2 gauge action at weak coupling. Details of improved simulation algorithms and preliminary results for the hadron spectrum and renormalised light and strange quark masses will be presented.


title: Description and comparison of Fat7-bar and HYP fat links

abstract: We study various methods of constructing fat links based upon the HYP (by Hasenfratz & Knechtli) and Fat7-bar (by W. Lee) algorithms. We present the minimum plaquette distribution for these fat links. This enables us to determine which algorithm is most effective at reducing the spread of plaquette values - a strong indicator of improved statistics for spectrum and static potential measurements, among other quantities.


title: Center Dominance in SU(2) Gauge-Higgs Theory

abstract: We study the SU(2) gauge-Higgs system in D=4 dimensions, and analyze the influence of the fundamental-representation Higgs field on the vortex content of the gauge field. It is shown that center projected Polyakov lines, at low temperature, are finite in the infinite volume limit, which means that the center vortex distribution is consistent with color screening. In addition we confirm and further investigate the presence of a "Kertesz-line" in the strong-coupling region of the phase diagram, which we relate to the percolation properties of center vortices. It is shown that this Kertesz-line separates the gauge-Higgs phase diagram into two regions: a confinement-like region, in which center vortices percolate, and a Higgs region, in which they do not. The free energy of the gauge-Higgs system, however, is analytic across the Kertesz line.


title: B_K from Two-flavor Dynamical Domain Wall Fermions

abstract: We report preliminary results from an ongoing calculation of $B_K$ for $N_f=2$ dynamical QCD with domain wall fermions. Simulations have been done with three dynamical quark masses on $16^3 \times 32$ volumes with $L_s = 12$, where the lattice spacing is $a^{-1} = 1.81(6)$ GeV. Using measurements on $\sim 70$ lattices for each dynamical mass and extrapolating $\mdyn=\mval$ to the kaon point, we find $B_K^{\rm \bar{MS}}(\mu=2{\rm GeV}) = 0.503(20)$.


title: Quarkonium at finite temperature

abstract: I discuss quarkonium spectral functions at finite temperature reconstructed using the Maximum Entropy Method. I show in particular that the J/psi survives in the deconfined phase up to 1.5Tc


title: Non-Perturbative Renormalization and the Fermilab Action

abstract: We discuss the application of the regularization independent (RI) scheme of Rome/Southampton to determine the normalization of heavy quark operators non-perturbatively using the Fermilab action.


title: Vector and chiral gauge theories on the lattice

abstract: We combine a pair of independent Weyl fermions to compose a Dirac fermion on the four-dimensional Euclidean lattice. The obtained Dirac operator is antihermitian and does not reproduce anomaly under the usual chiral transformation. To simulate the correct chiral anomaly, we modify the chiral transformation. We also show that chiral gauge theories can be constructed nonperturbatively with exact gauge invariance. The formulation is based on a doubler-free lattice derivative, which is a simple matrix defined as a discrete Fourier transform of momentum with antiperiodic boundary conditions. Long-range fermion hopping interactions are truncated using the Lanczos factor.


title: The Chiral Condensate of Strongly Coupled QCD in the 't Hooft Limit

abstract: Using the recently proposed generalization to an arbitrary number of colors of the strong coupling approach to lattice gauge theories\cite{Grignani:2003uv}, we compute the chiral condensate of massless QCD in the 't Hooft limit.


title: Heavy Quarks and Lattice QCD

abstract: This paper is a review of heavy quarks in lattice gauge theory, focusing on methodology. It includes a status report on some of the calculations that are relevant to heavy-quark spectroscopy and to flavor physics.


title: Lattice calculation of the lowest order hadronic contribution to the muon anomalous magnetic moment: an update with Kogut-Susskind fermions

abstract: I present a preliminary calculation of the hadronic vacuum polarization for 2+1 flavors of improved Kogut-Susskind quarks by utilizing a set of gauge configurations recently generated by the MILC collaboration. The polarization function $\Pi(q^2)$ is then used to calculate the lowest order (in $\alpha_{QED}$) hadronic contribution to the muon anomalous magnetic moment.


title: Pion Mass Dependence of Nucleon Magnetic Moments

abstract: The relevance of the pion mass, provenient from a term which explicitely breaks chiral symmetry in the Lagrangian, for nucleon magnetic moment in the frame of the Skyrmion model in two different versions: the usual Skyrme model and a modified one which includes a coupling to a light scalar meson field, the sigma $\sigma (\simeq 500-600$ MeV). The results are compared to other calculations. Our main motivation comes from usual extrapolations for values of low energy QCD observables obtained in lattices with large values of pion/quark masses toward realistic value of $m_{\pi}$. which do not allow it. We do a comparison with results from the Cloudy Bag Model and a chiral hadronic model from chiral perturbation theory. There are several resulting extrapolations from the region of large pion mass to the realistic value depending on the considered model for low energy QCD.


title: Monte Carlo Hamiltonian: Inverse Potential

abstract: The Monte Carlo Hamiltonian method developed recently allows to investigate ground state and low-lying excited states of a quantum system, using Monte Carlo algorithm with importance sampling. However, conventional MC algorithm has some difficulties when applying to inverse potentials. We propose to use effective potential and extrapolation method to solve the problem. We present examples from the hydrogen system.


title: Staggered versus overlap fermions: a study in the Schwinger model with $N_f=0,1,2$

abstract: We study the scalar condensate and the topological susceptibility for a continuous range of quark masses in the Schwinger model with $N_f=0,1,2$ dynamical flavors, using both the overlap and the staggered discretization. At finite lattice spacing the differences between the two formulations become rather dramatic near the chiral limit, but they get severely reduced, at the coupling considered, after a few smearing steps.


title: I=2 $\pi\pi$ scattering using G-parity boundary condition

abstract: To make the $\pi\pi$ state with non-zero relative momentum as the leading exponential, we impose anti-periodic boundary condition on the pion, which is implemented by imposing G-parity or H-parity on the quark fields at the boundary. With this, we calculate the I=2 $\pi\pi$ phase shift from lattice simulation by using L\"uscher's formula.


title: Further Improvements to Staggered Quarks

abstract: The improved staggered quark action has enabled breakthrough calculations of QCD with three light dynamical quarks. The largest remaining $a^2$ scaling violations come from "taste-changing" interactions, which are significantly reduced over naive quarks because of carefully chosen smearing of the gauge links. Here we examine further improvements to the quark action and show that additional smearing of the links can reduce the taste-changing errors by another factor of two, whilst retaining a fully $a^2$-improved action. The taste-changing errors can be understood within a perturbative calculation of the four-quark operators that appear at one-loop. We also discuss unquenching techniques with such a highly improved action.


title: Triviality of hierarchical O(N) spin model in four dimensions with large N

abstract: The renormalization group transformation for the hierarchical O(N) spin model in four dimensions is studied by means of characteristic functions of single-site measures, and convergence of the critical trajectory to the Gaussian fixed point is shown for a sufficiently large N. In the strong coupling regime, the trajectory is controlled by the help of the exactly solved O(\infty) trajectory, while, in the weak coupling regime, convergence to the Gaussian fixed point is shown by power decay of the effective coupling constant.


title: Electromagnetic fluxes, monopoles, and the order of the 4d compact U(1) phase transition

abstract: We consider the 4d compact U(1) gauge theory with extended action S=-beta sum_P cos theta_P -gamma sum_P cos 2 theta_P We give a full characterization of the phase diagram of this model using the notion of flux. The relation with the usual monopole picture is discussed. In analogy with the XY model we consider the helicity modulus \cite{Jose:1977gm} for this theory, and show that it is an order parameter. Analyzing the finite-size effects of the helicity modulus we conclude that the transition is first-order. The value of this order parameter is related to the renormalized coupling beta_R. We measure beta^c_R at the transition point and give a counterexample to its conjectured universal value \cite{Cardy:jg}.


title: The helicity modulus in gauge field theories

abstract: We consider the 4d compact U(1) theory with Wilson action and characterize its phase diagram using the notion of electromagnetic flux, instead of the more usual magnetic monopole. Taking inspiration from the flux picture, we consider the helicity modulus (h.m.) for this theory, and show that it is an order parameter for the confinement deconfinement phase transition. We extend the definition of the h.m. to an Abelian projected Yang-Mills theory, and discuss its behavior in SU(2).


title: Twisted-mass lattice QCD with mass non-degenerate quarks

abstract: The maximally twisted lattice QCD action of an $SU_f(2)$ doublet of mass degenerate Wilson quarks gives rise to a real positive fermion determinant and it is invariant under the product of standard parity times the change of sign of the coefficient of the Wilson term. The existence of this spurionic symmetry implies that O($a$) improvement is either automatic or achieved through simple linear combinations of quantities taken with opposite external three-momenta. We show that in the case of maximal twist all these nice results can be extended to the more interesting case of a mass non-degenerate quark pair.


title: A doubler-free lattice theory for QCD based on geometric fermions

abstract: We present doubler-free gauge-invariant lattice vector gauge action for some real representations of Wilson gauge fields on an octet of fermions. It is based on a geometric representation of the Dirac equation as an evolution equation on the three-dimensional exterior bundle /(R^3) for a single bispinor and of the bundle (/\x/)(R^3) for an octet. We find doubler-free lattice Dirac operators for above bundles. A gauge-invariant connection with Wilson lattice gauge fields is possible for some real representations of the gauge group. The QCD action of SU(3) is of this type. Application in lattice QCD seems useful: We don't have to waste time and memory for doublers as well as for correction terms to suppress them.


title: Dynamical overlap fermions, results with hybrid Monte-Carlo algorithm

abstract: We present first, exploratory results of a hybrid Monte-Carlo algorithm for dynamical, n_f=2, four-dimensional QCD with overlap fermions. As expected, the computational requirements are typically two orders of magnitude larger for the dynamical overlap formalism than for the more conventional (Wilson or staggered) formulations.


title: Recent development and perspectives of machines for lattice QCD

abstract: I highlight recent progress in cluster computer technology and assess status and prospects of cluster computers for lattice QCD with respect to the development of QCDOC and apeNEXT. Taking the LatFor test case, I specify a 512-processor QCD-cluster better than 1$/Mflops.


title: Axial Correlation Functions in the epsilon-Regime: a Numerical Study with Overlap Fermions

abstract: We present simulation results employing overlap fermions for the axial correlation functions in the epsilon-regime of chiral perturbation theory. In this regime, finite size effects and topology play a dominant role. Their description by quenched chiral perturbation theory is compared to our numerical results in quenched QCD. We show that lattices with a linear extent L > 1.1 fm are necessary to interpret the numerical data obtained in distinct topological sectors in terms of the epsilon-expansion. Such lattices are, however, still substantially smaller than the ones needed in standard chiral perturbation theory. However, we also observe severe difficulties at very low values of the quark mass, in particular in the topologically trivial sector.


title: Application of DWF to heavy-light mesons

abstract: We consider application of domain wall fermions to quarks with relatively heavy masses, aiming at precision calculations of charmed meson properties. Preliminary results for a few basic quantities are presented.


title: The Hamiltonian limit of (3+1)D SU(3) lattice gauge theory on anisotropic lattices

abstract: The extreme anisotropic limit of Euclidean SU(3) lattice gauge theory is examined to extract the Hamiltonian limit, using standard path integral Monte Carlo (PIMC) methods. We examine the mean plaquette and string tension and compare them to results obtained within the Hamiltonian framework of Kogut and Susskind. The results are a significant improvement upon previous Hamiltonian estimates, despite the extrapolation procedure necessary to extract observables. We conclude that the PIMC method is a reliable method of obtaining results for the Hamiltonian version of the theory. Our results also clearly demonstrate the universality between the Hamiltonian and Euclidean formulations of lattice gauge theory. It is particularly important to take into account the renormalization of both the anisotropy, and the Euclidean coupling $ \beta_E $, in obtaining these results.


title: The Study of Chiral Restoration using the Quark-Gluon Mixed Condensate g<\bar{q} \sigma G q> in Lattice QCD at Finite Temperature

abstract: The quark-gluon mixed condensate g<\bar{q} \sigma G q> is studied using SU(3)c lattice QCD with the Kogut-Susskind fermion at the quenched level. Using the lattices as \beta = 6.0 with 16^3 * N_t (N_t=16,12,10,8,6,4), \beta = 6.1 with 20^3 * N_t (N_t=20,12,10,8,6) and \beta = 6.2 with 24^3 * N_t (N_t=24,16,12,10,8) in high statistics of 100-1000 gauge configurations, we perform accurate measurement of the thermal effects on g<\bar{q} \sigma G q> as well as <\bar{q}q> in the chiral limit. We find that the thermal effects on both the condensates are very weak except for the vicinity of T_c, while both the condensates suddenly vanish around T_c \simeq 280 MeV, which indicates strong chiral restoration near T_c. We also find that the ratio m_0^2 = g<\bar{q} \sigma G q>/ <\bar{q}q> is almost independent of the temperature even in the very vicinity of T_c, which means that these two different condensates obey the same critical behavior. This nontrivial similarity between them would impose constraints on the chiral structure of the QCD vacuum near T_c.


title: Static potential, force, and flux-tube profile in 4D compact U(1) lattice gauge theory with the multi-level algorithm

abstract: The long range properties of four-dimensional compact U(1) lattice gauge theory with the Wilson action in the confinement phase is studied by using the multi-level algorithm. The static potential, force and flux-tube profile between two static charges are successfully measured from the correlation function involving the Polyakov loop. The universality of the coefficient of the 1/r correction to the static potential, known as the Luescher term, and the transversal width of the flux-tube profile as a function of its length are investigated. While the result supports the presence of the 1/r correction, the width of the flux tube shows an almost constant behavior at a large distance.


title: Structure functions and form factors close to the chiral limit from lattice QCD

abstract: Results for nucleon matrix elements (arising from moments of structure functions) and form factors from a mixture of runs using Wilson, clover and overlap fermions (both quenched and unquenched) are presented and compared in an effort to explore the size of the chiral `regime', lattice spacing errors and quenching artefacts. While no run covers this whole range of effects the partial results indicate a picture of small lattice spacing errors, small quenching effects and only reaching the chiral regime at rather light quark masses.


title: Analytic Smearing of SU(3) Link Variables in Lattice QCD

abstract: An analytic method of smearing link variables in lattice QCD is proposed and tested. The differentiability of the smearing scheme with respect to the link variables permits the use of modern Monte Carlo updating methods based on molecular dynamics evolution for gauge-field actions constructed using such smeared links. In examining the smeared mean plaquette and the static quark-antiquark potential, no degradation in effectiveness is observed as compared to link smearing methods currently in use, although an increased sensitivity to the smearing parameter is found.


title: Finite Density Lattice Gauge Theories with Positive Fermion Determinants

abstract: We perform simulations of (3-colour) QCD with 2 quark flavours at a finite chemical potential $\mu_I$ for isospin($I_3$), and of 2-colour QCD at a finite chemical potential $\mu$ for quark number. At zero temperature, QCD at finite $\mu_I$ has a mean-field phase transition at $\mu_I=m_\pi$ to a superfluid state with a charged pion condensate which spontaneously breaks $I_3$. We study the finite temperature transition as a function of $\mu_I$. For $\mu_I < m_\pi$, where this is closely related to the transition at finite $\mu$, this appears to be a crossover independent of quark mass, with no sign of the proposed critical endpoint. For $\mu_I > m_\pi$ this becomes a true phase transition where the pion condensate evaporates. For $\mu_I$ just above $m_\pi$ the transition seems to be second order, while for larger $\mu_I$ it appears to become first order. At zero temperature, 2-colour QCD also possesses a superfluid state with a diquark condensate. We study its spectrum of Goldstone and pseudo-Goldstone bosons associated with chiral and quark-number symmetry breaking.


title: Improved gauge action on an anisotropic lattice II - Anisotropy parameter in the medium coupling region -

abstract: The quantum correction of the anisotropy parameter, $\eta$, is calculated for $\xi=2$ and 3 in the $\beta$ region where numerical simulations such as hadron spectroscopy are currently carried out, for the improved actions composed of plaquette and rectangular 6-link loops. The $\beta$ dependences of $\eta$ for the renormalization group improved actions are quite different from those of the standard and Symanzik actions. In Iwasaki and DBW2 actions, $\eta$ stays almost constant in a wide range of $\beta$, which also differs from the one-loop perturbative result, while in the case of Symanzik action, it increases as $\beta$ decreases, which is qualitatively similar to the perturbative result, but the slope is steeper. In the calculation of the $\eta$ parameter close to and in the confined phase, we have applied the link integration method to suppress the fluctuation of the gauge fields. Some technical details are summarized.


title: A Lattice Formulation of Super Yang-Mills Theories with Exact Supersymmetry

abstract: We construct super Yang-Mills theories with extended supersymmetry on hypercubic lattices of various dimensions keeping one or two supercharges exactly. Gauge fields are represented by ordinary unitary link variables, and the exact supercharges are nilpotent up to gauge transformations. Among the models, we show that the desired continuum theories are obtained without any fine tuning of parameters for the cases ${\cal N}=2, 4, 8$ in two-dimensions.


title: Connecting short to long scales in the confining vacuum

abstract: We study approximate decimations in SU(N) LGT that connect the short to long distance regimes. Simple `bond-moving' decimations turn out to provide both upper and lower bounds on the exact partition function. This leads to a representation of the exact partition function in terms of successive decimations whose effective couplings flows are related to those of the easily computable bond-moving decimations. The implications for a derivation of confinement from first principles are discussed.


title: The pion mass in finite volume

abstract: We determine the relative pion mass shift $M_\pi(L)/M_\pi-1$ due to the finite spatial extent $L$ of the box by means of two-flavor chiral perturbation theory and the one-particle L\"uscher formula. We use as input the expression for the infinite volume $\pi\pi$ forward scattering amplitude up to next-to-next-to-leading order and can therefore control the convergence of the chiral series. A comparison to the full leading order chiral expression for the pion mass in finite volume allows us to check the size of subleading terms in the large-$L$ expansion.


title: Lattice calculation of gluon screening masses

abstract: We study SU(3) gluon electric and magnetic masses at finite temperatures using quenched lattice QCD on a $20^2 \times 32 \times 6$ lattice. We focus on temperature regions between $T=T_c$ and $6T_c$, which are realized in BNL Relativistic Heavy Ion Collider and CERN Large Hadron Collider experiments. Stochastic quantization with a gauge-fixing term is employed to calculate gluon propagators. The temperature dependence of the electric mass is found to be consistent with the hard-thermal-loop perturbation, and the magnetic mass has finite values in the temperature region of interest. Both screening masses have little gauge parameter dependence. The behavior of the gluon propagators is very different in confinement/deconfinement physics. The short distance magnetic part behaves like a confined propagator even in the deconfinement phase. A simulation with a larger lattice, $32^2 \times 48 \times 6 $, shows that the magnetic mass has a stronger finite size effect than the electric mass.


title: Numerical Methods for the QCD Overlap Operator: II. Optimal Krylov Subspace Methods

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


title: Matter degrees of freedom and string breaking in Abelian projected quenched SU(2) QCD

abstract: In the Abelian projection the Yang--Mills theory contains Abelian gauge fields (diagonal degrees of freedom) and the Abelian matter fields (off-diagonal degrees) described by a complicated action. The matter fields are essential for the breaking of the adjoint string. We obtain numerically the effective action of the Abelian gauge and the Abelian matter 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. We also observe the breaking of an Abelian analog of the adjoint string using Abelian operators. We show that the adjoint string tension is dominated by the Abelian and the monopole contributions similarly to the case of the fundamental particles. We conclude that the adjoint string breaking can successfully be described in the Abelian projection formalism.


abstract: A quick review of the features of FermiQCD, a C++ library for fast development of parallel Lattice QCD applications. FermiQCD is fully Object Oriented and is based on MPI but, it required no previsius knowledge of parallel programming. It includes parallel algorithms for Wilson, Clover, Staggered, Asqtad and Domain Wall fermions. Optimizations are based on M. Luescher's SSE2 macros.


title: Impact of the finite volume effects on the chiral behavior of fK and BK

abstract: We discuss the finite volume corrections to fK and BK by using the one-loop chiral perturbation theory in full, quenched, and partially quenched QCD. We show that the finite volume corrections to these quantities dominate the physical (infinite volume) chiral logarithms.


title: Center Vortex Model for the Infrared Sector of SU(3) Yang-Mills Theory - Confinement and Deconfinement

abstract: The center vortex model for the infrared sector of Yang-Mills theory, previously studied for the SU(2) gauge group, is extended to SU(3). This model is based on the assumption that vortex world-surfaces can be viewed as random surfaces in Euclidean space-time. The confining properties are investigated, with a particular emphasis on the finite-temperature deconfining phase transition. The model predicts a very weak first order transition, in agreement with SU(3) lattice Yang-Mills theory, and also reproduces a consistent behavior of the spatial string tension in the deconfined phase. The geometrical structure of the center vortices is studied, including vortex branchings, which are a new property of the SU(3) case.


title: Nonperturbative renormalisation of composite operators with overlap quarks

abstract: We compute non-perturbatively the renormalisation constants of composite operators on a $16^3 \times 32 $ lattice with lattice spacing $a$ = 0.093 fm for the overlap fermion action by using the regularisation independent (RI) scheme. The quenched gauge configurations are generated by tadpole improved plaquette plus rectangle action. We test the perturbative continuum relation $Z_A = Z_V$ and $ Z_S=Z_P$ and find that they agree well above $\mu$ = 1.6 GeV. We also perform a Renormalisation Group analysis at the next-to-next-to-leading order and convert the renormalisation constants to the $\bar{MS}$ scheme.


title: Discretization and Continuum Limit of Quantum Gravity on a Four-Dimensional Space-Time Lattice

abstract: The Regge Calculus is a powerful method to approximate a continuous manifold by a simplicial lattice, keeping the connectivities of the underlying lattice fixed and taking the edge lengths as degrees of freedom. The Discrete Regge Model limits the choice of the link lengths to a finite number. We examine the phase structure of Standard Regge Calculus in four dimensions and compare our Monte Carlo results with those of the $Z_2$-Regge Model as well as with another formulation of lattice gravity derived from group theoretical considerations. Within all of the three models of quantum gravity we find an extension of the well-defined phase to negative gravitational couplings and a new phase transition. We calculate two-point functions between geometrical quantities at the corresponding critical point and estimate the masses of the respective interaction particles. A main concern in lattice field theories is the existence of a continuum limit which requires the existence of a continuous phase transition. The recently conjectured second-order transition of the four-dimensional Regge skeleton at negative gravity coupling could be such a candidate. We examine this regime with Monte Carlo simulations and critically discuss its behavior.


title: Chiral perturbation theory for lattice QCD with a twisted mass term

abstract: Quantum Chromodynamics on a lattice with Wilson fermions and a chirally twisted mass term for two degenerate quark flavours is considered in the framework of chiral perturbation theory. The pion mass and decay constant are calculated in next-to-leading order including terms linear in the lattice spacing $a$.


title: The photon propagator in compact QED_{2+1}: the effect of wrapping Dirac strings

abstract: We discuss the influence of closed Dirac strings on the photon propagator in the Landau gauge emerging from a study of the compact U(1) gauge model in 2+1 dimensions. This gauge also minimizes the total length of the Dirac strings. Closed Dirac strings are stable against local gauge-fixing algorithms only due to the torus boundary conditions of the lattice. We demonstrate that these left-over Dirac strings are responsible for the previously observed unphysical behavior of the propagator of space-like photons (D_T) in the deconfinement (high temperature) phase. We show how one can monitor the number N_3 of thermal Dirac strings which allows to separate the propagator measurements into N_3 sectors. The propagator in N_3 \neq 0 sectors is characterized by a non--zero mass and an anomalous dimension similarly to the confinement phase. Both mass squared and anomalous dimension are found to be proportional to N_3. Consequently, in the N_3=0 sector the unphysical behavior of the D_T photon propagator is cured and the deviation from the free massless propagator disappears.


title: Minimalist's Linux Cluster

abstract: Using barebone PC components and NIC's, we construct a linux cluster which has 2-dimensional mesh structure. This cluster has smaller footprint, is less expensive, and use less power compared to conventional linux cluster. Here, we report our experience in building such a machine and discuss our current lattice project on the machine.


title: Two Particle States in an Asymmetric Box

abstract: The exact two-particle energy eigenstates in an asymmetric rectangular box with periodic boundary conditions in all three directions are studied. Their relation with the elastic scattering phases of the two particles in the continuum are obtained. These results can be viewed as a generalization of the corresponding formulae in a cubic box obtained by L\"uscher before. In particular, the s-wave scattering length is related to the energy shift in the finite box. Possible applications of these formulae are also discussed.


title: Teleparallel gravity on the lattice

abstract: We consider quantum gravity model with the squared curvature action. We construct lattice discretization of the model (both on hypercubic and simplicial lattices) starting from its teleparallel equivalent. The resulting lattice models have the actions that are bounded from below while Einstein equations (without matter) appear in their classical limit.


title: Supersymmetry on the lattice: where do predictions and results stand?

abstract: I summarize recent results in lattice supersymmetry with special attention to N=1 Super Yang-Mills (SYM) theory.


title: Chiral extrapolation and physical insights

abstract: It has recently been established that finite-range regularisation in chiral effective field theory enables the accurate extrapolation of modern lattice QCD results to the chiral regime. We review some of the highlights of extrapolations of quenched lattice QCD results, including spectroscopy and magnetic moments. The $\Delta$ resonance displays peculiar chiral features in the quenched theory which can be exploited to demonstrate the presence of significant chiral corrections.


title: Actions for dynamical fermion simulations: are we ready to go?

abstract: A critical review, playing devil's advocate, on present dynamical fermion simulations is given.


title: An Introduction to Lattice Chiral Fermions

abstract: This write-up starts by introducing lattice chirality to people possessing a fairly modern mathematical background, but little prior knowledge about modern physics. I then proceed to present two new and speculative ideas.


title: On practical problems to compute the ghost propagator in SU(2) lattice gauge theory

abstract: In SU(2) lattice pure gauge theory we study numerically the dependence of the ghost propagator G(p) on the choice of Gribov copies in Lorentz (or Landau) gauge. We find that the effect of Gribov copies is essential in the scaling window region, however, it tends to decrease with increasing beta. On the other hand, we find that at larger beta-values very strong fluctuations appear which can make problematic the calculation of the ghost propagator.


title: Lattice Sigma Models with Exact Supersymmetry

abstract: We show how to construct lattice sigma models in one, two and four dimensions which exhibit an exact fermionic symmetry. These models are discretized and {\it twisted} versions of conventional supersymmetric sigma models with N=2 supersymmetry. The fermionic symmetry corresponds to a scalar BRST charge built from the original supercharges. The lattice theories possess local actions and in many cases admit a Wilson term to suppress doubles. In the two and four dimensional theorie s we show that these lattice theories are invariant under additional discrete symmetries. We argue that the presence of these exact symmetries ensures that no fine tuning is required to achieve N=2 supersymmetry in the continuum limit. As a concrete example we show preliminary numerical results from a simulation of the O(3) supersymmetric sigma model in two dimensions.


title: SO(3) vs. SU(2) Yang-Mills theory on the lattice: an investigation at non-zero temperature

abstract: The adjoint SU(2) lattice gauge theory in 3+1 dimensions with the Wilson plaquette action modified by a Z(2) monopole suppression term is reinvestigated with special emphasis on the existence of a finite-temperature phase transition decoupling from the well-known bulk transitions.