title: The Schr\"odinger functional in lattice QCD with exact chiral symmetry

abstract: Similarly to the interaction lagrangian, the possible boundary conditions in quantum field theories on space-time manifolds with boundaries are strongly constrained by the symmetry and scaling properties of the theory. Based on this general insight, a lattice formulation of the QCD Schr\"odinger functional is proposed for the case where the lattice Dirac operator in the bulk of the lattice coincides with the Neuberger--Dirac operator. The construction satisfies all basic requirements (locality, symmetries, hermiticity) and is suitable for numerical simulations.


title: Sizes of Lightest Glueballs in SU(3) Lattice Gauge Theory

abstract: Standard Monte Carlo simulations have been performed on improved lattices to measure the wave functions and sizes of the scalar and tensor glueballs at four lattice spacings in the range $a= 0.05 - 0.145$ fm. Systematic errors from discretization and finite volume are studied. Our results in the continuum limit show that the size of the tensor state is approximately two times as large as that of the scalar glueball.


title: Two-flavor QCD Thermodynamics using Anisotropic Lattices

abstract: Numerical simulations of full QCD on anisotropic lattices provide a convenient way to study QCD thermodynamics with fixed physics scales and reduced lattice spacing errors. We report results from calculations with two flavors of dynamical staggered fermions, where all bare parameters and the renormalized anisotropy are kept constant and the temperature is changed in small steps by varying only the number of time slices. Including results from zero-temperature scale setting simulations, which determine the Karsch coefficients, allows for the calculation of the equation of state at finite temperatures.


title: Decay constants of P-wave heavy-light mesons from unquenched lattice QCD

abstract: We review some decays that require knowledge of the decay constants of $0^{+}$ heavy-light mesons. We compute the decay constants of P-wave heavy-light mesons from unquenched lattice QCD, with two degenerate flavours of sea quarks, at a single lattice spacing. The lightest sea quark mass used in the calculation is a third of the strange quark mass. For the charm-strange meson we obtain the decay constant: $f_{D_{s 0^+}} = 340(110)$ MeV using our normalisation conventions. We obtain the $f_{P_s}^{static}$ (static-strange P-wave) decay constant as $302(39)$ MeV.


title: Schroedinger functional formalism with domain-wall fermion

abstract: Finite volume renormalization scheme is one of the most fascinating scheme for non-perturbative renormalization on lattice. By using the step scaling function one can follow running of renormalized quantities with reasonable cost. It has been established the Schroedinger functional is very convenient to define a field theory in a finite volume for the renormalization scheme. The Schroedinger functional, which is characterized by a Dirichlet boundary condition in temporal direction, is well defined and works well for the Yang-Mills theory and QCD with the Wilson fermion. However one easily runs into difficulties if one sets the same sort of the Dirichlet boundary condition for the overlap Dirac operator or the domain-wall fermion. In this paper we propose an orbifolding projection procedure to impose the Schroedinger functional Dirichlet boundary condition on the domain-wall fermion.


title: Overlap lattice fermion in a gravitational field

abstract: We construct a lattice Dirac operator of overlap type that describes the propagation of a Dirac fermion in an external gravitational field. The local Lorentz symmetry is manifestly realized as a lattice gauge symmetry, while it is believed that the general coordinate invariance is restored only in the continuum limit. Our doubler-free Dirac operator satisfies the conventional Ginsparg-Wilson relation and possesses gamma_5 hermiticity with respect to the inner product, which is suggested by the general coordinate invariance. The lattice index theorem in the presence of a gravitational field holds, and the classical continuum limit of the index density reproduces the Dirac genus. Reduction to a single Majorana fermion is possible for 8k+2 and 8k+4 dimensions, but not for 8k dimensions, which is consistent with the existence of the global gravitational/gauge anomalies in 8k dimensions. Other Lorentz representations, such as the spinor-vector and the bi-spinor representations, can also be treated. Matter fields with a definite chirality (with respect to the lattice-modified chiral matrix) are briefly considered.


title: Deconfinement in dense 2-color QCD

abstract: We study SU(2) lattice gauge theory with two flavors of Wilson fermion at non-zero chemical potential mu and low temperature on a 8^3x16 system. We identify three regimes along the mu-axis. For mu<~m_pi/2 the system remains in the vacuum phase, and all physical observables considered remain essentially unchanged. The intermediate regime is characterised by a non-zero diquark condensate and an associated increase in the baryon density, consistent with what is expected for Bose-Einstein condensation of tightly bound diquarks. We also observe screening of the static quark potential here. In the high-density deconfined regime we find a non-zero Polyakov loop and a strong modification of the gluon propagator, including significant screening in the magnetic sector in the static limit, which must have a non-perturbative origin. The behaviour of thermodynamic observables and the superfluid order parameter are consistent with a Fermi surface disrupted by a BCS diquark condensate. The energy per baryon as a function of mu exhibits a minimum in the deconfined regime, implying that macroscopic objects such as stars formed in this theory are largely composed of quark matter.


title: Lattice fermions with complex mass

abstract: We present evidence in the Schwinger model that rooted staggered fermions may correctly describe the m<0 sector of a theory with an odd number of flavors. We point out that in QCD-type theories with a complex-valued quark mass every non-chiral action essentially "borrows" knowledge about the theta-transformation properties from the overlap action.


title: Non-perturbative mass spectrum of an extra-dimensional orbifold

abstract: We analyse non-perturbatively a five-dimensional SU(2) gauge theory compactified on the S^1/Z_2 orbifold. In particular, we present simulation results for the mass spectrum of the theory, which contains a Higgs and a photon. The Higgs mass is found to be free of divergences without fine-tuning. The photon mass is non-zero, thus providing us with the first lattice evidence for a Higgs mechanism derived from an extra dimension. Data from the static potential are consistent with dimensional reduction at low energies.


title: Analysis of systematic errors in the calculation of renormalization constants of the topological susceptibility on the lattice

abstract: A Ginsparg-Wilson based calibration of the topological charge is used to calculate the renormalization constants which appear in the field-theoretical determination of the topological susceptibility on the lattice. A systematic comparison is made with calculations based on cooling. The two methods agree within present statistical errors (3%-4%). We also discuss the independence of the multiplicative renormalization constant Z from the background topological charge used to determine it.


title: Pseudovector meson with strangeness and closed-charm

abstract: We investigate the mass spectrum of 1^{+} exotic mesons with quark content (\c\s\cbar\qbar) or (\c\q\cbar\sbar), using molecular and diquark-antidiquark operators, in quenched lattice QCD with exact chiral symmetry. For the molecular operator {(\qbar\gamma_i\c)(\cbar\gamma_5\s)-(\cbar\gamma_i\s)(\qbar\gamma_5\c)} and the diquark-antidiquark operator {(\q^T C \gamma_i \c)(\sbar C \gamma_5 \cbar^T)- (\q^T C \gamma_5 \c)(\sbar C \gamma_i^T \cbar^T)}, both detect a 1^{+} resonance with mass around 4010 +- 50 MeV in the limit m_q -> m_{u,d}.


title: Properties of light scalar mesons from lattice QCD

abstract: Lattice QCD with $N_f=2$ flavours of sea quark is used to explore the spectrum and decay of scalar mesons. We are able to determine the $b_1$ - $a_0$ mass difference and this leads to the conclusion that the lightest non-singlet scalar meson ($a_0$) has a mass of 1.01(4) GeV. We determine from the lattice the coupling strength to KK and $\pi \eta$. We compute the leptonic decay constant of the lightest non-singlet scalar meson. We discuss the impact of these lattice results on the interpretation of the $a_0(980)$ state. We also discuss $K^*_0$ states.


title: Validity of the Rooted Staggered Determinant in the continuum limit

abstract: We investigate the continuum limit of the rooted staggered determinant in the 2-dimensional Schwinger model. We match both the unrooted and rooted staggered determinant with an overlap fermion determinant of two (one) flavors and a local pure gauge effective action by fitting the coefficients of the effective action and the mass of the overlap operator. The residue of this fit measures the difference of the staggered and overlap fermion actions. We show that this residue scales at least as O(a^2), implying that any difference, be it local or non-local, between the staggered and overlap actions becomes irrelevant in the continuum limit. This observation justifies the rooting procedure.


title: Gauge invariant discretization of Poincare quantum gravity

abstract: In this paper we suggest gauge invariant discretization of Poincare quantum gravity. We generalize Regge calculus to the case of Riemann-Cartan space. The basic element of the constructed discretization is piecewize linear Riemann-Cartan space with flat pieces of hypercubic form. We consider the model with squared curvature action and calculate the correspondent lattice action. We construct local measure over the dynamical variables of the lattice model.


title: Lattice formulation of (2,2) supersymmetric gauge theories with matter fields

abstract: We construct lattice actions for a variety of (2,2) supersymmetric gauge theories in two dimensions with matter fields interacting via a superpotential.


title: The Gell-Mann -- Okubo mass relation among baryons from fully-dynamical mixed-action lattice QCD

abstract: We explore the Gell-Mann--Okubo mass relation among the octet baryons using fully-dynamical, mixed-action (domain-wall on rooted-staggered) lattice QCD calculations at a lattice spacing of b ~ 0.125 fm and pion masses of m_pi ~ 290 MeV, 350 MeV, 490 MeV and 590 MeV. Deviations from the Gell-Mann--Okubo mass relation are found to be small at each quark mass.


title: A strategy for implementing non-perturbative renormalisation of heavy-light four-quark operators in the static approximation

abstract: We discuss the renormalisation properties of the complete set of $\Delta B = 2$ four-quark operators with the heavy quark treated in the static approximation. We elucidate the role of heavy quark symmetry and other symmetry transformations in constraining their mixing under renormalisation. By employing the Schroedinger functional, a set of non-perturbative renormalisation conditions can be defined in terms of suitable correlation functions. As a first step in a fully non-perturbative determination of the scale-dependent renormalisation factors, we evaluate these conditions in lattice perturbation theory at one loop. Thereby we verify the expected mixing patterns and determine the anomalous dimensions of the operators at NLO in the Schroedinger functional scheme. Finally, by employing twisted-mass QCD it is shown how finite subtractions arising from explicit chiral symmetry breaking can be avoided completely.


title: Simulating an arbitrary number of flavors of dynamical overlap fermions

abstract: We present a set of related Hybrid Monte Carlo methods to simulate an arbitrary number of dynamical overlap fermions. Each fermion is represented by a chiral pseudo-fermion field. The new algorithm reduces critical slowing down in the chiral limit and for sectors of nontrivial topology.


title: Compact lattice formulation of Cho-Faddeev-Niemi decomposition: string tension from magnetic monopoles

abstract: In this paper we begin on a new lattice formulation of the non-linear change of variables called the Cho--Faddeev--Niemi decomposition in SU(2) Yang-Mills theory. This is a compact lattice formulation improving the non-compact lattice formulation proposed in our previous paper. Based on this formulation, we propose a new gauge-invariant definition of the magnetic monopole current which guarantees the magnetic charge quantization and reproduces the conventional magnetic-current density obtained in the Abelian projection based on the DeGrand--Toussaint method. Finally, we demonstrate the magnetic monopole dominance in the string tension in SU(2) Yang-Mills theory on a lattice. Our formulation enables one to reproduce in the gauge-invariant way remarkable results obtained so far only in the Maximally Abelian gauge.


title: Observations on staggered fermions at non-zero lattice spacing

abstract: We show that the use of the fourth-root trick in lattice QCD with staggered fermions corresponds to a non-local theory at non-zero lattice spacing, but argue that the non-local behavior is likely to go away in the continuum limit. We give examples of this non-local behavior in the free theory, and for the case of a fixed topologically non-trivial background gauge field. In both special cases, the non-local behavior indeed disappears in the continuum limit. Our results invalidate a recent claim that at non-zero lattice spacing an additive mass renormalization is needed because of taste-symmetry breaking.


title: Automatic O(a) improvement for twisted-mass QCD in the presence of spontaneous symmetry breaking

abstract: In this paper we present a proof for automatic O(a) improvement in twisted mass lattice QCD at maximal twist, which uses only the symmetries of the leading part in the Symanzik effective action. In the process of the proof we clarify that the twist angle is dynamically determined by vacuum expectation values in the Symanzik theory. For maximal twist according to this definition, we show that scaling violations of all quantities which have non-zero values in the continuum limit are even in a. In addition, using Wilson Chiral Perturbation Theory (WChPT), we investigate this definition for maximal twist and compare it to other definitions which were already employed in actual simulations.


title: Excited hadrons on the lattice: Baryons

abstract: We present results for masses of excited baryons from a quenched calculation with Chirally Improved quarks at pion masses down to 350 MeV. Our analysis of the correlators is based on the variational method. In order to provide a large basis set for spanning the physical states, we use interpolators with different Dirac structures and Jacobi smeared quark sources of different width. Our spectroscopy results for a wide range of ground state and excited baryons are discussed.


title: Overlap Dirac operator at nonzero chemical potential and random matrix theory

abstract: We show how to introduce a quark chemical potential in the overlap Dirac operator. The resulting operator satisfies a Ginsparg-Wilson relation and has exact zero modes. It is no longer gamma_5-hermitian, but its nonreal eigenvalues still occur in pairs. We compute the spectral density of the operator on the lattice and show that, for small eigenvalues, the data agree with analytical predictions of nonhermitian chiral random matrix theory for both trivial and nontrivial topology.


title: Dynamical QCD simulations on anisotropic lattices

abstract: The simulation of QCD on dynamical (Nf=2) anisotropic lattices is described. A method for nonperturbative renormalisation of the parameters in the anisotropic gauge and quark actions is presented. The precision with which this tuning can be carried out is determined in dynamical simulations on 8^3x48 and 8^3x80 lattices.


title: Precision Electromagnetic Structure of Octet Baryons in the Chiral Regime

abstract: The electromagnetic properties of the baryon octet are calculated in quenched QCD on a 20^3 x 40 lattice with a lattice spacing of 0.128 fm using the fat-link irrelevant clover (FLIC) fermion action. FLIC fermions enable simulations to be performed efficiently at quark masses as low as 300 MeV. By combining FLIC fermions with an improved-conserved vector current, we ensure that discretisation errors occur only at O(a^2) while maintaining current conservation. Magnetic moments and electric and magnetic radii are extracted from the electric and magnetic form factors for each individual quark sector. From these, the corresponding baryon properties are constructed. Our results are compared extensively with the predictions of quenched chiral perturbation theory. We detect substantial curvature and environment sensitivity of the quark contributions to electric charge radii and magnetic moments in the low quark mass region. Furthermore, our quenched QCD simulation results are in accord with the leading non-analytic behaviour of quenched chiral perturbation theory, suggesting that the sum of higher-order terms makes only a small contribution to chiral curvature.


title: Spin-3/2 baryons from an anisotropic lattice QCD action

abstract: The mass spectrum of baryons in the spin-3/2 sector is computed in quenched lattice QCD using a tadpole-improved anisotropic action. Both isospin 1/2 and 3/2 (the traditional decuplet) are considered, as well as members that contain strange quarks. States with positive and negative parities are isolated by parity projection, while states with spin-3/2 and spin-1/2 are separated by spin projection. The extent to which spin projection is needed is examined. The issue of optimal interpolating field is also investigated. The results are discussed in relation to previous calculations and experiment.


title: Glauber dynamics of phase transitions: SU(3) lattice gauge theory

abstract: Motivated by questions about the QCD deconfining phase transition, we studied in two previous papers Model A (Glauber) dynamics of 2D and 3D Potts models, focusing on structure factor evolution under heating (heating in the gauge theory notation, i.e., cooling of the spin systems). In the present paper we set for 3D Potts models (Ising and 3-state) the scale of the dynamical effects by comparing to equilibrium results at first and second order phase transition temperatures, obtained by re-weighting from a multicanonical ensemble. Our finding is that the dynamics entirely overwhelms the critical and non-critical equilibrium effects. In the second half of the paper we extend our results by investigating the Glauber dynamics of pure SU(3) lattice gauge on $N_{\tau} N_{\sigma}^3$ lattices directly under heating quenches from the confined into the deconfined regime. The exponential growth factors of the initial response are calculated, which give Debye screening mass estimates. The quench leads to competing vacuum domains of distinct $Z_3$ triality, which delay equilibration of pure gauge theory forever, while their role in full QCD remains a subtle question. As in spin systems we find for pure SU(3) gauge theory a dynamical growth of structure factors, reaching maxima which scale approximately with the volume of the system, before settling down to equilibrium. Their influence on various observables is studied and different lattice sizes are simulated to illustrate an approach to a finite volume continuum limit. Strong correlations are found during the dynamical process, but not in the deconfined phase at equilibrium.


title: Renormalisation of Composite Operators in Lattice Perturbation Theory with Clover Fermions: Non-forward Matrix Elements

abstract: We consider the renormalisation of lattice QCD operators with one and two covariant derivatives related to the first and second moments of generalised parton distributions and meson distribution amplitudes. Employing the clover fermion action we calculate their non-forward quark matrix elements in one-loop lattice perturbation theory. For some representations of the hypercubic group commonly used in simulations we determine the sets of all possible mixing operators and compute the matrices of renormalisation factors in one-loop approximation. We describe how tadpole improvement is applied to the results.


title: Chiral Lattice Gauge Theories from Warped Domain Walls and Ginsparg-Wilson Fermions

abstract: We propose a construction of a 2-dimensional lattice chiral gauge theory. The construction may be viewed as a particular limit of an infinite warped 3-dimensional theory. We also present a "single-site'' construction using Ginsparg-Wilson fermions which may avoid, in both 2 and 4 dimensions, the problems of waveguide-Yukawa models.


title: A deconstruction lattice description of the D1/D5 brane world-volume gauge theory

abstract: In this article, I generalize the deconstruction lattice formulation of Endres and Kaplan [hep-lat/0604012] to two-dimensional super-QCD with eight supercharges, denoted (4,4), and bifundamental matter. I specialize to a particularly interesting (4,4) gauge theory, with gauge group U(N_c) \times U(N_f), and U(N_f) weakly gauged. It describes the infrared limit of the D1/D5 brane system, which has been studied extensively as an example of the AdS_3/CFT_2 correspondence. The construction here preserves two supercharges exactly, and has a lattice structure quite similar to that which has previously appeared in the deconstruction approach; i.e., site, link and diagonal fields with both Bose and Fermi statistics. I remark on possible applications of the lattice theory that would test the AdS_3/CFT_2 correspondence, particularly one that would exploit the recent worldsheet instanton analysis of Chen and Tong [hep-th/0604090].


title: Study of Possible Proton-Antiproton Hexaquark State in Lattice QCD

abstract: We have used standard techniques of lattice quantum chromodynamics to look for evidence of the spin-zero six quark flavour singlet state ($J^{PC}=0^{-+}$) observed by BES Collaboration, and to determine the splitting between the mass of the possible proton-antiproton and the mass of two protons, its threshold. Using the various interpolating fields we find indications that for sufficiently light quarks proton-antiproton is slightly above the $2m_{p}$ threshold. The lattice resonance signal of binding observed near the physical and continuum regimes do not support the existence of proton-antiproton state as a spin-zero hexaquark state.


title: Exploring twisted mass Lattice QCD with the Clover term

abstract: It has been shown that in the twisted mass formulation of Lattice QCD at maximal twist large cutoff effects are generated when the quark mass becomes of O(a Lambda_QCD^2). In general, these effects can be suppressed in two ways: either by choosing the critical quark mass in an "optimal way", or by adding the Clover term to the twisted action. We investigate the second option by performing a quenched lattice QCD simulation with twisted Clover fermions and pion masses as low as 280 MeV. We show that the Clover term is indeed efficient in reducing the large cutoff effects. In particular, the so-called bending phenomenon observed in the determination of the pion decay constant is cured in this way. In addition, by using maximally twisted Clover fermions, we provide a non-perturbative determination of the vector current renormalization constant Z_V as well as of the non-perturbatively renormalized light quark masses. Finally, we calculate the connected contribution to the charged-neutral pseudoscalar meson mass splitting, finding that the introduction of the Clover term in the twisted action is also beneficial, in the quenched approximation, in reducing cutoff effects related to the isospin symmetry breaking at finite lattice spacing.


title: Leptonic widths of heavy quarkonia: S-Wave QCD/NRQCD matching coefficients for the electromagnetic vector annihilation current at O(alpha_s v^2)

abstract: We construct the S-wave part of the electromagnetic vector annihilation current to O(alpha_s v^2)$ on the lattice for heavy quarks whose dynamics are described by the NRQCD action, where v is the non-relativistic quark velocity. The lattice vector current for QQbar annihilation is expressed as a linear combination of lattice operators with quantum numbers L=0, J^P=1^-, and the coefficients are determined by matching this lattice current to the corresponding continuum current in QCD to O(v^2) to one-loop. The annihilation channel gives a complex amplitude and a proper choice for the contours of integration is needed; a simple Wick rotation is not possible. In this way, and with a careful choice of subtraction functions in the numerical integration, the Coulomb-exchange and infrared singularities appearing in the amplitudes are successfully treated. The matching coefficients are given as a function of the heavy quark mass Ma in lattice units. An automated vertex generation program written in Python is employed, allowing us to use a realistic NRQCD action and an improved gluon lattice action. A change in the definition of either action is easily accommodated in this procedure. The final result is applicable to the electromagnetic decays of heavy quarkonia, notably the Upsilon meson.


title: A Framework for Systematic Study of QCD Vacuum Structure I: Kolmogorov Entropy and the Principle of Chiral Ordering

abstract: In this series of articles we describe a systematic approach to studying QCD vacuum structure using the methods of lattice gauge theory. Our framework incorporates four major components. (i) The recently established existence of space-time order at all scales (fundamental structure) observed directly in typical configurations of regularized path-integral ensembles. (ii) The notion of scale-dependent vacuum structure (effective structure) providing the means for representing and quantifying the influence of fluctuations at various scales on physical observables (phenomena). (iii) The unified description of gauge and fermionic aspects of the theory which facilitates a high level of space-time order in the path-integral ensembles. (iv) The strict ``Bottom-Up'' approach wherein the process of identifying the vacuum structure proceeds inductively, using the information from valid lattice QCD ensembles as the only input. In this work we first elaborate on the meaning of the notion of space-time order in a given configuration which is conceptually at the heart of the path-integral approach to vacuum structure. It is argued that the algorithmic complexity of binary strings associated with coarse-grained descriptions of the configuration provides a relevant quantitative measure. The corresponding ensemble averages define the ranking of different lattice theories at given cutoff by the degree of space-time order generated via their dynamics. We then introduce the set of local transformations of a configuration, chiral orderings, in which the transformed gauge connection represents an effective matrix phase acquired by chiral fermion when hopping over a given link. It is proposed that chiral orderings facilitate the evolution in the set of actions which increases the degree of space-time order while preserving the physical content of the theory, and should thus be used in the search for the fundamental QCD vacuum structure. The relation to renormalization group ideas is discussed, and the first step in general formulation of effective lattice QCD realizing the notion of scale-dependent vacuum structure is given.


title: First lattice study of low-energy charmonium-hadron interaction

abstract: We study the scattering lengths of charmonia (J/psi and eta_c) with light hadrons (pi, rho and N) by the quenched lattice QCD simulations on 24x24x24x48, 32x32x32x48 and 48x48x48x48 lattices with the lattice spacing a = 0.068 fm. The scattering length is extracted by using the Luscher's phase-shift formula together with the measurement of the energy shift Delta E of two hadrons on the lattice. We find that there exist attractive interactions in all channels, J/psi(eta_c)-pi, J/psi(eta_c)-rho and J/psi(eta_c)-N: The s-wave J/psi-pi (eta_c-pi) scattering length is determined as 0.0119+-0.0039 fm (0.0113+-0.0035 fm) and the corresponding elastic cross section at the threshold becomes 0.018+0.013-0.010 mb (0.016+0.011-0.008 mb). Also, the J/psi-N (eta_c-N) spin-averaged scattering length is 0.71+-0.48 fm (0.70+-0.66 fm), which is at least an order of magnitude larger than the charmonium-pion scattering length. The volume dependence of the energy shifts is also investigated to check the expected 1/L^3 behavior of Delta E at a large spatial size L.


title: The nucleon axial-vector coupling beyond one loop

abstract: We analyze the nucleon axial-vector coupling to two loops in chiral perturbation theory. We show that chiral extrapolations based on this representation require lattice data with pion masses below 300 MeV.


title: Exploratory study of three-point Green's functions in Landau-gauge Yang-Mills theory

abstract: Green's functions are a central element in the attempt to understand non-perturbative phenomena in Yang-Mills theory. Besides the propagators, 3-point Green's functions play a significant role, since they permit access to the running coupling constant and are an important input in functional methods. Here we present numerical results for the two non-vanishing 3-point Green's functions in 3d pure SU(2) Yang-Mills theory in (minimal) Landau gauge, i.e. the three-gluon vertex and the ghost-gluon vertex, considering various kinematical regimes. In this exploratory investigation the lattice volumes are limited to 20^3 and 30^3 at beta=4.2 and beta=6.0. We also present results for the gluon and the ghost propagators, as well as for the eigenvalue spectrum of the Faddeev-Popov operator. Finally, we compare two different numerical methods for the evaluation of the inverse of the Faddeev-Popov matrix, the point-source and the plane-wave-source methods.


title: Phase Structure of Z(3)-Polyakov-Loop Models

abstract: We study effective lattice actions describing the Polyakov loop dynamics originating from finite-temperature Yang-Mills theory. Starting with a strong-coupling expansion the effective action is obtained as a series of Z(3)-invariant operators involving higher and higher powers of the Polyakov loop, each with its own coupling. Truncating to a subclass with two couplings we perform a detailed analysis of the statistical mechanics involved. To this end we employ a modified mean field approximation and Monte Carlo simulations based on a novel cluster algorithm. We find excellent agreement of both approaches concerning the phase structure of the theories. The phase diagram exhibits both first and second order transitions between symmetric, ferromagnetic and anti-ferromagnetic phases with phase boundaries merging at three tricritical points. The critical exponents nu and gamma at the continuous transition between symmetric and anti-ferromagnetic phases are the same as for the 3-state Potts model.


title: Overlap Hypercube Fermions in QCD Simulations Near the Chiral Limit

abstract: The overlap hypercube fermion is a variant of a chirally symmetric lattice fermion, which is endowed with a higher level of locality than the standard overlap fermion. We apply this formulation in quenched QCD simulations with light quarks. In the p-regime we evaluate the masses of light pseudoscalar and vector mesons, as well as the pion decay constant and the renormalisation constant Z_A. In the epsilon-regime we present results for the leading Low Energy Constants of the chiral Lagrangian, Sigma and F_pi. To this end, we perform fits to predictions by chiral Random Matrix Theory and by different versions of quenched Chiral Perturbation Theory, referring to distinct correlation functions. These results, along with an evaluation of the topological susceptibility, are also compared to the outcome based on the standard overlap operator.


title: Strong-Isospin Violation in the Neutron-Proton Mass Difference from Fully-Dynamical Lattice QCD and PQQCD

abstract: We determine the strong-isospin violating component of the neutron-proton mass difference from fully-dynamical lattice QCD and partially-quenched QCD calculations of the nucleon mass, constrained by partially-quenched chiral perturbation theory at one-loop level. The lattice calculations were performed with domain-wall valence quarks on MILC lattices with rooted staggered sea-quarks at a lattice spacing of b=0.125 fm, lattice spatial size of L=2.5 fm and pion masses ranging from m_pi ~ 290 MeV to ~ 350 MeV. At the physical value of the pion mass, we predict M_n - M_p |(d-u) = 2.26 +- 0.57 +- 0.42 +- 0.10 MeV where the first error is statistical, the second error is due to the uncertainty in the ratio of light-quark masses, eta=m_u/m_d, determined by MILC, and the third error is an estimate of the systematic due to chiral extrapolation.


title: Dynamic Critical Behavior of Percolation Observables in the 2d Ising Model

abstract: We present preliminary results of our numerical study of the critical dynamics of percolation observables for the two-dimensional Ising model. We consider the (Monte-Carlo) short-time evolution of the system obtained with a local heat-bath method and with the global Swendsen-Wang algorithm. In both cases, we find qualitatively different dynamic behaviors for the magnetization and Omega, the order parameter of the percolation transition. This may have implications for the recent attempts to describe the dynamics of the QCD phase transition using cluster observables.


title: K^0-\bar K^0 mixing beyond the standard model and CP-violating electroweak penguins in quenched QCD with exact chiral symmetry

abstract: We present results for the \Delta S=2 matrix elements which are required to study neutral kaon mixing in the standard model (SM) and beyond (BSM). We also provide leading chiral order results for the matrix elements of the electroweak penguin operators which give the dominant \Delta I=3/2 contribution to direct CP violation in K->\pi\pi decays. Our calculations were performed with Neuberger fermions on two sets of quenched Wilson gauge configurations at inverse lattice spacings of approximately 2.2 GeV and 1.5 GeV. All renormalizations were implemented non-perturbatively in the RI/MOM scheme, where we accounted for sub-leading operator product expansion corrections and discretization errors. We find ratios of non-SM to SM matrix elements which are roughly twice as large as in the only other dedicated lattice study of these amplitudes. On the other hand, our results for the electroweak penguin matrix elements are in good agreement with two recent domain-wall fermion calculations. As a by-product of our study, we determine the strange quark mass. Our main results are summarized and discussed in Sec. VII. Within our statistics, we find no evidence for scaling violations.


title: The nucleon electromagnetic form factors from Lattice QCD

abstract: We evaluate the isovector nucleon electromagnetic form factors in quenched and full QCD on the lattice using Wilson fermions. In the quenched theory we use a lattice of spatial size 3 fm at beta=6.0 enabling us to reach low momentum transfers and a lowest pion mass of about 400 MeV. In the full theory we use a lattice of spatial size 1.9 fm at beta=5.6 and lowest pion mass of about 380 MeV enabling comparison with the results obtained in the quenched theory. We compare our lattice results to the isovector part of the experimentally measured form factors.


title: Linking confinement to spectral properties of the Dirac operator

abstract: We represent Polyakov loops and their correlators as spectral sums of eigenvalues and eigenmodes of the lattice Dirac operator. The deconfinement transition of pure gauge theory is characterized as a change in the response of moments of eigenvalues to varying the boundary conditions of the Dirac operator. We argue that the potential between static quarks is linked to spatial correlations of Dirac eigenvectors.


title: Non-Perturbative U(1) Gauge Theory at Finite Temperature

abstract: For compact U(1) lattice gauge theory (LGT) we have performed a finite size scaling analysis on $N_{\tau} N_s^3$ lattices for $N_{\tau}$ fixed by extrapolating spatial volumes of size $N_s\le 18$ to $N_s\to\infty$. Within the numerical accuracy of the thus obtained fits we find for $N_{\tau}=4$, 5 and~6 second order critical exponents, which exhibit no obvious $N_{\tau}$ dependence. The exponents are consistent with 3d Gaussian values, but not with either first order transitions or the universality class of the 3d XY model. As the 3d Gaussian fixed point is known to be unstable, the scenario of a yet unidentified non-trivial fixed point close to the 3d Gaussian emerges as one of the possible explanations.


title: Improved Perturbation Theory for Improved Lattice Actions

abstract: We study a systematic improvement of perturbation theory for gauge fields on the lattice; the improvement entails resumming, to all orders in the coupling constant, a dominant subclass of tadpole diagrams. This method, originally proposed for the Wilson gluon action, is extended here to encompass all possible gluon actions made of closed Wilson loops; any fermion action can be employed as well. The effect of resummation is to replace various parameters in the action (coupling constant, Symanzik coefficients, clover coefficient) by ``dressed'' values; the latter are solutions to certain coupled integral equations, which are easy to solve numerically. Some positive features of this method are: a) It is gauge invariant, b) it can be systematically applied to improve (to all orders) results obtained at any given order in perturbation theory, c) it does indeed absorb in the dressed parameters the bulk of tadpole contributions. Two different applications are presented: The additive renormalization of fermion masses, and the multiplicative renormalization Z_V (Z_A) of the vector (axial) current. In many cases where non-perturbative estimates of renormalization functions are also available for comparison, the agreement with improved perturbative results is significantly better as compared to results from bare perturbation theory.


title: Discretization errors in the spectrum of the Hermitian Wilson-Dirac operator

abstract: I study the leading effects of discretization errors on the low energy part of the spectrum of the Hermitian Wilson-Dirac operator in infinite volume. The method generalizes that used to study the spectrum of the Dirac operator in the continuum, and uses partially quenched chiral perturbation theory for Wilson fermions. The leading-order corrections are proportional to a^2 (a being the lattice spacing). At this order I find that the method works only for one choice of sign of one of the three low energy constants describing discretization errors. If these constants have the relative magnitudes expected from large N_c arguments, then the method works if the theory has an Aoki phase for m of order a^2, but fails if there is a first-order transition. In the former case, the dependence of the gap and the spectral density on m and a^2 are determined. In particular, the gap is found to vanish more quickly as m_pi^2-> 0 than in the continuum. This reduces the region where simulations are safe from fluctuations in the gap.


title: Power-counting theorem for staggered fermions

abstract: Lattice power-counting is extended to QCD with staggered fermions. As preparation, the difficulties encountered by Reisz's original formulation of the lattice power-counting theorem are illustrated. One of the assumptions that is used in his proof does not hold for staggered fermions, as was pointed out long ago by Luscher. Finally, I generalize the power-counting theorem, and the methods of Reisz's proof, such that the difficulties posed by staggered fermions are overcome.


title: Exploring the HMC trajectory-length dependence of autocorrelation times in lattice QCD

abstract: We study autocorrelation times of physical observables in lattice QCD as a function of the molecular dynamics trajectory length in the hybrid Monte-Carlo algorithm. In an interval of trajectory lengths where energy and reversibility violations can be kept under control, we find a variation of the integrated autocorrelation times by a factor of about two in the quantities of interest. Trajectories longer than conventionally used are found to be superior both in the Nf=0 and Nf=2 examples considered here. We also provide evidence that they lead to faster thermalization of systems with light quarks.


title: Exotic Meson Decay Widths using Lattice QCD

abstract: A decay width calculation for a hybrid exotic meson h, with JPC=1-+, is presented for the channel h->pi+a1. This quenched lattice QCD simulation employs Luescher's finite box method. Operators coupling to the h and pi+a1 states are used at various levels of smearing and fuzzing, and at four quark masses. Eigenvalues of the corresponding correlation matrices yield energy spectra that determine scattering phase shifts for a discrete set of relative pi+a1 momenta. Although the phase shift data is sparse, fits to a Breit-Wigner model are attempted, resulting in a decay width of about 60 MeV when averaged over two lattice sizes.


title: Search for the Theta^+(1540) in lattice QCD

abstract: We report on a study of the pentaquark Theta^+(1540), using a variety of different interpolating fields. We use Chirally Improved fermions in combination with Jacobi smeared quark sources to improve the signal and get reliable results even for small quark masses. The results of our quenched calculations, which have been done on a 12^3 x 24 lattice with a lattice spacing of a = 0.148 fm, do not provide any evidence for the existence of a \Theta^+ with positive parity. We do observe, however, a signal compatible with nucleon-kaon scattering state. For the negative parity the results are inconclusive, due to the potential mixture with nucleon-kaon and N^*-kaon scattering states.


title: Quarkonia correlators and spectral functions from lattice QCD

abstract: I discuss recent progress in calculating quarkonia correlators and spectral functions on the lattice in relation with the problem of quarkonia dissolution at high temperatures and heavy quark transport in Quark Gluon Plasma.


title: Peculiarities in the Spectrum of the Adjoint Scalar Kinetic Operator in Yang-Mills Theory

abstract: We study the spectrum of low-lying eigenmodes of the kinetic operator for scalar particles, in the color adjoint representation of Yang-Mills theory. The kinetic operator is the covariant Laplacian, plus a constant which serves to renormalize mass. In the pure gauge theory, our data indicates that the interval between the lowest eigenvalue and the mobility edge tends to infinity in the continuum limit. On these grounds, it is suggested that the perturbative expression for the scalar propagator may be misleading even at distance scales that are small compared to the confinement scale. We also measure the density of low-lying eigenmodes, and find a possible connection to multi-critical matrix models of order m=1.


title: A gauge invariant cluster algorithm for the Ising spin glass

abstract: The frustrated Ising model in two dimensions is revisited. The frustration is quantified in terms of the number of non-trivial plaquettes which is invariant under the Nishimori gauge symmetry. The exact ground state energy is calculated using Edmond's algorithm. A novel cluster algorithm is designed which treats gauge equivalent spin glasses on equal footing and allows for efficient simulations near criticality. As a first application, the specific heat near criticality is investigated.


title: $Z_6$ symmetry, electroweak transition, and magnetic monopoles at high temperature

abstract: We consider the lattice realization of the Standard Model with an additional $Z_6$ symmetry. Numerical simulations were performed on the asymmetric lattice, which corresponds to the finite temperature theory. Our choice of parameters corresponds to large Higgs masses ($M_H > 90$ Gev). The phase diagram was investigated and has been found to be different from that of the usual lattice realization of the Standard Model. It has been found, that the confinement-deconfinement phase transition lines for the SU(2) and SU(3) fields coincide. The transition line between Higgs and symmetric deconfinement parts of the phase diagram and the confinement-deconfinement transition line meet in a triple point. The transition between Higgs and symmetric parts of the phase diagram corresponds to the finite temperature electroweak transition/crossover. We see for the first time evidence that Nambu monopoles are condensed at $T>T_c$ while at $T<T_c$ their condensate vanishes.


title: Numerical simulation of QCD with u, d, s and c quarks in the twisted-mass Wilson formulation

abstract: A first study of numerical Monte Carlo simulations with two quark doublets, a mass-degenerate one and a mass-split one, interpreted as u, d, s and c quarks, is carried out in the framework of the twisted mass Wilson lattice formulation. Tuning the bare parameters of this theory is explored on 12^3x24 and 16^3x32 lattices at lattice spacings a=0.20fm and a=0.15fm, respectively.


title: Moments of pseudoscalar meson distribution amplitudes from the lattice

abstract: Based on lattice simulations with two flavours of dynamical, O(a)-improved Wilson fermions we present results for the first two moments of the distribution amplitudes of pseudoscalar mesons at several values of the valence quark masses. By extrapolating our results to the physical masses of up/down and strange quarks, we find the first two moments of the K^+ distribution amplitude and the second moment of the pi^+ distribution amplitude. We use nonperturbatively determined renormalisation coefficients to obtain results in the MSbar scheme. At a scale of 4 GeV^2 we find a_2^pi=0.201(114) for the second Gegenbauer moment of the pion's distribution amplitude, while for the kaon, a_1^K=0.0453(9)(29) and a_2^K=0.175(18)(47).


title: Monte Carlo Simulation of a NC Gauge Theory on The Fuzzy Sphere

abstract: We find using Monte Carlo simulation the phase structure of noncommutative U(1) gauge theory in two dimensions with the fuzzy sphere S^2_N as a non-perturbative regulator. There are three phases of the model. i) A matrix phase where the theory is essentially SU(N) Yang-Mills reduced to zero dimension . ii) A weak coupling fuzzy sphere phase with constant specific heat and iii) A strong coupling fuzzy sphere phase with non-constant specific heat. The order parameter distinguishing the matrix phase from the sphere phase is the radius of the fuzzy sphere. The three phases meet at a triple point. We also give the theoretical one-loop and 1/N expansion predictions for the transition lines which are in good agreement with the numerical data. A Monte Carlo measurement of the triple point is also given.


title: Gap Domain Wall Fermions

abstract: I demonstrate that the chiral properties of Domain Wall Fermions (DWF) in the large to intermediate lattice spacing regime of QCD, 1 to 2 GeV, are significantly improved by adding to the action two standard Wilson fermions with supercritical mass equal to the negative DWF five dimensional mass. Using quenched DWF simulations I show that the eigenvalue spectrum of the transfer matrix Hamiltonian develops a substantial gap and that the residual mass decreases appreciatively. Furthermore, I confirm that topology changing remains active and that the hadron spectrum of the added Wilson fermions is above the lattice cutoff and therefore is irrelevant. I argue that this result should also hold for dynamical DWF and furthermore that it should improve the chiral properties of related fermion methods.


title: Dilute Liquid of Instanton and Its Topological Charge Dominate the QCD Vacuum

abstract: APE smearing and overlap-Dirac operator are combined to filter QCD vacuum configurations. The results obtained from overlap fermions and improved 5Li cooling are compared, both of them exhibit structures of dilute liquid of instanton. Finally the overlap fermions, improved 5Li cooling and APE smearing are combined to calculate the topological charge and identify the structure of QCD vacuum. The results suggest dilute liquid of instanton dominance of topological charge fluctuations in quenched lattice QCD.


title: Gluon flux-tube distribution and linear confinement in baryons

abstract: We have observed the formation of gluon flux-tubes within baryons using lattice QCD techniques. A high-statistics approach, based on translational and rotational symmetries of the four-dimensional lattice, enables us to observe correlations between vacuum action density and quark positions in a completely gauge independent manner. This contrasts with earlier studies which used gauge-dependent smoothing techniques. We used 200 O(a^2) improved quenched QCD gauge-field configurations on a 16^3x32 lattice with a lattice spacing of 0.123 fm. In the presence of static quarks flux tubes representing the suppression of gluon-field fluctuations are observed. We have analyzed 11 L-shapes and 8 T and Y shapes of varying sizes in order to explore a variety of flux-tube topologies, including the ground state. At large separations, Y-shape flux-tube formation is observed. T-shaped paths are observed to relax towards a Y-shaped topology, whereas L-shaped paths give rise to a large potential energy. We do not find any evidence for the formation of a Delta-shaped flux-tube (empty triangle) distribution. However, at small quark separations, we observe an expulsion of gluon-field fluctuations in the shape of a filled triangle with maximal expulsion at the centre of the triangle. Having identified the precise geometry of the flux distribution, we are able to perform quantitative comparison between the length of the flux-tube and the associated static quark potential. For every source configuration considered we find a universal string tension, and conclude that, for large quark separations, the ground state potential is that which minimizes the length of the flux-tube. The flux tube radius of the baryonic ground state potential is found to be 0.38 \pm 0.03 fm, with vacuum fluctuations suppressed by 7.2 \pm 0.6 %.


title: The eta mass and NNLO Three-Flavor Partially Quenched Chiral Perturbation Theory

abstract: We show how to resum neutral propagators to all orders in Partially Quenched Chiral Perturbation Theory. We calculate the relevant quantities to next-to-next-to-leading order (NNLO). Using these results we show how to extend the proposal by Sharpe and Shoresh for determining the parameters relevant for the eta mass from partially quenched lattice QCD calculations to NNLO.


title: Robustness of baryon-strangeness correlation and related ratios of susceptibilities

abstract: Using quenched lattice QCD simulations we investigate the continuum limit of baryon-strangeness correlation and other related conserved charge-flavour correlations for temperatures T_c<T\le2T_c. By working with lattices having large temporal extents (N_\tau=12, 10, 8, 4) we find that these quantities are almost independent of the lattice spacing, i.e, robust. We also find that these quantities have very mild dependence on the sea quark mass and acquire values which are very close to their respective ideal gas limits. Our results also confirm robustness of the Wroblewski parameter.


title: SU(3) gauge theory at finite temperature in 2 + 1 dimensions

abstract: In this article we will discuss numerical results on screening masses and thermodynamic quantities in 2 + 1 dimensional SU(3) gauge theory. We will also compare them to perturbation theory and the dimensionally reduced model.


title: QCD Thermodynamics with Almost Realistic Quark Masses

abstract: Ongoing calculations on the QCDOC supercomputer at Brookhaven National Laboratory and the APEnext installation at the University of Bielefeld aim to determine the critical temperature of the QCD phase transition as well as the equation of state with almost realistic quark masses. We will discuss preliminary results of the quark mass and cut-off dependence of order parameters, susceptibilities, static quark potentials and the critical temperature in (2+1)-flavor QCD. All these quantities are of immediate interest for heavy ion phenomenology.


title: Lattice Regularization and Symmetries

abstract: Finding the relation between the symmetry transformations in the continuum and on the lattice might be a nontrivial task as illustrated by the history of chiral symmetry. Lattice actions induced by a renormalization group procedure inherit all symmetries of the continuum theory. We give a general procedure which gives the corresponding symmetry transformations on the lattice.


title: Chiral violations in domain-wall QCD from one-loop perturbation theory at finite $N_s$

abstract: We present perturbative calculations made with domain-wall fermions which possess a finite number of points $N_s$ in the extra fifth dimension. We have derived the required propagator functions, investigated the one-loop properties of quark amplitudes at finite $N_s$ and evaluated three quantities that can provide insights on chirality-breaking effects from the perturbative side. First we have computed the residual mass for various choices of $N_s$ and of the domain-wall height $M$. We have found that this radiatively induced mass approaches zero reasonably fast with the extent of the fifth dimension, depending on $M$ and on a lesser extent on the coupling $g_0$. We have also computed the differences of the renormalization constants of the vector and axial-vector currents and of the scalar and pseudoscalar densities. Finally we have calculated the chirally-forbidden mixing (which at finite $N_s$ is suppressed only partially) of an operator which describes the lowest moment of the $g_2$ structure function. In general we see that at $M=1.8$, where simulations are usually performed, values of $N_s=20$ or larger would be desirable in order for chiral violations to be negligible. The quantities that we have studied turn out to lose gauge invariance when $N_s$ is not infinite. We have also found that anomalous dimensions of operators at finite $N_s$ generally depend on $N_s$ and $M$. In particular, the vector and axial-vector currents have in general a nonzero anomalous dimension at finite $N_s$.


title: f_K/f_pi in Full QCD with Domain Wall Valence Quarks

abstract: We compute the ratio of pseudoscalar decay constants f_K/f_pi using domain-wall valence quarks and rooted improved Kogut-Susskind sea quarks. By employing continuum chiral perturbation theory, we extract the Gasser-Leutwyler low-energy constant L_5, and extrapolate f_K/f_pi to the physical point. We find: f_K/f_pi = 1.218 (+- 0.002) (+0.011 -0.024) where the first error is statistical and the second error is an estimate of the systematic due to chiral extrapolation and fitting procedures. This value agrees within the uncertainties with the determination by the MILC collaboration, calculated using Kogut-Susskind valence quarks, indicating that systematic errors arising from the choice of lattice valence quark are small.


title: Extrapolations of Lattice Meson Form Factors

abstract: We use chiral perturbation theory to study the extrapolations necessary to make physical predictions from lattice QCD data for the electromagnetic form factors of pseudoscalar mesons. We focus on the quark mass, momentum, lattice spacing, and volume dependence and apply our results to simulations employing mixed actions of Ginsparg-Wilson valence quarks and staggered sea quarks. To determine charge radii at quark masses on the lattices currently used, we find that all extrapolations except the one to infinite volume make significant contributions to the systematic error.


title: Proton decay matrix elements with domain-wall fermions

abstract: Hadronic matrix elements of operators relevant to nucleon decay in grand unified theories are calculated numerically using lattice QCD. In this context, the domain-wall fermion formulation, combined with non-perturbative renormalization, is used for the first time. These techniques bring reduction of a large fraction of the systematic error from the finite lattice spacing. Our main effort is devoted to a calculation performed in the quenched approximation, where the direct calculation of the nucleon to pseudoscalar matrix elements, as well as the indirect estimate of them from the nucleon to vacuum matrix elements, are performed. First results, using two flavors of dynamical domain-wall quarks for the nucleon to vacuum matrix elements are also presented to address the systematic error of quenching, which appears to be small compared to the other errors. Our results suggest that the representative value for the low energy constants from the nucleon to vacuum matrix elements are given as |alpha| simeq |beta| simeq 0.01 GeV^3. For a more reliable estimate of the physical low energy matrix elements, it is better to use the relevant form factors calculated in the direct method. The direct method tends to give smaller value of the form factors, compared to the indirect one, thus enhancing the proton life-time; indeed for the pi^0 final state the difference between the two methods is quite appreciable.


title: Once more on the interrelation between Abelian monopoles and P-vortices in SU(2) LGT

abstract: We study the properties of configurations from which P-vortices on one hand or Abelian monopoles on the other hand have been removed. We confirm the loss of confinement in both cases and investigate in what respect the modified ensembles differ from the confining ones from the point of view of the complementary confinement scenario.


title: Multiparticle States and the Hadron Spectrum on the Lattice

abstract: The Clebsch-Gordan decomposition is calculated for direct products of the irreducible representations of the cubic space group. These results are used to identify multiparticle states which appear in the hadron spectrum on the lattice. Consideration of the cubic space group indicates how combinations of both zero momentum and non-zero momentum multiparticle states contribute to the spectrum.


title: The Effect of Reduced Spatial Symmetries on Lattice States: Results for Non-zero Linear Momentum

abstract: The irreducible representations of the cubic space group are described and used to determine the mapping of continuum states to lattice states with non-zero linear momentum. The Clebsch-Gordan decomposition is calculated from the character table for the cubic space group. These results are used to identify multiparticle states which appear in the hadron spectrum on the lattice.


title: Rationale for UV-filtered clover fermions

abstract: We study the contributions Sigma_0 and Sigma_1, proportional to a^0 and a^1, to the fermion self-energy in Wilson's formulation of lattice QCD with UV-filtering in the fermion action. We derive results for m_{crit} and the renormalization factors Z_S, Z_P, Z_V, Z_A to 1-loop order in perturbation theory for several filtering recipes (APE, HYP, EXP, HEX), both with and without a clover term. The perturbative series is much better behaved with filtering, in particular tadpole resummation proves irrelevant. Our non-perturbative data for m_{crit} and Z_A/(Z_m*Z_P) show that the combination of filtering and clover improvement efficiently reduces the amount of chiral symmetry breaking -- we find residual masses am_{res}=O(10^{-2}).


title: Renormalization-group analysis of the validity of staggered-fermion QCD with the fourth-root recipe

abstract: I develop a renormalization-group blocking framework for lattice QCD with staggered fermions. Under plausible, and testable, assumptions, I then argue that the fourth-root recipe used in numerical simulations is valid in the continuum limit. The taste-symmetry violating terms, which give rise to non-local effects in the fourth-root theory when the lattice spacing is non-zero, vanish in the continuum limit. A key role is played by reweighted theories that are local and renormalizable on the one hand, and that approximate the fourth-root theory better and better as the continuum limit is approached on the other hand.


title: Calculation of the nucleon axial charge in lattice QCD

abstract: Protons and neutrons have a rich structure in terms of their constituents, the quarks and gluons. Understanding this structure requires solving Quantum Chromodynamics (QCD). However QCD is extremely complicated, so we must numerically solve the equations of QCD using a method known as lattice QCD. Here we describe a typical lattice QCD calculation by examining our recent computation of the nucleon axial charge.


title: Nonperturbative determination of the QCD potential at O(1/m)

abstract: The relativistic correction to the QCD static inter-quark potential at O(1/m) is investigated nonperturbatively for the first time by using lattice Monte Carlo QCD simulations. The correction is found to be comparable with the Coulombic term of the static potential when applied to charmonium, and amounts to one-fourth of the Coulombic term for bottomonium.


title: Twisted supersymmetric sigma model on the lattice

abstract: In this paper we conduct a numerical study of the supersymmetric O(3) non-linear sigma model. The lattice formulation we employ was derived in \cite{sigma1} and corresponds to a discretization of a {\it twisted} form of the continuum action. The twisting process exposes a {\it nilpotent} supercharge $Q$ and allows the action to be rewritten in $Q$-exact form. These properties may be maintained on the lattice. We show how to deform the theory by the addition of potential terms which preserve the supersymmetry. A Wilson mass operator may be introduced in this way with a minimal breaking of supersymmetry. We additionally show how to rewrite the theory in the language of K\"{a}hler-Dirac fields and explain why this avenue does not provide a good route to discretization. Our numerical results provide strong evidence for a restoration of full supersymmetry in the continuum limit {\it without} fine tuning. We also observe a non-vanishing chiral condensate as expected from continuum instanton calculations.


title: Predictions with Lattice QCD

abstract: In recent years, we used lattice QCD to calculate some quantities that were unknown or poorly known. They are the $q^2$ dependence of the form factor in semileptonic $D\to Kl\nu$ decay, the leptonic decay constants of the $D^+$ and $D_s$ mesons, and the mass of the $B_c$ meson. In this paper, we summarize these calculations, with emphasis on their (subsequent) confirmation by measurements in $e^+e^-$, $\gamma p$ and $\bar{p}p$ collisions.


title: Meson spectral functions at nonzero momentum in hot QCD

abstract: We present preliminary results for meson spectral functions at nonzero momentum, obtained from quenched lattice QCD simulations at finite temperature using the Maximal Entropy Method. Twisted boundary conditions are used to have access to many momenta p~T. For light quarks, we observe a drastic modification when heating the system from below to above Tc. In particular, for the vector spectral density we find a nonzero spectral weight at all energies.


title: Spontaneous breaking of remnant gauge symmetries in zero-temperature SU(2) lattice gauge theory

abstract: The 4-d SU(2) lattice gauge theory is simulated in the minimal Coulomb gauge which aims to maximize the traces of all links in three directions. Fourth-direction links are interpreted as spins in a Heisenberg-like model with varying interactions. These spins magnetize in 3-d hyperlayers at weak coupling, breaking a remnant gauge symmetry, as well as the Polyakov-loop symmetry. They demagnetize at a phase transition around $\beta = 2.5$ on the infinite lattice, as determined by Binder cumulant crossings. Because $N$ symmetries are breaking on an $N^4$ lattice, the transition is unusually broad, encompassing most of the crossover region on typical lattices.


title: Dual superconductivity and vacuum properties in Yang--Mills theories

abstract: We address, within the dual superconductivity model for color confinement, the question whether the Yang-Mills vacuum behaves as a superconductor of type I or type II. In order to do that we compare, for the theory with gauge group SU(2), the determination of the field penetration depth $\lambda$ with that of the superconductor correlation length $\xi$. The latter is obtained by measuring the temporal correlator of a disorder parameter developed by the Pisa group to detect dual superconductivity. The comparison places the vacuum close to the border between type I and type II and marginally on the type II side. We also check our results against the study of directly measurable effects such as the interaction between two parallel flux tubes, obtaining consistent indications for a weak repulsive behaviour. Future strategies to improve our investigation are discussed.


title: Static forces in d=2+1 SU(N) gauge theories

abstract: Using a three-level algorithm we perform a high-precision lattice computation of the static force up to 1fm in the 2+1 dimensional SU(5) gauge theory. Discretization errors and the continuum limit are discussed in detail. By comparison with existing SU(2) and SU(3) data it is found that \sigma r_0^2=1.65-\pi/24 holds at an accuracy of 1% for all N>=2, where r_0 is the Sommer reference scale. The effective central charge c_{eff}(r) is obtained and an intermediate distance r_s is defined via the property c_{eff}(r_s)=\pi/24. It separates in a natural way the short-distance regime governed by perturbation theory from the long-distance regime described by an effective string theory. The ratio r_s/r_0 decreases significantly from SU(2) to SU(3) to SU(5), where r_s < r_0. We give a preliminary estimate of its value in the large-N limit. The static force in the smallest representation of N-ality 2, which tends to the k=2 string tension as r->oo, is also computed up to 0.7fm. The deviation from Casimir scaling is positive and grows from 0.1% to 1% in that range.


title: Applications of Chiral Perturbation theory to lattice QCD

abstract: These lectures describe the use of effective field theories to extrapolate results from the parameter region where numerical simulations of lattice QCD are possible to the physical parameters (physical quark masses, infinite volume, vanishing lattice spacing, etc.). After a brief introduction and overview, I discuss three topics: 1) Chiral perturbation theory in the continuum; 2) The inclusion of discretization effects into chiral perturbation theory, focusing on the application to Wilson and twisted-mass lattice fermions; 3) Extending chiral perturbation theory to describe partially quenched QCD.


title: The chiral critical line of N_f=2+1 QCD at zero and non-zero baryon density

abstract: We present numerical results for the location of the chiral critical line at finite temperature and zero and non-zero baryon density for QCD with N_f=2+1 flavours of staggered fermions on lattices with temporal extent N_t=4. For degenerate quark masses, we compare our results obtained with the exact RHMC algorithm with earlier, inexact R-algorithm results and find a reduction of 25% in the critical quark mass, for which the first order phase transition changes to a smooth crossover. Extending our analysis to non-degenerate quark masses, we map out the chiral critical line up to the neighbourhood of the physical point, which we confirm to be in the crossover region. Our data are consistent with a tricritical point at a strange quark mass of ~500 MeV. Finally, we investigate the shift of the critical line with finite baryon density, by simulating with an imaginary chemical potential for which there is no sign problem. We observe this shift to be very small or, conversely, the critical endpoint \mu^c(m_{u,d},m_s) to be extremely quark mass sensitive. Moreover, the sign of this shift is opposite to standard expectations. If confirmed on a finer lattice, it implies the absence of a critical endpoint or phase transition for chemical potentials \mu_B < 500 MeV. We thus argue that finer lattices are required to settle even the qualitative features of the QCD phase diagram.


title: A Lattice Computation of the First Moment of the Kaon's Distribution Amplitude

abstract: We present a lattice computation of the first moment of the kaon's leading-twist distribution amplitude. The results were computed using ensembles with 2+1 dynamical flavours with the domain wall fermion action and Iwasaki gauge action from the RBC and UKQCD joint dataset. The first moment is non-zero because of SU(3)-breaking effects, and we find that we are able to measure these effects very clearly. We observe the expected chiral behaviour and finally obtain <\xi>(2GeV) = 3/5*a_K^1(2GeV) = 0.032(3), which agrees very well with results obtained using sum-rules, but with a significantly smaller error.


title: Flavor Twisted Boundary Conditions and Isovector Form Factors

abstract: We use vector flavor symmetry to relate form factors of isospin changing operators to isovector form factors. Flavor twisted boundary conditions in lattice QCD thus allow isovector form factors of twist-two operators, e.g, to be computed at continuous values of the momentum transfer. These twisted boundary conditions, moreover, are implemented only in the valence sector. Effects of the finite volume must be addressed to extract isovector moments and radii at zero lattice momentum. As an example, we use chiral perturbation theory to assess the volume effects in extracting the isovector magnetic moment of the nucleon from simulations with twisted boundary conditions.


title: Lattice gauge action suppressing near-zero modes of H_W

abstract: We propose a lattice action including unphysical Wilson fermions with a negative mass m_0 of the order of the inverse lattice spacing. With this action, the exact zero mode of the hermitian Wilson-Dirac operator H_W(m_0) cannot appear and near-zero modes are strongly suppressed. By measuring the spectral density rho(lambda_W), we find a gap near lambda_W=0 on the configurations generated with the standard and improved gauge actions. This gap provides a necessary condition for the proof of the exponential locality of the overlap-Dirac operator by Hernandez, Jansen, and Luescher. Since the number of near-zero modes is small, the numerical cost to calculate the matrix sign function of H_W(m_0) is significantly reduced, and the simulation including dynamical overlap fermions becomes feasible. We also introduce a pair of twisted mass pseudo-fermions to cancel the unwanted higher mode effects of the Wilson fermions. The gauge coupling renormalization due to the additional fields is then minimized. The topological charge measured through the index of the overlap-Dirac operator is conserved during continuous evolutions of gauge field variables.


title: Massive Domain Wall Fermions on Four-dimensional Anisotropic Lattices

abstract: We formulate the massive domain wall fermions on anisotropic lattices. For the massive domain wall fermion, we find that the dispersion relation assumes the usual form in the low momentum region when the bare parameters are properly tuned. The quark self-energy and the quark field renormalization constants are calculated to one-loop in bare lattice perturbation theory. For light domain wall fermions, we verified that the chiral mode is stable against quantum fluctuations on anisotropic lattices. This calculation serves as a guidance for the tuning of the parameters in the quark action in future numerical simulations.


title: Determining F_pi from spectral sum rules

abstract: We derive spectral sum rules for a system with two quarks coupled to an imaginary isospin chemical potential in the \epsilon regime. The sum rules show an explicit dependence on the pion decay constant which should make it possible to measure F_pi from the eigenvalue spectrum of this particular Dirac operator.


title: Doubly Heavy Baryons and Quark-Diquark Symmetry in Quenched and Partially Quenched Chiral Perturbation Theory

abstract: We extend the chiral Lagrangian with heavy quark-diquark symmetry to quenched and partially quenched theories. These theories are used to derive formulae for the chiral extrapolation of masses and hyperfine splittings of doubly heavy baryons in lattice QCD simulations. A quark-diquark symmetry prediction for the hyperfine splittings of heavy mesons and doubly heavy baryons is rather insensitive to chiral corrections in both quenched and partially quenched QCD. Extrapolation formulae for the doubly heavy baryon electromagnetic transition moments are also determined for the partially quenched theory.


title: Hadron spectrum, quark masses and decay constants from light overlap fermions on large lattices

abstract: We present results from a simulation of quenched overlap fermions with L\"uscher-Weisz gauge field action on lattices up to $24^3 48$ and for pion masses down to $\approx 250$ MeV. Among the quantities we study are the pion, rho and nucleon masses, the light and strange quark masses, and the pion decay constant. The renormalization of the scalar and axial vector currents is done nonperturbatively in the $RI-MOM$ scheme. The simulations are performed at two different lattice spacings, $a \approx 0.1$ fm and $\approx 0.15$ fm, and on two different physical volumes, to test the scaling properties of our action and to study finite volume effects. We compare our results with the predictions of chiral perturbation theory and compute several of its low-energy constants. The pion mass is computed in sectors of fixed topology as well.


title: A new lattice measurement for potentials between static SU(3) sources

abstract: In this article, a new calculation of static potentials between sources of different representations in SU(3) gauge group is presented. The results of author's previous study \cite{Deld00} at the smallest lattice spacing $a_{s}\simeq0.11$~ fm are shown to have been affected by finite volume effects. Within statistical errors, the new results obtained here are still in agreement with both, Casimir scaling and flux tube counting. There is also no contradiction to the results obtained in Ref.~ \cite{Bali00} which however exclude flux counting.


title: Small Instantons in $CP^1$ and $CP^2$ Sigma Models

abstract: The anomalous scaling behavior of the topological susceptibility $\chi_t$ in two-dimensional $CP^{N-1}$ sigma models for $N\leq 3$ is studied using the overlap Dirac operator construction of the lattice topological charge density. The divergence of $\chi_t$ in these models is traced to the presence of small instantons with a radius of order $a$ (= lattice spacing), which are directly observed on the lattice. The observation of these small instantons provides detailed confirmation of L\"{u}scher's argument that such short-distance excitations, with quantized topological charge, should be the dominant topological fluctuations in $CP^1$ and $CP^2$, leading to a divergent topological susceptibility in the continuum limit. For the $\CP$ models with $N>3$ the topological susceptibility is observed to scale properly with the mass gap. These larger $N$ models are not dominated by instantons, but rather by coherent, one-dimensional regions of topological charge which can be interpreted as domain wall or Wilson line excitations and are analogous to D-brane or ``Wilson bag'' excitations in QCD. In Lorentz gauge, the small instantons and Wilson line excitations can be described, respectively, in terms of poles and cuts of an analytic gauge potential.


title: Probing the chiral weak Hamiltonian at finite volumes

abstract: Non-leptonic kaon decays are often described through an effective chiral weak Hamiltonian, whose couplings ("low-energy constants") encode all non-perturbative QCD physics. It has recently been suggested that these low-energy constants could be determined at finite volumes by matching the non-perturbatively measured three-point correlation functions between the weak Hamiltonian and two left-handed flavour currents, to analytic predictions following from chiral perturbation theory. Here we complete the analytic side in two respects: by inspecting how small ("epsilon-regime") and intermediate or large ("p-regime") quark masses connect to each other, and by including in the discussion the two leading Delta I = 1/2 operators. We show that the epsilon-regime offers a straightforward strategy for disentangling the coefficients of the Delta I = 1/2 operators, and that in the p-regime finite-volume effects are significant in these observables once the pseudoscalar mass M and the box length L are in the regime ML \lsim 5.0.


title: Non-perturbative renormalisation of left-left four-fermion operators with Neuberger fermions

abstract: We outline a general strategy for the non-perturbative renormalisation of composite operators in discretisations based on Neuberger fermions, via a matching to results obtained with Wilson-type fermions. As an application, we consider the renormalisation of the four-quark operators entering the Delta S=1 and Delta S=2 effective Hamiltonians. Our results are an essential ingredient for the determination of the low-energy constants governing non-leptonic kaon decays.


title: Domain decomposition improvement of quark propagator estimation

abstract: Applying domain decomposition to the lattice Dirac operator and the associated quark propagator, we arrive at expressions which, with the proper insertion of random sources therein, can provide improvement to the estimation of the propagator. Schemes are presented for both open and closed (or loop) propagators. In the end, our technique for improving open contributions is similar to the ``maximal variance reduction'' approach of Michael and Peisa, but contains the advantage, especially for improved actions, of dealing directly with the Dirac operator. Using these improved open propagators for the Chirally Improved operator, we present preliminary results for the static-light meson spectrum. The improvement of closed propagators is modest: on some configurations there are signs of significant noise reduction of disconnected correlators; on others, the improvement amounts to a smoothening of the same correlators.


title: The axial N to Delta transition form factors from Lattice QCD

abstract: We evaluate the N to Delta axial transition form factors in lattice QCD in the quenched theory, with two degenerate flavors of dynamical Wilson fermions and using domain wall valence fermions with staggered sea quark configurations. We predict the ratio $C_5^A(q^2)/C_3^V(q^2)$ relevant to the parity violating asymmetry and check the validity of the off-diagonal Goldberger-Treiman relation.


title: A Framework for Systematic Study of QCD Vacuum Structure II: Coherent Lattice QCD

abstract: We propose the formulation of lattice QCD wherein all elements of the theory (gauge action, fermionic action, theta-term, and all operators) are constructed from a single object, namely the lattice Dirac operator D with exact chiral symmetry. Several regularizations of this type are suggested via constructing scalar densities (gauge actions) that are explicit functions of D. The simplest of these is based on the proposition that classical limit of density associated with Tr D is (up to an additive constant) proportional to FF, while the corresponding operator is local. The possibilities of explicit interrelations between gauge and fermionic aspects of the theory are emphasized together with the utility of such formulations for exploring the QCD vacuum structure.


title: The decay constant of the first excited pion from lattice QCD

abstract: We review the theory that predicts that the decay constant of the excited light pseudo-scalar mesons are suppressed in the chiral limit relative to the pion decay constant. We compute the decay constant of the first excited pion ($\pi^\prime$) using unquenched QCD at a fixed lattice spacing with sea quarks of mass as low as a third the mass of the strange quark. The final result is very sensitive to the improvement conditions. We obtain $f_{\pi^\prime} / f_{\pi}$ = 0.078(93) in the chiral limit.


title: The 4-D Layer Phase as a Gauge Field Localization: Extensive Study of the 5-D Anisotropic U(1) Gauge Model on the Lattice

abstract: We study a 4+1 dimensional pure Abelian Gauge model on the lattice with two anisotropic couplings independent of each other and of the coordinates. A first exploration of the phase diagram using mean field approximation and monte carlo techniques has demonstrated the existence of a new phase, the so called Layer phase, in which the forces in the 4-D subspace are Coulomb-like while in the transverse direction (fifth dimension) the force is confining. This allows the possibility of a gauge field localization scheme. In this work the use of bigger lattice volumes and higher statistics confirms the existence of the Layer phase and furthermore clarifies the issue of the phase transitions' order. We show that the Layer phase is separated from the strongly coupled phase by a weak first order phase transition. Also we provide evidence that the Layer phase is separated by the five-dimensional Coulomb phase with a second order phase transition and we give a first estimation of the critical exponents.


title: Ergodic SO(3), monopole condensation and vortex free energy

abstract: We study the continuum limit of adjoint SU(2) LGT by means of a suppression term for Z2 monopoles. High barriers for tunnelling among different twist sectors are overcome through parallel tempering. Monopole condensation is used to study the deconfinement transition and the properties of the confined phase. Ergodicity in summing over all twist sectors allows an unbiased measure of the 't Hooft vortices free energy. Its behaviour in the SO(3) confined phase hints at differences from what conjectured for semi-integer discretizations.


title: Charm as a domain wall fermion in quenched lattice QCD

abstract: We report a study describing the charm quark by a domain-wall fermion (DWF) in lattice quantum chromodynamics (QCD). Our study uses a quenched gauge ensemble with the DBW2 rectangle-improved gauge action at a lattice cutoff of $a^{-1} \sim 3$ GeV. We calculate masses of heavy-light (charmed) and heavy-heavy (charmonium) mesons with spin-parity $J^P = 0^\mp$ and $1^\mp$, leptonic decay constants of the charmed pseudoscalar mesons ($D$ and $D_s$), and the $D^0$-$\bar{D^0}$ mixing parameter. The charm quark mass is found to be $m^{\bar{\rm MS}}_{c}(m_{c})=1.24(1)(18)$ GeV. The mass splittings in charmed-meson parity partners $\Delta_{q,J=0}$ and $\Delta_{q, J=1}$ are degenerate within statistical errors, in accord with experiment, and they satisfy a relation $\Delta_{q=ud, J} > \Delta_{q=s, J}$, also consistent with experiment. A C-odd axial vector charmonium state, $h_c), lies 22(11) MeV above the $\chi_{c1}$ meson, or $m_{h_{c}} = 3533(11)_{\rm stat.}$ MeV using the experimental $\chi_{c1}) mass. However, in this regard, we emphasize significant discrepancies in the calculation of hyperfine splittings on the lattice. The leptonic decay constants of $D$ and $D_s$ mesons are found to be $f_D=232(7)_{\rm stat.}(^{+6}_{-0})_{\rm chiral}(11)_{\rm syst.}$ MeV and $f_{D_s}/f_{D} = 1.05(2)_{\rm stat.}(^{+0}_{-2})_{\rm chiral}(2)_{\rm syst.}$, where the first error is statistical, the second a systematic due to chiral extrapolation and the third error combination of other known systematics. The $D^0$-$\bar{D^0}$ mixing bag parameter, which enters the $\Delta C = 2$ transition amplitude, is found to be $B_D(2{GeV})=0.845(24)_{\rm stat.}(^{+24}_{-6})_{\rm chiral}(105)_{\rm syst.}$.


title: Pi-K Scattering in Full QCD with Domain-Wall Valence Quarks

abstract: We calculate the pi+ K+ scattering length in fully-dynamical lattice QCD with domain-wall valence quarks on MILC lattices with rooted staggered sea-quarks at a lattice spacing of b=0.125 fm, lattice spatial size of L =2.5 fm and at pion masses of m_pi=290, 350, 490 and 600 MeV. The lattice data, analyzed at next-to-leading order in chiral perturbation theory, allows an extraction of the full pi K scattering amplitude at threshold. Extrapolating to the physical point gives m_pi a_3/2 = -0.0574 (+- 0.0016)(+0.0024 -0.0058) and m_pi a_1/2 = 0.1725 (+- 0.0017)(+0.0023 -0.0156) for the I=3/2 and I=1/2 scattering lengths, respectively, where the first error is statistical and the second error is an estimate of the systematic due to truncation of the chiral expansion.


title: Non-perturbative study of QCD correlators

abstract: This PhD dissertation is devoted to a non-perturbative study of QCD correlators. The main tool that we use is lattice QCD. We concentrated our efforts on the study of the main correlators of the pure Yang - Mills theory in the Landau gauge, namely the ghost and the gluon propagators. We are particularly interested in determining the $\Lqcd$ parameter. It is extracted by means of perturbative predictions available up to NNNLO. The related topic is the influence of non-perturbative effects that show up as appearance of power-corrections to the low-momentum behaviour of the Green functions. A new method of removing these power corrections allows a better estimate of $\Lqcd$. Our result is $\Lambda^{n_f=0}_{\ms} = 269(5)^{+12}_{-9}$ MeV. Another question that we address is the infrared behaviour of Green functions, at momenta of order and below $\Lqcd$. At low energy the momentum dependence of the propagators changes considerably, and this is probably related to confinement. The lattice approach allows to check the predictions of analytical methods because it gives access to non-perturbative correlators. According to our analysis the gluon propagator is finite and non-zero at vanishing momentum, and the power-law behaviour of the ghost propagator is the same as in the free case.


title: Algorithms for Dynamical Fermions

abstract: This is the write-up of three lectures on algorithms for dynamical fermions that were given at the ILFTN workshop 'Perspectives in Lattice QCD' in Nara during November 2005. The first lecture is on the fundamentals of Markov Chain Monte Carlo methods and introduces the Hybrid Monte Carlo (HMC) algorithm and symplectic integrators; the second lecture covers topics in approximation theory and thereby introduces the Rational Hybrid Monte Carlo (RHMC) algorithm and ways of evading integrator instabilities by means of multiple pseudofermion fields; the third lecture introduces on-shell chiral (Ginsparg-Wilson) lattice fermions and discusses five-dimensional formulations for computing fermion propagators for such fermions.


title: Chiral corrections to the Roper mass

abstract: We analyze the quark mass dependence of the Roper mass to one-loop order in relativistic baryon chiral perturbation theory. The loop integrals are evaluated using infrared regularization which preserves chiral symmetry and establishes a chiral counting scheme. The derived chiral expansion of the Roper mass may prove useful for chiral extrapolations of lattice data. For couplings of natural size the quark mass dependence of the Roper mass is similar to the one of the nucleon.


title: Convergence of the chiral expansion for the nucleon mass

abstract: A number of papers recently have used fourth-order chiral perturbation theory to extrapolate lattice data for the nucleon mass; the process seems surprisingly successful even for large pion masses. This paper shows that the inclusion of the fifth-order terms spoils the agreement.


title: Lattice Simulations of the Thermodynamics of Strongly Interacting Elementary Particles and the Exploration of New Phases of Matter in Relativistic Heavy Ion Collisions

abstract: At high temperatures or densities matter formed by strongly interacting elementary particles (hadronic matter) is expected to undergo a transition to a new form of matter - the quark gluon plasma - in which elementary particles (quarks and gluons) are no longer confined inside hadrons but are free to propagate in a thermal medium much larger in extent than the typical size of a hadron. The transition to this new form of matter as well as properties of the plasma phase are studied in large scale numerical calculations based on the theory of strong interactions - Quantum Chromo Dynamics (QCD). Experimentally properties of hot and dense elementary particle matter are studied in relativistic heavy ion collisions such as those currently performed at the relativistic heavy ion collider (RHIC) at BNL. We review here recent results from studies of thermodynamic properties of strongly interacting elementary particle matter performed on Teraflops-Computer. We present results on the QCD equation of state and discuss the status of studies of the phase diagram at non-vanishing baryon number density.


title: Charmonia from Lattice QCD

abstract: Recent lattice QCD results on charmonium properties are reviewed. I comment on molecules and hybrid states as well as on future studies of states near strong decay thresholds.


title: Non-perturbatively Determined Relativistic Heavy Quark Action

abstract: We present a method to non-perturbatively determine the parameters of the on-shell, $O(a)$-improved relativistic heavy quark action. These three parameters, $m_0$, $\zeta$, and $c_B=c_E$ are obtained by matching finite-volume, heavy-heavy and heavy-light meson masses to the exact relativistic spectrum through a finite-volume, step-scaling recursion procedure. We demonstrate that accuracy on the level of a few percent can be achieved by carrying out this matching on a pair of lattices with equal physical spatial volumes but quite different lattice spacings. A fine lattice with inverse lattice spacing $1/a=5.4$ GeV and $24^3 \times 48$ sites and a coarse, $1/a=3.6$ GeV, $16^3 \times 32$ lattice are used together with a heavy quark mass $m$ approximately that of the charm quark. This approach is unable to determine the initially expected, four heavy-quark parameters: $m_0$, $\zeta$, $c_B$ and $c_E$. This apparent non-uniqueness of these four parameters motivated the analytic result, presented in a companion paper, that this set is redundant and that the restriction $c_E=c_B$ is permitted through order $a|\vec p|$ and to all orders in $m a$ where $\vec p$ is the heavy quark's spatial momenta.


title: Relativistic Heavy Quark Effective Action

abstract: We study the fermion action needed to accurately describe the low energy physics of systems including heavy quarks in lattice QCD even when the heavy fermion mass $m$ is on the order of, or larger than, the inverse lattice spacing: $m \ge 1/a$. We carry out an expansion through first order in $|\vec p| a$ (where $\vec p$ is the heavy quark momentum) and all orders in $ma$, refining the analysis of the Fermilab and Tsukuba groups. We demonstrate that the spectrum of heavy quark bound states can be determined accurately through $|\vec p| a$ and $(ma)^n$ for arbitrary exponent $n$ by using a lattice action containing only three unknown coefficients: $m_0$, $\zeta$ and $c_P$ (a generalization of $c_{SW}$), which are functions of $ma$. In a companion paper, we show how these three coefficients can be precisely determined using non-perturbative techniques.


title: Anomalous superfluidity in 2+1 dimensional two-color lattice QCD

abstract: We study thermodynamics of strongly coupled lattice QCD with $two$ colors of staggered fermions in $(2+1)$ dimensions. The partition function of this model can be written elegantly as a statistical mechanics of dimers and baryonloops. The model is invariant under an $SO(3)\times U(1)$ symmetry. At low temperatures we find evidence for superfluidity in the U(1) symmetry sector while the SO(3) symmetry remains unbroken. The finite temperature phase transition appears to belong to the Kosterlitz-Thouless universality class, but the superfluid density jump $\rho_s(T_c)$ at the critical temperature $T_c$ is anomalously higher than the normal value of $2 T_c/\pi$. We show that by adding a small SO(3) symmetry breaking term to the model, the superfluid density jump returns to its normal value implying that the extra symmetry causes anomalous superfluid behavior. Our results may be of interest to researchers studying superfluidity in spin-1 systems.


title: Fixed twist dynamics of SO(3) gauge theory

abstract: We perform a throughout study of 3+1 dim. SO(3) LGT for any fixed-twist background. We concentrate in particular on the physically significant trivial and 1-twist sectors. Introducing a Z(2) monopole chemical potential the 1st order bulk transition is moved down in the strong coupling region and weakened to 2nd order in the 4-dim Ising model universality class. In this extended phase diagram we gain access to a confined phase in every fixed twist sector of the theory. The Pisa disorder operator is employed together with the Polyakov loop to study the confinement-deconfinement transition in each sector. Due to the specific properties of both operators, most results can be used to gain insight in the ergodic theory, where all twist sectors should be summed upon. An explicit mapping of each fixed twist theory to effective positive plaquette models with fixed twisted boundary conditions is applied to better establish their properties in the different phases.


title: Charmonium spectral functions in two-flavour QCD

abstract: We compute charmonium spectral functions in 2-flavour QCD using the maximum entropy method and anisotropic lattices. We find that the S-waves (J/psi and eta_c) survive up to temperatures close to 2T_c, while the P-waves (chi_c0 and chi_c1) melt away below 1.3T_c.


title: Topics in Effective Field Theory for Lattice QCD

abstract: In this work, we extend and apply effective field theory techniques to systematically understand a subset of lattice artifacts which pollute the lattice correlation functions for a few processes of physical interest. Where possible, we compare to existing lattice QCD calculations. In particular, we extend the heavy baryon Lagrangian to the next order in partially quenched chiral perturbation theory and use it to compute the masses of the lightest spin-1/2 and spin-3/2 baryons to next-to-next-to leading order. We then construct the twisted mass chiral Lagrangian for baryons and apply it to compute the lattice spacing corrections to the baryon masses simulated with twisted mass lattice QCD. We extend computations of the nucleon electromagnetic structure to account for finite volume effects, as these observables are particularly sensitive to the finite extent of the lattice. We resolve subtle peculiarities for lattice QCD simulations of polarizabilities and we show that using background field techniques, one can make predictions for the 4 spin-dependent nucleon polarizabilities, quantities which are difficult to access experimentally. We then discuss the two-pion system in finite volume, determining the exponentially small volume corrections necessary for lattice determinations of the scattering parameters. We also determine the lattice spacing artifacts that arise for a mixed-action lattice simulation of the two-pion system with Ginsparg-Wilson valence quarks and staggered sea quarks. We show that the isospin 2 scattering length has a near continuum like behavior, differing from the chiral perturbation theory calculation by a computable difference.


title: Calculating the hadronic vacuum polarization and leading hadronic contribution to the muon anomalous magnetic moment with improved staggered quarks

abstract: We present a lattice calculation of the hadronic vacuum polarization and the lowest-order hadronic contribution to the muon anomalous magnetic moment, a_\mu = (g-2)/2, using 2+1 flavors of improved staggered fermions. A precise fit to the low-q^2 region of the vacuum polarization is necessary to accurately extract the muon g-2. To obtain this fit, we use staggered chiral perturbation theory, including the vector particles as resonances, and compare these to polynomial fits to the lattice data. We discuss the fit results and associated systematic uncertainties, paying particular attention to the relative contributions of the pions and vector mesons. Using a single lattice spacing ensemble (a=0.086 fm), light quark masses as small as roughly one-tenth the strange quark mass, and volumes as large as (3.4 fm)^3, we find a_\mu^{HLO} = (713 \pm 15) \times 10^{-10} and (748 \pm 21) \times 10^{-10} where the error is statistical only and the two values correspond to linear and quadratic extrapolations in the light quark mass, respectively. Considering systematic uncertainties not eliminated in this study, we view this as agreement with the current best calculations using the experimental cross section for e^+e^- annihilation to hadrons, 692.4 (5.9) (2.4)\times 10^{-10}, and including the experimental decay rate of the tau lepton to hadrons, 711.0 (5.0) (0.8)(2.8)\times 10^{-10}. We discuss several ways to improve the current lattice calculation.


title: The spectrum of lattice QCD with staggered fermions at strong coupling

abstract: Using 4 flavors of staggered fermions at infinite gauge coupling, we compare various analytic results for the hadron spectrum with exact Monte Carlo simulations. Agreement with Ref. \cite{Martin_etal} is very good, at the level of a few percent. Our results give credence to a discrepancy between the baryon mass and the critical chemical potential, for which baryons fill the lattice at zero temperature and infinite gauge coupling. Independent determinations of the latter set it at about 30% less than the baryon mass. One possible explanation is that the nuclear attraction becomes strong at infinite gauge coupling.


title: The transition temperature in QCD

abstract: We present a detailed calculation of the transition temperature in QCD with two light and one heavier (strange) quark mass on lattices with temporal extent N_t =4 and 6. Calculations with improved staggered fermions have been performed for various light to strange quark mass ratios in the range, 0.05 <= m_l/m_s <= 0.5, and with a strange quark mass fixed close to its physical value. From a combined extrapolation to the chiral (m_l -> 0) and continuum (aT = 1/N_t -> 0) limits we find for the transition temperature at the physical point T_c r_0 = 0.457(7) where the scale is set by the Sommer-scale parameter r_0 defined as the distance in the static quark potential at which the slope takes on the value, (dV_qq(r)/dr)_r=r_0 = 1.65/r_0^2. Using the currently best known value for r_0 this translates to a transition temperature T_c = 192(7)(4)MeV. The transition temperature in the chiral limit is about 3% smaller. We discuss current ambiguities in the determination of T_c in physical units and also comment on the universal scaling behavior of thermodynamic quantities in the chiral limit.