How to compute Green's Functions for entire Mass Trajectories within Krylov Solvers

The availability of efficient Krylov subspace solvers play a vital role for the solution of a variety of numerical problems in computational science. Here we consider lattice field theory. We present a new general numerical method to compute many Green's functions for complex non-singular matrices within one iteration process. Our procedure applies to matrices of structure $A=D-m$, with $m$ proportional to the unit matrix, and can be integrated within any Krylov subspace solver. We can compute the derivatives $x^{(n)}$ of the solution vector $x$ with respect to the parameter $m$ and construct the Taylor expansion of $x$ around $m$. We demonstrate the advantages of our method using a minimal residual solver. Here the procedure requires $1$ intermediate vector for each Green's function to compute. As real life example, we determine a mass trajectory of the Wilson fermion matrix for lattice QCD. Here we find that we can obtain Green's functions at all masses $\geq m$ at the price of one inversion at mass $m$.Comment: 11 pages, 2 eps-figures, needs epsf.st

Light Spectrum and Decay Constants in Full QCD with Wilson Fermions

We present results from an analysis of the light spectrum and the decay constants f_{\pi} and f_V^{-1} in Full QCD with n_f=2 Wilson fermions at a coupling of beta=5.6 on a 16^3x32 lattice.Comment: 3 pages, LaTeX with 4 eps figures, Talk presented at LATTICE96(spectrum

An Estimate of alpha_S from Bottomonium in Unquenched QCD

We estimate the strong coupling constant from the perturbative expansion of the plaquette. The scale is set by the 2S-1S and 1P-1S splittings in bottomonium which are computed in NRQCD on dynamical gauge configurations with nf=2 degenerate Wilson quarks at intermediate masses. We have increased the statistics of our spectrum calculation in order to reliably extrapolate in the sea-quark mass. We find a value of alpha_MS(m_Z) = 0.1118(26) which is somewhat lower than previous estimates within NRQCD.Comment: LATTICE98(heavyqk

Flavour singlet pseudoscalar masses in N_f = 2 QCD

We perform a lattice mass analysis in the flavour singlet pseudoscalar channel on the SESAM and TXL full QCD vacuum configurations, with 2 active flavours of dynamical Wilson fermions at beta = 5.6. At our inverse lattice spacing, a^-1 = 2.3 GeV, we retrieve by a chiral extrapolation to the physical light quark masses the value m_eta' = 3.7(+8)(-4) m_pi. A crude extrapolation from (N_f = 3) phenomenology would suggest m_eta' \approx 5.1 m_pi for N_f = 2 QCD. we verify that the mass gap between the singlet state eta' and the pi flavour triplt state is due to gauge configurations with non-trivial topology.Comment: 8 pages, 10 figure

Glueballs and string breaking from full QCD

We present results on the static potential, and torelon and glueball masses from simulations of QCD with two flavours of dynamical Wilson fermions on $16^3\times 32$ and $24^3\times 40$ lattices at $\beta=5.6$.Comment: Talk presented by Gunnar Bali at International Symposium on Lattice Field Theories (Lattice 97), Edinburgh, July 1997, 3 pages LaTeX (epscrc2.sty) with 4 eps figure

To what distances do we know the confining potential?

We argue that asymptotically linear static potential is built in into the common procedure of extracting it from lattice Wilson loop measurements. To illustrate the point, we extract the potential by the standard lattice method in a model vacuum made of instantons. A beautiful infinitely rising linear potential is obtained in the case where the true potential is actually flattening. We argue that the flux tube formation might be also an artifact of the lattice procedure and not necessarily a measured physical effect. We conclude that at present the rising potential is known for sure up to no more than about 0.7 fm. It may explain why no screening has been clearly observed so far for adjoint sources and for fundamental sources but with dynamical fermions. Finally, we speculate on how confinement could be achieved even if the static potential in the pure glue theory is not infinitely rising.Comment: 16 pages, 5 figures. Additional arguments presented, a new figure and references adde

Flavor Singlet Meson Mass in the Continuum Limit in Two-Flavor Lattice QCD

We present results for the mass of the eta-prime meson in the continuum limit for two-flavor lattice QCD, calculated on the CP-PACS computer, using a renormalization-group improved gauge action, and Sheikholeslami and Wohlert's fermion action with tadpole-improved csw. Correlation functions are measured at three values of the coupling constant beta corresponding to the lattice spacing a approx. 0.22, 0.16, 0.11 fm and for four values of the quark mass parameter kappa corresponding to mpi over mrho approx. 0.8, 0.75, 0.7 and 0.6. For each beta, kappa pair, 400-800 gauge configurations are used. The two-loop diagrams are evaluated using a noisy source method. We calculate eta-prime propagators using local sources, and find that excited state contributions are much reduced by smearing. A full analysis for the smeared propagators gives metaprime=0.960(87)+0.036-0.248 GeV, in the continuum limit, where the second error represents the systematic uncertainty coming from varying the functional form for chiral and continuum extrapolations.Comment: 9 pages, 19 figures, 4 table