• ### Numerical determination of partial spectrum of Hermitian matrices using a Lanczos method with selective reorthogonalization(1211.1180)

Nov. 6, 2012 math.NA, hep-lat
We introduce a new algorithm for finding the eigenvalues and eigenvectors of Hermitian matrices within a specified region, based upon the LANSO algorithm of Parlett and Scott. It uses selective reorthogonalization to avoid the duplication of eigenpairs in finite-precision arithmetic, but uses a new bound to decide when such reorthogonalization is required, and only reorthogonalizes with respect to eigenpairs within the region of interest. We investigate its performance for the Hermitian Wilson--Dirac operator (\gamma_5D) in lattice quantum chromodynamics, and compare it with previous methods.
• ### Time symmetric Bloch-Nordsieck infrared cancellation(1006.2622)

July 26, 2010 hep-ph
It is usually assumed that the infrared divergences that occur for virtual photons in Feynman diagrams cancel against the corresponding divergences in phase space integrals for the emission of real soft photons. If this were the case then why do the corresponding phase space integrals for the absorption of real soft photons not diverge too, as one might expect on grounds of time-reversal symmetry, and if so what cancels these divergences? We show that this apparent time asymmetry is an artifact of the way that the diagrams are calculated, and we introduce a procedure for such calculations that is manifestly time symmetric. Moreover, we show that it is not mathematically meaningful to distinguish arbitrarily soft massless particles of positive and negative energy.
• ### Results from 2+1 flavours of SLiNC fermions(0910.2963)

Oct. 15, 2009 hep-lat
QCD results are presented for a 2+1 flavour fermion clover action (which we call the SLiNC action). A method of tuning the quark masses to their physical values is discussed. In this method the singlet quark mass is kept fixed, which solves the problem of different renormalisations (for singlet and non-singlet quark masses) occuring for non-chirally invariant lattice fermions. This procedure enables a wide range of quark masses to be probed, including the case with a heavy up-down quark mass and light strange quark mass. Preliminary results show the correct splittings for the baryon (octet and) decuplet spectrum.
• ### Physical Results from 2+1 Flavor Domain Wall QCD and SU(2) Chiral Perturbation Theory(0804.0473)

April 3, 2008 hep-lat
We have simulated QCD using 2+1 flavors of domain wall quarks on a $(2.74 {\rm fm})^3$ volume with an inverse lattice scale of $a^{-1} = 1.729(28)$ GeV. The up and down (light) quarks are degenerate in our calculations and we have used four values for the ratio of light quark masses to the strange (heavy) quark mass in our simulations: 0.217, 0.350, 0.617 and 0.884. We have measured pseudoscalar meson masses and decay constants, the kaon bag parameter $B_K$ and vector meson couplings. We have used SU(2) chiral perturbation theory, which assumes only the up and down quark masses are small, and SU(3) chiral perturbation theory to extrapolate to the physical values for the light quark masses. While next-to-leading order formulae from both approaches fit our data for light quarks, we find the higher order corrections for SU(3) very large, making such fits unreliable. We also find that SU(3) does not fit our data when the quark masses are near the physical strange quark mass. Thus, we rely on SU(2) chiral perturbation theory for accurate results. We use the masses of the $\Omega$ baryon, and the $\pi$ and $K$ mesons to set the lattice scale and determine the quark masses. We then find $f_\pi = 124.1(3.6)_{\rm stat}(6.9)_{\rm syst} {\rm MeV}$, $f_K = 149.6(3.6)_{\rm stat}(6.3)_{\rm syst} {\rm MeV}$ and $f_K/f_\pi = 1.205(0.018)_{\rm stat}(0.062)_{\rm syst}$. Using non-perturbative renormalization to relate lattice regularized quark masses to RI-MOM masses, and perturbation theory to relate these to $\bar{\rm MS}$ we find $m_{ud}^{\bar{\rm MS}}(2 {\rm GeV}) = 3.72(0.16)_{\rm stat}(0.33)_{\rm ren}(0.18)_{\rm syst} {\rm MeV}$ and $m_{s}^{\bar{\rm MS}}(2 {\rm GeV}) = 107.3(4.4)_{\rm stat}(9.7)_{\rm ren}(4.9)_{\rm syst} {\rm MeV}$.
• ### 2+1 flavor domain wall QCD on a (2 fm)^3 lattice: light meson spectroscopy with Ls = 16(hep-lat/0701013)

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

Dec. 8, 2005 hep-lat
We present results for light meson masses and psedoscalar meson decay constants in 2+1 flavour domain wall QCD with the DBW2 and Iwasaki gauge actions, using lattices with linear sizes in the range 1.6 to 2.2fm and $u$ and $d$ quark masses as low as one quarter of the strange quark mass. All data were generated on the QCDOC machines at the University of Edinburgh and Brookhaven National Laboratory. Despite large residual masses and a limited number of sea quark mass values with which to perform chiral extrapolations, our results agree with experiment and scale within errors.
• ### Baryons in 2+1 flavour domain wall QCD(hep-lat/0511011)

Nov. 4, 2005 hep-lat
We present results for some of the light baryon masses and their excited states in 2+1 flavour domain wall QCD. We considered several lattice spacings, with the DBW2 and Iwasaki gauge actions and different sea quark masses on a volume of $16^3\times32$ and a fifth dimension of size 8. All data were generated on the QCDOC machines. Despite large residual massses and a limited number of sea quark mass values with which to perform chiral extrapolations, our results are in reasonable agreement with experiment and scale within errors. Finite size effects on most ensembles appear to be small.
• ### Quark spectra and light hadron phenomenology from overlap fermions with improved gauge field action(hep-lat/0310028)

Jan. 8, 2004 hep-lat
We present first results from a simulation of quenched overlap fermions with improved gauge field action. Among the quantities we study are the spectral properties of the overlap operator, the chiral condensate and topological charge, quark and hadron masses, and selected nucleon matrix elements. To make contact with continuum physics, we compute the renormalization constants of quark bilinear operators in perturbation theory and beyond.
• ### QCD with dynamical Wilson fermions at $\beta=5.5$(hep-lat/9607043)

July 19, 1996 hep-lat
We study QCD with two flavors of dynamical Wilson fermions at $\beta = 5.5$ and three values of $\kappa$. The corresponding pion masses are 0.375, 0.324 and 0.262 in lattice units, with pion to rho mass ratios of 0.76, 0.71 and 0.62, respectively. We use the configurations to compute the heavy quark potential, leading to lattice spacings of 0.110, 0.105 and 0.099 fm, and to compute spectroscopy for several different valence quark $\kappa$'s.
• ### Instabilities and Non-Reversibility of Molecular Dynamics Trajectories(hep-lat/9606004)

June 13, 1996 hep-lat
The theoretical justification of the Hybrid Monte Carlo algorithm depends upon the molecular dynamics trajectories within it being exactly reversible. If computations were carried out with exact arithmetic then it would be easy to ensure such reversibility, but the use of approximate floating point arithmetic inevitably introduces violations of reversibility. In the absence of evidence to the contrary, we are usually prepared to accept that such rounding errors can be made small enough to be innocuous, but in certain circumstances they are exponentially amplified and lead to blatantly erroneous results. We show that there are two types of instability of the molecular dynamics trajectories which lead to this behavior, instabilities due to insufficiently accurate numerical integration of Hamilton's equations, and intrinsic chaos in the underlying continuous fictitious time equations of motion themselves. We analyze the former for free field theory, and show that it is essentially a finite volume effect. For the latter we propose a hypothesis as to how the Liapunov exponent describing the chaotic behavior of the fictitious time equations of motion for an asymptotically free quantum field theory behaves as the system is taken to its continuum limit, and explain why this means that instabilities in molecular dynamics trajectories are not a significant problem for Hybrid Monte Carlo computations. We present data for pure $SU(3)$ gauge theory and for QCD with dynamical fermions on small lattices to illustrate and confirm some of our results.
• ### QCD $\beta$ Function with Two Flavors of Dynamical Wilson Fermions(hep-lat/9602010)

May 6, 1996 hep-lat
We test the asymptotic scaling behavior of state-of-the-art simulations of QCD with two flavors of light Wilson fermions. This is done by matching $\pi$ and $\rho$ masses on lattices of size $16^3\times32$ and $8^3\times16$. We find that at $\beta=6/g^2=5.3$ matching is not possible over a range extending down to $\beta=3.5$. The large lattice data at $\beta=5.5$ matches the small lattice values at $\beta=4.9(1)$ leading to a shift $\Delta\beta=0.6(1)$, considerably larger than the perturbative prediction of $0.45$. In both cases we conclude that the simulations are very far from the asymptotic scaling region.
• ### Status of the 0.8 Teraflops Supercomputer at Columbia(hep-lat/9509075)

Sept. 21, 1995 hep-lat
The first stage in the construction of the 0.8 Teraflops Supercomputer at Columbia, a working, two node parallel computer, has been successfully completed. The next stage, a 512 node, 26 Gigaflops prototype, is in its final construction phase. A general description and current status of the hardware and software is presented.
• ### Cooling in QCD Spectroscopy(hep-lat/9412043)

Dec. 8, 1994 hep-lat
We test the cooling algorithm with gluonic and staggered hadronic spectroscopy on $SU(3)$ gauge field configurations generated with two flavors of staggered dynamical fermions. We find cooling is not reliable as the basis for improved hadronic operators. We also find that performing cooling sweeps to reveal more clearly the topological properties of the gauge fields eliminates the spin structure of the hadron spectrum.
• ### Toward the QCD $\beta$ Function with Dynamical Wilson Fermions(hep-lat/9411057)

Nov. 25, 1994 hep-lat
We present data for the scaling behavior of lattice QCD with two flavors of light Wilson fermions. This is done by matching $\pi$ n and $\rho$ masses at the two lattice sizes of $16^3\times32$ and $8^3\times16$. We find that at $\beta=6/g^2\equiv5.3$ the matching does not occur over a range extending down to $\beta=3.5$. For $\beta=5.5$ matching may be achieved at $\beta=4.8-4.9$, leading to a $\Delta \beta =0.6-0.7$ which is higher than the perturbative 2-loop value of 0.45. In both cases we conclude that the simulations are very far from the perturbative scaling region.
• ### The heavy quark potential in QCD with 2 flavors of dynamical quarks(hep-lat/9401025)

Jan. 20, 1994 hep-lat
We compute the heavy quark potential on configurations generated by the HEMCGC collaboration with dynamical staggered fermions at $6/g^2 = 5.6$ and with dynamical Wilson fermions at $6/g^2 = 5.3$. The computations are done on $16^3 \times 32$ lattices, corresponding to physical sizes of about 1.6 and 2.3 fm, respectively. Up to the distances probed no sign of string breaking is detectable. We also compute the recently proposed scale $r_0$ defined by $r_0^2 F(r_0) = 1.65$.
• ### Simple Hadronic Matrix Elements with Wilson Valence Quarks and Dynamical Staggered Fermions at ${\bf 6/g^2=5.6}$(hep-lat/9301018)

Jan. 26, 1993 hep-lat
We have measured some simple matrix elements for pseudoscalar and vector mesons made of Wilson valence quarks and staggered sea quarks at $\beta=5.6$ at sea quark masses $am_q=0.01$ and 0.025. Our measurements include the decay constants of pseudoscalars (including $f_D$), the wave function at the origin (or decay constant) of vector mesons, and the calculation of quark masses from current algebra. The effects of sea quarks on the simulations are small. We make comparisons to quenched simulations at similar values of the lattice spacing ($1/a \simeq 2$ GeV).
• ### On the Dynamics of Light Quarks in QCD(hep-lat/9212034)

Dec. 24, 1992 hep-lat
We describe recent results concerning the behavior of lattice QCD with light dynamical Wilson and Staggered quarks. We show that it is possible to reach regions of parameter space with light pions $m_\pi\approx 0.2/a$ using Wilson fermions. If the Hybrid Molecular Dynamics (HMD) algorithm is used with the same parameters it gives incorrect results. We also present preliminary results using a higher-order integration scheme.