Continuing Progress on a Lattice QCD Software Infrastructure

We report on the progress of the software effort in the QCD Application Area of SciDAC. In particular, we discuss how the software developed under SciDAC enabled the aggressive exploitation of leadership computers, and we report on progress in the area of QCD software for multi-core architectures.Comment: 5 Pages, to appear in the Proceedings of SciDAC 2008 conference, (Seattle, July 13-17, 2008), Conference Poster Presentation Proceeding

Application of heavy-quark effective theory to lattice QCD: III. Radiative corrections to heavy-heavy currents

We apply heavy-quark effective theory (HQET) to separate long- and short-distance effects of heavy quarks in lattice gauge theory. In this paper we focus on flavor-changing currents that mediate transitions from one heavy flavor to another. We stress differences in the formalism for heavy-light currents, which are discussed in a companion paper, showing how HQET provides a systematic matching procedure. We obtain one-loop results for the matching factors of lattice currents, needed for heavy-quark phenomenology, such as the calculation of zero-recoil form factors for the semileptonic decays $B\to D^{(*)}l\nu$. Results for the Brodsky-Lepage-Mackenzie scale $q^*$ are also given.Comment: 35 pages, 17 figures. Program LatHQ2QCD to compute matching one-loop coefficients available at http://theory.fnal.gov/people/kronfeld/LatHQ2QCD

Dual Superconductor Mechanism of Confinement on the Lattice

We investigate the dual superconductor mechanism of confinement for pure SU(2) lattice gauge theory in the maximally abelian gauge. We focus on the the dual Meissner effect. We find that the transverse distribution of the longitudinal chromoelectric field due to a static quark-antiquark pair satisfies the dual London equation. Moreover we show that the size of the flux tube scales according to asymptotic freedom.Comment: LaTeX, 9 pages, 6 figures available as a PostScript file from L. Cosmai, BARI - TH 110/9

Charmonium mass splittings at the physical point

We present results from an ongoing study of mass splittings of the lowest lying states in the charmonium system. We use clover valence charm quarks in the Fermilab interpretation, an improved staggered (asqtad) action for sea quarks, and the one-loop, tadpole-improved gauge action for gluons. This study includes five lattice spacings, 0.15, 0.12, 0.09, 0.06, and 0.045 fm, with two sets of degenerate up- and down-quark masses for most spacings. We use an enlarged set of interpolation operators and a variational analysis that permits study of various low-lying excited states. The masses of the sea quarks and charm valence quark are adjusted to their physical values. This large set of gauge configurations allows us to extrapolate results to the continuum physical point and test the methodology.Comment: 7 pp, 6 figs, Lattice 201

Creutz Ratios From Color-Truncated Lattice Configurations

We investigate whether information about Creutz ratios is encoded, separately, in each gluon color component of numerically generated lattice configurations. Working in SU(2) lattice gauge theory in Landau gauge, we set two of the three gluon color components to zero, and compensate for the loss of two-thirds of the fluctuation by simply rescaling the remaining component by a factor of $\sqrt{3}$. Creutz ratios are then computed with this "abelianized" configuration. We find that the Creutz ratios of loops constructed from abelianized links converge to the usual Creutz ratios in the scaling regime.Comment: 5 pages, plain LaTeX, one figur

Heavy and Light Quarks with Lattice Chiral Fermions

The feasibility of using lattice chiral fermions which are free of $O(a)$ errors for both the heavy and light quarks is examined. The fact that the effective quark propagators in these fermions have the same form as that in the continuum with the quark mass being only an additive parameter to a chirally symmetric antihermitian Dirac operator is highlighted. This implies that there is no distinction between the heavy and light quarks and no mass dependent tuning of the action or operators as long as the discretization error $O(m^2 a^2)$ is negligible. Using the overlap fermion, we find that the $O(m^2a^2)$ (and $O(ma^2)$) errors in the dispersion relations of the pseudoscalar and vector mesons and the renormalization of the axial-vector current and scalar density are small. This suggests that the applicable range of $ma$ may be extended to $\sim 0.56$ with only 5% error, which is a factor of $\sim 2.4$ larger than that of the improved Wilson action. We show that the generalized Gell-Mann-Oakes-Renner relation with unequal masses can be utilized to determine the finite $ma$ errors in the renormalization of the matrix elements for the heavy-light decay constants and semileptonic decay constants of the B/D meson.Comment: final version to appear in Int. Jou. Mod. Phys.

Low lying charmonium states at the physical point

We present results for the mass splittings of low-lying charmonium states from a calculation with Wilson clover valence quarks with the Fermilab interpretation on an asqtad sea. We use five lattice spacings and two values of the light sea quark mass to extrapolate our results to the physical point. Sources of systematic uncertainty in our calculation are discussed and we compare our results for the 1S hyperfine splitting, the 1P-1S splitting and the P-wave spin orbit and tensor splittings to experiment.Comment: For the Fermilab Lattice and MILC Collaborations; 7 pages, 6 figures; Contribution to the 32nd International Symposium on Lattice Field Theory, 23-28 June, 2014, Columbia University New York, N

Remarks on abelian dominance

We used a renormalisation group based smoothing to address two questions related to abelian dominance. Smoothing drastically reduces short distance fluctuations but it preserves the long distance physical properties of the SU(2) configurations. This enabled us to extract the abelian heavy-quark potential from time-like Wilson loops on Polyakov gauge projected configurations. We obtained a very small string tension which is inconsistent with the string tension extracted from Polyakov loop correlators. This shows that the Polyakov gauge projected abelian configurations do not have a consistent physical meaning. We also applied the smoothing on SU(2) configurations to test how sensitive abelian dominance in the maximal abelian gauge is to the short distance fluctuations. We found that on smoothed SU(2) configurations the abelian string tension was about 30% smaller than the SU(2) string tension which was unaffected by smoothing. This suggests that the approximate abelian dominance found with the Wilson action is probably an accident and it has no fundamental physical relevance.Comment: 13 pages, LaTeX, 3 eps figure