title: Lattice QCD at finite density: imaginary chemical potential

abstract: I describe the results for the critical line and the thermodynamics of different phases of QCD which have been obtained by lattice simulations with an imaginary chemical potential, and review motivations and merits of the different strategies -- Taylor expansion, Fourier analysis and Pade' approximants -- for analytic continuation from imaginary to real chemical potential. I consider phenomenological models which can be easily extended to the complex chemical potential plane, thus affording a direct comparison with lattice data at imaginary $\mu$: the hadronic phase and the high temperature limit are amenable to a simple description, while a rather subtle interplay between thermodynamics and critical behaviour emerges in the hot phase close to $T_c$


title: Analytic continuation from imaginary to real chemical potential in two-color QCD

abstract: The method of analytic continuation from imaginary to real chemical potential is one of the most powerful tools to circumvent the sign problem in lattice QCD. Here we test this method in a theory, 2-color QCD, which is free from the sign problem. We find that the method gives reliable results, within appropriate ranges of the chemical potential, and that a considerable improvement can be achieved if suitable functions are used to interpolate data with imaginary chemical potential.


title: Bound States of (Anti-)Scalar-Quarks in SU(3)_c Lattice QCD

abstract: Light scalar-quarks \phi (colored scalar particles or idealized diquarks) and their color-singlet hadronic states are studied with quenched SU(3)_c lattice QCD in terms of mass generation. We investigate ``scalar-quark mesons'' \phi^\dagger \phi and ``scalar-quark baryons'' \phi\phi\phi as the bound states of scalar-quarks \phi. We also investigate the bound states of scalar-quarks \phi and quarks \psi, i.e., \phi^\dagger \psi, \psi\psi\phi and \phi\phi\psi, which we name ``chimera hadrons''. All the new-type hadrons including \phi are found to have a large mass due to large quantum corrections by gluons, even for zero bare scalar-quark mass m_\phi=0 at a^{-1}\sim 1{\rm GeV}. We conjecture that all colored particles generally acquire a large effective mass due to dressed gluon effects.


title: Complete spectra of the Dirac operator and their relation to confinement

abstract: We compute complete spectra of the staggered lattice Dirac operator for quenched SU(3) gauge configurations below and above the critical temperature. The confined and the deconfined phase are characterized by a different response of the Dirac eigenvalues to a change of the fermionic boundary conditions. We analyze the role of the eigenvalues in recently developed spectral sums representing the Polyakov loop. We show that the Polyakov loop gets its main contributions from the UV end of the spectrum.


title: Precision study of the SU(3) topological susceptibility in the continuum

abstract: We determine the topological susceptibility in the SU(3) pure gauge theory. We perform a series of high-statistics lattice studies and take the combined continuum and infinite volume limit. We find chi_{top}r_0^4=0.0524(7)(6) which translates into chi_{top}^{1/4}=193(1)(8)MeV with the second error exclusively due to the intrinsic scale ambiguity.


title: Perturbative renormalization for static and domain-wall bilinears and four-fermion operators with improved gauge actions

abstract: We calculate one-loop renormalization factors for heavy-light bilinears as well as four-fermion operators relevant for $B^{0} - \bar{B}^{0}$ mixing calculations on the lattice. We use the static approximation for heavy quarks and the domain-wall formulation for light quarks. We present results for different choices of improved gauge action.


title: Hot QCD and warm dark matter

abstract: One of the possible explanations for the dark matter needed in the standard cosmological model is so-called warm dark matter, in the form of right-handed ("sterile") neutrinos with a mass in the keV range. I describe how various properties of QCD at temperatures of a few hundred MeV play an important role in the theoretical computations that are needed for consolidating or falsifying this scenario. In particular the points where lattice QCD could help are underlined.


title: Quantitative comparison of filtering methods in lattice QCD

abstract: We systematically compare filtering methods used to extract topological excitations (like instantons, calorons, monopoles and vortices) from lattice gauge configurations, namely APE-smearing and spectral decompositions based on lattice Dirac and Laplace operators. Each of these techniques introduces ambiguities, which can invalidate the interpretation of the results. We show, however, that all these methods, when handled with care, reveal very similar topological structures. Hence, these common structures are free of ambiguities and faithfully represent infrared degrees of freedom in the QCD vacuum. As an application we discuss an interesting power-law for the clusters of filtered topological charge.


title: The spectrum of charmonium-like vector mesons in lattice QCD

abstract: We present the first lattice results of the spectrum of exotic vector mesons extracted from the molecular and diquark-antidiquark operators, with quark fields c-q-cbar-qbar, and c-s-cbar-qbar (c-q-cbar-sbar) respectively, in lattice QCD with exact chiral symmetry. Our results suggest that X(3872) and Y(4260) are in the spectrum of QCD, with J^{PC} = 1^{++} and 1^{--} respectively. Moreover, we obtain the spectrum of heavier exotic mesons with c-s-cbar-ubar (c-u-cbar-sbar), c-s-cbar-dbar (c-d-cbar-sbar), c-s-cbar-sbar, and c-c-cbar-cbar, as the first theoretical predictions from lattice QCD.


title: Hyperon-Nucleon Scattering from Fully-Dynamical Lattice QCD

abstract: We present results of the first fully-dynamical lattice QCD determination of hyperon-nucleon scattering. One s-wave phase shift was determined for n Lambda scattering in both spin-channels at pion masses of 350, 490, and 590 MeV, and for n Sigma- scattering in both spin channels at pion masses of 490, and 590 MeV. The calculations were performed with domain-wall valence quarks on dynamical, staggered gauge configurations with a lattice spacing of b ~ 0.125 fm.


title: Generalised parton distributions of the pion in partially-quenched chiral perturbation theory

abstract: We consider the pion matrix elements of the isoscalar and isovector combinations of the vector and tensor twist-two operators that determine the moments of the various pion generalised parton distributions. Our analysis is performed using partially-quenched chiral perturbation theory. We work in the SU(2) and SU(4|2) theories and present our results at infinite volume and also at finite volume where some subtleties arise. These results are useful for extrapolations of lattice calculations of these matrix elements at small momentum transfer to the physical regime.


title: Matrix elements of the complete set of \Delta B = 2 and \Delta C = 2 operators in heavy meson chiral perturbation theory

abstract: Using heavy meson chiral perturbation theory, we consider the light quark-mass and spatial volume dependence of the matrix elements of \Delta B=2 and \Delta C=2 four-quark operators relevant for B^{0}_{(s)}{-}\bar{B}^{0}_{(s)} and D^{0}{-}\bar{D}^{0} mixing, and the B_{s} meson width difference. Our results for these matrix elements are obtained in the N_{f}=2+1 partially quenched theory, which becomes full QCD in the limit where sea and valence quark masses become equal. They can be used in extrapolation of lattice calculations of these matrix elements to the physical light quark masses and to infinite volume. An important conclusion of this paper is that the chiral extrapolations for matrix elements of heavy-light meson mixing beyond the Standard Model, and those relevant for the B_s width difference are more complicated than that for the Standard Model mixing matrix elements.


title: Precision study of B^* B\pi coupling for the static heavy-light meson

abstract: We compute the B^*B\pi coupling \hat{g}_{\infty} for static heavy-light meson using all-to-all propagators. It is shown that low-mode averaging with 100 low-lying eigenmodes indeed improves the signal for the 2-point and 3-point functions for heavy-light meson significantly. Our study suggests that the all-to-all propagator will be a very efficient method for high precision computation of the B^*B\pi coupling especially in unquenched QCD where the number of configurations is limited.


title: Flavor symmetry breaking and scaling for improved staggered actions in quenched QCD

abstract: We present a study of the flavor symmetry breaking in the pion spectrum for various improved staggered fermion actions. To study the effects of link fattening and tadpole improvement, we use three different variants of the p4 action - p4fat3, p4fat7, and p4fat7tad. These are compared to Asqtad and also to naive staggered. To study the pattern of symmetry breaking, we measure all 15 meson masses in the 4-flavor staggered theory. The measurements are done on a quenched gauge background, generated using a one-loop improved Symanzik action with $\beta=10/g^2 = 7.40, 7.75,$ and 8.00, corresponding to lattice spacings of approximately a = .31 fm., .21 fm., and .14 fm. We also study how the lattice scale set by the $\rho$ mass on each of these ensembles compares to one set by the static quark potential.


title: Three dimensional four-fermion models - A Monte Carlo study

abstract: We present results from numerical simulations of three different 3d four-fermion models that exhibit Z_2, U(1), and SU(2) x SU(2) chiral symmetries, respectively. We performed the simulations by using the hybrid Monte Carlo algorithm. We employed finite size scaling methods on lattices ranging from 8^3 to 40^3 to study the properties of the second order chiral phase transition in each model. The corresponding critical coupling defines an ultraviolet fixed point of the renormalization group. In our high precision simulations, we detected next-to-leading order corrections for various critical exponents and we found them to be in good agreement with existing analytical large-N_f calculations.


title: Transverse spin structure of the nucleon from lattice QCD simulations

abstract: We present the first calculation in lattice QCD of the lowest two moments of transverse spin densities of quarks in the nucleon. They encode correlations between quark spin and orbital angular momentum. Our dynamical simulations are based on two flavors of clover-improved Wilson fermions and Wilson gluons. We find significant contributions from certain quark helicity flip generalized parton distributions, leading to strongly distorted densities of transversely polarized quarks in the nucleon. In particular, based on our results and recent arguments by Burkardt [Phys. Rev. D 72 (2005) 094020], we predict that the Boer-Mulders-function $h_1^\perp$, describing correlations of transverse quark spin and intrinsic transverse momentum of quarks, is large and negative for both up and down quarks.


title: Couplings of hybrid operators to ground and excited states of bottomonia

abstract: We analyze the overlap of local color-octet meson operators with the $\Upsilon$ and the $\eta_b$ and their low-lying excited states, especially the first radial excitations. Our analysis is based on NRQCD and includes all terms up to order $v^4$. We use a variety of source and sink operators as a basis for the variational method, which enables us to clearly separate the mass eigenstates and hence to extract the desired amplitudes. The results show the usefulness of the variational method for determining couplings to excited hadronic states.


title: QCD phase diagram: an overview

abstract: The aim of this review is to summarize the contemporary understanding of the QCD phase diagram as a function of temperature and baryo-chemical potential. The focus is on recent theoretical developments due to lattice simulations of the phase diagram.


title: Understanding nucleon structure using lattice simulations

abstract: This review focuses on the discussion of three key results of nucleon structure calculations on the lattice. These three results are the quark contribution to the nucleon spin, J_q, the nucleon-Delta transition form factors, and the nucleon axial coupling, g_A. The importance for phenomenology and experiment is discussed and the requirements for future simulations are pointed out.


title: Chiral Lattice Gauge Theories and The Strong Coupling Dynamics of a Yukawa-Higgs Model with Ginsparg-Wilson Fermions

abstract: The Yukawa-Higgs/Ginsparg-Wilson-fermion construction of chiral lattice gauge theories described in hep-lat/0605003 uses exact lattice chirality to decouple the massless chiral fermions from a mirror sector, whose strong dynamics is conjectured to give cutoff-scale mass to the mirror fermions without breaking the chiral gauge symmetry. In this paper, we study the mirror sector dynamics of a two-dimensional chiral gauge theory in the limitof strong Yukawa and vanishing gauge couplings, in which case it reduces to an XY model coupled to Ginsparg-Wilson fermions. For the mirror fermions to acquire cutoff-scale mass it is believed to be important that the XY model remain in its "high temperature" phase, where there is no algebraic ordering--a conjecture supported by the results of our work. We use analytic and Monte-Carlo methods with dynamical fermions to study the scalar and fermion susceptibilities, and the mirror fermion spectrum. Our results provide convincing evidence that the strong dynamics does not "break" the chiral symmetry (more precisely, that the mirror fermions do not induce algebraic ordering in two-dimensions), and that the mirror fermions decouple from the infrared physics.


title: Constant contribution in meson correlators at finite temperature

abstract: We discuss a constant contribution to meson correlators at finite temperature. In the deconfinement phase of QCD, a colored single quark state is allowed as a finite energy state, which yields to a contribution of wraparound quark propagation to temporal meson correlators. We investigate the effects in the free quark case and quenched QCD at finite temperature. The ``scattering'' contribution causes a constant mode in meson correlators with zero spatial momentum and degenerate quark masses, which can dominate the correlators in the region of large imaginary times. In the free spectral function, the contribution yields a term proportional to $\omega\delta(\omega)$. Therefore this contribution is related to transport phenomena in the quark gluon plasma. It is possible to distinguish the constant contribution from the other part using several analysis methods proposed in this paper. As a result of the analyses, we find that drastic changes in charmonium correlators for $\chi_c$ states just above the deconfinement transition are due to the constant contribution. The other differences in the $\chi_c$ states are small. It may indicate the survival of $\chi_c$ states after the deconfinement transition until, at least, $1.4T_c$.


title: Advances and applications of lattice supersymmetry

abstract: I review motivations for the study of supersymmetric field theories by lattice techniques. In particular, some of the more interesting potential applications are described. These are models of quantum gravity, that rely on the AdS/CFT correspondence. The dual gauge theory is an ideal place for lattice studies to make a contribution. I also survey some of the recent advances in lattice formulations of supersymmetric theories.


title: Deconfining Phase Transition on Lattices with Boundaries at Low Temperature

abstract: In lattice gauge theory (LGT) equilibrium simulations of QCD are usually performed with periodic boundary conditions (BCs). In contrast to that deconfined regions created in heavy ion collisions are bordered by the confined phase. Here we discuss BCs in LGT, which model a cold exterior of the lattice volume. Subsequently we perform Monte Carlo (MC) simulations of pure SU(3) LGT with a thus inspired simple change of BCs using volumes of a size comparable to those typically encountered in the BNL relativistic heavy ion collider (RHIC) experiment. Corrections to the usual LGT results survive in the finite volume continuum limit and we estimate them as function of the volume size. In magnitude they are found comparable to those of including quarks. As observables we use a pseudocritical temperature, which rises opposite to the effect of quarks, and the width of the transition, which broadens similar to the effect of quarks.


title: Lattice measurement of $B_{B_s}$ with a chiral light quark action

abstract: The computation on the lattice of the bag parameter $B_{B_s}$ associated to the $B_s - \bar{B_s}$ mixing amplitude in the Standard Model is presented. The estimation has been made by combining the static limit of HQET and the Neuberger light quark action which preserves the chiral symmetry on the lattice. We find $B^{\bar{MS}}{stat}}_{B_s}(m_b)=0.92(3)$.


title: QCD with light Wilson quarks on fine lattices (II): DD-HMC simulations and data analysis

abstract: In this second report on our recent numerical simulations of two-flavour QCD, we provide further technical details on the simulations and describe the methods we used to extract the meson masses and decay constants from the generated ensembles of gauge fields. Among the topics covered are the choice of the DD-HMC parameters, the issue of stability, autocorrelations and the statistical error analysis. Extensive data tables are included as well as a short discussion of the quark-mass dependence in partially quenched QCD, supplementing the physics analysis that was presented in the first paper in this series.


title: Equilibrium Thermodynamics of Lattice QCD

abstract: Lattice QCD allows us to simulate QCD at non-zero temperature and/or densities. Such equilibrium thermodynamics calculations are relevant to the physics of relativistic heavy-ion collisions. I give a brief review of the field with emphasis on our work.


title: Infrared-suppressed gluon propagator in 4d Yang-Mills theory in a Landau-like gauge

abstract: The infrared behavior of the gluon propagator is directly related to confinement in QCD. Indeed, the Gribov-Zwanziger scenario of confinement predicts an infrared vanishing (transverse) gluon propagator in Landau-like gauges, implying violation of reflection positivity and gluon confinement. Finite-volume effects make it very difficult to observe (in the minimal Landau gauge) an infrared suppressed gluon propagator in lattice simulations of the four-dimensional case. Here we report results for the SU(2) gluon propagator in a gauge that interpolates between the minimal Landau gauge (for gauge parameter lambda equal to 1) and the minimal Coulomb gauge (corresponding to lambda = 0). For small values of lambda we find that the spatially-transverse gluon propagator D^tr(0,|\vec p|), considered as a function of the spatial momenta |\vec p|, is clearly infrared suppressed. This result is in agreement with the Gribov-Zwanziger scenario and with previous numerical results in the minimal Coulomb gauge. We also discuss the nature of the limit lambda -> 0 (complete Coulomb gauge) and its relation to the standard Coulomb gauge (lambda = 0). Our findings are corroborated by similar results in the three-dimensional case, where the infrared suppression is observed for all considered values of lambda.


title: Dynamical Twisted Mass Fermions with Light Quarks

abstract: We present results of dynamical simulations with 2 flavours of degenerate Wilson twisted mass quarks at maximal twist in the range of pseudo scalar masses from 300 to 550 MeV. The simulations are performed at one value of the lattice spacing a \lesssim 0.1 fm. In order to have O(a) improvement and aiming at small residual cutoff effects, the theory is tuned to maximal twist by requiring the vanishing of the untwisted quark mass. Precise results for the pseudo scalar decay constant and the pseudo scalar mass are confronted with chiral perturbation theory predictions and the low energy constants F, \bar{l}_3 and \bar{l}_4 are evaluated with small statistical errors.


title: 2+1 flavor domain wall QCD on a (2 fm)^3 lattice: light meson spectroscopy with Ls = 16

abstract: We present results for light meson masses and pseudoscalar decay constants from the first of a series of lattice calculations with 2+1 dynamical flavors of domain wall fermions and the Iwasaki gauge action. The work reported here was done at a fixed lattice spacing of about 0.12 fm on a 16^3\times32 lattice, which amounts to a spatial volume of (2 fm)^3 in physical units. The number of sites in the fifth dimension is 16, which gives m_{res} = 0.00308(4) in these simulations. Three values of input light sea quark masses, m_l^{sea} \approx 0.85 m_s, 0.59 m_s and 0.33 m_s were used to allow for extrapolations to the physical light quark limit, whilst the heavier sea quark mass was fixed to approximately the physical strange quark mass m_s. The exact rational hybrid Monte Carlo algorithm was used to evaluate the fractional powers of the fermion determinants in the ensemble generation. We have found that f_\pi = 127(4) MeV, f_K = 157(5) MeV and f_K/f_\pi = 1.24(2), where the errors are statistical only, which are in good agreement with the experimental values.


title: Renormalisation group evolution for the $\Delta S = 1$ effective Hamiltonian with $N_f=2+1$

abstract: We discuss the renormalisation group (RG) evolution for the $\Delta S = 1$ operators in unquenched QCD with $N_f = 3$ ($m_u=m_d=m_s$) or, more generally, $N_f = 2+1$ ($m_u=m_d \ne m_s$) flavors. In particular, we focus on the specific problem of how to treat the singularities which show up only for $N_f=3$ or $N_f = 2+1$ in the original solution of Buras {\it et al.} for the RG evolution matrix at next-to-leading order. On top of Buras {\it et al.}'s original treatment, we use a new method of analytic continuation to obtain the correct solution in this case. It is free of singularities and can therefore be used in numerical analysis of data sets calculated in lattice QCD.


title: Decay constants of charm and beauty pseudoscalar heavy-light mesons on fine lattices

abstract: We compute decay constants of heavy-light mesons in quenched lattice QCD with a lattice spacing of a ~ 0.04 fm using non-perturbatively O(a) improved Wilson fermions and O(a) improved currents. We obtain f_{D_s} = 220(6)(5)(11) MeV, f_D = 206(6)(3)(22) MeV, f_{B_s} = 205(7)(26)(17) MeV and f_B = 190(8)(23)(25) MeV, using the Sommer parameter r_0 = 0.5 fm to set the scale. The first error is statistical, the second systematic and the third from assuming a +-10% uncertainty in the experimental value of r_0. A detailed discussion is given in the text. We also present results for the meson decay constants f_K and f_\pi and the \rho meson mass.


title: Critical Flavor Number in the Three Dimensional Thirring Model

abstract: We present results of a Monte Carlo simulation of the three dimensional Thirring model with the number of fermion flavors N_f varied between 2 and 18. By identifying the lattice coupling at which the chiral condensate peaks, simulations are be performed at couplings g^2(N_f) corresponding to the strong coupling limit of the continuum theory. The chiral symmetry restoring phase transition is studied as N_f is increased, and the critical number of flavors estimated as N_{fc}=6.6(1). The critical exponents measured at the transition do not agree with self-consistent solutions of the Schwinger-Dyson equations; in particular there is no evidence for the transition being of infinite order. Implications for the critical flavor number in QED_3 are briefly discussed.


title: Heavy quark thermodynamics in full QCD

abstract: We analyze the large-distance behaviour of static quark-anti-quark pair correlations in QCD. The singlet free energy is calculated and the entropy contribution to it is identified allowing us to calculate the excess internal energy. The free energy has a sharp drop in the critical region, leading to sharp peaks in both excess entropy and internal energy.


title: Chiral anomalies and rooted staggered fermions

abstract: A popular approximation in lattice gauge theory is an extrapolation in the number of fermion species away from the four fold degeneracy natural with the staggered fermion formulation. I show that the extrapolation procedure mutilates the expected continuum holomorphic behavior in the quark masses. The conventional resolution proposes canceling the unphysical singularities with a plethora of extra states appearing at finite lattice spacing. This unproven conjecture requires an explicit loss of unitarity and locality. Even if correct, the approach implies large cutoff effects in the low-energy flavor-neutral sector.


title: QCD thermodynamics at zero and non-zero density

abstract: We present recent results on thermodynamics of QCD with almost physical light quark masses and a physical strange quark mass value. These calculations have been performed with an improved staggerd action especially designed for finite temperature lattice QCD. In detail we present a calculation of the transition temperature, using a combined chiral and continuum extrapolation. Furthermore we present preliminary results on the interaction measure and energy density at almost realistic quark masses. Finally we disscuss the response of the pressure to a finite quark chemical potential. Within the Taylor expansion formalism we calculate quark number susceptibilities and leading order corrections to finite chemical potential. This is particularly usefull for mapping out the critical region in the QCD phase diagram.


title: The Stefan-Boltzmann law in a small box and the pressure deficit in hot SU(N) lattice gauge theory

abstract: The blackbody radiation in a box L^3 with periodic boundary conditions in thermal equilibrium at a temperature T is affected by finite-size effects. These bring about modifications of the thermodynamic functions which can be expressed in a closed form in terms of the dimensionless parameter LT. For instance, when LT~4 - corresponding to the value where the most reliable SU(N) gauge lattice simulations have been performed above the deconfining temperature T_c - the deviation of the free energy density from its thermodynamic limit is about 5%. This may account for almost half of the pressure deficit observed in lattice simulations at T~ 4 T_c.

categories: hep-lat astro-ph cond-mat.stat-mech hep-ph hep-th


title: Radially Excited States of 1P Charmonia and X(3872)

abstract: The excited states of charmonia are numerically investigated in quenched lattice QCD with improved gauge and Wilson fermion actions formulated on anisotropic lattices. Through a constrained curve fitting algorithm, the masses of the first excited states in $0^{++}$, $1^{++}$, and $1^{+-}$ channels are determined to be 3.825(88), 3.853(57), and 3.858(70) GeV, respectively. Furthormore, a node structure is also observed in the Bethe-Salpeter amplitude of the $1^{++}$ first excited state. These observations indicate that X(3872) could be the first radial excitation of $\chi_{c1}$.


title: The density of states method at non-zero chemical potential

abstract: We study the QCD phase diagram by first principle lattice calculations at so far unreached high densities. For this purpose we employ the density of states method. Unimproved staggered fermions, which describe four quark flavors in the continuum are used in this analysis. Though the method is quite expensive, small lattices show an indication for a triple-point connecting three different phases on the phase diagram.


title: Diquark correlations in baryons on the lattice with overlap quarks

abstract: We evaluate baryon wave functions in both the Coulomb and Landau gauge in lattice QCD. These are constructed from quark propagators calculated with the overlap Dirac operator on quenched gauge configurations at beta = 6. By comparing baryon states that differ in their diquark content, we find evidence for enhanced correlation in the scalar diquark channel, as favored by quark models. We also summarize earlier results for diquark masses in the Landau gauge, casting them in a form more easily compared with subsequent studies.


title: Phase diagram at finite temperature and quark density in the strong coupling region of lattice QCD for color SU(3)

abstract: We study the phase diagram of quark matter at finite temperature (T) and chemical potential (mu) in the strong coupling region of lattice QCD for color SU(3). Baryon has effects to extend the hadron phase to a larger mu direction relative to Tc at low temperatures in the strong coupling limit. With the 1/g^2 corrections, Tc is found to decrease rapidly as g decreases, and the shape of the phase diagram becomes closer to that expected in the real world.


title: A semi-variational approach to QCD at finite temperature and baryon density

abstract: Recently a new bosonization method has been used to derive, at zero fermion density, an effective action for relativistic field theories whose partition function is dominated by fermionic composites, chiral mesons in the case of QCD. This approach shares two important features with variational methods: the restriction to the subspace of the composites, and the determination of their structure functions by a variational calculation. But unlike standard variational methods it treats excited states at the same time and on the same footing as the ground state. I extend this method including states of nonvanishing fermion (baryon) number and derive an effective action for QCD at finite temperature and baryon density. I test the result on a four-fermion interaction model.


title: Infrared behavior of the Faddeev-Popov operator in Coulomb gauge QCD

abstract: We calculate the eigenvalue distribution of the Faddeev-Popov operator in Coulomb gauge QCD using quenched SU(3) lattice simulation. In the confinement phase, the density of the low-lying eigenvalues increases with lattice volume, and the confinement criterion is satisfied. Moreover, even in the deconfinement phase, the behavior of the FP eigenvalue density is qualitatively the same as in the confinement phase. This is consistent with the fact that the color-Coulomb potential is not screened in the deconfined phase.


title: Two-flavor lattice QCD simulation in the epsilon-regime with exact chiral symmetry

abstract: We perform lattice simulations of two-flavor QCD using Neuberger's overlap fermion, with which the exact chiral symmetry is realized at finite lattice spacings. The epsilon-regime is reached by decreasing the light quark mass down to 3 MeV on a 16^3 32 lattice with a lattice spacing \sim 0.11 fm. We find a good agreement of the low-lying Dirac eigenvalue spectrum with the analytical predictions of the chiral random matrix theory, which reduces to the chiral perturbation theory in the epsilon-regime. The chiral condensate is extracted as \Sigma(2 GeV) = (251(7)(11) MeV)^3, where the errors are statistical and an estimate of the higher order effects in the epsilon-expansion.


title: Heavy-Quark Free Energy, Debye Mass, and Spatial String Tension at Finite Temperature in Two Flavor Lattice QCD with Wilson Quark Action

abstract: We study Polyakov loop correlations and spatial Wilson loop at finite Temperature in two-flavor QCD simulations with the RG-improved gluon action and the clover-improved Wilson quark action on a $ 16^3 \times 4$ lattice. From the line of constant physics at $m_{\rm PS}/m_{\rm V}=0.65$ and 0.80, we extract the heavy-quark free energies, the effective running coupling $g_{\rm eff}(T)$ and the Debye screening mass $m_D(T)$ for various color channels of heavy quark--quark and quark--anti-quark pairs above the critical temperature. The free energies are well approximated by the screened Coulomb form with the appropriate Casimir factors at high temperature. The magnitude and the temperature dependence of the Debye mass are compared to those of the next-to-leading order thermal perturbation theory and to a phenomenological formula in terms of $g_{\rm eff}(T)$. We make a comparison between our results with the Wilson quark action and the previous results with the staggered quark action. The spatial string tension is also studied in the high temperature phase and is compared to the next-to-next-leading order prediction in an effective theory with dimensional reduction.


title: Thermodynamics of two-flavor lattice QCD with an improved Wilson quark action at non-zero temperature and density

abstract: We report the current status of our systematic studies of the QCD thermodynamics by lattice QCD simulations with two flavors of improved Wilson quarks. We evaluate the critical temperature of two flavor QCD in the chiral limit at zero chemical potential and show the preliminary result. Also we discuss fluctuations at none-zero temperature and density by calculating the quark number and isospin susceptibilities and their derivatives with respect to chemical potential.


title: Energies and radial distributions of B_s mesons on the lattice

abstract: This is a follow-up to our earlier work for the energies and the charge (vector) and matter (scalar) distributions for S-wave states in a heavy-light meson, where the heavy quark is static and the light quark has a mass about that of the strange quark. We now study excited states of these mesons with higher angular momentum and with radial nodes. The calculation is carried out with dynamical fermions on a 16^3 times 32 lattice with a lattice spacing approximately 0.10 fm. The lattice configurations were generated by the UKQCD Collaboration. Attempts are now being made to understand these results in terms of the Dirac equation. In nature the closest equivalent of this heavy-light system is the B_s meson, which allows us to compare our lattice calculations to experimental results (where available) or give a prediction where the P-wave states should lie. We pay particular attention to the spin-orbit splitting, to see which one of the states (for a given angular momentum L) has the lower energy.


title: Loop Approach to Lattice Gauge Theories

abstract: We solve the Gauss law and the corresponding Mandelstam constraints in the loop Hilbert space ${\cal H}^{L}$ using the prepotential formulation of $(d+1)$ dimensional SU(2) lattice gauge theory. The resulting orthonormal and complete loop basis, explicitly constructed in terms of the $d(2d-1)$ prepotential intertwining operators, is used to transcribe the gauge dynamics directly in ${\cal H}^{L}$ without any redundant gauge and loop degrees of freedom. Using generalized Wigner-Eckart theorem and Biedenharn -Elliot identity in ${\cal H}^L$, we show that the loop dynamics for pure SU(2) lattice gauge theory in arbitrary dimension, is given by the real symmetric $3nj$ symbols of first kind (e.g., n=6, 10 for d=2, 3 respectively). The corresponding "ribbon diagrams" representing SU(2) loop dynamics are constructed. The prepotential techniques are trivially extended to include fundamental matter fields leading to a description in terms of loops and strings. The SU(N) gauge group is briefly discussed.


title: Lattice QCD with a chiral twist

abstract: In these lectures I explain how chiral symmetry of continuum QCD naturally leads to a class of lattice regularisations known as twisted mass QCD (tmQCD). As compared to standard Wilson quarks, its advantages are the absence of unphysical zero modes, the possibility to circumvent lattice renormalisation problems and automatic O(a) improvement. On the other hand, the physical parity and flavour symmetries are explicitly broken. I discuss these aspects and then turn to the theory in a finite space-time volume with Schroedinger functional boundary conditions. Again, chiral transformations of the continuum theory may be used as a guide to formulate an alternative lattice regularisation of the Schr\"odinger functional, with interesting applications to renormalization problems in QCD.


title: Color Screening and Quark-Quark Interactions in Finite Temperature QCD

abstract: We analyze the screening of static diquark sources in 2-flavor QCD and compare results with the screening of static quark-antiquark pairs. We show that a two quark system in a fixed color representations is screened at short distances like a single quark source in the same color representation whereas at large distances the two quarks are screened independently. At high temperatures we observe that the relative strength of the interaction in diquark and quark-antiquark systems, respectively, obeys Casimir scaling. We use this result to examine the possible existence of heavy quark-quark bound states in the high temperature phase of QCD. We find support for the existence of $bb$ states up to about $2T_c$ while $cc$ states are unlikely to be formed above $T_c$.


title: Radial and orbital excitations of static-light mesons

abstract: We present results for the spectrum of static-light mesons from Nf=2 lattice QCD. These results were obtained using all-to-all light quark propagators on an anisotropic lattice, yielding an improved signal resolution when compared to more conventional lattice techniques. With a light quark mass close to the strange quark, we have measured the splittings between the ground-state S-wave static-light meson and higher excitations. We attempt to identify the quantum numbers of the excited states in the context of the reduced spatial symmetries of the lattice.


title: The QCD Sign Problem for Small Chemical Potential

abstract: The expectation value of the complex phase factor of the fermion determinant is computed in the microscopic domain of QCD at nonzero chemical potential. We find that the average phase factor is non-vanishing below a critical value of the chemical potential equal to the half the pion mass and vanishes exponentially in the volume for larger values of the chemical potential. This holds for QCD with dynamical quarks as well as for quenched and phase quenched QCD. The average phase factor has an essential singularity for zero chemical potential and cannot be obtained by analytic continuation from imaginary chemical potential or by means of a Taylor expansion. The leading order correction in the $p$-expansion of the chiral Lagrangian is calculated as well.


title: The Delta-resonance in a finite volume

abstract: We study the extraction of Delta-resonance parameters from lattice data for small quark masses, corresponding to the case of an unstable Delta. To this end, we calculate the spectrum of the correlator of two Delta-fields in a finite Euclidian box up-to-and-including O(epsilon^3) in the small scale expansion using infrared regularization. On the basis of our numerical study, we argue that the extraction of the parameters of the Delta-resonance (in particular, of the mass and the pion-nucleon-delta coupling constant) from the measured volume dependence of the lowest energy levels should be feasible.


title: Spontaneous chiral symmetry breaking in QCD:a finite-size scaling study on the lattice

abstract: Spontaneous chiral symmetry breaking in QCD with massless quarks at infinite volume can be seen in a finite box by studying, for instance, the dependence of the chiral condensate from the volume and the quark mass. We perform a feasibility study of this program by computing the quark condensate on the lattice in the quenched approximation of QCD at small quark masses. We carry out simulations in various topological sectors of the theory at several volumes, quark masses and lattice spacings by employing fermions with an exact chiral symmetry, and we focus on observables which are infrared stable and free from mass-dependent ultraviolet divergences. The numerical calculation is carried out with an exact variance-reduction technique, which is designed to be particularly efficient when spontaneous symmetry breaking is at work in generating a few very small low-lying eigenvalues of the Dirac operator. The finite-size scaling behaviour of the condensate in the topological sectors considered agrees, within our statistical accuracy, with the expectations of the chiral effective theory. Close to the chiral limit we observe a detailed agreement with the first Leutwyler-Smilga sum rule. By comparing the mass, the volume and the topology dependence of our results with the predictions of the chiral effective theory, we extract the corresponding low-energy constant.


title: Light dynamical fermions on the lattice: toward the chiral regime of QCD

abstract: Algorithmic and technical progress achieved over the last few years makes QCD simulations with light dynamical quarks much faster than before. As a result lattices with pions as light as 250--300 MeV can be simulated with the present generation of computers. I review recent conceptual and numerical progress in this field, with particular emphasis on results obtained and difficulties encountered in simulations with significantly smaller quark masses with respect to previous computations. I also attempt to compare physical results for pion masses and decay constants available to date in the two-flavour theory with expectations from chiral perturbation theory.


title: Renormalization Group Therapy

abstract: We point out a general problem with the procedures commonly used to obtain improved actions from MCRG decimated configurations. Straightforward measurement of the couplings from the decimated configurations, by one of the known methods, can result into actions that do not correctly reproduce the physics on the undecimated lattice. This is because the decimated configurations are generally not representative of the equilibrium configurations of the assumed form of the effective action at the measured couplings. Curing this involves fine-tuning of the chosen MCRG decimation procedure, which is also dependent on the form assumed for the effective action. We illustrate this in decimation studies of the SU(2) LGT using Swendsen and Double Smeared Blocking decimation procedures. A single-plaquette improved action involving five group representations and free of this pathology is given.


title: Lattice Field Theory with the Sign Problem and the Maximum Entropy Method

abstract: Although numerical simulation in lattice field theory is one of the most effective tools to study non-perturbative properties of field theories, it faces serious obstacles coming from the sign problem in some theories such as finite density QCD and lattice field theory with the $\theta$ term. We reconsider this problem from the point of view of the maximum entropy method.


title: Flavour symmetry restoration and kaon weak matrix elements in quenched twisted mass QCD

abstract: We simulate two variants of quenched twisted mass QCD (tmQCD), with degenerate Wilson quarks of masses equal to or heavier than half the strange quark mass. We use Ward identities in order to measure the twist angles of the theory and thus check the quality of the tuning of mass parameters to a physics condition which stays constant as the lattice spacing is varied. Flavour symmetry breaking in tmQCD is studied in a framework of two fully twisted and two standard Wilson quark flavours, tuned to be degenerate in the continuum. Comparing pseudoscalar masses, obtained from connected quark diagrams made of tmQCD and/or standard Wilson quark propagators, we confirm that flavour symmetry breaking effects, which are at most 5%, decrease as we approach the continuum limit. We also compute the pseudoscalar decay constant in the continuum limit, with reduced systematics. As a consequence of improved tuning of the mass parameters at $\beta = 6.1$, we reanalyse our previous $B_K$ results. Our main phenomenological findings are $r_0 f_K = 0.421(7)$ and $\hat B_K = 0.735(71)$.


title: Normalized entropy density of the 3D 3-state Potts model

abstract: Using a multicanonical Metropolis algorithm we have performed Monte Carlo simulations of the 3D 3-state Potts model on $L^3$ lattices with L=20, 30, 40, 50. Covering a range of inverse temperatures from $\beta_{\min}=0$ to $\beta_{\max}=0.33$ we calculated the infinite volume limit of the entropy density $s(\beta)$ with its normalization obtained from $s(0)=\ln 3$. At the transition temperature the entropy and energy endpoints in the ordered and disordered phase are estimated employing a novel reweighting procedure. We also evaluate the transition temperature and the order-disorder interface tension. The latter estimate increases when capillary waves are taken into account.


title: Continuum limit of string formation in 3-d SU(2) LGT

abstract: We study the continuum limit of the string-like behaviour of flux tubes formed between static quarks and anti-quarks in three dimensional SU(2) lattice gauge theory. We compare our simulation data with the predictions of both effective string models as well as perturbation theory. On the string side we obtain clear evidence for convergence of data to predictions of Nambu-Goto theory. We comment on the scales at which the static potential starts departing from one loop perturbation theory and then again being well described by effective string theories. We also estimate the leading corrections to the one-loop perturbative potential as well as the Nambu-Goto effective string. In the intermediate regions we find that a modified Lennard-Jones type potential gives surprisingly good fits.


title: The Monte Carlo method in quantum field theory

abstract: This series of six lectures is an introduction to using the Monte Carlo method to carry out nonperturbative studies in quantum field theories. Path integrals in quantum field theory are reviewed, and their evaluation by the Monte Carlo method with Markov-chain based importance sampling is presented. Properties of Markov chains are discussed in detail and several proofs are presented, culminating in the fundamental limit theorem for irreducible Markov chains. The example of a real scalar field theory is used to illustrate the Metropolis-Hastings method and to demonstrate the effectiveness of an action-preserving (microcanonical) local updating algorithm in reducing autocorrelations. The goal of these lectures is to provide the beginner with the basic skills needed to start carrying out Monte Carlo studies in quantum field theories, as well as to present the underlying theoretical foundations of the method.


title: Fitting two nucleons inside a box: exponentially suppressed corrections to the Luscher's formula

abstract: Scattering observables can be computed in lattice field theory by measuring the volume dependence of energy levels of two particle states. The dominant volume dependence, proportional to inverse powers of the volume, is determined by the phase shifts. This universal relation (\Lu's formula) between energy levels and phase shifts is distorted by corrections which, in the large volume limit, are exponentially suppressed. They may be sizable, however, for the volumes used in practice and they set a limit on how small the lattice can be in these studies. We estimate these corrections, mostly in the case of two nucleons. Qualitatively, we find that the exponentially suppressed corrections are proportional to the {\it square} of the potential (or to terms suppressed in the chiral expansion) and the effect due to pions going ``around the world'' vanishes. Quantitatively, the size of the lattice should be greater than $\approx(5 {fm})^3$ in order to keep finite volume corrections to the phase less than $1^\circ$ for realistic pion mass.


title: Infrared properties of propagators in Landau-gauge pure Yang-Mills theory at finite temperature

abstract: The finite-temperature behavior of gluon and of Faddeev-Popov-ghost propagators is investigated for pure SU(2) Yang-Mills theory in Landau gauge. We present nonperturbative results, obtained using lattice simulations and Dyson-Schwinger equations. Possible limitations of these two approaches, such as finite-volume effects and truncation artifacts, are extensively discussed. Both methods suggest a very different temperature dependence for the magnetic sector when compared to the electric one. In particular, a clear thermodynamic transition seems to affect only the electric sector. These results imply in particular the confinement of transverse gluons at all temperatures and they can be understood inside the framework of the so-called Gribov-Zwanziger scenario of confinement.


title: Lattice Study of Low-lying Nonet Scalar Mesons in Quenched Approximation

abstract: Using lattice QCD simulation in the quenched approximation, we study the $\kappa$ meson, which is ^3P_0 in the quark model, and compare experimental and other lattice data. The $\kappa$ is the lowest scalar meson with strangeness and constitutes the scalar nonet. The obtained mass is much higher than the recent experimental value, and therefore the $\kappa(800)$ is difficult to consider as a simple two-body constituent-quark structure, and may have another unconventional structure.


title: The chiral limit in lattice QCD

abstract: It has been a big challenge for lattice QCD to simulate dynamical quarks near the chiral limit. Theoretically, it is well-known that the naive chiral symmetry cannot be realized on the lattice (the Nielsen-Ninomiya theorem). Also practically, the computational cost rapidly grows as the quark mass is reduced. The JLQCD collaboration started a project to perform simulations with exact but modified chiral symmetry using the the overlap-Dirac operator and the topology conserving action. The latter is helpful to reduce the numerical cost of the dynamical quarks. Our simulation of two-flavor QCD has been successful to reduce the sea quark mass down to a few MeV.


title: Phase diagram of QCD with 2+1 flavors of Wilson quarks at finite temperature and chemical potential

abstract: The first results on lattice at finite temperature $T$ and chemical potential $\mu$ with 2+1 flavors of Wilson quarks are presented. We have studied the dependence of chiral phase transition and deconfinement phase transition on quark mass. Finite volume size analysis and Binder cumulants are used to determine the properties of phase transition. Phase diagram of QCD with 2+1 flavors of Wilson quarks are presented.


title: K_l3 form factor with N_f=2+1 dynamical domain wall fermions: A progress report

abstract: We present the latest results from the UKQCD/RBC collaborations for the K_l3 form factor with 2+1 flavours of dynamical domain wall quarks. Simulations are performed on 16^3x32x16 and 24^3x64x16 lattices with three values of the light quark mass, allowing for an extrapolation to the chiral limit. After interpolating to zero momentum transfer, we obtain the preliminary result f_+^{K\pi}(0)=0.9609(51) (or \Delta f = -0.0161(51)), which is in agreement with the result of Leutwyler & Roos.


title: Improving the improved action

abstract: We investigate the construction of improved actions by the Monte Carlo Renormalization Group method in the context of SU(2) gauge theory utilizing different decimation procedures and effective actions. We demonstrate that the basic self-consistency requirement for correct application of MCRG, i.e. that the decimated configurations are equilibrium configurations of the adopted form of the effective action, can only be achieved by careful fine-tuning of the choice of decimation prescription and/or action.


title: Hypercubic Smeared Links for Dynamical Fermions

abstract: We investigate a variant of hypercubic gauge link smearing where the SU(3) projection is replaced with a normalization to the corresponding unitary group. This smearing is differentiable and thus suitable for use in dynamical fermion simulations using molecular dynamics type algorithms. We show that this smearing is as efficient as projected hypercubic smearing in removing ultraviolet noise from the gauge fields. We test the normalized hypercubic smearing in dynamical improved (clover) Wilson and valence overlap simulations.


title: Can we study Quark Matter in the Quenched Approximation?

abstract: We study a quenched SU(2) lattice gauge theory in which, in an attempt to distinguish between timelike and spacelike gauge fields, the gauge ensemble {U_mu} is generated from a 3 dimensional gauge-Higgs model, the timelike link variables being "reconstructed" from the Higgs fields. The resulting ensemble is used to study quenched quark propagation with non-zero chemical potential mu; in particular, the quark density, chiral and superfluid condensates, meson, baryon and gauge-fixed quark propagators are all studied as functions of mu. While it proves possible to alter the strength of the inter-quark interaction by changing the parameters of the dimensionally reduced model, there is no evidence for any region of parameter space where quarks exhibit deconfined behaviour or thermodynamic observables scale as if there were a Fermi surface.


title: Abelian dominance and the dual Meissner effect in local unitary gauges in SU(2) gluodynamics

abstract: Performing highly precise Monte-Carlo simulations of SU(2) gluodynamics, we observe for the first time Abelian dominance in the confining part of the static potential in local unitary gauges such as the F12 gauge. We also study the flux-tube profile between the quark and antiquark in these local unitary gauges and find a clear signal of the dual Meissner effect. The Abelian electric field is found to be squeezed into a flux tube by the monopole supercurrent. This feature is the same as that observed in the non-local maximally Abelian gauge. These results suggest that the Abelian confinement scenario is gauge independent. Observing the important role of space-like monopoles in the Polyakov gauge also indicates that the monopoles defined on the lattice do not necessarily correspond to those proposed by 't Hooft in the context of Abelian projection.


title: Abelian dominance in local unitary gauges and without gauge-fixing in pure SU(2) QCD

abstract: We perform lattice Monte-Carlo simulations of pure SU(2) QCD using the multi-level method. We find Abelian dominance in local unitary gauges such as those diagonalizing a plaquette. A static potential described by Abelian link fields alone gives us the same string tension as that of a non-Abelian potential. Abelian dominance of the string tension and Abelian flux tube profiles are observed also without gauge-fixing, i.e., without any Abelian projection. On the basis of these results, we propose a simple gauge-independent Abelian confinement scenario without any Abelian projection. All color components of the non-Abelian field strength become Abelian dominant in the infrared region. The Abelian dual Meissner effect works in any color direction. Abelian neutral states in any color directions which are just non-Abelian color-singlet can exist as a physical state. In this way, the non-Abelian color confinement could be understood in the framework of the Abelian dual Meissner effect.


title: The Sign Problem via Imaginary Chemical Potential

abstract: We calculate an analogue of the average phase factor of the staggered fermion determinant at imaginary chemical potential. Our results from the lattice agree well with the analytical predictions in the microscopic regime for both quenched and phase-quenched QCD. We demonstrate that the average phase factor in the microscopic domain is dominated by the lowest-lying Dirac eigenvalues.


title: Hadronic form factors in Lattice QCD at small and vanishing momentum transfer

abstract: The introduction of partially twisted boundary conditions allows weak and electromagnetic form factors to be evaluated at specified values of the hadronic momenta (and hence momentum transfers) in lattice simulations. We present and demonstrate this technique for the computation of the K->pi semileptonic form factor at zero momentum transfer and for the electromagnetic form factor of the pion at arbitrarily small momentum transfers. These exploratory computations are carried out in full QCD with 3 flavours of sea quarks, but with only two values of m_u=m_d which limits our ability to perform the chiral extrapolations. The results should therefore be viewed primarily as a demonstration of the feasibility of the method. For the K-> pi form factor we compare the new technique to the conventional approach and for the pion form factor we assess our results for very small momentum transfer with the help of chiral perturbation theory.


title: Non-perturbative improvement of the axial current with three dynamical flavors and the Iwasaki gauge action

abstract: We perform a non-perturbative determination of the improvement coefficient c_A to remove O(a) discretization errors in the axial vector current in three-flavor lattice QCD with the Iwasaki gauge action and the standard O$(a)$-improved Wilson quark action. An improvement condition with a good sensitivity to c_A is imposed at constant physics. Combining our results with the perturbative expansion, c_A is now known rather precisely for 1/a \gtrsim 1.6 GeV.


title: Confining string and its widening in HP1 embedding approach

abstract: Structure of confining string in terms of the topological charge density and the action density is studied in SU(2) Yang-Mills theory on the lattice using HP1 sigma-model embedding approach. We find that the confining flux tube noticeably suppresses both the topological charge and the action densities. Beyond the string formation length the string cross section in terms of these quantities is well described by a Gaussian profile. In both cases the squared string width is found to be a logarithmic function of the string length confirming the Luscher widening of the chromoelectric string. Characteristic string scales in terms of the topological and action densities are estimated as well.


title: Spectral functions at small energies and the electrical conductivity in hot, quenched lattice QCD

abstract: In lattice QCD, the Maximum Entropy Method can be used to reconstruct spectral functions from euclidean correlators obtained in numerical simulations. We show that at finite temperature the most commonly used algorithm, employing Bryan's method, is inherently unstable at small energies and give a modification that avoids this. We demonstrate this approach using the vector current-current correlator obtained in quenched QCD at finite temperature. Our first results indicate a small electrical conductivity above the deconfinement transition.

categories: hep-lat cond-mat.other hep-ph hep-th nucl-th


title: BB Potentials in Quenched Lattice QCD

abstract: The potentials between two B-mesons are computed in the heavy-quark limit using quenched lattice QCD at $m_\pi\sim 400~{\rm MeV}$. Non-zero central potentials are clearly evident in all four spin-isospin channels, (I,s_l) = (0,0) , (0,1) , (1,0) , (1,1), where s_l is the total spin of the light degrees of freedom. At short distance, we find repulsion in the $I\ne s_l$ channels and attraction in the I=s_l channels. Linear combinations of these potentials that have well-defined spin and isospin in the t-channel are found, in three of the four cases, to have substantially smaller uncertainties than the potentials defined with the s-channel (I,s_l), and allow quenching artifacts from single hairpin exchange to be isolated. The BB*\pi coupling extracted from the long-distance behavior of the finite-volume t-channel potential is found to be consistent with quenched calculations of the matrix element of the isovector axial-current. The tensor potentials in both of the s_l = 1 channels are found to be consistent with zero within calculational uncertainties.


title: Gauge Field Strength Tensor from the Overlap Dirac Operator

abstract: We derive the classical continuum limit of the operator tr$_s \sigma_{\mu\nu} D^{ov}(x,x)$ with $D^{ov}$ being the overlap Dirac operator and show that it corresponds to the gauge field strength tensor $F_{\mu\nu}(x)$.


title: Energies and radial distributions of B_s mesons - the effect of hypercubic blocking

abstract: This is a follow-up to our earlier work for the energies and the charge (vector) and matter (scalar) distributions for S-wave states in a heavy-light meson, where the heavy quark is static and the light quark has a mass about that of the strange quark. We study the radial distributions of higher angular momentum states, namely P- and D-wave states, using a "fuzzy" static quark. A new improvement is the use of hypercubic blocking in the time direction, which effectively constrains the heavy quark to move within a 2a hypercube (a is the lattice spacing). The calculation is carried out with dynamical fermions on a 16^3 times 32 lattice with a lattice spacing approximately 0.10 fm generated using the non-perturbatively improved clover action. The configurations were generated by the UKQCD Collaboration using lattice action parameters beta = 5.2, c_SW = 2.0171 and kappa = 0.1350. In nature the closest equivalent of this heavy-light system is the B_s meson. Attempts are now being made to understand these results in terms of the Dirac equation.


title: Mixed Action Lattice Spacing Effects on the Nucleon Axial Charge

abstract: We study the nucleon axial charge in the chiral perturbation theory for a mixed lattice action of Ginsparg-Wilson valence quarks and staggered sea quarks. In particular, we investigate the lattice spacing $a^2$-dependence of the neutron to proton axial transition matrix element. By using the known low-energy constants (LEC's) and an estimated value of a new LEC which appears in the calculation, we see a large lattice spacing effect on $g_A$.


title: Domain wall filters

abstract: We propose using the extra dimension separating the domain walls carrying lattice quarks of opposite handedness to gradually filter out the ultraviolet fluctuations of the gauge fields that are felt by the fermionic excitations living in the bulk. This generalization of the homogeneous domain wall construction has some theoretical features that seem nontrivial.


title: Pseudoscalar and vector meson form factors from lattice QCD

abstract: We present a study of the pseudoscalar and vector meson form factors, calculated using the Fat-Link Irrelevant Clover (FLIC) action in the framework of Quenched Lattice QCD. Of particular interest is the determination of a negative quadrupole moment, indicating that the $\rho$ meson is not spherically symmetric.


title: Hadron Scattering in an Asymmetric Box

abstract: We propose to study hadron-hadron scattering using lattice QCD in an asymmetric box which allows one to access more non-degenerate low-momentum modes for a given volume. The conventional L\"{u}scher's formula applicable in a symmetric box is modified accordingly. To illustrate the feasibility of this approach, pion-pion elastic scattering phase shifts in the I=2, J=0 channel are calculated within quenched approximation using improved gauge and Wilson fermion actions on anisotropic lattices in an asymmetric box. After the chiral and continuum extrapolation, we find that our quenched results for the scattering phase shifts in this channel are consistent with the experimental data when the three-momentum of the pion is below 300MeV. Agreement is also found when compared with previous theoretical results from lattice and other means. Moreover, with the usage of asymmetric volume, we are able to compute the scattering phases in the low-momentum range (pion three momentum less than about 350MeV in the center of mass frame) for over a dozen values of the pion three-momenta, much more than using the conventional symmetric box with comparable volume.


title: Phase diagram of the three dimensional Thirring model - A Monte Carlo study

abstract: Certain approximate solutions of the continuum Schwinger-Dyson Equations (SDEs) predict chiral symmetry breaking in the 3d Thirring model when the number of fermion flavors N_f<4.32 whereas others predict symmetry breaking for all N_f. Our results from Monte Carlo simulations with N_f=6, predict a second order chiral phase transition. The critical coupling in this case corresponds to an ultra-violet fixed point of the renormalization group defining a non-trivial continuum limit. Further, our numerical simulations provide an estimate for the critical number of fermion flavors, N_fc \approx 6.5.


title: Simulating Dense Matter

abstract: I review the Sign Problem hindering lattice QCD simulations of dense baryonic matter, focussing where possible on its physical relevance. The possibility of avoiding the Sign Problem via a duality transformation is also briefly considered. Finally, I review evidence for deconfinement at non-zero quark density in recent simulations of Two Color QCD.


title: Spectral sums of the Dirac-Wilson Operator and their relation to the Polyakov loop

abstract: We investigate and compute spectral sums of the Wilson lattice Dirac operator for quenched SU(3) gauge theory. It is demonstrated that there exist sums which serve as order parameters for the confinement-deconfinement phase transition and get their main contribution from the IR end of the spectrum. They are approximately proportional to the Polyakov loop. In contrast to earlier studied spectral sums some of them are expected to have a well-defined continuum limit.


title: Scalar-Quark Systems and Chimera Hadrons in SU(3)_c Lattice QCD

abstract: Light scalar-quarks \phi (colored scalar particles or idealized diquarks) and their color-singlet hadronic states are studied with quenched SU(3)_c lattice QCD in terms of mass generation in strong interaction without chiral symmetry breaking. We investigate ``scalar-quark mesons'' \phi^\dagger \phi and ``scalar-quark baryons'' \phi\phi\phi which are the bound states of scalar-quarks \phi. We also investigate the bound states of scalar-quarks \phi and quarks \psi, i.e., \phi^\dagger \psi, \psi\psi\phi and \phi\phi\psi, which we name ``chimera hadrons''. All the new-type hadrons including \phi are found to have a large mass even for zero bare scalar-quark mass m_\phi=0 at a^{-1}\simeq 1GeV. We find that the constituent scalar-quark and quark picture is satisfied for all the new-type hadrons. Namely, the mass of the new-type hadron composed of m \phi's and n \psi's, M_{{m}\phi+{n}\psi}, satisfies M_{{m}\phi+{n}\psi}\simeq {m} M_\phi +{n} M_\psi, where M_\phi and M_\psi are the constituent scalar-quark and quark mass, respectively. M_\phi at m_\phi=0 estimated from these new-type hadrons is 1.5-1.6GeV, which is larger than that of light quarks, M_\psi\simeq 400{\rm MeV}. Therefore, in the systems of scalar-quark hadrons and chimera hadrons, scalar-quarks acquire large mass due to large quantum corrections by gluons. Together with other evidences of mass generations of glueballs and charmonia, we conjecture that all colored particles generally acquire a large effective mass due to dressed gluon effects.


title: Infra-red divergences in the large-N expansion

abstract: We investigate a vectorial O(N) model with a generic nearest-neighbor interaction W(\bsigma_i\cdot \bsigma_j) (depending on {\cal N} tunable parameters), a Yukawa (and Gross Neveu) model with N_f fermions at finite temperature and the vectorial \phi^6 model, in the large N (N_f) limit. All this models exhibit a Mean Field critical point for N=\infinity. When 1/N fluctuations are included, infra red divergences appear near the critical point. In the framework of a generalized 1/N expansion we show that these divergences are related to a universal crossover mechanism between the Mean Field universality class (N=\infinity) and the nonclassical one for N<\infinity. For the generic nearest-neighbor interaction multicritical points are also investigated.


title: Speeding up Domain Wall Fermion Algorithms using QCDLAB

abstract: Simulating lattice QCD with chiral fermions and indeed using Domain Wall Fermions continues to be challenging project however large are concurrent computers. One obvious bottleneck is the slow pace of prototyping using the low level coding which prevails in most, if not all, lattice projects. Recently, we came up with a new proposal, namely QCDLAB, a high level language interface, which we believe will boost our endeavours to rapidly code lattice prototype applications in lattice QCD using MATLAB/OCTAVE language and environment. The first version of the software, QCDLAB 1.0 offers the general framework on how to achieve this goal by simulating set of the lattice Schwinger model {\tt http://phys.fshn.edu.al/qcdlab.html}. In this talk we introduce QCDLAB 1.1, which extends QCDLAB 1.0 capabilities for real world lattice computations with Wilson and Domain Wall fermions.


title: Scaling behavior and positivity violation of the gluon propagator in full QCD

abstract: The Landau-gauge gluon propagator is studied using the coarse and fine dynamical MILC configurations. The effects of dynamical quarks are clearly visible and lead to a reduction of the nonperturbative infrared enhancement relative to the quenched case. Lattice spacing effects are studied and found to be small. The gluon spectral function is shown to clearly violate positivity in both quenched and full QCD.


title: Low-lying Dirac eigenmodes and monopoles in 4D compact QED

abstract: We study the properties of low-lying Dirac modes in quenched compact QED at $\beta$=0.99, 1.01 and 1.03, employing $12^3\times 12$ lattices and the overlap formalism for the fermion action. We pay special attention to the spatial distributions of the low-lying Dirac modes. Near-zero modes are found to have universal anti-correlations with monopole currents below/above the critical $\beta$. We also study the nearest-neighbor level spacing distribution of Dirac eigenvalues and find a signal of a Wigner-Poisson transition. We make a possible speculation on the chiral phase transition in 4D compact QED.


title: How to Put a Heavier Higgs on the Lattice

abstract: Lattice work, exploring the Higgs mass triviality bound, seems to indicate that a strongly interacting scalar sector in the minimal standard model cannot exist while low energy QCD phenomenology seems to indicate that it could. We attack this puzzle using the 1/N expansion and discover a simple criterion for selecting a lattice action that is more likely to produce a heavy Higgs particle. Our large $N$ calculation suggests that the Higgs mass bound might be around $850 GeV$, which is about 30% higher than previously obtained.


title: Spectral Density Study of the SU(3) Deconfining Phase Transition

abstract: We present spectral density reweighting techniques adapted to the analysis of a time series of data with a continuous range of allowed values. In a first application we analyze action and Polyakov line data from a Monte Carlo simulation on $L_t L^3 (L_t=2,4)$ lattices for the SU(3) deconfining phase transition. We calculate partition function zeros, as well as maxima of the specific heat and of the order parameter susceptibility. Details and warnings are given concerning i) autocorrelations in computer time and ii) a reliable extraction of partition function zeros. The finite size scaling analysis of these data leads to precise results for the critical couplings $\beta_c$, for the critical exponent $\nu$ and for the latent heat $\triangle s$. In both cases ($L_t=2$ and 4), the first order nature of the transition is substantiated.


title: Vectorized Cluster Search

abstract: Contrary to conventional wisdom, the construction of clusters on a lattice can easily be vectorized, namely over each ``generation'' in a breadth first search. This applies directly to, e.g., the {\it single cluster} variant of the Swendsen-Wang algorithm. On a Cray Y-MP, total CPU time was reduced by a factor 3.5 -- 7 in actual applications.


title: Multi-Grid Monte Carlo III. Two-Dimensional O(4)-Symmetric Nonlinear $\sigma$-Model

abstract: We study the dynamic critical behavior of the multi-grid Monte Carlo (MGMC) algorithm with piecewise-constant interpolation applied to the two-dimensional O(4)-symmetric nonlinear $\sigma$-model [= SU(2) principal chiral model], on lattices up to $256 \times 256$. We find a dynamic critical exponent $z_{int,{\cal M}^2} = 0.60 \pm 0.07$ for the W-cycle and $z_{int,{\cal M}^2} = 1.13 \pm 0.11$ for the V-cycle, compared to $z_{int,{\cal M}^2} = 2.0 \pm 0.15$ for the single-site heat-bath algorithm (subjective 68% confidence intervals). Thus, for this asymptotically free model, critical slowing-down is greatly reduced compared to local algorithms, but not completely eliminated. For a $256 \times 256$ lattice, W-cycle MGMC is about 35 times as efficient as a single-site heat-bath algorithm.


title: Monte Carlo Simulations of Higgs-Fermion Systems

abstract: To gain understanding of the Higgs-fermion sector of the standard model, we study the one-component $Z_2$ symmetric and the four-component O(4) symmetric scalar models coupled to staggered fermions using the hybrid Monte Carlo algorithm. We map out the phase diagrams, and show that the $Z_2$ model has a tree level perturbative behaviour at all points in the broken phase. The O(4) model on the other hand is shown to have two characteristically different behaviours; one for large Yukawa couplings where the fermions get infinitely heavy and decouple in the continuum limit, and one for small Yukawa couplings where the fermions remain light. For very small Yukawa couplings the fermions show the expected tree level perturbative behaviour and for larger values the influence of the fermions becomes substantial. After estimating the finite size effects at small Yukawa couplings we make relatively accurate measurements of the scalar mass and wave function renormalization constants at the point $\kappa=0.0$ and $y=0.85-0.95$. Even though this is not the largest value possible for the Yukawa coupling we are able to show that the bound of the Higgs mass will move up significantly, from around $600 GeV$ to around $900 GeV$, by including fermions in the model. Likewise we show that a bound can be put on the fermion mass, around $200 GeV$. The largest value of the bare Yukawa coupling is obtained at rather large negative $\kappa$. Due to bad convergence rates in the inversion of the fermion matrix, which is needed in the updating procedure, this region has not been possible to investigate.


title: A Study of Symmetry Restoration at Finite Temperature in the O(4) Model Using Anisotropic Lattices

abstract: Results of investigations of the O(4) spin model at finite temperature using anisotropic lattices are presented. In both the large $N$ approximation and the numerical simulations using the Wolff cluster algorithm we find that the ratio of the symmetry restoration temperature $T_{\rm SR}$ to the Higgs mass $m_{\rm H}$ is independent of the anisotropy. We obtain a lower bound of $0.59 \pm 0.04$ for the ratio, $T_{\rm SR}/m_{\rm H}$, at $m_{\rm H}a \simeq 0.5$, which is lowered further by about 10% at $m_{\rm H}a \simeq 1.$


title: Dynamic Critical Behaviour of Wolff's Algorithm for $RP^N$ $\sigma$-Models

abstract: We study the performance of a Wolff-type embedding algorithm for $RP^N$ $\sigma$-models. We find that the algorithm in which we update the embedded Ising model \`a la Swendsen-Wang has critical slowing-down as $z_\chi \approx 1$. If instead we update the Ising spins with a perfect algorithm which at every iteration produces a new independent configuration, we obtain $z_\chi \approx 0$. This shows that the Ising embedding encodes well the collective modes of the system, and that the behaviour of the first algorithm is connected to the poor performance of the Swendsen-Wang algorithm in dealing with a frustrated Ising model.


abstract: The cutoff dependence of the Scalar Sector of the Minimal Standard Model can result in an increase of the existing triviality bound estimates of the Higgs mass. We present a large $N$ calculation and some preliminary N=4 results that suggest that the increase can be as large as 30%, resulting to a bound of about 850 G eV.


title: QCD Hadron Spectroscopy with Staggered Dynamical Quarks at $\beta = 5.6$

abstract: We present preliminary results from the 1991 HEMCGC simulations with staggered dynamical fermions on a $16^3 \times 32$ lattice at $\beta = 5.6$ with sea quark masses $am_q = 0.025$ and 0.01. The spectroscopy was done both for staggered valence quarks with mass equal to the sea quark masses and for Wilson valence quarks at six different values for $\kappa$, 0.1320, 0.1410, 0.1525, 0.1565, 0.1585, and 0.1600. In addition to the measurements performed in our earlier work, we also measured the $\Delta$ and other `extended' hadrons for staggered valence quarks and pseudo-scalar decay constants and vector meson matrix elements, the wave function at the origin, for Wilson valence quarks.


title: Path Integrals and Voronin's Theorem on the Universality of the Riemann Zeta Function

abstract: We explore a new approach to the path integral for a latticized quantum theory. This talk is based on work with N. Khuri and H. Ren.


title: On the Dynamics of Light Wilson Quarks

abstract: We describe recent results obtained as part of the High Energy Monte Carlo Grand Challenge (HEMCGC) project concerning the behaviour of lattice QCD with light dynamical Wilson quarks. We show that it is possible to reach regions of parameter space with light pions $m_\pi\ll0.2/a$, but that the equilibration time for such a system is at least of the order of 1,000 unit-length Hybrid Monte Carlo (HMC) trajectories (about a Gigaflop/sec-year). If the Hybrid Molecular Dynamics (HMD) algorithm is used with the same parameters it gives incorrect results.


title: Monte Carlo Calculation of the Surface Free Energy for the 2d 7-state Potts Model, and an Estimate for 4d SU(3) Gauge Theory

abstract: Using the recently proposed multicanonical ensemble, we perform Monte Carlo simulation for the 2d 7-state Potts model and calculate its surface free energy density (surface tension) to be $2 f^s = 0.0241 \pm 0.0010$. This is an order of magnitude smaller than other estimates in the recent literature. Relying on existing Monte Carlo data, we also give a preliminary estimate for the surface tension of 4d SU(3) lattice gauge theory with $L_t=2$.


title: Random Walks and the Correlation Length Critical Exponent in Scalar Quantum Field Theory

abstract: The distance scale for a quantum field theory is the correlation length $\xi$, which diverges with exponent $\nu$ as the bare mass approaches a critical value. If $t=m^{2}-m_{c}^{2}$, then $\xi=m_{P}^{-1} \sim t^{-\nu}$ as $t \to 0$. The two-point function of a scalar field has a random walk representation. The walk takes place in a background of fluctuations (closed walks) of the field itself. We describe the connection between properties of the walk and of the two-point function. Using the known behavior of the two point function, we deduce that the dimension of the walk is $d_{w}=\phi / \nu$ and that there is a singular relation between $t$ and the energy per unit length of the walk $\theta \sim t^{\phi}$ that is due to the singular behavior of the background at $t=0$. ($\phi$ is a computable crossover exponent.)


title: Thermodynamics of Lattice QCD with 2 Light Dynamical (Staggered) Quark Flavours on a $16^3\times 8$ Lattice

abstract: Lattice QCD with 2 light staggered quark flavours is being simulated on a $16^3\times8$ lattice to study the transition from hadronic matter to a quark gluon plasma. We have completed runs at $m_q=0.0125$ and are extending this to $m_q=0.00625$. We also examine the addition of a non-dynamical "strange" quark. Thermodynamic order parameters are being measured across the transition and further into the plasma phase, as are various screening lengths. No evidence for a first order transition is seen, and we estimate the transition temperature to be $TY_c=143(7) MeV$.


title: Multicanonical Ensemble: A New Approach to Simulate First-order Phase Transitions

abstract: Relying on the recently proposed multicanonical algorithm, we present a numerical simulation of the first order phase transition in the 2d 10-state Potts model on lattices up to sizes $100\times100$. It is demonstrated that the new algorithm $lacks$ an exponentially fast increase of the tunneling time between metastable states as a function of the linear size $L$ of the system. Instead, the tunneling time diverges approximately proportional to $L^{2.65}$. Thus the computational effort as counted per degree of freedom for generating an independent configuration in the unstable region of the model rises proportional to $V^{2.3}$, where $V$ is the volume of the system. On our largest lattice we gain more than two orders of magnitude as compared to a standard heat bath algorithm. As a first physical application we report a high precision computation of the interfacial tension.


title: Is there a $\rho$ in the O(4) $\lambda \phi^4_4$ theory?

abstract: A Monte Carlo simulation of the O(4) $\lambda \phi^4$ theory in the broken phase is performed on a hypercubic lattice in search of an I=1, J=1 resonance. The region of the cutoff theory where the interaction is strong is investigated since it is there that a resonance would be expected to have a better chance to form. In that region the presence of an I=1, J=1 resonance with mass below the cutoff is excluded.


title: Multimagnetical Simulations

abstract: We modified the recently proposed multicanonical MC algorithm for the case of a magnetic field driven order--order phase transition. We test this {\it multimagnetic} Monte Carlo algorithm for the D=2 Ising model at $\beta=0.5$ and simulate square lattices up to size $100 \times 100$. On these lattices with periodic boundary conditions it is possible to enhance the appearance of order-order interfaces during the simulation by many orders of magnitude as compared to the standard Monte Carlo simulation.


title: A New Approach to Spin Glass Simulations

abstract: We present a recursive procedure to calculate the parameters of the recently introduced multicanonical ensemble and explore the approach for spin glasses. Temperature dependence of the energy, the entropy and other physical quantities are easily calculable and we report results for the zero temperature limit. Our data provide evidence that the large $L$ increase of the ergodicity time is greatly improved. The multicanonical ensemble seems to open new horizons for simulations of spin glasses and other systems which have to cope with conflicting constraints.


title: A Numerical Test of KPZ Scaling: Potts Models Coupled to Two-Dimensional Quantum Gravity

abstract: We perform Monte Carlo simulations using the Wolff cluster algorithm of the q=2 (Ising), 3, 4 and q=10 Potts models on dynamical phi-cubed graphs of spherical topology with up to 5000 nodes. We find that the measured critical exponents are in reasonable agreement with those from the exact solution of the Ising model and with those calculated from KPZ scaling for q=3,4 where no exact solution is available. Using Binder's cumulant we find that the q=10 Potts model displays a first order phase transition on a dynamical graph, as it does on a fixed lattice. We also examine the internal geometry of the graphs generated in the simulation, finding a linear relationship between ring length probabilities and the central charge of the Potts model


title: Multiple Potts Models Coupled to Two-Dimensional Quantum Gravity

abstract: We perform Monte Carlo simulations using the Wolff cluster algorithm of {\it multiple} $q=2,3,4$ state Potts models on dynamical phi-cubed graphs of spherical topology in order to investigate the $c>1$ region of two-dimensional quantum gravity. Contrary to naive expectation we find no obvious signs of pathological behaviour for $c>1$. We discuss the results in the light of suggestions that have been made for a modified DDK ansatz for $c>1$.


title: Critical Behavior of Dynamically Triangulated Quantum Gravity in Four Dimensions

abstract: We performed detailed study of the phase transition region in Four Dimensional Simplicial Quantum Gravity, using the dynamical triangulation approach. The phase transition between the Gravity and Antigravity phases turned out to be asymmetrical, so that we observed the scaling laws only when the Newton constant approached the critical value from perturbative side. The curvature susceptibility diverges with the scaling index $-.6$. The physical (i.e. measured with heavy particle propagation) Hausdorff dimension of the manifolds, which is 2.3 in the Gravity phase and 4.6 in the Antigravity phase, turned out to be 4 at the critical point, within the measurement accuracy. These facts indicate the existence of the continuum limit in Four Dimensional Euclidean Quantum Gravity.


title: Microcanonical cluster algorithms

abstract: I propose a numerical simulation algorithm for statistical systems which combines a microcanonical transfer of energy with global changes in clusters of spins. The advantages of the cluster approach near a critical point augment the speed increases associated with multi-spin coding in the microcanonical approach. The method also provides a limited ability to tune the average cluster size.


title: Point-Split Lattice Operators for B Decays

abstract: The matrix element which determines the B meson decay constant can be measured on the lattice using an effective field theory for heavy quarks. Various discretizations of the heavy-light bilinears which appear in this and other B decay matrix elements are possible. The heavy-light bilinear currently used for the determination of the B meson decay constant on the lattice suffers a substantial one-loop renormalization. In this paper, we compute the one-loop renormalizations of the discretizations in which the heavy and light fields in the bilinear are separated by one lattice spacing, and discuss their application. Readers of this paper may also be interested in our paper on the application of Symanzik's improvement program to heavy-light currents (paper number 9203221 on hep-ph).


title: Chiral Perturbation Theory for the Quenched Approximation of QCD

abstract: [This version is a minor revision of a previously submitted preprint. Only references have been changed.] We describe a technique for constructing the effective chiral theory for quenched QCD. The effective theory which results is a lagrangian one, with a graded symmetry group which mixes Goldstone bosons and fermions, and with a definite (though slightly peculiar) set of Feynman rules. The straightforward application of these rules gives automatic cancellation of diagrams which would arise from virtual quark loops. The techniques are used to calculate chiral logarithms in $f_K/f_\pi$, $m_\pi$, $m_K$, and the ratio of $\langle{\bar s}s\rangle$ to $\langle{\bar u}u\rangle$. The leading finite-volume corrections to these quantities are also computed. Problems for future study are described.


title: Hadron Spectrum in QCD with Valence Wilson Fermions and Dynamical Staggered Fermions at $6/g^2=5.6

abstract: We present an analysis of hadronic spectroscopy for Wilson valence quarks with dynamical staggered fermions at lattice coupling $6/g^2 = \beta=5.6$ at sea quark mass $am_q=0.01$ and 0.025, and of Wilson valence quarks in quenched approximation at $\beta=5.85$ and 5.95, both on $16^3 \times 32$ lattices. We make comparisons with our previous results with dynamical staggered fermions at the same parameter values but on $16^4$ lattices doubled in the temporal direction.


title: Spectrum of Relativistic Fermions in a 2d Doped Lattice

abstract: Motivated by some previous work on fermions on random lattices and by suggestions that impurities could trigger parity breaking in 2d crystals, we have analyzed the spectrum of the Dirac equation on a two dimensional square lattice where sites have been removed randomly --- a doped lattice. We have found that the system is well described by a sine-Gordon action. The solitons of this model are the lattice fermions, which pick a quartic interaction due to the doping and become Thirring fermions. They also get an effective mass different from the lagrangian mass. The system seems to exhibit spontaneous symmetry breaking, exactly as it happens for a randomly triangulated lattice. The associated ``Goldstone boson" is the sine-Gordon scalar. We argue, however, that the peculiar behaviour of the chiral condensate is due to finite size effects.


title: Influence of the Measure on Simplicial Quantum Gravity in Four Dimensions

abstract: We investigate the influence of the measure in the path integral for Euclidean quantum gravity in four dimensions within the Regge calculus. The action is bounded without additional terms by fixing the average lattice spacing. We set the length scale by a parameter $\beta$ and consider a scale invariant and a uniform measure. In the low $\beta$ region we observe a phase with negative curvature and a homogeneous distribution of the link lengths independent of the measure. The large $\beta$ region is characterized by inhomogeneous link lengths distributions with spikes and positive curvature depending on the measure.


title: Phenomenology with Wilson fermions using smeared sources

abstract: We investigate the use of two types of non-local (``smeared'') sources for quark propagators in quenched lattice QCD at $\beta=6.0$ using Wilson fermions at $\kappa=0.154$ and $0.155$. We present results for the hadron mass spectrum, meson decay constants, quark masses, the chiral condensate and the quark distribution amplitude of the pion. The use of smeared sources leads to a considerable improvement over previous results. We find a disturbing discrepancy between the baryon spectra obtained using Wuppertal and wall sources. We find good signals in the ratio of correlators used to calculate the quark mass and the chiral condensate and show that the extrapolation to the chiral limit is smooth.


title: Scalar and axial matrix elements of the nucleon: sea quark content

abstract: Sea quark contributions to the scalar density and the axial current matrix elements of the nucleon are studied in lattice qcd with two flavours of dynamical wilson fermions. the results are compared to trends in heavy quark mass expansions, and contrasted with the numbers obtained using dynamical staggered fermions.


title: The Free Energy of Spherical Bubbles in Lattice SU(3) Gauge Theory

abstract: We study the coefficients of the expansion $F(R) = 1/3 c_3 R^3 + 1/2 c_2 R^2 + c_1 R$ of the free energy of spherical bubbles at $T=T_c$ in pure glue QCD using lattice Monte Carlo techniques. The coefficient $c_3$ vanishes at $T=T_c$ and our results suggest that the sign and the order of magnitude of $c_1$ is in agreement with the value $c_1=\pm 32\pi T_c^2/9$ (- for hadronic bubbles in quark phase, + for quark bubbles in hadronic phase) computed by Mardor and Svetitsky from the MIT bag model. The parameter $c_2$ is small in agreement with earlier determinations.


title: Multigrid Methods in Lattice Field Computations

abstract: The multigrid methodology is reviewed. By integrating numerical processes at all scales of a problem, it seeks to perform various computational tasks at a cost that rises as slowly as possible as a function of $n$, the number of degrees of freedom in the problem. Current and potential benefits for lattice field computations are outlined. They include: $O(n)$ solution of Dirac equations; just $O(1)$ operations in updating the solution (upon any local change of data, including the gauge field); similar efficiency in gauge fixing and updating; $O(1)$ operations in updating the inverse matrix and in calculating the change in the logarithm of its determinant; $O(n)$ operations per producing each independent configuration in statistical simulations (eliminating CSD), and, more important, effectively just $O(1)$ operations per each independent measurement (eliminating the volume factor as well). These potential capabilities have been demonstrated on simple model problems. Extensions to real life are explored.


title: Blockspin Cluster Algorithms for Quantum Spin Systems

abstract: Cluster algorithms are developed for simulating quantum spin systems like the one- and two-dimensional Heisenberg ferro- and anti-ferromagnets. The corresponding two- and three-dimensional classical spin models with four-spin couplings are maped to blockspin models with two-blockspin interactions. Clusters of blockspins are updated collectively. The efficiency of the method is investigated in detail for one-dimensional spin chains. Then in most cases the new algorithms solve the problems of slowing down from which standard algorithms are suffering.


title: Acceptance Rates in Multigrid Monte Carlo

abstract: An approximation formula is derived for acceptance rates of nonlocal Metropolis updates in simulations of lattice field theories. The predictions of the formula agree quite well with Monte Carlo simulations of 2-dimensional Sine Gordon, XY and phi**4 models. The results are consistent with the following rule: For a critical model with a fundamental Hamiltonian H(phi) sufficiently high acceptance rates for a complete elimination of critical slowing down can only be expected if the expansion of < H(phi+psi) > in terms of the shift psi contains no relevant term (mass term).


title: Decoupling a Fermion Whose Mass Comes from a Yukawa Coupling: Nonperturbative Considerations

abstract: Perturbative analyses seem to suggest that fermions whose mass comes solely from a Yukawa coupling to a scalar field can be made arbitrarily heavy, while the scalar remains light. The effects of the fermion can be summarized by a local effective Lagrangian for the light degrees of freedom. Using weak coupling and large N techniques, we present a variety of models in which this conclusion is shown to be false when nonperturbative variations of the scalar field are considered. The heavy fermions contribute nonlocal terms to the effective action for light degrees of freedom. This resolves paradoxes about anomalous and nonanomalous symmetry violation in these models. Application of these results to lattice gauge theory imply that attempts to decouple lattice fermion doubles by the method of Swift and Smit cannot succeed, a result already suggested by lattice calculations.


title: Scaling topological charge in the CP^3 spin model

abstract: The CP^3 spin model is simulated at large correlation lengths in two dimensions. An overrelaxation algorithm is employed which yields reduced critical slowing down with dynamical exponents z around unity. We compare our results with recent multigrid data on the massgap m and the spin susceptibility and confirm the absence of asymptotic scaling. As a new result we find scaling for the universal topological susceptibility with values extrapolating to chi_t / m^2 = 0.156(2) in the continuum limit.


title: Chiral Invariance and Species Doublers in Generic Fermion Models on the Lattice

abstract: Discussions are made on the structures of chirally invariant lattice actions without any restriction of hermiticity. With the help of the Ward-Takahashi identity a general conclusion can be derived that there must be species doublers in any chirally invariant model provided that the model is chosen as well-regularized, that is, there is no singularity in the propagator after introducing fermion mass on the lattice. Various examples are discussed to pick up better models defined in the sense that the number of species doubler is smaller than that of the naive Dirac action.


title: Non--perturbative bounds on the Higgs mass in the minimal standard model

abstract: This is a slightly extended version of the talk delivered at the Topical Workshop ``Non perturbative aspects of chiral gauge theories'', Accademia Nazionale dei Lincei, Roma, 9-11 March, 1992. Abstract: The Higgs mass in the minimal standard model is bounded by triviality and vacuum stability in the range 50--100 $GeV$ to 700--900 $GeV$. Recent results will be presented in brief and directions for future work will be proposed.


title: Spontaneous symmetry breaking on the lattice generated by Yukawa interaction

abstract: We study by numerical simulation a lattice Yukawa model with naive fermions at intermediate values of the Yukawa coupling $y$ when the nearest neighbour coupling $\kp$ of the scalar field $\Phi$ is very weakly ferromagnetic ($\kp \approx 0$) or even antiferromagnetic ($\kappa < 0$) and the nonvanishing value of $\vev$ is generated by the Yukawa interaction. The renormalized Yukawa coupling $y_R$ achieves here its maximal value and this $y$-region is thus of particular importance for lattice investigations of strong Yukawa interaction. However, here the scalar field propagators have a very complex structure caused by fermion loop corrections and by the proximity of phases with antiferromagnetic properties. We develop methods for analyzing these propagators and for extracting the physical observables. We find that going into the negative $\kp$ region, the scalar field renormalization constant becomes small and $y_R$ does not seem to exceed the unitarity bound, making the existence of a nontrivial fixed point in the investigated Yukawa model quite unlikely.


title: Wolff-Type Embedding Algorithms for General Nonlinear $\sigma$-Models

abstract: We study a class of Monte Carlo algorithms for the nonlinear $\sigma$-model, based on a Wolff-type embedding of Ising spins into the target manifold $M$. We argue heuristically that, at least for an asymptotically free model, such an algorithm can have dynamic critical exponent $z \ll 2$ only if the embedding is based on an (involutive) isometry of $M$ whose fixed-point manifold has codimension 1. Such an isometry exists only if the manifold is a discrete quotient of a product of spheres. Numerical simulations of the idealized codimension-2 algorithm for the two-dimensional $O(4)$-symmetric $\sigma$-model yield $z_{int,{\cal M}^2} = 1.5 \pm 0.5$ (subjective 68\% confidence interval), in agreement with our heuristic argument.


title: Rigorous QCD Predictions for Decays of P-Wave Quarkonia

abstract: Rigorous QCD predictions for decay rates of the P-wave states of heavy quarkonia are presented. They are based on a new factorization theorem which is valid to leading order in the heavy quark velocity and to all orders in the running coupling constant of QCD. The decay rates for all four P states into light hadronic or electromagnetic final states are expressed in terms of two phenomenological parameters, whose coefficients are perturbatively calculable. Logarithms of the binding energy encountered in previous perturbative calculations of P-wave decays are factored into a phenomenological parameter that is related to the probability for the heavy quark-antiquark pair to be in a color-octet S-wave state. Applying these predictions to charmonium, we use measured decay rates for the $\chione$ and $\chitwo$ to predict the decay rates of the $\chizero$ and $h_c$.