title: Lattice simulations of 2-colour QCD with Wilson fermions

abstract: We report on the status of our simulations of two-colour QCD with two flavours of Wilson quarks, at nonzero chemical potential and diquark source. Preliminary results are presented for the static quark potential, gluon propagator, and meson and diquark spectrum.


title: Baryon operators and spectroscopy in lattice QCD

abstract: The construction of the operators and correlators required to determine the excited baryon spectrum is presented, with the aim of exploring the spatial and spin structure of the states while minimizing the number of propagator inversions. The method used to construct operators that transform irreducibly under the symmetries of the lattice is detailed, and the properties of example operators are studied using domain-wall fermion valence propagators computed on MILC asqtad dynamical lattices.


title: Two Nucleons on a Lattice

abstract: The two-nucleon sector is near an infrared fixed point of QCD and as a result the S-wave scattering lengths are unnaturally large compared to the effective ranges and shape parameters. It is usually assumed that a lattice QCD simulation of the two-nucleon sector will require a lattice that is much larger than the scattering lengths in order to extract quantitative information. In this paper we point out that this does not have to be the case: lattice QCD simulations on much smaller lattices will produce rigorous results for nuclear physics.


title: Effective String Picture for Confining Gauge Theories at Finite Temperature

abstract: We study the ``effective string picture'' of confinement, deriving theoretical predictions for the interquark potential at finite temperature. At low temperatures, the leading string correction to the linear confining potential between a heavy quark-antiquark pair is the "L\"uscher term''. Assuming a Nambu-Goto effective string action, subleading contributions can be worked out in an analytical way. We also discuss the contribution given by a possible ``boundary term'' in the effective action, and compare these predictions with high precision results from simulations of lattice Z(2) gauge theory in three dimensions, obtained with an algorithm that exploits the duality of the Z(2) gauge model with the Ising spin model.


title: Study of QCD thermodynamics at finite density by Taylor expansion

abstract: We discuss the phase structure and the equation of state for 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 for 2-flavor QCD, enabling estimates of the pressure, quark number density and associated susceptibilities as functions of $\mu_q$ via a Taylor series expansion. Also, the phase transition line for 2 and 3-flavor QCD and the critical endpoint in the $(T, \mu_q)$ plane are investigated in the low density regime.


title: Excited B mesons from the lattice

abstract: We determine the energies of the excited states of a heavy-light meson $Q\bar{q}$, with a static heavy quark and light quark with mass approximately that of the strange quark from both quenched lattices and with dynamical fermions. We are able to explore the energies of orbital excitations up to L=3, the spin-orbit splitting up to L=2 and the first radial excitation. These $b \bar{s}$ mesons will be very narrow if their mass is less than 5775 MeV -- the $BK$ threshold. We investigate this in detail and present evidence that the scalar meson (L=1) will be very narrow and that as many as 6 $b \bar{s}$ excited states will have energies close to the $BK$ threshold and should also be relatively narrow.


title: Mesonic Form Factors

abstract: We have started a program to compute the electromagnetic form factors of mesons. We discuss the techniques used to compute the pion form factor and present preliminary results computed with domain wall valence fermions on MILC asqtad lattices, as well as Wilson fermions on quenched lattices. These methods can easily be extended to rho-to-gamma-pi transition form factors.


title: Anisotropic lattices for precision computations in heavy flavor physics

abstract: We study the anisotropic lattice QCD for precision computations of heavy-light matrix elements. Our previous study in which the lattices are calibrated with a few percent accuracy has already given results comparable to the existing calculations. This suggests that even higher precision may be achieved by a more precise calibration of anisotropic lattices. We describe our strategy to tune the gauge and quark parameters with accuracies much less than 1 % in the quenched approximation.


title: Numerical study of staggered fermion on anisotropic lattices

abstract: We study calibration procedures of the staggered quarks on anisotropic lattices in the quenched approximation and in $N_f = 2$ dynamical simulations. For the calibration conditions we adopt the hadronic radii $r_0$ and the meson masses in the temporal and spatial directions. On the quenched lattice, we calibrate the quark field and compare the result with the result determined using the meson dispersion relation. In dynamical simulations, we determine the anisotropy parameters $\gamma_G$ and $\gamma_F$ simultaneously within 1% accuracy at renormalized anisotropy $\xi = a_s/a_t = 4$.


title: Non-perturbative calculation of $Z_V$ and $Z_A$ in domain-wall QCD on a finite box

abstract: We report on a non-perturbative evaluation of the renormalization factors for the vector and axial-vector currents, $Z_V$ and $Z_A$, in the quenched domain-wall QCD (DWQCD) with plaquette and renormalization group improved gauge actions. We take the Dirichlet boundary condition for both gauge and domain-wall fermion fields on the finite box, and introduce the flavor-chiral Ward-Takahashi identities to calculate the renormalization factors. As a test of the method, we numerically confirm the expected relation that $Z_V \simeq Z_A$ in DWQCD. Employing two different box sizes for the numerical simulations at several values of the gauge coupling constant $g^2$ and the domain-wall height $M$, we extrapolate $Z_V$ to the infinite volume to remove $a/L$ errors. We finally give the interpolation formula of $Z_V$ in the infinite volume as a function of $g^2$ and $M$.


title: Low-energy couplings of QCD from topological zero-mode wave functions

abstract: By matching 1/m^2 divergences in finite-volume two-point correlation functions of the scalar or pseudoscalar densities with those obtained in chiral perturbation theory, we derive a relation between the Dirac operator zero-mode eigenfunctions at fixed non-trivial topology and the low-energy constants of QCD. We investigate the feasibility of using this relation to extract the pion decay constant, by computing the zero-mode correlation functions on the lattice in the quenched approximation and comparing them with the corresponding expressions in quenched chiral perturbation theory.


title: Scaling test for Wilson twisted mass QCD

abstract: We present a first scaling test of twisted mass QCD with pure Wilson quarks for a twisting angle of pi/2. We have computed the vector meson mass and the pseudoscalar decay constant for different values of beta at fixed value of r_0 m_PS. The results obtained in the quenched approximation are compared with data for pure Wilson and non-perturbatively O(a) improved Wilson computations. We show that our results from Wilson twisted mass QCD show clearly reduced lattice spacing errors, consistent with O(a) improvement and without the need of any improvement terms added. These results thus provide numerical evidence of the prediction in ref. [1].


title: Transverse Structure of Nucleon Parton Distributions from Lattice QCD

abstract: This work presents the first calculation in lattice QCD of three moments of spin-averaged and spin-polarized generalized parton distributions in the proton. It is shown that the slope of the associated generalized form factors decreases significantly as the moment increases, indicating that the transverse size of the light-cone quark distribution decreases as the momentum fraction of the struck parton increases.


title: Heavy quark free energies and the renormalized Polyakov loop in full QCD

abstract: We study the renormalized free energy of a heavy quark anti-quark pair in the different colour channels in full QCD at finite temperature. Similarities and differences to the quenched case are discussed and the temperature dependence as well as their short distance behavior are analyzed. The asymptotic large distance behavior of the free energy is used to define the non-perturbatively renormalized Polyakov loop which is well behaved in the continuum limit.


title: Effective heavy-light meson energies in small-volume quenched QCD

abstract: We study effective energies of heavy-light meson correlation functions in lattice QCD and a small volume of (0.2 fm)^4 to non-perturbatively calculate their dependence on the heavy quark mass in the continuum limit. Our quenched results obtained here constitute an essential intermediate step of a first fully non-perturbative computation of the b-quark's mass in the static approximation that has recently been presented as an application of a new proposal to non-perturbatively renormalize the Heavy Quark Effective Theory. The renormalization constant and the improvement coefficients relating the renormalized current and subtracted quark mass are determined in the relevant parameter region at weak couplings, which allows to perform the numerical simulations at several, precisely fixed values of the renormalization group invariant heavy quark mass in a range from 3 GeV to 15 GeV.


title: One-Loop Matching of the Heavy-Light A_0 and V_0 Currents with NRQCD Heavy and Improved Naive Light Quarks

abstract: One-loop matching of heavy-light currents is carried out for a highly improved lattice action, including the effects of dimension 4 O(1/M) and O(a) operators. We use the NRQCD action for heavy quarks, the Asqtad improved naive action for light quarks, and the Symanzik improved glue action. As part of the matching procedure we also present results for the NRQCD self energy and for massless Asqtad quark wavefunction renormalization with improved glue.


title: Spontaneous violation of CP symmetry in the strong interactions

abstract: Some time ago Dashen pointed out that spontaneous CP violation can occur in the strong interactions. I show how a simple effective Lagrangian exposes the remarkably large domain of quark mass parameters for which this occurs. I close with some warnings for lattice simulations.


title: Excitations of the static quark-antiquark system in several gauge theories

abstract: The spectrum of gluons in the presence of a static quark-antiquark pair is studied using Monte Carlo simulations on anisotropic space-time lattices. For very small quark-antiquark separations R, the level orderings and approximate degeneracies disagree with the expectations from an effective string theory. As the quark-antiquark separation R increases, a dramatic rearrangement of the energies occurs, and above 2 fm, all of the levels studied show behavior consistent with an effective string description. The energy spacings are nearly pi/R, but a tantalizing fine structure remains. In addition to 4-dimensional SU(3) gauge theory, results from 3-dimensional SU(2) and compact U(1) gauge theories are also presented.


title: Deconstruction, 2d lattice Yang-Mills, and the dynamical lattice spacing

abstract: We study expectation values related to the dynamical lattice spacing that occurs in the recent supersymmetric 2d lattice Yang-Mills constructions of Cohen et al. [hep-lat/0307012]. For the purposes of this preliminary analysis, we restrict our attention to the bosonic part of that theory. That is, we compute observables in the fully quenched ensemble, equivalent to non-supersymmetric 2d lattice Yang-Mills with a dynamical lattice spacing and adjoint scalars. Our numerical simulations indicate difficulties with the proposed continuum limit. We find that expectation values tend to those of the undeformed ``daughter theory,'' in spite of the deformation suggested by Cohen et al. In an effort to understand these results, we examine the zero action configurations, with and without the deformation. Based on these considerations, we are able to interpret the simulation results in terms of entropic effects.


title: General chiral gauge theories

abstract: Only requiring that Dirac operators decribing massless fermions on the lattice decompose into Weyl operators we arrive at a large class of them. After deriving general relations from spectral representations we study correlation functions of Weyl fermions for any value of the index, stressing the related conditions for basis transformations and getting the precise behaviors under gauge and CP transformations. Using the detailed structure of the chiral projections we also obtain a form of the correlation functions with a determinant in the general case.


title: The Deconfinement Phase Transition of Sp(2) and Sp(3) Yang-Mills Theories in 2+1 and 3+1 Dimensions

abstract: Some time ago, Svetitsky and Yaffe have argued that -- if the deconfinement phase transition of a (d+1)-dimensional Yang-Mills theory with gauge group G is second order -- it should be in the universality class of a d-dimensional spin model symmetric under the center of G. For d=3 these arguments have been confirmed numerically only in the SU(2) case with center Z(2), simply because all SU(N) Yang-Mills theories with N>=3 seem to have non-universal first order phase transitions. The symplectic groups Sp(N) also have the center Z(2) and provide another extension of SU(2) = Sp(1) to general N. Using lattice simulations, we find that the deconfinement phase transition of Sp(2) Yang-Mills theory is first order in 3+1 dimensions, while in 2+1 dimensions stronger fluctuations induce a second order transition. In agreement with the Svetitsky-Yaffe conjecture, for (2+1)-d Sp(2) Yang-Mills theory we find the universal critical behavior of the 2-d Ising model. For Sp(3) Yang-Mills theory the transition is first order both in 2+1 and in 3+1 dimensions. This suggests that the size of the gauge group -- and not the center symmetry -- determines the order of the deconfinement phase transition.


title: Pion structure from improved lattice QCD: form factor and charge radius at low masses

abstract: The charge form factor of the pion is calculated in lattice QCD. The non-perturbatively improved Sheikholeslami-Wohlert action is used together with the $\mathcal{O}(a)$ improved vector current. Other choices for the current are examined. The form factor is extracted for pion masses from 970 MeV down to 360 MeV and for momentum transfers $Q^2 \leq 2 \mathrm{GeV}^2$. The mean square charge radius is extracted, compared to previous determinations and its extrapolation to lower masses discussed.


title: A note on the power divergence in lattice calculations of $\Delta I = 1/2$ $K\to\pi\pi$ amplitudes at $M_{K}=M_{\pi}$

abstract: In this note, we clarify a point concerning the power divergence in lattice calculations of $\Delta I = 1/2$ $K\to\pi\pi$ decay amplitudes. There have been worries that this divergence might show up in the Minkowski amplitudes at $M_{K}=M_{\pi}$ with all the mesons at rest. Here we demonstrate, via an explicit calculation in leading-order Chiral Perturbation Theory, that the power divergence is absent at the above kinematic point, as predicted by CPS symmetry.


abstract: A simplified test of universality in Lattice QCD is performed by analytically evaluating the continuous Euclidean time limits of various lattice fermion determinants, both with and without a Wilson term to lift the fermion doubling on the Euclidean time axis, and comparing them with each other and with the zeta-regularised fermion determinant in the continuous time--lattice space setting. The determinant relations expected from universality considerations are found to be violated by a certain gauge field-dependent factor, i.e. we uncover a "universality anomaly". The physical significance, or lack thereof, of this factor is a delicate question which remains to be settled.


title: The Calculation of Critical Parameters in SU(2) Gauge Theory with Kouvel-Fisher Method

abstract: We calculate the critical coupling $4/g_c^2$ and critical exponent $\beta$ for the order parameter in SU(2) lattice gauge theory by applying of the finite size scaling technique and the method proposed by Kouvel and Fisher for analysis of experimental data. In contrast to the standard finite size scaling approach, this method allows to determine simultaneously both $g_c^2/4$ and $\beta$ as two parameters of the linear fit to the Monte-Carlo data.


title: A Finite Baryon Density Algorithm

abstract: I will review the progress toward a finite baryon density algorithm in the canonical ensemble approach which entails particle number projection from the fermion determinant. These include an efficient Pad\'{e}-Z$_2$ stochastic estimator of the $Tr \log$ of the fermion matrix and a Noisy Monte Carlo update to accommodate unbiased estimate of the probability. Finally, I will propose a Hybrid Noisy Monte Carlo algorithm to reduce the large fluctuation in the estimated $Tr \log$ due to the gauge field which should improve the acceptance rate. Other application such as treating $u$ and $d$ as two separate flavors is discussed.


title: Large N QCD -- Continuum reduction

abstract: Numerical evidence combined with Eguchi-Kawai reduction indicate that there are no finite volumes effects in the large N limit of QCD as long as the linear extent of the four-torus is bigger than a critical size. This is referred to as continuum reduction. Since fermions in the fundamental representation are naturally quenched in the large N limit, as long as we only have a finite number of flavors, continuum reduction provides us with the exciting possibility to numerically solve large N QCD using chiral fermions and present day computers.


title: Calculation of Finite Size Effects on the Nucleon Mass in Unquenched QCD using Chiral Perturbation Theory

abstract: The finite size effects on nucleon masses are calculated in relativistic chiral perturbation theory. Results are compared with two-flavor lattice results.


title: The nucleon mass in N_f=2 lattice QCD: finite size effects from chiral perturbation theory

abstract: In the framework of relativistic SU(2)_f baryon chiral perturbation theory we calculate the volume dependence of the nucleon mass up to and including O(p^4). Since the parameters in the resulting finite size formulae are fixed from the pion mass dependence of the large volume nucleon masses and from phenomenology, we obtain a parameter-free prediction of the finite size effects. We present mass data from the recent N_f=2 simulations of the UKQCD and QCDSF collaborations and compare these data as well as published mass values from the dynamical simulations of the CP-PACS and JLQCD collaborations with the theoretical expectations. Remarkable agreement between the lattice data and the predictions of chiral perturbation theory in a finite volume is found.


title: Detailed Lattice-QCD Study for the Three-Quark Potential and Y-type Flux-Tube Formation

abstract: We study the ground-state three-quark (3Q) potential $V_{\rm 3Q}^{\rm g.s.}$ and the excited-state 3Q potential $V_{\rm 3Q}^{\rm e.s.}$ using SU(3) lattice Q CD. From the accurate and thorough calculation for more than 300 different patterns of 3Q systems, the static ground-state 3Q potential $V_{\rm 3Q}^{\rm g.s.}$ is found to be well described by the Coulomb plus Y-type linear potential, i.e., Y-ansatz, within 1%-level deviation. With lattice QCD, we calculate also the excited-state potential in the 3Q system, and find the gluonic excitation energy, $V_{\rm 3Q}^{\rm e.s.}-V_{\rm 3Q}^{\rm g.s.}$, to be about 1 GeV. This large gluonic-excitation energy would play an essential role to the success of the quark model for the low-lying hadrons in terms of the absence of the gluonic mode.


title: Vacuum Polarisation and Hadronic Contribution to muon g-2 from Lattice QCD

abstract: We compute the vacuum polarisation on the lattice in quenched QCD using non-perturbatively improved Wilson fermions. Above Q^2 of about 2 GeV^2 the results are very close to the predictions of perturbative QCD. Below this scale we see signs of non-perturbative effects which we can describe by the use of dispersion relations. We use our results to estimate the light quark contribution to the muon's anomalous magnetic moment. We find the result 446(23) x 10^{-10}, where the error only includes statistical uncertainties. Finally we make some comments on the applicability of the Operator Product Expansion to our data.


title: Negative-Parity Baryons in Quenched Anisotropic Lattice QCD

abstract: We study negative-parity baryon spectra in quenched anisotropic lattice QCD. The negative-parity baryons are measured as the parity partner of the ground-state baryons. In addition to the flavor octet and decuplet baryons, we pay much attention to the flavor-singlet negative-parity baryon as a three-quark state and compare it with the Lambda(1405) baryon. Numerical results of the flavor octet and decuplet negative-parity baryon masses are close to experimental values of lowest-lying negative-parity baryons, while the flavor-singlet baryon is much heavier than Lambda(1405). This indicates that the Lambda(1405) would be a multi-quark state such as the N-Kbar molecule rather than the flavor-singlet 3 quark state.


title: The Yang-Mills Spectrum from a 2-level Algorithm

abstract: We investigate in detail a 2-level algorithm for the computation of 2-point functions of fuzzy Wilson loops in lattice gauge theory. Its performance and the optimization of its parameters are described in the context of 2+1D SU(2) gluodynamics. In realistic calculations of glueball masses, it is found that the reduction in CPU time for given error bars on the correlator at time-separation ~0.2fm, where a mass-plateau sets in, varies between 1.5 and 7 for the lightest glueballs in the non-trivial symmetry channels; only for the lightest glueball is the 2-level algorithm not helpful. For the heavier states, or for larger time-separations, the gain increases as expected exponentially in (mt). We present further physics applications in 2+1 and 3+1 dimensions and for different gauge groups that confirm these conclusions.


title: Topological Charge Evolution in the Markov-Chain of QCD

abstract: The topological charge is studied on lattices of large physical volume and fine lattice spacing. We illustrate how a parity transformation on the SU(3) link-variables of lattice gauge configurations reverses the sign of the topological charge and leaves the action invariant. Random applications of the parity transformation are proposed to traverse from one topological charge sign to the other. The transformation provides an improved unbiased estimator of the ensemble average and is essential in improving the ergodicity of the Markov chain process.


title: Propagators and running coupling from SU(2) lattice gauge theory

abstract: We perform numerical studies of the running coupling constant alpha_R(p^2) and of the gluon and ghost propagators for pure SU(2) lattice gauge theory in the minimal Landau gauge. Different definitions of the gauge fields and different gauge-fixing procedures are used respectively for gaining better control over the approach to the continuum limit and for a better understanding of Gribov-copy effects. We find that the ghost-ghost-gluon-vertex renormalization constant is finite in the continuum limit, confirming earlier results by all-order perturbation theory. In the low momentum regime, the gluon form factor is suppressed while the ghost form factor is divergent. Correspondingly, the ghost propagator diverges faster than 1/p^2 and the gluon propagator appears to be finite. Precision data for the running coupling alpha_R(p^2) are obtained. These data are consistent with an IR fixed point given by lim_{p \to 0} alpha_R(p^2) = 5(1).


title: Behavior of Charmonium Systems after Deconfinement

abstract: We present a study of charmonia in hot gluonic plasma, for temperatures upto three times the deconfinement transition temperature Tc. q \bar{q} systems with quark masses close to the charm mass and different spin-parity quantum numbers were studied on very fine isotropic lattices. The analysis of temporal correlators, and spectral functions constructed from them, shows that the J/psi and eta_c survive up to quite high temperatures, with little observable change up to 1.5 Tc, and then gradually weakening and disappearing by 3 Tc. For the scalar and axial vector channels, serious modifications are induced by the hot medium already close to Tc, possibly dissociating the mesons by 1.1 Tc.


title: Density of States Method at Finite Isospin Density

abstract: The density of states method is applied for lattice QCD at a finite isospin density. The advantage of this method is that one can easily obtain results for various values of parameters (quark mass, coupling constant and the number of flavors). We compare results for the chiral condensate and the quark number density with those from the R-algorithm and find that they are in good agreement. By calculating the chiral condensate we obtain information on the phase structure for various quark flavors and isospin chemical potentials. We also show results for the chiral condensate at two different quark masses and at two different isospin densities which are not easily obtainable in the conventional Monte Carlo method.


title: A new scheme for the Klein-Gordon and Dirac fields on the lattice with axial anomaly

abstract: Using the method of finite differences a scheme is proposed to solve exactly the Klein-Gordon and Dirac free field equations, in a (1+1)-dimensional lattice. The hamiltonian of the Dirac field is translational invariant, hermitian, avoids fermion doubling, and, for the massless case, preserves global chiral symmetry. Coupling the fermion field to the electromagnetic vector potential we construct a gauge invariant vector current leading to the correct axial anomaly.


title: Properties of Thermal Glueballs

abstract: We study the properties of the 0++ glueball at finite temperature using SU(3) quenched lattice QCD. We find a significant thermal effects near T_c. We perform the \chi^2 fit analyses adopting two Ansaetze for the spectral function, i.e., the conventional narrow-peak Ansatz and an advanced Breit-Wigner Ansatz. The latter is an extension of the former, taking account of the appearance of the thermal width at T>0. We also perform the MEM analysis. These analyses indicate that the thermal effect on the glueball is a significant thermal-width broadening \Gamma(T_c) \sim 300 MeV together with a modest reduction in the peak center \Delta\omega_0(T_c) \sim 100 MeV.


title: A relativistic invariant scheme for the quantum Klein-Gordon and Dirac fields on the lattice

abstract: We explore the connection between the symmetry transformations and conservation laws for the Klein-Gordon and Dirac fields on the lattice. The generators of the space time translations and Lorentz boost (defined on the lattice) are constants of motion and satisfy the standard commutation relations.


title: Representations of classical groups on the lattice and its application to the field theory on discrete space-time

abstract: We explore the mathematical consequences of the assumption of a discrete space-time. The fundamental laws of physics have to be translated into the language of discrete mathematics. We find integral transformations that leave the lattice of any dimension invariant and apply these transformations to field equations.


title: Primordial function and ambiguity in its determination

abstract: We study the possibility to reconstruct the primordial function for some periodic function. The procedure includes an analytical continuation of a discrete function for Fourier coefficients computation, that introduces an ambiguity. To establish conditions under which no ambiguity appears, the modification of Carlson theorem is suggested. In a case when no subsidiary condition is imposed, ambiguity in the definition of primordial function does not go beyond the functions which belong to space $\Omega^{\prime}$.


title: A realistic interpretation of lattice gauge theories

abstract: Following recent assumptions to unify quantum mechanics and general relativity, the structure of spacetime is suppose to be a consequence of the relations among some fundamental objects, and its concept can be formulated without the reference to the intuition. As physical consequences the continuous laws should be translated in to difference equations and the lattice field theories should be interpreted as a realistic model.


title: Discrete Differential Geometry and Lattice Field Theory

abstract: We develope a difference calculus analogous to the differential geometry by translating the forms and exterior derivatives to similar expressions with difference operators, and apply the results to fields theory on the lattice [Ref. 1]. Our approach has the advantage with respect to other attempts [Ref. 2-6] that the Lorentz invariance is automatically preserved as it can be seen explicitely in the Maxwell, Klein-Gordon and Dirac equations on the lattice.


title: Y-type Flux-Tube Formation in Baryons

abstract: For more than 300 different patterns of the 3Q systems, the ground-state 3Q potential $V_{\rm 3Q}^{\rm g.s.}$ is investigated using SU(3) lattice QCD with $12^3\times 24$ at $\beta=5.7$ and $16^3\times 32$ at $\beta=5.8, 6.0$ at the quenched level. As a result of the detailed analyses, we find that the ground-state potential $V_{\rm 3Q}^{\rm g.s.}$ is well described with so-called Y-ansatz as $V_{\rm 3Q}=-A_{\rm 3Q}\sum_{i<j}\frac1{|{\bf r}_i-{\bf r}_j|} +\sigma_{\rm 3Q} L_{\rm min}+C_{\rm 3Q}$, with the accuracy better than 1%. Here, $L_{\rm min}$ denotes the minimal value of total flux-tube length. We also studythe excited-state potential $V_{\rm 3Q}^{\rm e.s.}$ using lattice QCD with $16^3\times 32$ at $\beta=5.8, 6.0$ for more than 100 patterns of the 3Q systems. The energy gap between $V_{\rm 3Q}^{\rm g.s.}$ and $V_{\rm 3Q}^{\rm e.s.}$, which physically means the gluonic excitation energy, is found to be about 1 GeV in the typical hadronic scale. Finally, we suggest a possible scenario which connects the success of the quark model to QCD.


title: Connecting Lattice QCD with Chiral Perturbation Theory at Strong Coupling

abstract: We study the difficulties associated with detecting chiral singularities predicted by chiral perturbation theory (ChPT) in lattice QCD. We focus on the physics of the remnant O(2) chiral symmetry of staggered fermions in the strong coupling limit using the recently discovered directed path algorithm. Since it is easier to look for power-like singularities as compared to logarithmic ones, our calculations are performed at a fixed finite temperature in the chirally broken phase. We show that the behavior of the chiral condensate, the pion mass and the pion decay constant, for small masses, are all consistent with the predictions of ChPT. However, the values of the quark masses that we need to demonstrate this are much smaller than those being used in dynamical QCD simulations. We also need to use higher order terms in the chiral expansion to fit our data.


title: Coulomb Energy, Remnant Symmetry, and the Phases of Non-Abelian Gauge Theories

abstract: We show that the confining property of the one-gluon propagator, in Coulomb gauge, is linked to the unbroken realization of a remnant gauge symmetry which exists in this gauge. An order parameter for the remnant gauge symmetry is introduced, and its behavior is investigated in a variety of models via numerical simulations. We find that the color-Coulomb potential, associated with the gluon propagator, grows linearly with distance both in the confined and - surprisingly - in the high-temperature deconfined phase of pure Yang-Mills theory. We also find a remnant symmetry-breaking transition in SU(2) gauge-Higgs theory which completely isolates the Higgs from the (pseudo)confinement region of the phase diagram. This transition exists despite the absence, pointed out long ago by Fradkin and Shenker, of a genuine thermodynamic phase transition separating the two regions.


title: A partially quenched analysis of the eta-eta' system in N_f=2 QCD

abstract: We report on a first, comprehensive partially quenched study of the eta-eta' problem, based on SESAM configurations on a 16^3x32 lattice at beta=5.6 QCD with two (mass degenerate) active sea quark flavours. By means of the spectral approximation of the two-loop (hairpin) diagrams, we find clear plateau formation in the effective masses which enables us both to determine the eta-eta' mass matrix and the alpha-parameter in the effective chiral Lagrangian for the flavour singlet sector, alpha=0.028 +- 0.013.


title: Computing the eta and eta' Mesons in Lattice QCD

abstract: It has been known for a long time that the large experimental singlet-octet mass gap in the pseudoscalar meson mass spectrum originates from the anomaly of the axial vector current, i.e. from nonperturbative effects and the nontrivial topological structure of the QCD vacuum. In the N_colour -> infinity limit of the theory, this connection elucidates in the famous Witten-Veneziano relation between the eta'-mass and the topological susceptibility of the quenched QCD vacuum.While lattice QCD has by now produced impressive high precision results on the flavour nonsinglet hadron spectrum, the determination of the pseudoscalar singlet mesons from direct correlator studies is markedly lagging behind, due to the computational complexity in handling observables that include OZI-rule violating diagrams, like the eta' propagator. In this article, we report on some recent progress in dealing with the numerical bottleneck problem.


title: Discrete reflection groups and induced representations of Poincare group on the lattice

abstract: We continue the program, presented in previous Symposia, of discretizing physical models. In particular we calculate the integral Lorentz transformations with the help of discrete reflection groups, and use them for the covariance of Klein-Gordon and Dirac wave equation on the lattice. Finally we define the unitary representation of Poincar group on discrete momentum and configuration space, induced by integral representations of its closed subgroup.


title: Least-Squared Optimized Polynomials for Smeared Link Actions

abstract: We introduce a numerical method for generating the approximating polynomials used in fermionic calculations with smeared link actions. We investigate the stability of the algorithm and determine the optimal weight function and the optimal type of discretization. The achievable order of polynomial approximation reaches several thousands allowing fermionic calculations using the Hypercubic Smeared Link action even with physical quark masses.


title: Induced representations of Poincare group on the lattice: spin 1/2 and 1 case

abstract: Following standard methods we explore the construction of the discrete Poincare group, the semidirect product of discrete translations and integral Lorentz transformations, using the Wigner-Mackey construction restricted to the momentum and position space on the lattice. The orbit condition, irreducibility and assimptotic limit are discussed.


title: Numerical study of staggered quark action on quenched anisotropic lattices

abstract: The staggered quark action on anisotropic lattices is studied. We carry out numerical simulations in the quenched approximation at three values of lattice spacing ($a_{\sigma}^{-1}=1-2$ GeV) with the anisotropy $\xi= a_{\sigma}/a_{\tau}=4$, where $a_{\sigma}$ and $a_{\tau}$ are the spatial and temporal lattice spacings, respectively. The bare anisotropy $\gamma_F$ in the quark action is numerically tuned through the ratio of meson masses in the fine and coarse directions, and through the dispersion relation of a meson, so that the renormalized fermionic anisotropy coincides with that of the gauge field. The discrepancy between these two calibration schemes provides an estimate of the finite lattice artifact, which is found to be sizable in the range of cutoff explored in this work. We also compute the meson masses using correlators with the wall source at the tuned anisotropy parameter. The flavor symmetry breaking effect smoothly decreases as $\beta$ increases. The effect of uncertainty in $\gamma_F$ on the meson masses are examined. We also discuss a perspective on dynamical simulations.


title: Charmonium correlators at finite temperature in quenched lattice QCD

abstract: We study charmonium correlators at finite temperature using quenched lattice QCD simulations. Two analysis procedures are applied to extract information on the spectral function: the maximum entropy method, and the $\chi^2$ fit analyses including the constrained curve fitting. We focus on the low energy structure of the spectral function by applying the smearing technique. We first discuss the applicability of these methods to the finite temperatures by analyzing the data at T=0 with restricted numbers of degrees of freedom. Then we apply these methods to the correlators at $T>0$. We find no indication of mass shift and finite width for the charmonium states at $T\simeq 0.9T_c$. The results at $T\simeq 1.1 T_c$ imply that bound-state-like structures may survive even above $T_c$.


title: Improving meson two-point functions in lattice QCD

abstract: We describe and test a method to compute Euclidean meson two-point functions in lattice QCD. The contribution from the low-lying eigenmodes of the Dirac operator is averaged over all positions of the quark sources. The contribution from the higher modes is estimated in the traditional way with one or a few source points per lattice. In some channels, we observe a significant improvement in the two-point functions for small quark masses.


title: Remarks on the multi-parameter reweighting method for the study of lattice QCD at non-zero temperature and density

abstract: We comment on the reweighting method for the study of finite density lattice QCD. We discuss the applicable parameter range of the reweighting method for models which have more than one simulation parameter. The applicability range is determined by the fluctuations of the modification factor of the Boltzmann weight. In some models having a first order phase transition, the fluctuations are minimized along the phase transition line if we assume that the pressure in the hot and the cold phase is balanced at the first order phase transition point. This suggests that the reweighting method with two parameters is applicable in a wide range for the purpose of tracing out the phase transition line in the parameter space. To confirm the usefulness of the reweighting method for 2 flavor QCD, the fluctuations of the reweighting factor are measured by numerical simulations for the cases of reweighting in the quark mass and chemical potential directions. The relation with the phase transition line is discussed. Moreover, the sign problem caused by the complex phase fluctuations is studied.


title: Exact solutions of the Dirac equation and induced representations of the Poincare group on the lattice

abstract: We deduce the structure of the Dirac field on the lattice from the discrete version of differential geometry and from the representation of the integral Lorentz transformations. The analysis of the induced representations of the Poincare group on the lattice reveals that they are reducible, a result that can be considered a group theoretical approach to the problem of fermion doubling.


title: Finite Temperature QCD with Two Flavors of Non-perturbatively Improved Wilson Fermions

abstract: We study QCD with two flavors of non-perturbatively improved Wilson fermions at finite temperature on the $16^3 8$ lattice. We determine the transition temperature at lattice spacings as small as $a \sim 0.12$ fm, and study string breaking below the finite temperature transition. We find that the static potential can be fitted by a two-state ansatz, including a string state and a two-meson state. We investigate the role of Abelian monopoles at finite temperature.


title: Discretisation effects in the topological susceptibility in lattice QCD

abstract: We study the topological susceptibility, chi, in QCD with two quark flavours using lattice field configurations that have been produced with an O(a)-improved clover quark action. We find clear evidence for the expected suppression at small quark mass, m_q, and examine the variation of $\chi$ with this mass and the lattice spacing, a. A joint continuum and chiral extrapolation yields good agreement with theoretical expectations as (a, m_q) -> 0. A moderate increase in autocorrelation is observed on the more chiral ensembles, but within large statistical errors. Finite volume effects are negligible for Leutwyler-Smilga parameter x_LS > 10, and no evidence for a nearby phase transition is observed.


title: Equation of state at finite temperature and chemical potential, lattice QCD results

abstract: We present an N_t=4 lattice study for the equation of state of 2+1 flavour staggered, dynamical QCD at finite temperature and chemical potential. We use the overlap improving multi-parameter reweighting technique to extend the equation of state for non-vanishing chemical potentials. The results are obtained on the line of constant physics and our physical parameters extend in temperature and baryon chemical potential upto \approx 500-600 MeV.


title: Super Yang-Mills Theories on the Two-Dimensional Lattice with Exact Supersymmetry

abstract: We construct super Yang-Mills theories with ${\cal N}=2, 4$ supersymmetries on the two-dimensional square lattice keeping one or two supercharges exactly. Along the same line as the previous paper \cite{sugino}, the construction is based on topological field theory formulation. In order to resolve the problem of degenerate classical vacua encountered in the previous paper, we consider two kinds of modifications of the action. For one of them, any supersymmetry breaking terms do not need to be introduced, and the formulations exactly realize some of supersymmetries at the lattice level. Our lattice actions flow to the desired continuum theories without any fine-tuning of parameters.


title: Numerical Portrait of a Relativistic BCS Gapped Superfluid

abstract: We present results of numerical simulations of the 3+1 dimensional Nambu - Jona-Lasinio (NJL) model with a non-zero baryon density enforced via the introduction of a chemical potential mu not equal to 0. The triviality of the model with a number of dimensions d>=4 is dealt with by fitting low energy constants, calculated analytically in the large number of colors (Hartree) limit, to phenomenological values. Non-perturbative measurements of local order parameters for superfluidity and their related susceptibilities show that, in contrast to the 2+1 dimensional model, the ground-state at high chemical potential and low temperature is that of a traditional BCS superfluid. This conclusion is supported by the direct observation of a gap in the dispersion relation for 0.5<=(mu a)<=0.85, which at (mu a)=0.8 is found to be roughly 15% the size of the vacuum fermion mass. We also present results of an initial investigation of the stability of the BCS phase against thermal fluctuations. Finally, we discuss the effect of splitting the Fermi surfaces of the pairing partners by the introduction of a non-zero isospin chemical potential.


title: Representations of the discrete inhomogeneous Lorentz group and Dirac wave equation on the lattice

abstract: We propose the fundamental and two dimensional representation of the Lorentz groups on a (3+1)-dimensional hypercubic lattice, from which representations of higher dimensions can be constructed. For the unitary representation of the discrete translation group we use the kernel of the Fourier transform. From the Dirac representation of the Lorentz group (including reflections) we derive in a natural way the wave equation on the lattice for spin 1/2 particles. Finally the induced representation of the discrete inhomogeneous Lorentz group is constructed by standard methods and its connection with the continuous case is discussed.


title: The deconfining phase transition in full QCD with two dynamical flavors

abstract: We investigate the deconfining phase transition in SU(3) pure gauge theory and in full QCD with two flavors of staggered fermions. The phase transition is detected by measuring the free energy in presence of an abelian monopole background field. In the pure gauge case our finite size scaling analysis is in agreement with the well known presence of a weak first order phase transition. In the case of 2 flavors full QCD we find, using the standard pure gauge and staggered fermion actions, that the phase transition is consistent with weak first order, contrary to the expectation of a crossover for not too large quark masses and in agreement with results obtained by the Pisa group.


title: Lattice QCD at finite density

abstract: QCD at finite density presents specific challenges to lattice gauge theory. Nonetheless, a region of the QCD phase diagram up to moderately large baryon chemical potentials has been successfully explored on the lattice and new results and idea are continuously emerging. I will outline the lattice formulation of QCD, introduce the calculational schemes currently used to treat a nonzero baryon density, and mention lattice methods alternative to MonteCarlo, including the strong coupling expansion which might give access to the the superconducting phase of QCD. The results for the critical line, and the different phases will be discussed highlighting the strength of the different methods, as well as the possible comparisons with phenomenological models.


title: The QCD equation of state at finite T/\mu on the lattice

abstract: We present N_t=4 lattice results for the equation of state of 2+1 flavour staggered, dynamical QCD at finite temperature and chemical potential. We use the overlap improving multi-parameter reweighting technique to extend the equation of state for non-vanishing chemical potentials. The results are obtained along the line of constant physics. Our physical parameters extend in temperature and baryon chemical potential upto \approx 500-600 MeV.


title: Finite T/\mu lattice QCD and the critical point

abstract: We present results on lattice QCD at finite temperature (T) and chemical potential (\mu). We apply the overlap-improving multi-parameter reweighting technique for n_f=4 and 2+1 staggered QCD. We use semi-realistic masses on L_t=4 lattices. The critical endpoint (E) of 2+1 flavour QCD on the \mu-T plane is located. Our results are based on {\cal{O}}(10^3-10^4) configurations.


title: Real and imaginary chemical potential in 2-color QCD

abstract: In this paper we study the finite temperature SU(2) gauge theory with staggered fermions for non-zero imaginary and real chemical potential. The method of analytical continuation of Monte Carlo results from imaginary to real chemical potential is tested by comparison with simulations performed {\em directly} for real chemical potential. We discuss the applicability of the method in the different regions of the phase diagram in the temperature -- imaginary chemical potential plane.


title: A numerical solution to the local cohomology problem in U(1) chiral gauge theories

abstract: We consider a numerical method to solve the local cohomology problem related to the gauge anomaly cancellation in U(1) chiral gauge theories. In the cohomological analysis of the chiral anomaly, it is required to carry out the differentiation and the integration of the anomaly with respect to the continuous parameter for the interpolation of the admissible gauge fields. In our numerical approach, the differentiation is evaluated explicitly through the rational approximation of the overlap Dirac operator with Zolotarev optimization. The integration is performed with a Gaussian Quadrature formula, which turns out to show rather good convergence. The Poincare lemma is reformulated for the finite lattice and is implemented numerically. We compute the current associated with the cohomologically trivial part of the chiral anomaly in two-dimensions and check its locality properties.


title: Baryonic Flux in Quenched and Two-Flavor Dynamical QCD after Abelian projection

abstract: We study the distribution of color electric flux of the three-quark system in quenched and full QCD (with N_f = 2 flavors of dynamical quarks) at zero and finite temperature. To reduce ultra-violet fluctuations, the calculations are done in the abelian projected theory fixed to the maximally abelian gauge. In the confined phase we find clear evidence for a Y--shape flux tube surrounded and formed by the solenoidal monopole current, in accordance with the dual superconductor picture of confinement. In the deconfined, high temperature phase monopoles cease to condense, and the distribution of the color electric field becomes Coulomb--like.


title: Profiles of the broken string in two-flavor QCD below and above the finite temperature transition

abstract: We study the Abelian flux tube profile of mesonic and the baryonic configurations of static quarks below and above the finite-temperature transition in $N_f=2$ full QCD. To reduce effects of ultra-violet fluctuations we measure Abelian distributions of the action density, the color-electric field and the monopole current after fixing to the maximally Abelian gauge. Changing the distance $R$ between quarks for fixed $T/T_c<1$, one can see a clear signal of string breaking for large $R$. The electric field becomes Coulomb-like and the circulating monopole current disappears. We study also the temperature dependence of the flux profile for fixed $R \approx 0.7$ fm. The disappearance of the squeezed flux is observed clearly above $T_c$. Similar behaviors are observed both in the mesonic and the baryonic systems.


title: Topology of SU(N) gauge theories at T=0 and T=Tc

abstract: We calculate the topological charge density of SU(N) lattice gauge fields for values of N up to N=8. Our T=0 topological susceptibility appears to approach a finite non-zero limit at N=infinity that is consistent with earlier extrapolations from smaller values of N. Near the deconfining temperature Tc we are able to investigate separately the confined and deconfined phases, since the transition is quite strongly first order. We find that the topological susceptibility of the confined phase is always very similar to that at T=0. By contrast, in the deconfined phase at larger N there are no topological fluctuations except for rare, isolated and small instantons. This shows that as N->infinity the large-T suppression of large instantons and the large-N suppression of small instantons overlap, even at T=Tc, so as to suppress all topological fluctuations in the deconfined phase. In the confined phase by contrast, the size distribution is much the same at all T, becoming more peaked as N grows, suggesting that D(rho) is proportional to a delta function at N=infinity, centered on rho close to 1/Tc.


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: Perturbative Determination of Mass Dependent $O(a)$ Improvement Coefficients for the Vector and Axial Vector Currents with a Relativistic Heavy Quark Action

abstract: We carry out a perturbative determination of mass dependent renormalization factors and $O(a)$ improvement coefficients for the vector and axial vector currents with 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}$ with $m_Q$ the heavy quark mass. We discuss what kind of improvement operators are required for the heavy-heavy and the heavy-light cases under the condition that the Euclidean rotational symmetry is not retained anymore because of the $m_Qa$ corrections. Our calculation is performed employing the ordinary perturbation theory with the fictitious gluon mass as an infrared regulator. We show that all the improvement coefficients are determined free from infrared divergences. Results of the renormalization factors and the improvement coefficients are presented as a function of $m_Q a$ for various improved gauge actions as well as the plaquette action.


title: Lattice QCD at non-zero chemical potential and the resonance gas model

abstract: We present results from lattice calculations on the thermodynamics of QCD at non-zero temperature and baryon chemical potential and discuss the role of resonances for the occurrence of the transition to the quark-gluon plasma in hot and dense matter. Properties of a hadronic resonance gas are compared to lattice results on the equation of state at zero as well as non-zero baryon chemical potential. Furthermore, it is shown that the quark mass dependence of the transition temperature can be understood in terms of lines of constant energy density in a resonance gas.


title: QCD String formation and the Casimir Energy

abstract: Three distinct scales are identified in the excitation spectrum of the gluon field around a static quark-antiquark pair as the color source separation R is varied. The spectrum, with string-like excitations on the largest length scales of 2-3 fm, provides clues in its rich fine structure for developing an effective bosonic string description. New results are reported from the three-dimensional Z(2) and SU(2) gauge models, providing further insight into the mechanism of bosonic string formation. The precocious onset of string-like behavior in the Casimir energy of the static quark-antiquark ground state is observed below R=1 fm where most of the string eigenmodes do not exist and the few stable excitations above the ground state are displaced. We find no firm theoretical foundation for the widely held view of discovering string formation from high precision ground state properties below the 1 fm scale.


title: Renormalization Constants of Quark Operators for the Non-Perturbatively Improved Wilson Action

abstract: We present the results of an extensive lattice calculation of the renormalization constants of bilinear and four-quark operators for the non-perturbatively O(a)-improved Wilson action. The results are obtained in the quenched approximation at four values of the lattice coupling by using the non-perturbative RI/MOM renormalization method. Several sources of systematic uncertainties, including discretization errors and final volume effects, are examined. The contribution of the Goldstone pole, which in some cases may affect the extrapolation of the renormalization constants to the chiral limit, is non-perturbatively subtracted. The scale independent renormalization constants of bilinear quark operators have been also computed by using the lattice chiral Ward identities approach and compared with those obtained with the RI-MOM method. For those renormalization constants the non-perturbative estimates of which have been already presented in the literature we find an agreement which is typically at the level of 1%.


title: Nf=1 QCD simulation with improved gauge action at finite temperature

abstract: We study the deconfinement transition of $N_f=1$ QCD by the hybrid Monte Carlo algorithm with Wilson fermions. We calculate the Polyakov loop, its susceptibility and Binder cumulant and use the $\chi^2$ method to locate the phase transition point. Our results are similar to the previous results obtained by the multiboson algorithm.


title: Low-energy couplings of QCD from current correlators near the chiral limit

abstract: We investigate a new numerical procedure to compute fermionic correlation functions at very small quark masses. Large statistical fluctuations, due to the presence of local ``bumps'' in the wave functions associated with the low-lying eigenmodes of the Dirac operator, are reduced by an exact low-mode averaging. To demonstrate the feasibility of the technique, we compute the two-point correlator of the left-handed vector current with Neuberger fermions in the quenched approximation, for lattices with a linear extent of L~1.5 fm, a lattice spacing a~0.09 fm, and quark masses down to the epsilon-regime. By matching the results with the corresponding (quenched) chiral perturbation theory expressions, an estimate of (quenched) low-energy constants can be obtained. We find agreement between the quenched values of F extrapolated from the p-regime and extracted in the epsilon-regime.


title: Chiral perturbation theory for partially quenched twisted mass lattice QCD

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


title: An additional symmetry in the Weinberg - Salam model

abstract: An additional $Z_6$ symmetry hidden in the fermion and Higgs sectors of the Standard Model has been found recently. It has a singular nature and is connected to the centers of the SU(3) and SU(2) subgroups of the gauge group. A lattice regularization of the Standard Model was constructed that possesses this symmetry. In this paper we report our results on the numerical simulation of its Electroweak sector.


title: Thermal effects on quark-gluon mixed condensate g<\bar{q} \sigma G q> from lattice QCD

abstract: We present the first study of the thermal effects on the quark-gluon mixed condensate g<\bar{q} \sigma G q>, which is another chiral order parameter, in SU(3)c lattice QCD with the Kogut-Susskind fermion at the quenched level. Using the lattices at \beta=6.0, 6.1, 6.2 in high statistics, we calculate the mixed condensate as well as the quark condensate in the chiral limit for the temperature 0MeV <= T <= 500MeV. Except for the sharp decrease of both the condensates around T_c \simeq 280MeV, the thermal effects are found to be weak below 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. This result indicates nontrivial similarity in the chiral behaviors of the two different condensates.


title: Critical point of QCD at finite T and \mu, lattice results for physical quark masses

abstract: A critical point (E) is expected in QCD on the temperature (T) versus baryonic chemical potential (\mu) plane. Using a recently proposed lattice method for \mu \neq 0 we study dynamical QCD with n_f=2+1 staggered quarks of physical masses on L_t=4 lattices. Our result for the critical point is T_E=162 \pm 2 MeV and \mu_E= 360 \pm 40 MeV. For the critical temperature at \mu=0 we obtained T_c=164 \pm 2 MeV. This work extends our previous study [Z. Fodor and S.D.Katz, JHEP 0203 (2002) 014] by two means. It decreases the light quark masses (m_{u,d}) by a factor of three down to their physical values. Furthermore, in order to approach the thermodynamical limit we increase our largest volume by a factor of three. As expected, decreasing m_{u,d} decreased \mu_E. Note, that the continuum extrapolation is still missing


title: Supersymmetry breaking in two dimensions: the lattice N=1 Wess-Zumino model

abstract: We study dynamical supersymmetry breaking by non perturbative lattice techniques in a class of two-dimensional N=1 Wess-Zumino models. We work in the Hamiltonian formalism and analyze the phase diagram by analytical strong-coupling expansions and explicit numerical simulations with Green Function Monte Carlo methods.


title: Topological lumps and Dirac zero modes in SU(3) lattice gauge theory on the torus

abstract: We compute eigenmodes of the lattice Dirac operator for quenched SU(3) gauge configurations on the 4-torus with topological charge 1. We find a strong dependence of the zero modes on the boundary conditions which we impose for the Dirac operator. The lumps seen by the eigenmodes often change their position when changing the boundary conditions, while the local chirality of the lumps remains the same. Our results show that the zero mode of a charge 1 configuration can couple to more than one object. We address the question whether these objects could be fractionally charged lumps.


title: An alternate smearing method for Wilson loops in lattice QCD

abstract: A gauge field link smearing method developed for calculations with staggered fermions, namely the use of unitarized fat7 links, is applied to mesonic and baryonic Wilson loop calculations. This method is found to be very effective for reducing statistical fluctuations for large Wilson loops. Examination of chromo-electric field distributions shows that self-interactions of the static sources are reduced when unitarized fat7 smearing is used but long-distance inter-quark effects are unchanged.


title: Recombination of dyons into calorons in SU(2) lattice fields at low temperatures

abstract: By cooling of equilibrium lattice fields at finite temperature in SU(2) gauge theory it has been shown that topological objects (calorons) observed on the lattice in the confined phase possess a dyonic substructure which becomes visible under certain circumstances. Here we show that, with decreasing temperature of the equilibrium ensemble, the distribution in the caloron parameter space is modified such that the calorons appear non-dissociated into constituent dyons. Still the calorons have nontrivial holonomy which is demonstrated by the Polyakov line behaviour for these configurations. At vanishing temperature (on a symmetric lattice) topological lumps obtained by cooling possess rotational symmetry in 4D and a characteristic double peak structure of Polyakov lines (defined with respect to temporal and spatial directions) with non-trivial asymptotics.


title: Lattice formulation of chiral gauge theories

abstract: We present a general formulation of chiral gauge theories, which admits Dirac operators with more general spectra, reveals considerably more possibilities for the structure of the chiral projections, and nevertheless allows appropriate realizations. In our analyses we use two forms of the correlation functions which both also apply in the presence of zero modes and for any value of the index. To account properly for the conditions on the bases the concept of equivalence classes of pairs of them is introduced. The behaviors under gauge transformations and under CP transformations are unambiguously derived.


title: An estimate of the flavour singlet contributions to the hyperfine splitting in charmonium

abstract: We explore the splitting between flavour singlet and non-singlet mesons in charmonium. This has implications for the hyperfine splitting in charmonium.


title: Monopoles in Real Time for Classical U(1) Gauge Field Theory

abstract: U(1) gauge fields are decomposed into a monopole and photon part across the phase transition from the confinement to the Coulomb phase. We analyze the leading Lyapunov exponents of such gauge field configurations on the lattice which are initialized by quantum Monte Carlo simulations. We observe that the monopole field carries the same Lyapunov exponent as the original U(1) field. As a long awaited result, we show that monopoles are created and annihilated in pairs as a function of real time in excess to a fixed average monopole number.


title: Full Spectrum of Lyapunov Exponents in Gauge Field Theory

abstract: We analyze the Lyapunov exponents of U(1) gauge fields across the phase transition from the confinement to the Coulomb phase on the lattice which are initialized by quantum Monte Carlo simulations. We observe all features of a strange attractor with a tendency to regularity towards the continuum limit. Results are also displayed for the full spectrum of Lyapunov exponents of the SU(2) gauge system.


title: Quantum chaos in supersymmetric QCD at finite density

abstract: We investigate the distribution of the spacings of adjacent eigenvalues of the lattice Dirac operator. At zero chemical potential $\mu$, the nearest-neighbor spacing distribution $P(s)$ follows the Wigner surmise of random matrix theory both in the confinement and in the deconfinement phase. This is indicative of quantum chaos. At nonzero chemical potential, the eigenvalues of the Dirac operator become complex and we discuss how $P(s)$ can be defined in the complex plane. Numerical results from an SU(2) simulation with staggered fermions in fundamental and adjoint representations are compared with predictions from non-hermitian random matrix theory, and agreement with the Ginibre ensemble is found for $\mu\approx 0.5$.


title: Hybrid Meson Spectrum from the FLIC action

abstract: The spectral properties of hybrid meson interpolating fields are investigated. The quantum numbers of the meson are carried by smeared-source fermion operators and highly-improved chromo-electric and -magnetic field operators composed with APE-smeared links. The effective masses of standard and hybrid operators indicate that the ground state meson is effectively isolated using both standard and hybrid interpolating fields. Focus is placed on interpolating fields in which the large spinor components of the quark and antiquark fields are merged. In particular, the effective mass of the exotic $1^{-+}$ meson is reported. Further, we report some values for excited mesonic states using a variational process.


title: Properties of P-vortex and monopole clusters in lattice SU(2) gauge theory

abstract: We study the action and geometry of P-vortices, discriminating between the percolating and finite clusters. We also discuss the interrelation of the monopoles and P-vortices. To define P-vortices we use both the direct maximal center projection and indirect maximal center projection. We find, in particular, that the action density of the P-vortices in short clusters is substantially higher than in the percolating cluster. The surface of the percolating cluster appears random at short distances, with action density depending on the shape.


title: A novel probe of the vacuum of the lattice gluodynamics

abstract: We introduce a notion of minimal number of negative links on the lattice for a given original configuration of SU(2) fields. Negative links correspond to a large potential, not necessarily large action. The idea is that the minimal number of negative links is a gauge invariant notion. To check this hypothesis we measure correlator of two negative links, averaged over all the directions, as function of the distance between the links. The inverse correlation length coincides within the error bars with the lightest glueball mass.


title: Dynamic critical behavior of the Swendsen--Wang Algorithm for the three-dimensional Ising model

abstract: We have performed a high-precision Monte Carlo study of the dynamic critical behavior of the Swendsen-Wang algorithm for the three-dimensional Ising model at the critical point. For the dynamic critical exponents associated to the integrated autocorrelation times of the "energy-like" observables, we find z_{int,N} = z_{int,E} = z_{int,E'} = 0.459 +- 0.005 +- 0.025, where the first error bar represents statistical error (68% confidence interval) and the second error bar represents possible systematic error due to corrections to scaling (68% subjective confidence interval). For the "susceptibility-like" observables, we find z_{int,M^2} = z_{int,S_2} = 0.443 +- 0.005 +- 0.030. For the dynamic critical exponent associated to the exponential autocorrelation time, we find z_{exp} \approx 0.481. Our data are consistent with the Coddington-Baillie conjecture z_{SW} = \beta/\nu \approx 0.5183, especially if it is interpreted as referring to z_{exp}.

categories: hep-lat cond-mat.stat-mech math-ph math.MP


title: On the stability of Dirac sheet configurations

abstract: Using cooling for SU(2) lattice configurations, purely Abelian constant magnetic field configurations were left over after the annihilation of constituents that formed metastable Q=0 configurations. These so-called Dirac sheet configurations were found to be stable if emerging from the confined phase, close to the deconfinement phase transition, provided their Polyakov loop was sufficiently non-trivial. Here we show how this is related to the notion of marginal stability of the appropriate constant magnetic field configurations. We find a perfect agreement between the analytic prediction for the dependence of stability on the value of the Polyakov loop (the holonomy) in a finite volume and the numerical results studied on a finite lattice in the context of the Dirac sheet configurations.


title: Duality and scaling in 3-dimensional scalar electrodynamics

abstract: Three-dimensional scalar electrodynamics, with a local U(1) gauge symmetry, is believed to be dual to a scalar theory with a global U(1) symmetry, near the phase transition point. The conjectured duality leads to definite predictions for the scaling exponents of the gauge theory transition in the type II region, and allows thus to be scrutinized empirically. We review these predictions, and carry out numerical lattice Monte Carlo measurements to test them: a number of exponents, characterising the two phases as well as the transition point, are found to agree with expectations, supporting the conjecture. We explain why some others, like the exponent characterising the photon correlation length, appear to disagree with expectations, unless very large system sizes and the extreme vicinity of the transition point are considered. Finally, we remark that in the type I region the duality implies an interesting quantitative relationship between a magnetic flux tube and a 2-dimensional non-topological soliton.


title: The continuum limit of the quark mass step scaling function in quenched lattice QCD

abstract: The renormalisation group running of the quark mass is determined non-perturbatively for a large range of scales, by computing the step scaling function in the Schroedinger Functional formalism of quenched lattice QCD both with and without O(a) improvement. A one-loop perturbative calculation of the discretisation effects has been carried out for both the Wilson and the Clover-improved actions and for a large number of lattice resolutions. The non-perturbative computation yields continuum results which are regularisation independent, thus providing convincing evidence for the uniqueness of the continuum limit. As a byproduct, the ratio of the renormalisation group invariant quark mass to the quark mass, renormalised at a hadronic scale, is obtained with very high accuracy.


title: B to rho ell nu form factors in lattice QCD

abstract: We present results from quenched lattice QCD for the form factors for the decay $B\to\rho l\nu$. The calculations are performed using a nonperturbatively improved action and operators at two values of the lattice spacing. The bottom quark mass is reached by extrapolation from simulations performed with heavy quark masses around the charm mass. Our primary result is for the partially integrated decay rate $\Gamma_{PI}$ over the range $12.7\gev^2 < q^2 < 18.2\gev^2$: \Gamma_PI = (4.9^{+12}_{-10}{+0}_{-14}) X 10^12 s^-1


title: From the Abelian projected flux tube to the dual Abelian Higgs model

abstract: A recent detailed study of the Q-anti-Q flux tube and its DAHM analysis is reported.


title: I=2 $\pi\pi$ Scattering Phase Shift with two Flavors of $O(a)$ Improved Dynamical Quarks

abstract: We present a lattice QCD calculation of phase shift including the chiral and continuum extrapolations in two-flavor QCD. The calculation is carried out for I=2 S-wave $\pi\pi$ scattering. The phase shift is evaluated for two momentum systems, the center of mass and laboratory systems, by using the finite volume method proposed by L\"uscher in the center of mass system and its extension to general systems by Rummukainen and Gottlieb. The measurements are made at three different bare couplings $\beta = 1.80$, 1.95 and 2.10 using a renormalization group improved gauge and a tadpole improved clover fermion action, and employing a set of configurations generated for hadron spectroscopy in our previous work. The illustrative values we obtain for the phase shift in the continuum limit are $\delta$(deg.) $= - 3.50(64)$, $ - 9.5(30)$ and $ - 16.9(64)$ for $\sqrt{s}({\rm GeV})$ $=0.4$, $ 0.6$ and $ 0.8$, which are consistent with experiment.


title: Lattice QCD calculation of the proton decay matrix element in the continuum limit

abstract: We present a quenched lattice QCD calculation of the \alpha and \beta parameters of the proton decay matrix element. The simulation is carried out using the Wilson quark action at three values of the lattice spacing in the range a\approx 0.1-0.064 fm to study the scaling violation effect. We find only mild scaling violation when the lattice scale is determined by the nucleon mass. We obtain in the continuum limit, |\alpha(NDR,2GeV)|=0.0090(09)(^{+5}_{-19})GeV^3 and |\beta(NDR,2GeV)|=0.0096(09)(^{+6}_{-20})GeV^3 with \alpha and \beta in a relatively opposite sign, where the first error is statistical and the second is due to the uncertainty in the determination of the physical scale.


title: Topological susceptibility in full QCD with Ginsparg-Wilson fermions

abstract: We show that, if the formula for the topological charge density operator suggested by fermions obeying the Ginsparg-Wilson relation is employed, it is possible to give a precise and unambiguous definition of the topological susceptibility in full QCD, chi_tL, for finite quark masses on the lattice. The lattice expression of chi_tL looks like the formal continuum one, in the sense that no power divergent subtractions are needed for its proper definition. As a consequence, the small mass behaviour of chi_tL leads directly to a multiplicative renormalizable definition of the chiral condensate that does not require any power divergent subtraction.


title: Lattice representation of vector and chiral gauge theories

abstract: A lattice derivative is defined as a discrete Fourier transform of momentum on a finite lattice. Species doublers are removed with anti-periodic boundary conditions. U(1) chiral transformation is modified to reproduce chiral anomaly. Chiral gauge theories can be constructed on the lattice using a single Weyl fermion as a building block.


title: Hadronic Electromagnetic Properties at Finite Lattice Spacing

abstract: Electromagnetic properties of the octet mesons as well as the octet and decuplet baryons are augmented in quenched and partially quenched chiral perturbation theory to include O(a) corrections due to lattice discretization. We present the results for the SU(3) flavor group in the isospin limit as well as the results for SU(2) flavor with non-degenerate quarks. These corrections will be useful for extrapolation of lattice calculations using Wilson valence and sea quarks, as well as calculations using Wilson sea quarks and Ginsparg-Wilson valence quarks.


title: Light hadrons with improved staggered quarks: approaching the continuum limit

abstract: We have extended our program of QCD simulations with an improved Kogut-Susskind quark action to a smaller lattice spacing, approximately 0.09 fm. Also, the simulations with a approximately 0.12 fm have been extended to smaller quark masses. In this paper we describe the new simulations and computations of the static quark potential and light hadron spectrum. These results give information about the remaining dependences on the lattice spacing. We examine the dependence of computed quantities on the spatial size of the lattice, on the numerical precision in the computations, and on the step size used in the numerical integrations. We examine the effects of autocorrelations in "simulation time" on the potential and spectrum. We see effects of decays, or coupling to two-meson states, in the 0++, 1+, and 0- meson propagators, and we make a preliminary mass computation for a radially excited 0- meson.


title: Non-Commutativity of the Zero Chemical Potential Limit and the Thermodynamic Limit in Finite Density Systems

abstract: Monte Carlo simulations of finite density systems are often plagued by the complex action problem. We point out that there exists certain non-commutativity in the zero chemical potential limit and the thermodynamic limit when one tries to study such systems by reweighting techniques. This is demonstrated by explicit calculations in a Random Matrix Theory, which is thought to be a simple qualitative model for finite density QCD. The factorization method allows us to understand how the non-commutativity, which appears at the intermediate steps, cancels in the end results for physical observables.


title: Unquenched Gluon Propagator in Landau Gauge

abstract: Using lattice quantum chromodynamics (QCD) we perform an unquenched calculation of the gluon propagator in Landau gauge. We use configurations generated with the AsqTad quark action by the MILC collaboration for the dynamical quarks and compare the gluon propagator of quenched QCD (i.e., the pure Yang-Mills gluon propagator) with that of 2+1 flavor QCD. The effects of the dynamical quarks are clearly visible and lead to a significant reduction of the nonperturbative infrared enhancement relative to the quenched case.


title: Perturbation theory vs. simulation for tadpole improvement factors in pure gauge theories

abstract: We calculate the mean link in Landau gauge for Wilson and improved SU(3) anisotropic gauge actions, using two loop perturbation theory and Monte Carlo simulation employing an accelerated Langevin algorithm. Twisted boundary conditions are employed, with a twist in all four lattice directions considerably improving the (Fourier accelerated) convergence to an improved lattice Landau gauge. Two loop perturbation theory is seen to predict the mean link extremely well even into the region of commonly simulated gauge couplings and so can be used remove the need for numerical tuning of self-consistent tadpole improvement factors. A three loop perturbative coefficient is inferred from the simulations and is found to be small. We show that finite size effects are small and argue likewise for (lattice) Gribov copies and double Dirac sheets.


title: Wess-Zumino model with exact supersymmetry on the lattice

abstract: A lattice formulation of the four dimensional Wess-Zumino model that uses Ginsparg-Wilson fermions and keeps exact supersymmetry is presented. The supersymmetry transformation that leaves invariant the action at finite lattice spacing is determined by performing an iterative procedure in the coupling constant. The closure of the algebra, generated by this transformation is also showed.


title: Computational methods for the fermion determinant and the link between overlap and domain wall fermions

abstract: This paper reviews the most popular methods which are used in lattice QCD to compute the determinant of the lattice Dirac operator: Gaussian integral representation and noisy methods. Both of them lead naturally to matrix function problems. We review the most recent development in Krylov subspace evaluation of matrix functions. The second part of the paper reviews the formal relationship and algebraic structure of domain wall and overlap fermions. We review the multigrid algorithm to invert the overlap operator. It is described here as a preconditioned Jacobi iteration where the preconditioner is the Schur complement of a certain block of the truncated overlap matrix.


title: Dynamical fat link 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. Frequently, fat link actions make 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: Approximation Theory for Matrices

abstract: We review the theory of optimal polynomial and rational Chebyshev approximations, and Zolotarev's formula for the sign function over the range (\epsilon \leq |z| \leq1). We explain how rational approximations can be applied to large sparse matrices efficiently by making use of partial fraction expansions and multi-shift Krylov space solvers.


title: Fast Evaluation of Zolotarev Coefficients

abstract: We review the theory of elliptic functions leading to Zolotarev's formula for the sign function over the range (\epsilon \leq|x| \leq1). We show how Gauss' arithmetico-geometric mean allows us to evaluate elliptic functions cheaply, and thus to compute Zolotarev coefficients ``on the fly'' as a function of (\epsilon). This in turn allows us to calculate the matrix functions (\sgn H), (\sqrt H), and (1/\sqrt H) both quickly and accurately for any Hermitian matrix (H) whose spectrum lies in the specified range.


title: Critical slowing down of topological modes

abstract: We investigate the critical slowing down of the topological modes using local updating algorithms in lattice 2-d CP^(N-1) models. We show that the topological modes experience a critical slowing down that is much more severe than the one of the quasi-Gaussian modes relevant to the magnetic susceptibility, which is characterized by $\tau_{\rm mag} \sim \xi^z$ with $z\approx 2$. We argue that this may be a general feature of Monte Carlo simulations of lattice theories with non-trivial topological properties, such as QCD, as also suggested by recent Monte Carlo simulations of 4-d SU(N) lattice gauge theories.


title: Infrared and ultraviolet properties of the Landau gauge quark propagator

abstract: We present a current summary of a program to study the quark propagator using lattice QCD. We use the Overlap and ``Asqtad'' quark actions on a number of lattice ensembles to assess systematic errors. We comment on the place of this work amongst studies of QCD Green's functions in other formulations. A preliminary calculation of the running quark mass is presented.


title: The overlap Dirac operator as a continued fraction

abstract: We use a continued fraction expansion of the sign-function in order to obtain a five dimensional formulation of the overlap lattice Dirac operator. Within this formulation the inverse of the overlap operator can be calculated by a single Krylov space method and nested conjugate gradient procedures are avoided. We point out that the five dimensional linear system can be made well conditioned using equivalence transformations on the continued fractions.


title: Short distance behaviour of the effective string

abstract: We study the Polyakov loop correlator in the (2+1) dimensional Z_2 gauge model. An algorithm that we have presented recently, allows us to reach high precision results for a large range of distances and temperatures, giving us the opportunity to test predictions of the effective Nambu-Goto string model. Here we focus on the regime of low temperatures and small distances. In contrast to the high temperature, large distance regime, we find that our numerical results are not well described by the two loop-prediction of the Nambu-Goto model. In addition we compare our data with those for the SU(2) and SU(3) gauge models in (2+1) dimensions obtained by other authors. We generalize the result of L\"uscher and Weisz for a boundary term in the interquark potential to the finite temperature case.


title: Electroweak Matrix Elements in the Two-Nucleon Sector from Lattice QCD

abstract: We demonstrate how to make rigorous predictions for electroweak matrix elements in nuclear systems directly from QCD. More precisely, we show how to determine the short-distance contributions to low-momentum transfer electroweak matrix elements in the two-nucleon sector from lattice QCD. In potential model descriptions of multi-nucleon systems, this is equivalent to uniquely determining the meson-exchange currents, while in the context of nuclear effective field theory, this translates into determining the coefficients of local, gauge-invariant, multi-nucleon-electroweak current operators. The energies of the lowest-lying states of two nucleons on a finite volume lattice with periodic boundary conditions in the presence of a background magnetic field are sufficient to determine the local four-nucleon operators that contribute to the deuteron magnetic moment and to the threshold cross-section of n + p -> d + gamma. Similarly, the energy-levels of two nucleons immersed in a background isovector axial weak field can be used to determine the coefficient of the leading local four-nucleon operator contributing to the neutral- and charged-current break-up of the deuteron. This is required for the extraction of solar neutrino fluxes at SNO and future neutrino experiments.


title: Effective Field Theories and Quantum Chromodynamics on the Lattice

abstract: We give a selection of results on spectrum and decay constants of light and heavy-light hadrons. Effective field theories relevant for their lattice calculation, namely non-relativistic QCD (NRQCD) for heavy quarks on the lattice and Chiral Perturbation Theory for light quarks, are briefly discussed.


title: Improved Wilson QCD simulations with light quark masses

abstract: We present results from simulations using 2 flavours of O(a)-improved Wilson quarks whose masses are about 1/3 of the physical strange quark mass. We present new data on the mass of the singlet pseudoscalar meson and evidence of the onset of chiral logarithms in the pion decay constant. The previously observed suppression of the topological susceptibility at lighter quark masses is confirmed. We report on the performance of the hybrid Monte Carlo algorithm at light quark masses.


title: Density matrix renormalization group in a two-dimensional $\lambda\phi^4$ Hamiltonian lattice model

abstract: Density matrix renormalization group (DMRG) is applied to a (1+1)-dimensional $\lambda\phi^4$ model. Spontaneous breakdown of discrete $Z_2$ symmetry is studied numerically using vacuum wavefunctions. We obtain the critical coupling $(\lambda/\mu^2)_{\rm c}=59.89\pm 0.01$ and the critical exponent $\beta=0.1264\pm 0.0073$, which are consistent with the Monte Carlo and the exact results, respectively. The results are based on extrapolation to the continuum limit with lattice sizes $L=250,500$, and 1000. We show that the lattice size L=500 is sufficiently close to the the limit $L\to\infty$.


title: Finite size effects of a pion matrix element

abstract: We investigate finite size effects of the pion matrix element of the non-singlet, twist-2 operator corresponding to the average momentum of non-singlet quark densities. Using the quenched approximation, they come out to be surprisingly large when compared to the finite size effects of the pion mass. As a consequence, simulations of corresponding nucleon matrix elements could be affected by finite size effects even stronger which could lead to serious systematic uncertainties in their evaluation.


title: Systematic errors due to linear congruential random-number generators with the Swendsen-Wang algorithm: A warning

abstract: We show that linear congruential pseudo-random-number generators can cause systematic errors in Monte Carlo simulations using the Swendsen-Wang algorithm, if the lattice size is a multiple of a very large power of 2 and one random number is used per bond. These systematic errors arise from correlations within a single bond-update half-sweep. The errors can be eliminated (or at least radically reduced) by updating the bonds in a random order or in an aperiodic manner. It also helps to use a generator of large modulus (e.g. 60 or more bits).


title: Signals of confinement in Green functions of SU(2) Yang-Mills theory

abstract: The vortex picture of confinement is employed to explore the signals of confinement in Yang-Mills Green functions. By using SU(2) lattice gauge theory, it has been well established that the removal of the center vortices from the lattice configurations results in the loss of confinement. The running coupling constant, the gluon and the ghost form factors are studied in Landau gauge for both cases, the full and the vortex removed theory. In the latter case, a strong suppression of the running coupling constant and the gluon form factor at low momenta is observed. At the same time, the singularity of the ghost form factor at vanishing momentum disappears. This observation establishes an intimate correlation between the ghost singularity and confinement. The result also shows that a removal of the vortices generates a theory for which Zwanziger's horizon condition for confinement is no longer satisfied.


title: Heavy meson chiral perturbation theory in finite volume

abstract: We study finite volume effects in heavy quark systems in the framework of heavy meson chiral perturbation theory for full, quenched, and partially quenched QCD. A novel feature of this investigation is the role played by the scales Delta and delta_s, where Delta is the mass difference between the heavy-light vector and pseudoscalar mesons of the same quark content, and delta_s is the mass difference due to light flavour SU(3) breaking. The primary conclusion of this work is that finite volume effects arising from the propagation of Goldstone particles in the effective theory can be altered by the presence of these scales. Since Delta varies significantly with the heavy quark mass, these volume effects can be amplified in both heavy and light quark mass extrapolations (interpolations). As an explicit example, we present results for B parameters of neutral B meson mixing matrix elements and heavy-light decay constants to one-loop order in finite volume heavy meson chiral perturbation theory for full, quenched, and N_f=2+1 partially quenched QCD. Our calculation shows that for high-precision determinations of the phenomenologically interesting SU(3) breaking ratios, finite volume effects are significant in quenched and not negligible in partially quenched QCD, although they are generally small in full QCD.


title: On projection (in)dependence of monopole condensate in lattice SU(2) gauge theory

abstract: We study the temperature dependence of the monopole condensate in different Abelian projections of the SU(2) lattice gauge theory. Using the Frohlich-Marchetti monopole creation operator we show numerically that the monopole condensate depends on the choice of the Abelian projection. Contrary to the claims in the current literature we observe that in the Abelian Polyakov gauge and in the field strength gauge the monopole condensate does not vanish at the critical temperature and thus is not an order parameter.


title: Quark mass dependence of pseudoscalar masses and decay constants on a lattice

abstract: Our previous calculations of the sea- and valence-quark mass dependence of the pseudoscalar meson masses and decay constants is repeated on a 16^3x32 lattice which allows for a better determination of the quantities in question. The conclusions are similar as before on the 16^4 lattice [1]. The two light dynamical quark flavours we simulate have masses in the range 1/4 m_s < m_{u,d} < 2/3 m_s. The sea quark mass dependence of f_pi and m_pi^2/m_q is well described by the next-to-leading order (NLO) Chiral Perturbation Theory (ChPT) formulas and clearly shows the presence of chiral logarithms. The valence quark mass dependence requires the presence of NNLO contributions in Partially Quenched ChPT (PQChPT) -- in addition to the NLO terms. The O(a) lattice artifacts in these quantities turn out to be small.


title: Nucleon Masses and Magnetic Moments in a Finite Volume

abstract: We compute finite-size corrections to nucleon masses and magnetic moments in a periodic, spatial box of size L, both in QCD and in partially-quenched QCD.


title: News from Lattice QCD on Heavy Quark Potentials and Spectral Functions of Heavy Quark States

abstract: We discuss recent lattice results on in-medium properties of hadrons and focus on thermal properties of heavy quark bound states. We will clarify the relation between heavy quark free energies and potentials used to analyze the melting of heavy quark bound states. Furthermore, we present calculations of meson spectral functions which indicate that the charmonium ground states, J/psi and eta_c, persist in the quark gluon plasma as well defined resonances with no significant change of their zero temperature masses at least up to T ~ 1.5 T_c. We also briefly comment on the current status of lattice calculations at non-vanishing baryon number density.


title: Charmonium systems after the deconfinement transition

abstract: The behavior of charmonia after the deconfinement transition is investigated on quenched lattices. Analysis of temporal correlators on fine lattices at temperatures upto 3 T_c show that the J/psi and eta_c survive the deconfinement transition with little significant changes, and survive as bound states at least upto 2.25 T_c. The spatially excited chi_c states suffer serious system modifications, maybe dissolution, already a little above T_c.


title: A novel improved action for SU(3) lattice gauge theory

abstract: SU(3) lattice gauge theory is studied by means of an improved action where a $2 \times 2$ Wilson loop is supplemented to the standard plaquette term. By contrast to earlier studies using a tree level improvement, the prefactor of the $2 \times 2$ Wilson term is determined by minimizing the breaking of rotational symmetry detected from the static quark-antiquark potential. On coarse lattices, the novel action is superior to the Iwasaki action and comparable with DBW2 action. The scaling behavior of the novel action is studied by using the static quark potential and the ratio of the deconfinement temperature and the string tension.


title: Hybrid Monte Carlo with Fat Link Fermion Actions

abstract: The use of APE smearing or other blocking techniques in lattice 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. The necessary equations of motion for FLIC fermions are derived, and some initial simulation results are presented. The technique is more general however, and is straightforwardly applicable to other smearing techniques or fat link actions.


title: Dimensional reduction in QCD: Lessons from lower dimensions

abstract: In this contribution we present the results of a series of investigations of dimensional reduction, applied to SU(3) gauge theory in 2 + 1 dimensions. We review earlier results, present a new reduced model with Z(3) symmetry, and discuss the results of numerical simulations of this model.


title: What the Gribov copy tells on the confinment and the theory of dynamical chiral symmetry breaking

abstract: We performed lattice Landau gauge QCD simulation on \beta=6.0, 16^4, 24^4, 32^4 and \beta=6.4, 32^4, 48^4 and 56^4 by adopting the gauge fixing that minimizes the norm of the gauge field, and measured the running coupling by using the gluon propagator and the ghost propagator. In view of ambiguity in the vertex renormalization factor \tilde Z_1 in the lattice, we adjust the normalization of the running coupling by the perturbative QCD results near the highest momentum point. It has a maximum \alpha_s(q)~ 2.1(3) at around q=0.5 GeV and decreases as q approaches 0, and the Kugo-Ojima parameter reached -0.83(2). The infrared exponent of the ghost propagator at 0.4GeV region is \alpha_G=0.20 but there is an exceptional Gribov copy with \alpha_G=0.27. The features of the exceptional Gribov copy are investigated by measuring four one-dimensional Fourier transform(1-d FT) of the gluon propagator transverse to each lattice axis. We observe, in general, correlation between absolute value of the Kugo-Ojima parameter and the degree of reflection positivity violation in the 1-d FT of the gluon propagator. The 1-d FT of the exceptional Gribov copy has an axis whose gluon propagator manifestly violates reflection positivity, and the average of the Cartan subalgebra components of the Kugo-Ojima parameter along this axis is consistent to -1. The running coupling of the enemble average shows a suppression at 0 momentum, but when the ghost propagator of the exceptional Gribov copy is adopted, the suppression disappears and the data implies presence of the infrared fixed point \alpha_s(0)~ 2.5(5) and \kappa=0.5 suggested by the Dyson-Schwinger approach in the multiplicative renormalizable scheme. Comparison with the SU(2) QCD and N_f=2 unquenched SU(3) QCD are also made.