title: Polyakov loop in different representations of SU(3) at finite temperature

abstract: We investigate the Polyakov loop in different representations of SU(3) in pure gauge at finite $T$. We discuss Casimir scaling for the Polyakov loop in the deconfined phase and test and generalize the renormalization procedure for the Polyakov loop from \cite{Kaczmarek:2002mc} to arbitrary representations. In the confined phase we extract the renormalized adjoint Polyakov loop, which is finite in the thermodynamic limit. For our numerical calculations we used the tree-level improved Symanzik action on $32^3\times 4,6,8$ lattices.


title: Accelerating Dynamical Fermion Computations using the Rational Hybrid Monte Carlo (RHMC) Algorithm with Multiple Pseudofermion Fields

abstract: There has been much recent progress in the understanding and reduction of the computational cost of the Hybrid Monte Carlo algorithm for Lattice QCD as the quark mass parameter is reduced. In this letter we present a new solution to this problem, where we represent the fermionic determinant using (n) pseudofermion fields, each with an (\nth) root kernel. We implement this within the framework of the Rational Hybrid Monte Carlo (RHMC) algorithm. We compare this algorithm with other recent methods in this area and find it is competitive with them.


title: Testing Algorithms for Finite Temperature Lattice QCD

abstract: We discuss recent algorithmic improvements in simulating finite temperature QCD on a lattice. In particular, the Rational Hybrid Monte Carlo(RHMC) algorithm is employed to generate lattice configurations for 2+1 flavor QCD. Unlike the Hybrid R Algorithm, RHMC is reversible, admitting a Metropolis accept/reject step that eliminates the $\mathcal{O}(\delta t^2)$ errors inherent in the R Algorithm. We also employ several algorithmic speed-ups, including multiple time scales, the use of a more efficient numerical integrator, and Hasenbusch pre-conditioning of the fermion force.


title: The RHMC algorithm for theories with unknown spectral bounds

abstract: The Rational Hybrid Monte Carlo (RHMC) algorithm extends the Hybrid Monte Carlo algorithm for lattice QCD simulations to situations involving fractional powers of the determinant of the quadratic Dirac operator. This avoids the updating increment ($dt$) dependence of observables which plagues the Hybrid Molecular-dynamics (HMD) method. The RHMC algorithm uses rational approximations to fractional powers of the quadratic Dirac operator. Such approximations are only available when positive upper and lower bounds to the operator's spectrum are known. We apply the RHMC algorithm to simulations of 2 theories for which a positive lower spectral bound is unknown: lattice QCD with staggered quarks at finite isospin chemical potential and lattice QCD with massless staggered quarks and chiral 4-fermion interactions ($\chi$QCD). A choice of lower bound is made in each case, and the properties of the RHMC simulations these define are studied. Justification of our choices of lower bounds is made by comparing measurements with those from HMD simulations, and by comparing different choices of lower bounds.


title: Precision lattice QCD calculations and predictions of fundamental physics in heavy quark systems

abstract: I describe the recent success in performing accurate calculations of the effects of the strong force on particles containing bottom and charm quarks. Since quarks are never seen in isolation, and so cannot be studied directly, numerical simulations are key to understanding the properties of these particles and extracting information about the quarks. The results have direct impact on the worldwide experimental programme that is aiming to determine the parameters of the Standard Model of particle physics precisely and thereby uncover or constrain the possibilities for physics beyond the Standard Model. The numerical simulation of the strong force is a huge computational task and the recent success is the result of international collaboration in developing techniques that are fast enough to do the calculations on powerful supercomputers.


title: Quark condensate in two-flavor QCD

abstract: We compute the condensate in QCD with two flavors of dynamical fermions using numerical simulation. The simulations use overlap fermions, and the condensate is extracted by fitting the distribution of low lying eigenvalues of the Dirac operator in sectors of fixed topological charge to the predictions of Random Matrix Theory.


title: Diseases with rooted staggered quarks

abstract: Calculations using staggered quarks augmented with a root of the fermion determinant to reduce doubling give a qualitatively incorrect behavior in the small quark mass region. Attempts to circumvent this problem for the continuum limit involve an unproven combination of unphysical states, a loss of unitarity, and a rather peculiar non-commutation of limits.


title: The pion form factor from lattice QCD with two dynamical flavours

abstract: We compute the electromagnetic form factor of the pion using non-perturbatively O(a) improved Wilson fermions. The calculations are done for a wide range of pion masses and lattice spacings. We check for finite size effects by repeating some of the measurements on smaller lattices. The large number of lattice parameters we use allows us to extrapolate to the physical point. For the square of the charge radius we find <r^2>=0.441(20) fm^2, in good agreement with experiment.


title: The critical end point in QCD

abstract: In this talk I present the logic behind, and examine the reliability of, estimates of the critical end point (CEP) of QCD using the Taylor expansion method.


title: Strong-Isospin Breaking in the Neutron-Proton Mass Difference

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


title: String-like behaviour of 4d SU(3) Yang-Mills flux tubes

abstract: We present here results on the fine structure of the static q\bar q potential in d=4 SU(3) Yang-Mills theory. The potential is obtained from Polyakov loop correlators having separations between 0.3 and 1.2 fermi. Measurements were carried out on lattices of spatial extents of about 4 and 5.4 fermi. The temporal extent was 5.4 fermi in both cases. The results are analyzed in terms of the force between a q\bar q pair as well as in terms of a scaled second derivative of the potential. The data is accurate enough to distinguish between different effective string models and it seems to favour the expression for ground state energy of a Nambu-Goto string.


title: Hadron Deformation and Form Factors from Lattice QCD

abstract: We review the current status of lattice QCD studies of the nucleon system. In particular, we focus on the determination of the shape of the nucleon by probing its wave function as well as by evaluating the N to Delta transition form factors.


title: A lattice study of the masses of singlet 0++ mesons

abstract: We compute the masses of the flavour singlet 0++ mesons using (n_f=2) unquenched lattice QCD with the Iwasaki and Wilson gauge actions. Both fermionic and glueball interpolating operators are used to create the states. The mass of the lightest 0++ meson is suppressed relative to the mass of the 0++ glueball in quenched QCD at an equivalent lattice spacing. We discuss two possible physical reasons for this.


title: Probing Nucleon Structure on the Lattice

abstract: The QCDSF/UKQCD collaboration has an ongoing program to calculate nucleon matrix elements with two flavours of dynamical O(a) improved Wilson fermions. Here we present recent results on the electromagnetic form factors, the quark momentum fraction <x> and the first three moments of the nucleon's spin-averaged and spin-dependent generalised parton distributions, including preliminary results with pion masses as low as 320 MeV.


title: Leptonic widths of heavy quarkonia: QCD/NRQCD matching for the electromagnetic current at O(\alpha_s v^2)

abstract: We construct the S-wave part of the electromagnetic vector annihilation current to $O(\alpha_s v^2)$, where $v$ is the non-relativistic quark velocity, for heavy quarks whose dynamics are described by the NRQCD action on the lattice. The NRQCD vector current for $Q\bar{Q}$ annihilation is expressed as a linear combination of lattice operators with quantum numbers L=0, $J^P=1^-$, and the coefficients are determined by matching to the corresponding continuum current in QCD to$O(v^2)$ at one-loop. The annihilation channel gives a complex amplitude with Coulomb-exchange and infrared singularities, making a careful choice for the contours of integration and infrared subtraction functions in the numerical integration necessary. An automated vertex generation program written in Python is employed, allowing us to use a realistic NRQCD action and an improved gluon lattice action; a change in the definition of either action is easily accommodated in this procedure. The final result is applicable to simulations of electromagnetic decays of heavy quarkonia, notably the $\Upsilon$ meson.


title: Dynamical overlap fermion simulations with a preconditioned Hybrid Monte Carlo force

abstract: We present simulation results for the 2-flavour Schwinger model with dynamical Ginsparg-Wilson fermions. Our Dirac operator is constructed by inserting an approximately chiral hypercube operator into the overlap formula, which yields the overlap hypercube operator. Due to the similarity with the hypercubic kernel, a low polynomial of this kernel can be used as a numerically cheap way to evaluate the fermionic part of the Hybrid Monte Carlo force. We verify algorithmic requirements like area conservation and reversibility, and we discuss the viability of this approach in view of the acceptance rate. Next we confirm a high level of locality for this formulation. Finally we evaluate the chiral condensate at light fermion masses, based on the density of low lying Dirac eigenvalues in different topological sectors. The results represent one of the first measurements with dynamical overlap fermions, and they agree very well with analytic predictions at weak coupling.


title: Evidence for diquarks in lattice QCD

abstract: Diquarks may play an important role in hadron spectroscopy, baryon decays and color superconductivity. We investigate the existence of diquark correlations in lattice QCD by considering systematically all the lowest energy diquark channels in a color gauge-invariant setup. We measure mass differences between the various channels and show that the positive parity scalar diquark is the lightest. Quark-quark correlations inside the diquark are clearly seen in this channel, and yield a diquark size of order 1 fm.


title: The non-perturbative part of the plaquette in quenched QCD

abstract: We define the non-perturbative part of a quantity as the difference between its numerical value and the perturbative series truncated by dropping the order of minimal contribution and the higher orders. For the anharmonic oscillator, the double-well potential and the single plaquette gauge theory, the non-perturbative part can be parametrized as A (lambda)^B exp{-C/lambda} and the coefficients can be calculated analytically. For lattice QCD in the quenched approximation, the perturbative series for the average plaquette is dominated at low order by a singularity in the complex coupling plane and the asymptotic behavior can only be reached by using extrapolations of the existing series. We discuss two extrapolations that provide a consistent description of the series up to order 20-25. These extrapolations favor the idea that the non-perturbative part scales like (a/r_0)^4 with a/r_0 defined with the force method. We discuss the large uncertainties associated with this statement. We propose a parametrization of ln((a/r_0)) as the two-loop universal terms plus a constant and exponential corrections. These corrections are consistent with a_{1-loop}^2 and play an important role when beta<6. We briefly discuss the possibility of calculating them semi-classically at large beta.


title: Critical Exponents for U(1) Lattice Gauge Theory at Finite Temperature

abstract: For compact U(1) lattice gauge theory (LGT) we have performed a finite size scaling analysis on $N_{\tau} N_s^3$ lattices for $N_{\tau}$ fixed and $N_s\to\infty$, approaching the phase transition from the confined phase. For $N_{\tau}=4$, 5 and 6 our data contradict the expected scenario that this transition is either first order or in the universality class of the 3d XY model. If there are no conceptional flaws in applying the argument that the Gaussian fixed point in 3d is unstable to our systems, estimates of the critical exponents $\alpha/\nu$, $\gamma/\nu$, $(1-\beta)/\nu$ and $2-\eta$ indicate the existence of a new, non-trivial renormalization group fixed point for second order phase transitions in 3d. Such a fixed point would be of importance for renormalization group theory and statistical physics.


title: Sum Rules and Cutoff Effects in Wilson Lattice QCD

abstract: We use the transfer matrix formalism to derive non-perturbative sum rules in Wilson's lattice QCD with N_f flavours of quarks. The discretization errors on these identities are treated in detail. As an application, it is shown how the sum rules can be exploited to give improved estimates of the continuum spectrum and static potential.


title: Hadronic decays from the lattice

abstract: I review the lattice QCD approach to determining hadronic decay transitions. Examples considered include rho to pi pi; b_1 to pi omega; hybrid meson decays and scalar meson decays. I discuss what lattices can provide to help understand the composition of hadrons.


title: The Kaon B-parameter in Mixed Action Chiral Perturbation Theory

abstract: We calculate the kaon B-parameter, B_K, in chiral perturbation theory for a partially quenched, mixed action theory with Ginsparg-Wilson valence quarks and staggered sea quarks. We find that the resulting expression is similar to that in the continuum, and in fact has only two additional unknown parameters. At one-loop order, taste-symmetry violations in the staggered sea sector only contribute to flavor-disconnected diagrams by generating an O(a^2) shift to the masses of taste-singlet sea-sea mesons. Lattice discretization errors also give rise to an analytic term which shifts the tree-level value of B_K by an amount of O(a^2). This term, however, is not strictly due to taste-breaking, and is therefore also present in the expression for B_K for pure G-W lattice fermions. We also present a numerical study of the mixed B_K expression in order to demonstrate that both discretization errors and finite volume effects are small and under control on the MILC improved staggered lattices.


title: Mass spectrum of 1-+ exotic mesons from lattice QCD

abstract: Time correlation functions of a hybrid exotic meson operator, with JPC=1-+, generated in quenched lattice QCD are subjected to a (Bayesian) maximum entropy analysis. Five distinct spectral levels are uncovered. Their extrapolation into the physical pion mass region suggests a possible relationship to experimentally known states pi1(1400) and pi1(1600), and also to a state in the 2 GeV region carrying the same quantum numbers.


title: Improvement of quark propagator estimation through domain decomposition

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


title: Status Report on ILDG activities

abstract: A status report about the International Lattice Data Grid (ILDG) is given. Different countries participating in the ILDG have created regional lattice data grid solutions that are implemented, working and used. The remaining task and the focus of present activities is the development of the interoperability of these regional grids. A first, successful step in this direction is a metadata catalogue service which is already working interoperable.


title: Quarkonium Spectral Function from Anisotropic Lattice

abstract: We discuss the behavior of charmonia and bottomonia correlators and spectral functions above the deconfinement temperature and determine melting temperatures for different mesonic states.


title: Excited mesons from $N_f=2$ dynamical Clover Wilson lattices

abstract: We study mesons on the lattice with a special focus on excited states. For that purpose we construct several quark sources with different spatial smearings, including p-waves. These quark sources are then combined with the appropiate Dirac structures to form meson interpolators of definite spin. We use these operators to construct a cross correlation matrix from which we extract ground and excited meson states using the variational method. For the calculations we use gauge configurations with $N_f=2$ dynamical Clover Wilson fermions provided by the CP-PACS collaboration. We show preliminary results for pseudoscalar, scalar, vector and pseudovector mesons.


title: Plaquette expectation value and lattice free energy of three-dimensional SU(N) gauge theory

abstract: We use high precision lattice simulations to calculate the plaquette expectation value in three-dimensional SU(N) gauge theory for N=2,3,4,5,8. Using these results, we study the N-dependence of the first non-perturbative coefficient in the weak-coupling expansion of hot QCD. We demonstrate that, in the limit of large N, the functional form of the plaquette expectation value with ultraviolet divergences subtracted is 15.9(2)-44(2)/N^2.


title: The infrared behavior of lattice QCD Green's functions

abstract: We investigate different aspects of lattice QCD in Landau gauge using Monte Carlo simulations. In particular, we focus on the low momentum behavior of gluon and ghost propagators. The gauge group is SU(3). Different systematic effects on the gluon and ghost propagators are studied, e.g. the dependence on the choice of Gribov copies or the influence of dynamical Wilson fermions. We compare our data with results from studies of Dyson-Schwinger equations for the gluon and ghost propagators. We demonstrate that the infrared behavior of both propagators, as found in this thesis, is consistent with different criteria for confinement. However, the running coupling constant, given as a renormalization-group-invariant combination of the gluon and ghost dressing functions, does not expose a finite infrared fixed point. We also report on a first nonperturbative computation of the SU(3) ghost-gluon-vertex renormalization constant and on an investigation of the spectral properties of the Faddeev-Popov operator.


title: Topology, chiral and screening transitions at finite density in two colour QCD

abstract: The behaviour of the topological susceptibility in QCD with two colours and 8 flavours of quarks is studied at nonzero temperature on the lattice across the finite density transition. It is shown that its signal drops at a (pseudo-)critical chemical potential mu_c. The Polyakov loop and the chiral condensate undergo their transitions at the same value. Pauli blocking supervenes at a value of the chemical potential larger than mu_c.


title: Direct numerical computation of disorder parameters

abstract: In the framework of various statistical models as well as of mechanisms for color confinement, disorder parameters can be developed which are generally expressed as ratios of partition functions and whose numerical determination is usually challenging. We develop an efficient method for their computation and apply it to the study of dual superconductivity in 4d compact U(1) gauge theory.


title: Quantum Operator Design for Lattice Baryon Spectroscopy

abstract: A previously-proposed method of constructing spatially-extended gauge-invariant three-quark operators for use in Monte Carlo lattice QCD calculations is tested, and a methodology for using these operators to extract the energies of a large number of baryon states is developed. This work is part of a long-term project undertaken by the Lattice Hadron Physics Collaboration to carry out a first-principles calculation of the low-lying spectrum of QCD. These techniques are then applied in the construction of nucleon operators. Correlation matrix elements between these operators are estimated using 200 configurations on a $12^3 \times 48$ anisotropic lattice in the quenched approximation with unphysically heavy u, d quark masses (the pion mass is approximately 700 MeV). After a change of basis operators using a variational method is applied, the energies of up to eight states are extracted in each symmetry channel. Although comparison with experiment is not justified, the pattern of levels obtained qualitatively agrees with the observed spectrum. A comparison with quark model predictions is also made; the quark model predicts more low-lying even-parity states than this study yields, but both the quark model and this study predict more odd-parity states near 2 GeV than currently observed in experiments.


title: Nonzero chemical potential in the overlap Dirac operator and comparison to random matrix theory

abstract: In this talk we present the results published recently in Ref. [1], where we showed how to introduce a quark chemical potential in the overlap Dirac operator. The resulting operator satisfies a Ginsparg-Wilson relation and has exact zero modes. It is no longer gamma_5-Hermitian, but its nonreal eigenvalues still occur in pairs. We compute the spectral density of the operator on the lattice and show that, for small eigenvalues, the data agree with analytical predictions of non-Hermitian chiral random matrix theory for both trivial and nontrivial topology. We also explain an observed change in the number of zero modes as a function of chemical potential.


title: Trajectory Length and Autocorrelation Times -- Nf=2 simulations in the Schroedinger functional

abstract: A status report is presented on the large-volume simulations in the Schroedinger functional with two flavours of O(a) improved Wilson quarks performed by the ALPHA collaboration. The physics goal is to set the scale for the computation of the fundamental parameters of QCD. In this talk the emphasis is on aspects of the Hybrid Monte-Carlo algorithm, which we use with (symmetric) even-odd and Hasenbusch preconditioning. We study the dependence of aucorrelation times on the trajectory length. The latter is found to be significant for fermionic correlators, the trajectories longer than unity performing better than the shorter ones.


title: Gross-Neveu model as a laboratory for fermion discretization

abstract: We introduce a finite volume renormalization scheme for the N-Majorana-component O(N) invariant Gross-Neveu model. Universal observables are defined that are accessible to precise numerical simulation in various discretizations and allow for an extrapolation to the continuum limit. Here first numerical results with Wilson fermions are reported. For N=2 they reproduce exact finite volume continuum results in the massless Thirring model. Our N=8 data are ready for comparison for instance with staggered results in the future.


title: Iterative methods for overlap and twisted mass fermions

abstract: We present a comparison of a number of iterative solvers of linear systems of equations for obtaining the fermion propagator in lattice QCD. In particular, we consider chirally invariant overlap and chirally improved Wilson (maximally) twisted mass fermions. The comparison of both formulations of lattice QCD is performed at four fixed values of the pion mass between 230MeV and 720MeV. For overlap fermions we address adaptive precision and low mode preconditioning while for twisted mass fermions we discuss even/odd preconditioning. Taking the best available algorithms in each case we find that calculations with the overlap operator are by a factor of 30-120 more expensive than with the twisted mass operator.


title: Correlation of the ghost and the quark in the lattice Landau gauge QCD

abstract: Effects of the quark field on the ghost propagator of the lattice Landau gauge are investigated by using the quenched gauge configuration of SU(2) first copy and the parallel tempering (PT) gauge fixed samples, quenched SU(3) $56^4$ configuration and unquenched SU(3) configurations produced by the MILC collaboration. We measure color symmetric and color antisymmetric ghost propagator and the Binder cumulant of $l^1$ norm and $l^2$ norm of color antisymmetric ghost propagator and investigate deviation from that of Gaussian distribution. The dynamical quarks make color antisymmetric ghost propagator closer to Gaussian distribution and Kugo-Ojima color confinement parameter $c$ closer to 1.


title: Rotational Symmetry and A Light Mode in Two-Dimensional Staggered Fermions

abstract: To obtain a light mode in two-dimensional staggered fermions, we introduce four new local operators keeping the rotational invariance for a staggered Dirac operator. To split masses of tastes, three cases are considered. The mass matrix and the propagator for free theories are analyzed. We find that one of three cases is a good candidate for taking a single mode by the mass splitting. In the case, it is possible that a heavy mode obtains infinite mass on even sites or odd sites.


title: Perturbative chiral violations for domain-wall QCD with improved gauge actions

abstract: We investigate, in the framework of perturbation theory at finite $N_s$, the effectiveness of improved gauge actions in suppressing the chiral violations of domain-wall fermions. Our calculations show substantial reductions of the residual mass when it is compared at the same value of the gauge coupling, the largest suppression being obtained when the DBW2 action is used. Similar effects can also be observed for a power-divergent mixing coefficient which is chirally suppressed. No significant reduction instead can be seen in the case of the difference between the vector and axial-vector renormalization constants when improved gauge actions are used in place of the plaquette action. We also find that one-loop perturbation theory is not an adequate tool to carry out comparisons at the same energy scale (of about 2 GeV), and in fact in this case even an enhancement of the chiral violations is frequently obtained.


title: Infrared Gluon and Ghost Propagators from Lattice QCD. Results from large asymmetric lattices

abstract: We report on the infrared limit of the quenched lattice Landau gauge gluon and ghost propagators as well as the strong coupling constant computed from large asymmetric lattices. The infrared lattice propagators are compared with the pure power law solutions from Dyson-Schwinger equations (DSE). For the gluon propagator, the lattice data is compatible with the DSE solution. The preferred measured gluon exponent being $\sim 0.52$, favouring a null zero momentum propagator. The lattice ghost propagator shows finite volume effects and, for the volumes considered, the propagator does not follow a pure power law. Furthermore, the strong coupling constant is computed and its infrared behaviour investigated.


title: Investigation of the overlap of excited bottomonium states with hybrid operators

abstract: We analyze the overlap of color-octet meson operators with the $\Upsilon$ and the $\eta_b$ and their 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: Dynamical simulations of QCD at finite temperature with a truncated perfect action

abstract: The Hypercube operator determines a variant of the approximate, truncated perfect fermion action. In this pilot study we are going to report on first experiences in dynamical QCD simulations with the Hypercube fermions. We apply this formulation in an investigation of the finite temperature transition for two flavours. On lattices of size $8^3\times 4$ we explore the phase diagram. Physical scales are estimated from pseudoscalar and vector meson masses obtained on $8^3\times 16$ lattices. We observe the presence of a metastability region but do not find evidence for an Aoki phase. The Hypercube operator allows us to simulate at ratios of pseudoscalar to vector meson masses at least as small as 0.8 at the thermal crossover at $N_t=4$, which renders this formulation cheaper than the Wilson like fermions.


title: The rooted staggered determinant in the Schwinger model

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


title: The one-flavor quark condensate and related problems

abstract: I describe a recent calculation (by me, Hoffmann, Liu, and Schaefer) of the chiral condensate in one-flavor QCD using numerical simulations with overlap fermions. The condensate is extracted by fitting the distribution of low lying eigenmodes of the Dirac operator in sectors of fixed topological charge to the predictions of Random Matrix Theory. Our results are in excellent agreement with estimates from the orientifold large-N expansion. Much interesting physics surrounds this calculation, which I will highlight.


title: Thermodynamics of (2+1)-flavor QCD

abstract: We report on the status of our QCD thermodynamics project. It is performed on the QCDOC machine at Brookhaven National Laboratory and the APEnext machine at Bielefeld University. Using a 2+1 flavor formulation of QCD at almost realistic quark masses we calculated several thermodynamical quantities. In this proceeding we show the susceptibilites of the chiral condensate and the Polyakov loop, the static quark potential and the spatial string tension.


title: Simulation of Lattice Gauge Action from the Overlap Operator

abstract: We show the result of the lattice gauge field tensor which is derived from the classical continuum limit of the overlap Dirac operator. By analogous construction, it was recently proposed that the gauge action can be obtained from the overlap operator as well which is proportional to Tr a D_{ov}. We discuss how to carry out Monte Carlo simulations of such a gauge action together with the overlap fermion.


title: On the Locality and Scaling of Overlap Fermions at Coarse Lattice Spacings

abstract: The overlap fermion offers the considerable advantage of exact chiral symmetry on the lattice, but is numerically intensive. This can be made affordable while still providing large lattice volumes, by using coarse lattice spacing, given that good scaling and localization properties are established. Here, using overlap fermions on quenched Iwasaki gauge configurations, we demonstrate directly that, with appropriate choice of negative Wilson's mass, the overlap Dirac operator's range is comfortably small in lattice units for each of the lattice spacings 0.20 fm, 0.17 fm, and 0.13 fm (and scales to zero in physical units in the continuum limit). In particular, our direct results contradict recent speculation that an inverse lattice spacing of 1 GeV is too low to have satisfactory localization. Furthermore, hadronic masses (available on the two coarser lattices) scale very well.


title: Lattice Results for the QCD Phase Transition

abstract: We discuss recent results for the phase transition in finite-temperature QCD from numerical (Monte Carlo) simulations of the lattice-regularized theory. Emphasis is given to the case of two degenerate light-quark flavors. The order of the transition in this case, which could have cosmological implications, has not yet been established.


title: Exploring the infrared gluon and ghost propagators using large asymmetric lattices

abstract: We report on the infrared limit of the quenched lattice Landau gauge gluon propagator computed from large asymmetric lattices. In particular, the compatibility of the pure power law infrared solution $(q^2)^{2\kappa}$ of the Dyson-Schwinger equations is investigated and the exponent $\kappa$ is measured. The lattice data favours $\kappa \sim 0.52$, which would imply a vanishing zero momentum gluon propagator as predicted by the Kugo-Ojima confinement mechanism and the Zwanziger horizon condition. Results for the ghost propagator and for the running coupling constant are shown.


title: Avoiding the boson sign problem at finite chemical potential

abstract: The usual path integral formulation for scalar particles at finite density involves a sign problem, making numerical simulation impractical. We present alternative methods free of this difficulty. We apply these approaches to phi^4 theory in 2+1 dimensions, presenting preliminary numerical results for physical quantities in the vicinity of the phase transition in the T-mu plane.


title: Baryon correlators containing different diquarks from lattice simulations

abstract: Point to point vacuum correlators containing diquarks in the color anti-triplet representation are computed both in the quenched approximation and dynamical overlap simulations with two flavors. The scalar, pseudoscalar and axial vector diquarks are combined with light quarks to form color singlets. The scalar ("good") diquark channel shows a stronger attraction than the axial vector ("bad") channel in the quenched data set. The pseudoscalar diquark channel shows a finite volume zero mode artifact: the correlator becomes negative at large distance when the quark mass is small. By separating configurations without zero modes from those with zero modes, we found that the zero modes have an important contribution in both the attraction in the scalar channel and the repulsion in the pseudoscalar channel. In the axial vector diquark channel, we did not find apparent zero mode effects.


title: Overlap fermions on a twisted mass sea

abstract: We present first results of a mixed action project. We analyze gauge configurations generated with two flavors of dynamical twisted mass fermions. Neuberger's overlap Dirac operator is used for the valence sector. The various choices in the setup of the simulation are discussed. We employ chiral perturbation theory to describe the effects of using different actions in the sea and valence sector at non-zero lattice spacing.


title: Lattice QCD at finite temperature

abstract: I discuss recent developments in finite temperature lattice QCD, including the calculation of the transition temperature, equation of state, color screening and meson spectral functions.


title: Searching for the elusive critical endpoint at finite temperature and isospin density

abstract: We consider 3-flavour lattice QCD with a finite chemical potential mu_I for isospin, close to the finite temperature transition from hadronic matter to a quark-gluon plasma. In this region one can argue that the position and probably the nature of this transition mimic those at finite quark-number chemical potential mu. The quark mass is chosen to be close to the critical mass at zero chemical potentials. Since the Binder cumulants used to determine the nature of this transition in HMD(R) simulations are very sensitive to the updating increment dt, we have switched to the newer exact RHMC algorithm for our simulations. Preliminary results indicate that there is no critical endpoint in the small mu_I regime, at least none connected with the critical point at zero chemical potentials.


title: Multi-Step stochastic correction in dynamical fermion updating algorithms

abstract: The advantages of using Multi-Step corrections for simulations of lattice gauge theories with dynamical fermions will be discussed. This technique is suited for algorithms based on the Multi-Boson representation of the dynamical fermions as well as for the Hybrid Monte-Carlo (HMC) algorithm and variants of the latter, like the Polynomial-HMC. Especially the latter has the power to deal with an odd number of fermion fields--an essential feature necessary for realistic QCD-simulations with up-, down-, and strange-quarks. In particular, we will discuss the application of the multi-step (actually two-step) correction technique to a PHMC updating algorithm for twisted-mass Wilson fermions with non-degenerate fermion masses, as it was used in recent dynamical simulations for N_f=2+1+1 fermion flavors.


title: Landau Gauge Gluon and Ghost Propagators from Lattice QCD

abstract: We report on recent numerical computations of the Landau gauge gluon and ghost propagators as well as of the ghost-gluon vertex function in pure SU(3) Yang-Mills theory and in full QCD on the lattice. Special emphasis is paid to the low momentum region. In particular, we present new data for the gluon propagator at momenta below 300 MeV. We also discuss different systematic effects as there are finite-size, lattice discretization and Gribov copy but also unquenching effects. A MOM-scheme running coupling \alpha_s(q^2) based on the ghost-gluon vertex is calculated and found to decrease for momenta below 550 MeV, even though the renormalization constant of the vertex deviates only weakly from being constant.


title: Exploring the chiral singularities of two color lattice QCD at strong coupling

abstract: We use a modified directed path algorithm to study the finite temperature chiral singularities of two color lattice QCD with staggered fermions at strong coupling. Our lattice calculations are done at a fixed finite temperature in the broken phase with a variety of different quark masses. We find that to the lowest order, the behavior of our observables, namely the condensates and the decay constants are all consistent with the predictions of the low energy effective theory of our model. We further notice that in order to see consistency, the quark masses used in our lattice calculations need to be quite small than typically used.


title: Lattice Gauge Theory Approach to Spontaneous Symmetry Breaking from an Extra Dimension

abstract: We present lattice simulation results corresponding to an SU(2) pure gauge theory defined on the orbifold space E_4 x I_1, where E_4 is the four-dimensional Euclidean space and I_1 is an interval, with the gauge symmetry broken to a U(1) subgroup at the two ends of the interval by appropriate boundary conditions. We demonstrate that the U(1) gauge boson acquires a mass from a Higgs mechanism. The mechanism is driven by two of the extra-dimensional components of the five-dimensional gauge field which play respectively the role of the longitudinal component of the gauge boson and a massive real physical scalar, the Higgs particle. Despite the non-renormalizable nature of the theory, we observe only a mild cut-off dependence of the physical observables. We also show evidence that there is a region in the parameter space where the system behaves in a way consistent with dimensional reduction.


title: GW Meson Scattering on a Staggered Sea

abstract: We discuss the structure of the NLO corrections to the chiral formulae for mesonic scattering processes in mixed action simulations using Ginsparg-Wilson valence quarks and staggered sea quarks. In particular, we show that the analytic contribution of the NLO chiral Lagrangian is the same as in QCD. We also comment on how this result restricts the dependence of the amplitudes on the unknown parameter C_Mix appearing in the chiral theory appropriate for these systems. We conclude with some comments on the explicit scattering lengths.


title: RG decimation study of SU(2) gauge theory

abstract: We report on numerical studies of RG decimations in SU(2) gauge theory. We study in particular a class of plaquette actions involving sums of group representations. We measure a number of observables representative of different length scales in order to investigate the transformation of the system under different choices of spin blocking, and examine the flow of the effective action couplings. A need for a projection to some class of effective actions on the configurations resulting from the adopted numerical decimation procedures is demonstrated. A numerical decimation procedure resulting in the effective single plaquette Lagrangian tailored to reproduce different medium/large scale physics is devised.


title: Renormalization of infrared contributions to the QCD pressure

abstract: Thanks to dimensional reduction, the infrared contributions to the QCD pressure can be obtained from two different three-dimensional effective field theories, called the Electrostatic QCD (Yang-Mills plus adjoint Higgs) and the Magnetostatic QCD (pure Yang-Mills theory). Lattice measurements have been carried out within these theories, but a proper interpretation of the results requires renormalization, and in some cases also improvement, i.e. the removal of terms of O(a) or O(a^2). We discuss how these computations can be implemented and carried out up to 4-loop level with the help of Numerical Stochastic Perturbation Theory.


title: Monopole condensation in two-flavour Adjoint QCD

abstract: Two distinct phase transitions occur at different temperatures in QCD with adjoint fermions (aQCD): deconfinement and chiral symmetry restoration. In this model, quarks do no explicitely break the center Z(3) symmetry and therefore the Polyakov loop is a good order parameter for the deconfinement transition. We study monopole condensation by inspecting the expectation value of an operator which creates a monopole. Such a quantity is expected to be an order parameter for the deconfinement transition as in the case of fundamental fermions.


title: Color Screening, Casimir Scaling, and Domain Structure in G(2) and SU(N) Gauge Theories

abstract: We argue that screening of higher-representation color charges by gluons implies a domain structure in the vacuum state of non-abelian gauge theories, with the color magnetic flux in each domain quantized in units corresponding to the gauge group center. Casimir scaling of string tensions at intermediate distances results from random spatial variations in the color magnetic flux within each domain. The exceptional G(2) gauge group is an example rather than an exception to this picture, although for G(2) there is only one type of vacuum domain, corresponding to the single element of the gauge group center. We present some numerical results for G(2) intermediate string tensions and Polyakov lines, as well as results for certain gauge-dependent projected quantities. In this context, we discuss critically the idea of projecting link variables to a subgroup of the gauge group. It is argued that such projections are useful only when the representation-dependence of the string tension, at some distance scale, is given by the representation of the subgroup.


title: Why (staggered fermions)^{1/4} fail at finite density

abstract: Because the staggered fermion determinant is complex at nonzero mu, taking its fourth root leads to phase ambiguities. These unphysical effects cause the measure to become discontinuous; the problem becomes acute when Re mu exceeds approximately half the pion mass (when T>0 this rough bound probably moves towards larger mu). We show how to overcome the problem, but only very close to the continuum limit. This regime may be beyond reach with current resources.


title: Lattice QCD determination of states with spin 5/2 or higher in the spectrum of nucleons

abstract: Energies for excited isospin 1/2 states that include the nucleon are computed using quenched, anisotropic lattices. Baryon interpolating field operators that are used include nonlocal operators that provide $G_2$ irreducible representations of the octahedral group. The decomposition of spin 5/2 or higher states is realized for the first time in a lattice QCD calculation. We observe patterns of degenerate energies in the irreducible representations of the octahedral group that correspond to the subduction of the continuum spin 5/2 or higher.


title: Update on the physics of light pseudoscalar mesons

abstract: We present an update of the MILC investigation of the properties of light pseudoscalar mesons using three flavors of improved staggered quarks. Results are presented for the $\pi$ and $K$ leptonic decay constants, the CKM matrix element $V_{us}$, the up, down and strange quark masses, and the coefficients of the $O(p^4)$ chiral lagrangian. We have new data for lattice spacing $a \approx 0.15$ fm with several values of the light quark mass down to one-tenth the strange quark mass, higher statistics for $a \approx 0.09$ fm with the light quark mass equal to one-tenth the strange quark mass, and initial results for our smallest lattice spacing, $a \approx 0.06$ fm with light quark mass two-fifths of the strange quark mass.


title: Field theoretic study of a cold Fermi gas in the unitary limit

abstract: Trapped and cooled gases of alkali atoms can be manipulated to exhibit a variety of interesting phenomena. For example, dilute gases of fermionic atoms, in 2 hyperfine states, can be cooled to temperatures where they become superfluid. An external field can be applied to tune the scattering length a. When |a| exceeds the interparticle spacing, nonperturbative tools are needed to study the system theoretically. The unitary limit, |a|\to\infty, is particularly interesting due to its universality and symmetry. Lattice field theory and effective field theory can be used to systematically calculate properties of this system. Results are presented for the finite temperature phase transition and for behavior near zero temperature.


title: Quantum fluctuations of k-strings: A case study

abstract: K strings in Yang-Mills theory can be considered as bound states of k elementary confining strings carrying one unit of colour flux. Current estimates of k-string tension sigma_k are very sensitive to the leading corrections due to quantum fluctuations of the string. In this study we address this problem by comparing Polyakov-Polyakov correlators in the fundamental representation (k=1) with the corresponding ones with k=2 in the confining phase of a Z_4 gauge theory in three dimensions. Highly efficient simulation techniques are available in this case. Although the k=1 Polyakov-Polyakov correlator matches nicely with the expected bosonic string effects up to the Next-to-Leading-Order, the k=2 Polyakov-Polyakov correlators show large deviations. This is an important source of potential systematic errors in the current estimates of sigma_k.


title: An Analysis of the Hadronic Spectrum from Lattice QCD

abstract: I begin by discussing the basic ideas of quantum field theory (QFT). I provide a review of symmetries in physics and then move on to discuss the quark model. I then review lattice gauge theory with particular attention paid to lattice QCD and some of the associated problems. I discuss gauge fields on the lattice along with free lattice fermions. I follow this with a definition of the lattice QCD action and discuss how to reproduce the correct continuum physics. I include a discussion of the basic numerical techniques employed in lattice simulations. I review methods for putting particles onto the lattice and describe how to fit the resulting data. The symmetries of the QCD Lagrangian are reviewed along with various forms of symmetry breaking in physics, the PCAC relation, the Goldberger-Treiman relation and the spontaneous breakdown of the axial symmetry. I move on to discuss sigma models and finally arrive at a basic chiral perturbation theory. I present research completed with my supervisor C. Allton and collaborators A.W. Thomas, D.B. Leinweber and R. Young in the final three chapters. This work involves making lattice predictions for the Hadronic mass spectrum using extrapolation techniques based on the predictions of chiral perturbation theory which have been developed by the Adelaide group.


title: Confining properties of 2-color QCD at finite density

abstract: We study the confining properties of QCD with two colors across the finite density phase transition. A disorder parameter detecting dual superconductivity of the QCD vacuum is used as a probe for the confinement/deconfinement phase transition.


title: Lattice simulations of real-time quantum fields

abstract: We investigate lattice simulations of scalar and nonabelian gauge fields in Minkowski space-time. For SU(2) gauge-theory expectation values of link variables in 3+1 dimensions are constructed by a stochastic process in an additional (5th) ``Langevin-time''. A sufficiently small Langevin step size and the use of a tilted real-time contour leads to converging results in general. All fixed point solutions are shown to fulfil the infinite hierarchy of Dyson-Schwinger identities, however, they are not unique without further constraints. For the nonabelian gauge theory the thermal equilibrium fixed point is only approached at intermediate Langevin-times. It becomes more stable if the complex time path is deformed towards Euclidean space-time. We analyze this behavior further using the real-time evolution of a quantum anharmonic oscillator, which is alternatively solved by diagonalizing its Hamiltonian. Without further optimization stochastic quantization can give accurate descriptions if the real-time extend of the lattice is small on the scale of the inverse temperature.


title: Testing UV-filtered ("fat-link") clover fermions

abstract: We investigate filtered clover fermions, built from fat gauge links, both in one-loop perturbation theory and in numerical simulations. We use a variety of filtering recipes (APE, HYP, EXP, HEX), some of which are suitable for a HMC with dynamical fermions. A generic filtering together with a (fat-link) clover term yields fermions with much reduced chiral symmetry breaking.


title: Broadening of the confining flux tube joining higher-representation sources

abstract: In the context of confining gauge theories we study the flux tube generated by a pair of static sources belonging to higher rank representations of the gauge group. Using a simple geometric approach based on minimal surfaces describing the world-sheet of the underlying k-string, we infer logarithmic broadening of the flux tube as a function of source separation. It turns out that the coefficient of the logarithm does not depend on the specific representation of the source nor on its N-ality and is universal. Numerical results on a Z_4 gauge theory in three dimensions with static sources carrying two units of elementary flux compare very favourably with these predictions.


title: Dual Superconductivity in G2 group

abstract: We investigate the dual superconductivity mechanism in the exceptional group $G_2$. This is a centerless group (no 't Hooft flux vortices are allowed) and we check for the presence of a magnetic monopole condensate in the confined phase by measuring on the lattice a disorder parameter related to the vacuum expectation value of an operator carrying magnetic charge. The behaviour of the disorder parameter is consistent with the dual superconductor picture. A first step of an analysis on the thermodynamical properties of the theory is conducted by mean of this operator.


title: Chiral Symmetry Restoration in Anisotropic QED(3)

abstract: We present results from a Monte Carlo simulation of non-compact lattice QED in 3 dimensions in which an explicit anisotropy $\kappa$ between $x$ and $y$ hopping terms has been introduced into the action. Using a parameter set corresponding to broken chiral symmetry in the isotropic limit $\kappa=1$, we study the chiral condensate on $16^3$, $20^3$, and $24^3$ lattices as $\kappa$ is varied, and fit the data to an equation of state which incorporates anisotropic volume corrections. The value $\kappa_c$ at which chiral symmetry ispparently restored is strongly volume-dependent, suggesting that the transition may be a crossover rather than a true phase transition. In addition we present results on $16^3$ lattices for the scalar meson propagator, and for the Landau gauge-fixed fermion propagator. The scalar mass approaches the pion mass at large $\kappa$, consistent with chiral symmetry restoration, but the fermion remains massive at all values of $\kappa$ studied, suggesting that strong infra-red fluctuations persist into the chirally symmetric regime. Implications for models of high-$T_c$ superconductivity based on anisotropic QED$_3$ are discussed.


title: Algorithms for dynamical overlap fermions

abstract: An overview of the current status of algorithmic approaches to dynamical overlap fermions is given. In particular the issue of changing the topological sector is discussed.


title: Lattice Field Theory Methods in Modern Biophysics

abstract: An effective field theory exists describing a very large class of biophysically interesting Coulomb gas systems: the lowest order (mean-field) version of this theory takes the form of a generalized Poisson-Boltzmann theory. Interaction terms depend on details (finite-size effects, multipole properties, etc). Convergence of the loop expansion holds only if mutual interactions of mobile charges are small compared to their interaction with the fixed-charge environment, which is frequently not the case. Problems with the strongly- coupled effective theory can be circumvented with an alternative local lattice formulation, with real positive action. In realistic situations, with variable dielectric, a determinant of the Poisson operator must be inserted to generate correct electrostatics. Methods adopted from unquenched lattice QCD do this very efficiently.


title: m_b and f_Bs from a combination of HQET and QCD

abstract: We compute the mass of the b-quark and the B_s meson decay constant in quenched lattice QCD using a combination of HQET and the standard relativistic QCD Lagrangian. We start from a small volume, where one can directly deal with the b-quark, and compute the evolution to a big volume, where the finite size effects are negligible through step scaling functions which give the change of the observables when L is changed to 2L. In all steps we extrapolate to the continuum limit, separately in HQET and in QCD for masses below m_b. The point m_b is then reached by an interpolation of the continuum results. With r_0=0.5 fm and the experimental B_s and K masses we find f_Bs=191(6) MeV and the renormalization group invariant mass M_b=6.89(11) GeV, translating into m_b(m_b)=4.42(7) GeV in the MS-bar scheme.


title: Mass Splitting of Staggered Fermion and $SO(2D)$ Clifford Algebra

abstract: We present a new method to introduce rotationally invariant terms in staggered fermions which is based on an $SO(2D)$ Clifford algebra formulation, where $D$ means the number of space-time dimensions. We have four candidates for improved mass terms that can split the degenerate mass of staggered fermions. Among them, we analyze three types of combinations and find only one case that can identify with the light single Dirac mode.


title: Mixed Action Simulations on Staggered Background; Interpretation and Result for the 2-flavor QCD Chiral Condensate

abstract: Growing evidence indicates that in the continuum limit the rooted staggered action is in the correct QCD universality class, the non-local terms arising from taste breaking can be viewed as lattice artifacts. In this paper we consider the 2-flavor Asqtad staggered action at a lattice spacing of approximately 0.13 fm and probe the properties of the staggered configurations by an overlap valence Dirac operator. By comparing the distribution of the overlap eigenmodes to continuum QCD predictions we investigate if/when the lattice artifacts are small as a function of the staggered quark mass. We define a matching overlap quark mass where the lattice corrections are minimal for the topological susceptibility and from the eigenmode distribution we predict the 2-flavor chiral condensate. Our results indicate that the staggered configurations are consistent with 2-flavor continuum QCD up to small lattice artifacts, and predict a consistent value for the infinite volume chiral condensate.


title: The QCD transition temperature: results with physical masses in the continuum limit

abstract: The transition temperature ($T_c$) of QCD is determined by Symanzik improved gauge and stout-link improved staggered fermionic lattice simulations. We use physical masses both for the light quarks ($m_{ud}$) and for the strange quark ($m_s$). Four sets of lattice spacings ($N_t$=4,6,8 and 10) were used to carry out a continuum extrapolation. It turned out that only $N_t$=6,8 and 10 can be used for a controlled extrapolation, $N_t$=4 is out of the scaling region. Since the QCD transition is a non-singular cross-over there is no unique $T_c$. Thus, different observables lead to different numerical $T_c$ values even in the continuum and thermodynamic limit. The peak of the renormalized chiral susceptibility predicts $T_c$=151(3)(3) MeV, wheres $T_c$-s based on the strange quark number susceptibility and Polyakov loops result in 24(4) MeV and 25(4) MeV larger values, respectively. Another consequence of the cross-over is the non-vanishing width of the peaks even in the thermodynamic limit, which we also determine. These numbers are attempted to be the full result for the $T$$\neq$0 transition, though other lattice fermion formulations (e.g. Wilson) are needed to cross-check them.


title: Exploring the infrared Landau gauge propagators using large asymmetric lattices

abstract: We report on the infrared limit of the quenched lattice Landau gauge gluon and ghost propagators computed from large asymmetric lattices. In particular, the compatibility of the pure power law infrared solutions of the Dyson-Schwinger equations with the lattice data is investigated and the exponent $\kappa$ is measured. The gluon lattice data favour $\kappa \sim 0.52$, which would imply a vanishing zero momentum gluon propagator. For the subset of lattices where the ghost propagator was computed, the data is not compatible with a pure power law. Our data also show a decreasing running coupling in the infrared region. Furthermore, positivity violation for the gluon propagator is also verified.


title: The chiral condensate on 2-flavor staggered configurations from an overlap operator

abstract: We measure the low lying eigenmodes of an overlap Dirac operator on 2--flavor staggered configurations. By comparing the eigenmode distribution to the predictions of Random Matrix Theory we test to what accuracy staggered configurations describe continuum QCD. The agreement between the numerical data and RMT implies that at our quark mass values the lattice artifacts of the staggered configurations are comparable to overlap configurations. We identify the overlap valence mass that best matches the staggered sea quarks and predict the value of the infinite volume 2-flavor chiral condensate.


title: Toward Group Theory Operators for tmLQCD Hadrons

abstract: Extraction of the mass spectrum from twisted mass lattice QCD is facilitated by respecting the octahedral group of rotations and accommodating the broken parity and flavor symmetries of the theory. In this work, group theory meson operators adapted to these constraints are constructed for the special case of quark and antiquark fields at a common lattice site, connected by extended gauge field paths.


title: Light baryon spectrum using improved interpolating operators

abstract: Energies for excited light baryons are computed in quenched QCD with a pion mass of 490 MeV. Operators used in the simulations include local operators and the simplest nonlocal operators that have nontrivial orbital structures. All operators are designed with the use of Clebsch-Gordan coefficients of the octahedral group so that they transform irreducibly under the group rotations. Matrices of correlation functions are computed for each irreducible representation, and then the variational method is applied to separate mass eigenstates. We obtained 17 states for isospin 1/2 and 11 states for isospin 3/2 in various spin-parity channels including $J^P=5/2^\pm$. The pattern of the lowest-lying energies from each irrep is discussed. We use anisotropic lattices of volume $24^3\times 64$ with temporal lattice spacing $a_t^{-1}=6.05$ GeV with renormalized anisotropy $\xi=3.0$.


title: Mobility edge and locality of the overlap-Dirac operator with and without dynamical overlap fermions

abstract: We perform a systematic study of low-lying eigenmodes of $H_w$ with various gauge actions to find the optimal choice for dynamical overlap fermion simulations, with which one may achieve lower numerical cost for HMC and better locality property of the overlap kernel. For this purpose, our study is made with emphasis on the distribution of low-lying eigenvalues and the mobility edge with and without dynamical overlap fermions.


title: Towards continuum limit of screening lengths with chiral Fermions

abstract: We investigate mesonic screening correlators at T=2T_c using the overlap Fermions in the quenched approximation, where T_c is the QCD phase transition temperature. Using lattices with temporal extent up to 8, we found that both pseudoscalar and vector correlators exhibit a nice $cosh$ behaviour, leading to a plateau behaviour in the local screening masses as a function of distance. The rho and pi masses so determined show very little variation with the lattice spacing a. This augurs well for the use of chiral Fermions, and further suggests the small deviations of these masses from the ideal gas values are genuine effects of interactions.


title: Equation of state for two-flavor QCD with an improved Wilson quark action at non-zero chemical potential

abstract: The QCD thermodynamics on the lattice provides fundamental theoretical grounds to analyze the various experimental data in relativistic heavy ion collisions. So far, most of the numerical simulations on the lattice have been performed by using the staggered-type fermion actions. Therefore it is important to carry out studies using different fermion formulations to test the uncertainties of the lattice QCD results. For this purpose, we perform systematic simulations of two-flavor QCD with an improved Wilson quark action to investigate the equation of state. We report the current status of our project and show the preliminary results of the Taylor expansion coefficients of the thermodynamic grand partition function in terms of chemical potential.


title: Phase of the Fermion Determinant at Nonzero Chemical Potential

abstract: We show that in the microscopic domain of QCD (also known as the $\epsilon$-domain) at nonzero chemical potential the average phase factor of the fermion determinant is nonzero for $\mu < m_\pi/2$ and is exponentially suppressed for larger values of the chemical potential. This follows from the chiral Lagrangian that describes the low-energy limit of the expectation value of the phase factor. Explicit expressions for the average phase factor are derived using a random matrix formulation of the zero momentum limit of this chiral Lagrangian.


title: Renormalization constants for Lattice QCD: new results from Numerical Stochastic Perturbation Theory

abstract: By making use of Numerical Stochastic Perturbation Theory (NSPT) we can compute renormalization constants for Lattice QCD to high orders, e.g. three or four loops for quark bilinears. We report on the status of our computations, which provide several results for Wilson quarks and in particular (values and/or ratios of) Z_V, Z_A, Z_S, Z_P. Results are given for various number of flavors (n_f = 0, 2, 3, 4). While we recall the care which is due for the computation of quantities for which an anomalous dimension is in place, we point out that our computational framework is well suited to a variety of other calculations and we briefly discuss the application of NSPT to other regularizations (in particular the Clover action).


title: Spin-dependent potentials from lattice QCD

abstract: The spin-dependent corrections to the static inter-quark potential are phenomenologically relevant to describing the fine and hyperfine spin splitting of the heavy quarkonium spectra. We investigate these corrections, which are represented as the field strength correlators on the quark-antiquark source, in SU(3) lattice gauge theory. We use the Polyakov loop correlation function as the quark-antiquark source, and by employing the multi-level algorithm, we obtain remarkably clean signals for these corrections up to intermediate distances of around 0.6 fm. Our observation suggests several new features of the corrections.


title: Analysis of the finite temperature transition for 3 flavor QCD using p4-improved staggered fermions

abstract: We present a calculation of the transition temperature for 3 flavor QCD using p4-improved staggered fermions with two different variants of fattened links. We examine various susceptibilities at two different values for the temporal extent - $N_t = 4$ and $N_t = 6$ - in the vicinity of the 3 flavor transition at vanishing chemical potential. For $N_t = 4$, we study bare quark masses in the range $m_qa = 0.005$ to $m_qa = 0.1$, and for $N_t = 6$, we use $m_qa = 0.02$ to $m_qa = 0.2$. We also calculate hadron masses and the static quark potential at zero temperature in order to set the scale for the transition temperature and to study quark mass and cut-off dependent scaling properties of the transition temperature. A comparison of the R Algorithm and the RHMC algorithm for finite-temperature simulations is also made


title: Nuclear Force from Lattice QCD

abstract: The first lattice QCD result on the nuclear force (the NN potential) is presented in the quenched level. The standard Wilson gauge action and the standard Wilson quark action are employed on the lattice of the size 16^3\times 24 with the gauge coupling beta=5.7 and the hopping parameter kappa=0.1665. To obtain the NN potential, we adopt a method recently proposed by CP-PACS collaboration to study the pi pi scattering phase shift. It turns out that this method provides the NN potentials which are faithful to those obtained in the analysis of NN scattering data. By identifying the equal-time Bethe-Salpeter wave function with the Schroedinger wave function for the two nucleon system, the NN potential is reconstructed so that the wave function satisfies the time-independent Schroedinger equation. In this report, we restrict ourselves to the J^P=0^+ and I=1 channel, which enables us to pick up unambiguously the ``central'' NN potential V_{central}(r). The resulting potential is seen to posses a clear repulsive core of about 500 MeV at short distance (r < 0.5 fm). Although the attraction in the intermediate and long distance regions is still missing in the present lattice set-up, our method is appeared to be quite promising in reconstructing the NN potential with lattice QCD.


title: Regularizing QCD with staggered fermions and the fourth root trick

abstract: We investigate the properties of staggered-fermion lattice QCD in which the fourth root of the fermion determinant is taken. We show that this theory is non-local at non-zero lattice spacing $a$, and that the non-locality is caused by the breaking of taste symmetry at $a\ne 0$. We then present a renormalization-group based argument that the theory restores taste symmetry in the continuum limit. As a consequence the theory is local in that limit, and falls into the correct universality class. Finally, we argue that the correct effective theory for the physics of Goldstone bosons at $a\ne 0$ is given by staggered chiral perturbation theory with the replica trick.


title: The Spectrum of tmLQCD with Quark and Link Smearing

abstract: The effect of using smeared sink operators on the hadron spectrum is studied for quenched twisted mass lattice QCD with up, down, and strange quarks. Gaussian smearing is used for quark fields, and stout link smearing for gauge fields. Smeared correlators are found to be dominated by the ground state with a small contribution from excited states, leading to an improved determination of some ground state masses.


title: Quenched penguins, the Delta I=1/2 rule, and epsilon'/epsilon

abstract: The transformation properties of strong penguin operators under the action of the flavor group change when they are considered as operators in (partially) quenched QCD instead of the unquenched theory. An ambiguity arises, which is parameterized by new low-energy constants in the effective theory describing non-leptonic kaon decays in the (partially) quenched setting. Here we summarize results of the analysis for the complete set of three-flavor strong penguin operators, consisting of products of two left-handed flavor currents, or a left- and a right-handed current. Our results imply that (partially) quenched lattice computations of the Delta I=1/2 rule and epsilon'/epsilon are both affected by ambiguities intrinsic to the use of the quenched approximation at leading order in the chiral expansion. The only exception is the partially quenched case with three light sea quarks, consistent with general expectations. We also address the issue of quenched ambiguities in the case of an active charm, correcting and extending that in Phys. Rev. D 74, 014509 (2006).


title: On the infrared behavior of Green's functions in Yang-Mills theory

abstract: Non-perturbative properties of QCD, such as color confinement, are encoded in the infrared behavior of correlation functions, e.g. propagators and vertices. Various analytic predictions have been suggested for these quantities in various gauges. Here we numerically test these predictions using lattice gauge theory. In particular, we present results for the 2- and 3-point functions for SU(2) Landau-gauge Yang-Mills theory in three and in four dimensions. Special attention is paid to systematic finite-volume effects. The gluon and ghost propagators are also evaluated in the so-called interpolating gauge (between the Landau and the Coulomb gauge), in order to study their gauge-dependence. Finally, we consider these propagators in Landau gauge at finite temperature, with the aim of understanding the effect of the deconfinement phase transition on their infrared behavior. All our results are compatible with the so-called Gribov-Zwanziger confinement scenario.


title: Nucleon structure in the chiral regime with domain wall fermions on an improved staggered sea

abstract: Moments of unpolarized, helicity, and transversity distributions, electromagnetic form factors, and generalized form factors of the nucleon are presented from a preliminary analysis of lattice results using pion masses down to 359 MeV. The twist two matrix elements are calculated using a mixed action of domain wall valence quarks and asqtad staggered sea quarks and are renormalized perturbatively. Several observables are extrapolated to the physical limit using chiral perturbation theory. Results are compared with experimental moments of quark distributions and electromagnetic form factors and phenomenologically determined generalized form factors, and the implications on the transverse structure and spin content of the nucleon are discussed.


title: Definition and parametrization of non-perturbative effects in quenched QCD

abstract: The notion of a non-perturbative effect is ambiguous if it requires the subtraction of a perturbative part defined by a diverging series. A common procedure consists in dropping the order of minimal contribution and the higher orders. This allows us to isolate very accurately the one-instanton effect for the double-well potential. For the one plaquette gauge theory, an exact analytical expression can be written for the non-perturbative part. We report recent attempts to extend this approach to the average plaquette of quenched QCD. Our goal is to express the non-perturbative effects in terms of expressions of the form beta^B exp(-A beta) calculable semi-classically. The situation is complicated by zeroes of the partition function in the complex $\beta$ plane (presumably near 5.75 pm i 0.2). We discuss two methods to describe the intermediate and large order behavior of the perturbative series. One is inspired by mean field theory (logarithmic specific heat) and reproduces accurately the known perturbative series with only two free parameters. A diagrammatic interpretation of this fact is still lacking. The other is based on infra-red renormalons with a factorial growth showing up at order larger than 20 and a possible effective theory interpretation. These extrapolations are compatible with the non-perturbative part of the plaquette being proportional to a^4. We propose an exponential parametrization to the corrections to the universal part of the beta function and find results compatible with the suggestion of a^2 corrections made by C. Allton.


title: Heavy baryon mass spectrum from lattice QCD with 2+1 flavors

abstract: We study the heavy baryon mass spectrum on gauge configurations that include 2+1 flavors of dynamical improved staggered quarks. A valence clover heavy quark is combined with two improved staggered light quark propagators to form baryons with different flavors. We are using MILC coarse gauge configurations with a lattice spacing of about 0.12 fm. In this preliminary investigation, we explore the chiral limit by studying two light sea quark masses, three different strange valence quark masses, and nine different light valence quark masses ranging from 0.1 to 0.4 times the nominal strange quark mass.


title: Quark propagators at finite temperature with the clover action

abstract: We study properties of the finite temperature quark propagator by using the SU(3) quenched lattice simulation in the Landau gauge and report numerical results of the standard Wilson quark case as well as the improved clover one. The mass function in the deconfinement phase is different from that of the confinement phase, especially at low momentum regions.


title: Dynamical overlap fermion at fixed topology

abstract: We launched a project to perform dymanical fermion simulations using the overlap fermion formulation for sea quarks. In order to avoid the appearace of near-zero modes of the hermitian Wilson-Dirac operator $H_W$, we introduce a pair of extra Wilson fermions with a large negative mass term. Crossing of the topological boundary is then strictly prohibited, and the topological charge is conserved during simulations. It makes the simulations substantially faster compared to the algorithms which allow the topology change. We discuss on the finite volume effects due to the fixed global topology.


title: On the phase structure of a chiral invariant Higgs-Yukawa model

abstract: In the past the construction of Higgs-Yukawa models on the lattice was blocked by the lack of a consistent definition of a chiral invariant Yukawa coupling term. Here, we consider a chiral invariant Higgs-Yukawa model based on the overlap operator, realized by the Neuberger-Dirac operator. As a first step towards a numerical examination of this model we study its phase diagram by means of an analytic 1/N-expansion, which is possible for small and for large values of the Yukawa coupling constant. In the case of strong Yukawa couplings the model effectively becomes an O(4)-symmetric non-linear sigma-model.


title: Static quark free energies at finite temperature with two flavors of improved Wilson quarks

abstract: Polyakov loop correlations at finite temperature in two-flavor QCD are studied in lattice simulations with the RG-improved gluon action and the clover-improved Wilson quark action. From the simulations on a $16^3 \times 4$ lattice, we extract the 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. 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 given in terms of $g_{\rm eff}(T)$. Also we made a comparison between our results with the Wilson quark and those with the staggered quark previously reported.


title: Full QCD in external chromomagnetic field

abstract: We investigate the deconfining phase transition in full QCD with two flavors of staggered fermions in presence of a constant abelian chromomagnetic field. We find that the deconfinement temperature decreases and eventually goes to zero by increasing the strength of the chromomagnetic field. Moreover our results suggest that the chiral transition coincides with the deconfinement transition and therefore even the chiral critical temperature depends on the applied chromomagnetic field. We also find that the chiral condensate increases with the strength of the chromomagnetic field.


title: First results with two light flavours of quarks with maximally twisted mass

abstract: We report on first results of an ongoing effort to simulate lattice QCD with two degenerate flavours of quarks by means of the twisted mass formulation tuned to maximal twist. By utilising recent improvements of the HMC algorithm, pseudo-scalar masses well below 300 MeV are simulated on volumes with spatial size larger than 2 fm at values of the lattice spacing similar or smaller than 0.1 fm. We present first evidence that scaling violations in the pseudo-scalar decay constant are small and well compatible with O(a) improvement. Additionally, exploratory results for the case of N_f=2+1+1 flavours are discussed.


title: Partially quenched QCD with a chemical potential

abstract: Using a chiral random matrix theory we can now derive the low energy partition functions and Dirac eigenvalue correlations of QCD with different chemical potentials for the dynamical and valence quarks. The results can also be extended to complex (and purely imaginary) chemical potential. We also discuss possible applications such as fitting to low energy constants and understanding the phase diagrams of the full and partially quenched theories.


title: The QCD equation of state with asqtad staggered fermions

abstract: We report on our result for the equation of state (EOS) with a Symanzik improved gauge action and the asqtad improved staggered fermion action at $N_t=4$ and 6. In our dynamical simulations with 2+1 flavors we use the inexact R algorithm and here we estimate the finite step-size systematic error on the EOS. Finally we discuss the non-zero chemical potential extension of the EOS and give some preliminary results.


title: Electromagnetic and spin polarisabilities in lattice QCD

abstract: We discuss the extraction of the electromagnetic and spin polarisabilities of nucleons from lattice QCD. We show that the external field method can be used to measure all the electromagnetic and spin polarisabilities including those of charged particles. We then turn to the extrapolations required to connect such calculations to experiment in the context of chiral perturbation theory, finding a strong dependence on the lattice volume and quark masses.


title: QCD thermodynamics with $N_f=2+1$ near the continuum limit at realistic quark masses

abstract: We report on our study of QCD thermodynamics with 2+1 flavors of dynamical quarks. In this proceeding we present several thermodynamic quantities and our recent calculation of the critical temperature. In order to investigate the thermodynamic properties of QCD near the continuum limit we adopt improved staggered (p4) quarks coupled with tree-level Symanzik improved glue on $N_t=4$ and 6 lattices. The simulations are performed with a physical value of the strange quark mass and light quark masses which are in the range of $m_q/m_s=0.05-0.4$. The lightest quark mass corresponds to a pion mass of about 150 MeV.


title: $\rho$ meson decay from the lattice

abstract: We present preliminary results on the $\rho$ meson decay width estimated from the scattering phase shift of the I=1 two-pion system. The phase shift is calculated by the finite size formula for non-zero total momentum frame (the moving frame) derived by Rummukainen and Gottlieb, using the $N_f=2$ improved Wilson fermion action at $m_\pi/m_\rho=0.41$ and $L=2.53 {\rm fm}$.


title: An estimate of the eta and eta-prime meson masses in Nf=2+1 lattice QCD

abstract: Masses of the eta and eta-prime mesons are estimated in Nf=2+1 lattice QCD with the non-perturbatively O(a) improved Wilson quark action and the Iwasaki RG-improved gluon action, using CP-PACS/JLQCD configurations on a 16^3 x 32 lattice at beta=1.83 (lattice spacing is 0.122 fm). We apply a stochastic noise estimator technique combined with smearing method to evaluate correlators among flavor SU(2) singlet pseudoscalar operators and strange pseudoscalar operators for 10 combinations of up/down and strange quark masses. The correlator matrix is then diagonalized to identify signals for mass eigenstates. Masses of the ground state and the first excited state extrapolated to the physical point are m_eta= 0.545(16) GeV and m_eta-prime= 0.871(46) GeV, being close to the experimental values of the eta and eta-prime masses.


title: Calculation of Neutron EDM in quenched and full QCD

abstract: We report on a direct lattice calculation of the neutron EDM(NEDM) using the external electric field method in both quenched and full QCD. In quenched QCD, we use a $24^3\times 32$ lattice at $\beta=2.6$ with the Iwasaki gauge action and the clover fermion action to examine the viability of this method. In particular we investigate possible effects of violation of the periodic boundary condition of the external electric field on the NEDM signal. We also study the quark mass dependence of NEDM in quenched QCD, and observe that NEDM seems to remain non-zero toward the chiral limit because of the quenched artifact. In 2-flavor full QCD we employ configurations generated by the CP-PACS collaboration on a $24^3\times 48$ lattice at $\beta=2.1$ with the same gluon and quark actions as in the quenched case. Since the number of configurations is limited, we employ 8 different source points per one configuration and take an average over them. Our preliminary result at three quark masses ($m_{PS}/m_V\simeq 0.81, 0.76, 0.69$) indicates that non-zero value for NEDM can be obtained in full QCD. Statistical errors, however, are still too large to show the theoretically expected behaviour for NEDM in full QCD that it vanishes in the chiral limit.


title: Gluon mass generation and infrared Abelian dominance in Yang-Mills theory

abstract: The dual superconductivity is believed to be a promising mechanism for quark confinement. Indeed, what this picture is true has been confirmed in the maximal Abelian (MA) gauge. However, it is not yet confirmed in any other gauge and the MA gauge explicitly breaks color symmetry. To remedy this defect, we propose to use our compact formulation of a non-linear change of variables on a lattice. This formulation has succeeded to extract the magnetic monopole with integer-valued magnetic charge in the gauge-invariant way. In this talk, we present measurements of various correlation functions for the operators constructed from the CFN variables in SU(2) Yang-Mills theory. Some of our results reproduce previous results obtained in MA gauge, e.g., DeGrant-Toussaint monopole, infrared Abelian dominance and off-diagonal gluon mass generation. These studies preserve color symmetry, in sharp contrast to the conventional MA gauge. We argue the gauge fixing independence of these results and the implications to quark confinement


title: Dynamical overlap fermions in the epsilon-regime

abstract: We report on the two-flavor QCD simulation in the epsilon-regime using the overlap fermion formulation. Sea quark mass is reduced to ~ 2 MeV on a 16^3x32 lattice with the lattice spacing ~ 0.11fm. Topological charge is fixed at Q=0. We compare the eigenvalue distribution of the overlap-Dirac operator with the prediction of the chiral random matrix theory. Preliminary results on meson correlators are also reported.


title: The First Moment of the Kaon Distribution Amplitude from N_f=2+1 Domain Wall Fermions

abstract: We present a lattice computation of the first moment of the kaon's leading-twist distribution amplitude. We use ensembles with 2+1 dynamical flavours of domain wall fermions and the Iwasaki gauge action from the RBC and UKQCD joint dataset. We observe the expected chiral behaviour and obtain <\xi>(2GeV)=3/5a_K^1(2GeV)=0.032(3), which agrees very well with other results obtained using QCD sum-rules and the recent lattice result from the UKQCD/QCDSF collaboration.


title: Improvement of algorithms for dynamical overlap fermions

abstract: We investigate the algorithms for dynamical overlap fermions aiming at improving the performance for large-scale simulations. We look for the best combination of Hybrid Monte Carlo options and iterative quark solvers with respect to the numerical costs. Our main target is a $N_f=2$ simulation with overlap fermion on a $16^3\times 32$ lattice at lattice spacing around 0.12 fm.


title: A Strategy to Study Confinement in QCD

abstract: The order and the universality class of the deconfining phase transition can provide insight into the mechanism of color confinement, in particular for N_f=2. The mechanism of confinement by monopole condensation is reviewed.


title: Nucleon structure in lattice QCD with dynamical domain-wall fermions quarks

abstract: We report RBC and RBC/UKQCD lattice QCD numerical calculations of nucleon electroweak matrix elements with dynamical domain-wall fermions (DWF) quarks. The first, RBC, set of dynamical DWF ensembles employs two degenerate flavors of DWF quarks and the DBW2 gauge action. Three sea quark mass values of 0.04, 0.03 and 0.02 in lattice units are used with about 200 gauge configurations each. The lattice cutoff is about 1.7 GeV and the spatial volume is about (1.9 fm)^3. Despite the small volume, the ratio of the isovector vector and axial charges g_A/g_V and that of structure function moments <x>_{u-d}/<x>_{Delta u - Delta d} are in agreement with experiment, and show only very mild quark mass dependence. The second, RBC/UK, set of ensembles employs one strange and two degenerate (up and down) dynamical DWF quarks and Iwasaki gauge action. The strange quark mass is set at 0.04, and three up/down mass values of 0.03, 0.02 and 0.01 in lattice units are used. The lattice cutoff is about 1.6 GeV and the spatial volume is about (3.0 fm)^3. Even with preliminary statistics of 25-30 gauge configurations, the ratios g_A/g_V and <x>_{u-d}/<x>_{Delta u - Delta d} are consistent with experiment and show only very mild quark mass dependence. Another structure function moment, d_1, though yet to be renormalized, appears small in both sets.


title: Method for simulating O(N) lattice models at finite density

abstract: We present a method for simulating relativistic and nonrelativistic scalar field theories at finite density, with matter transforming in the fundamental representation of the global symmetry group O(N). The method avoids the problem of complex probability weights which is present in conventional path integral Monte Carlo algorithms. To verify our approach, we simulate the free and interacting relativistic U(1)=O(2) theory in 2+1 dimensions. We compute the two-point correlation function and charge density as a function of chemical potential in the free theory. At weak phi^4 coupling and zero temperature we map the m^2-mu phase diagram and compare our numerical results with perturbative calculations. Finally, we compute properties of theT-mu phase diagram in the vicinity of the phase transition and at bare self-couplings large compared to the temperature and chemical potential.


title: Quarkonium at finite temperature

abstract: Lattice QCD studies on charmonium at finite temperature are presented After a discussion about problems for the Maximum Entropy Method applied to finite temperature lattice QCD, I show several results on charmonium spectral functions. The ``wave function'' of charmonium is also discussed to study the spatial correlation between quark and anti-quark in deconfinement phase.


title: Taste breaking effects in scalar meson correlators

abstract: As a consistency check of the staggered-fermion fourth-root approximation, we analyze the a_0 and f_0 correlators, including the effects of two-meson taste-multiplet intermediate states. Rooted staggered chiral perturbation theory describes the contributions from the pseudoscalar taste multiplets in terms of only a few low energy constants, which have all been previously determined by the MILC collaboration. In previous work one of us (Prelovsek) showed that the two-meson ``bubble'' contributions could explain the observed anomalies in the lattice data for the isovector a_0 channel. In the present work we extend this analysis to the f_0 channel. On a MILC collaboration lattice ensemble at 0.12 fm with 2+1 flavors of Asqtad-improved staggered fermions, we have made new measurements of correlators in both channels for a variety of momenta. A fit to these correlators gives low energy constants that are reasonably consistent with previous determinations by the MILC collaboration.}


title: Quark confinement and gauge invariant monopoles in SU(2) YM

abstract: We give a short review of recently obtained results on a new lattice formulation of the non-linear change of variables which was once called the Cho--Faddeev--Niemi decomposition in SU(2) Yang-Mills theory. Based on this formulation, we proposed a new gauge-invariant definition of the magnetic monopole current which guarantees the magnetic charge quantization. We also demonstrated the magnetic monopole dominance in the string tension in SU(2) Yang-Mills theory on a lattice. Our formulation enables one to reproduce in the gauge-invariant way remarkable results obtained so far only in the Maximally Abelian gauge.


title: Overlap Fermion in External Gravity

abstract: On a lattice, we construct an overlap Dirac operator which describes the propagation of a Dirac fermion in external gravity. The local Lorentz symmetry is manifestly realized as a lattice gauge symmetry, while the general coordinate invariance is expected to be restored only in the continuum limit. The lattice index density in the presence of a gravitational field is calculated.


title: About Orientifold Planar Equivalence on the Lattice

abstract: The orientifold planar equivalence is the equivalence in the large-N limit of the bosonic sectors of the super Yang-Mills and the QCD with a quark in the antisymmetric representation. I give a sketch of the proof of the orientifold planar equivalence in the strong-coupling and large-mass phase on the lattice. It is still matter of discussion, if its validity extends also in the continuum limit.


title: String tensions of SU(N) gauge theories in 2+1 dimensions

abstract: We calculate the energy spectrum of closed strings in SU(N) gauge theories with N=2,3,4,6,8 in 2+1 dimensions to a high accuracy. We attempt to control all systematic errors, and this allows us to perform a precise comparison with different theoretical predictions. When we study the dependence of the string mass on its length L we find that the Nambu-Goto prediction is a very good approximation down to relatively short lengths, where the Luscher term alone is insufficient. We then isolate the corrections to the Luscher term, and compare them to recent theoretical predictions, which indeed seem to be mildly preferred by the data. When we take these corrections into account and extract string tensions from the string masses, we find that their continuum limit is lower by 2%-1% from the predictions of Karabli, Kim, and Nair. The discrepancy decreases with N, but when we extrapolate our results to N=oo we still find a discrepancy of 0.88% which is a 4.5 sigma effect.


title: JLQCD's dynamical overlap project

abstract: We present an overview of our project of large-scale simulations with dynamical overlap fermions. The first production run in two-flavor QCD is on-going using the Iwasaki gauge action on a 16^3 \times 32 lattice at the lattice spacing of 0.12fm with six sea quark masses down to m_{s,phys}/6, where m_{s,phys} is the physical strange quark mass. We briefly introduce our choice of the lattice action and simulation algorithm, and describe the present status of the production run. Preliminary results on the light meson masses and the static quark potential are also reported.


title: An application of the UV-filtering preconditioner to the Polynomial Hybrid Monte Carlo algorithm

abstract: We apply the UV-filtering preconditioner, previously used to improve the Multi-Boson algorithm, to the Polynomial Hybrid Monte Carlo (UV-PHMC) algorithm. The performance test for the algorithm is given for the plaquette gauge action and the $O(a)$-improved Wilson action at $\beta=5.2, c_{\mathrm{sw}}=2.02, M_{\pi}/M_{\rho}\sim 0.8$ and 0.7 on a $16^3\times 48$ lattice. We find that the UV-filtering reduces the magnitude of the molecular dynamics force from the pseudo fermion by a factor 3 by tuning the UV-filter parameter. Combining with the multi-time scale molecular dynamics integrator we achieve a factor 2 improvement.


title: Finite temperature phase transition of two-flavor QCD with an improved Wilson quark action

abstract: We study the phase structure of QCD at finite temperatures with two flavors of dynamical quarks on a lattice with the size $N_s^3 \times N_t=16^3 \times 4$, using a renormalization group improved gauge action and a clover improved Wilson quark action. The simulations are made along the lines of constant physics determined in terms of $m_{\rm PS}/m_{\rm V}$ at zero-temperature. We show preliminary results for the spatial string tension in the high temperature phase.


title: Hadron form factors using density-density correlators

abstract: Gauge invariant density-density correlators yield detailed information on hadron structure. Hadron deformation and form factors can be extracted for momentum transfers up to about 6 GeV$^2$. We use stochastic techniques and dilution to compute the all to all propagator required for the exact evaluation of density-density correlators. We present first results for the pion form factor.


title: Low energy scattering parameters from the solutions of the non-relativistic Yukawa model on a 3-dimensional lattice

abstract: The numerical solutions of the non-relativistic Yukawa model on a 3-dimensional size lattice with periodic boundary conditions are obtained. The possibility to extract the corresponding -- infinite space -- low energy parameters and bound state binding energies from eigensates computed at finite lattice size is discussed.


title: Singlet Free Energies and Renormalized Polyakov Loop in full QCD

abstract: We calculate the free energy of a static quark anti-quark pair and the renormalized Polyakov loop in 2+1- and 3- flavor QCD using $16^3 \times 4$ and $16^3 \times 6$ lattices and improved staggered p4 action. We also compare the renormalized Polyakov loop with the results of our earlier studies.


title: Exploring Residual Gauge Symmetry Breaking

abstract: Simulations of pure-gauge SU(2) lattice gauge theory are performed in the minimal Coulomb gauge. This leaves a residual or remnant gauge symmetry still active which is global in three directions but still local in one. Using averaged fourth-dimension pointing links as a spin-like order parameter, the remnant symmetry appears to undergo spontaneous symmetry breaking at around $\beta = 2.6$. Both the Binder cumulant and the magnetization itself exhibit crossings in this region using lattices up to $20^4$, and a susceptibility peak is also observed. Finite size scaling indicates a weak first-order transition. The symmetry breaking is also observed to take place in the fundamental-adjoint plane, and is coincident with the strong first-order transition that exists there at large $\beta_{\rm{adjoint}}$. This provides confirmation that this phase transition is a symmetry-breaking transition. A well-known theorem concerning the instantaneous Coulomb potential has previously proven that a transition where such a Coulomb-gauge remnant symmetry breaks is necessarily deconfining.


title: Generalized Potts-Models and their Relevance for Gauge Theories

abstract: We study the Polyakov loop dynamics originating from finite-temperature Yang-Mills theory. The effective actions contain center-symmetric terms involving powers of the Polyakov loop, each with its own coupling. For a subclass with two couplings we perform a detailed analysis of the statistical mechanics involved. To this end we employ a modified mean field approximation and Monte Carlo simulations based on a novel cluster algorithm. We find excellent agreement of both approaches. The phase diagram exhibits both first and second order transitions between symmetric, ferromagnetic and antiferromagnetic phases with phase boundaries merging at three tricritical points. The critical exponents $\nu$ and $\gamma$ at the continuous transition between symmetric and antiferromagnetic phases are the same as for the 3-state spin Potts model.