• In this document we address an error discovered in the ensemble generation for our calculation of the $I=0$ $K\to\pi\pi$ amplitude (Phys. Rev. Lett. 115, 212001 (2015), arXiv:1505.07863) whereby the same random numbers were used for the two independent quark flavors, resulting in small but measurable correlations between gauge observables separated by 12 units in the y-direction. We conclude that the effects of this error are negligible compared to the overall errors on our calculation.
  • We report the first lattice QCD calculation of the complex kaon decay amplitude $A_0$ with physical kinematics, using a $32^3\times 64$ lattice volume and a single lattice spacing $a$, with $1/a= 1.3784(68)$ GeV. We find Re$(A_0) = 4.66(1.00)(1.26) \times 10^{-7}$ GeV and Im$(A_0) = -1.90(1.23)(1.08) \times 10^{-11}$ GeV, where the first error is statistical and the second systematic. The first value is in approximate agreement with the experimental result: Re$(A_0) = 3.3201(18) \times 10^{-7}$ GeV while the second can be used to compute the direct CP violating ratio Re$(\varepsilon'/\varepsilon)=1.38(5.15)(4.59)\times 10^{-4}$, which is $2.1\sigma$ below the experimental value $16.6(2.3)\times 10^{-4}$. The real part of $A_0$ is CP conserving and serves as a test of our method while the result for Re$(\varepsilon'/\varepsilon)$ provides a new test of the standard-model theory of CP violation, one which can be made more accurate with increasing computer capability.
  • We present new results for the amplitude $A_2$ for a kaon to decay into two pions with isospin $I=2$: Re$A_2 = 1.50(4)_\mathrm{stat}(14)_\mathrm{syst}\times 10^{-8}$ GeV; Im$A_2 = -6.99(20)_\mathrm{stat}(84)_\mathrm{syst}\times 10^{-13}$ GeV. These results were obtained from two ensembles generated at physical quark masses (in the isospin limit) with inverse lattice spacings $a^{-1}=1.728(4)$ GeV and $2.358(7)$ GeV. We are therefore able to perform a continuum extrapolation and hence largely to remove the dominant systematic uncertainty from our earlier results, that due to lattice artefacts. The only previous lattice computation of $K\to\pi\pi$ decays at physical kinematics was performed using an ensemble at a single, rather coarse, value of the lattice spacing ($a^{-1}\simeq 1.37(1)$ GeV). We confirm the observation that there is a significant cancellation between the two dominant contributions to Re$A_2$ which we suggest is an important ingredient in understanding the $\Delta I=1/2$ rule, Re$A_0$/Re$A_2\simeq 22.5$, where the subscript denotes the total isospin of the two-pion final state. Our result for $A_2$ implies that the electroweak penguin contribution to $\epsilon^\prime/\epsilon$ is Re($\epsilon^\prime/\epsilon)_\textrm{EWP}=-(6.6\pm 1.0)\times 10^{-4}$.
  • We present physical results for a variety of light hadronic quantities obtained via a combined analysis of three 2+1 flavour domain wall fermion ensemble sets. For two of our ensemble sets we used the Iwasaki gauge action with beta=2.13 (a^-1=1.75(4) GeV) and beta=2.25 (a^-1=2.31(4) GeV) and lattice sizes of 24^3 x 64 and 32^3 x 64 respectively, with unitary pion masses in the range 293(5)-417(10) MeV. The extent L_s for the 5^th dimension of the domain wall fermion formulation is L_s=16 in these ensembles. In this analysis we include a third ensemble set that makes use of the novel Iwasaki+DSDR (Dislocation Suppressing Determinant Ratio) gauge action at beta = 1.75 (a^-1=1.37(1) GeV) with a lattice size of 32^3 x 64 and L_s=32 to reach down to partially-quenched pion masses as low as 143(1) MeV and a unitary pion mass of 171(1) MeV, while retaining good chiral symmetry and topological tunneling. We demonstrate a significant improvement in our control over the chiral extrapolation, resulting in much improved continuum predictions for the above quantities. The main results of this analysis include the pion and kaon decay constants, f_\pi=127(3)_{stat}(3)_{sys} MeV and f_K = 152(3)_{stat}(2)_{sys} MeV respectively (f_K/f_\pi = 1.199(12)_{stat}(14)_{sys}); the average up/down quark mass and the strange-quark mass in the MSbar-scheme at 3 GeV, m_{ud}(MSbar, 3 GeV) = 3.05(8)_{stat}(6)_{sys} MeV and m_s(MSbar, 3 GeV) = 83.5(1.7)_{stat}(1.1)_{sys}; the neutral kaon mixing parameter in the MSbar-scheme at 3 GeV, B_K(MSbar,3 GeV) = 0.535(8)_{stat}(13)_{sys}, and in the RGI scheme, \hat B_K = 0.758(11)_{stat}(19)_{sys}; and the Sommer scales r_1 = 0.323(8)_{stat}(4)_{sys} fm and r_0 = 0.480(10)_{stat}(4)_{sys} (r_1/r_0 = 0.673(11)_{stat}(3)_{sys}). We also obtain values for the SU(2) ChPT effective couplings, \bar{l_3} = 2.91(23)_{stat}(7)_{sys}$ and \bar{l_4} = 3.99(16)_{stat}(9)_{sys}.
  • We develop and demonstrate techniques needed to compute the long distance contribution to the $K_{L}$-$K_{S}$ mass difference, $\Delta M_K$, in lattice QCD and carry out a first, exploratory calculation of this fundamental quantity. The calculation is performed on 2+1 flavor, domain wall fermion, $16^3\times32$ configurations with a 421 MeV pion mass and an inverse lattice spacing $1/a=1.73$ GeV. We include only current-current operators and drop all disconnected and double penguin diagrams. The short distance part of the mass difference in a 2+1 flavor calculation contains a quadratic divergence cut off by the lattice spacing. Here, this quadratic divergence is eliminated through the GIM mechanism by introducing a valence charm quark. The inclusion of the charm quark makes the complete calculation accessible to lattice methods provided the discretization errors associated with the charm quark can be controlled. The long distance effects are discussed for each parity channel separately. While we can see a clear signal in the parity odd channel, the signal to noise ratio in the parity even channel is exponentially decreasing as the separation between the two weak operators increases. We obtain a mass difference $\Delta M_K$ which ranges from $6.58(30)\times 10^{-12}$ MeV to $11.89(81)\times 10^{-12}$ MeV for kaon masses varying from 563 MeV to 839 MeV. Extensions of these methods are proposed which promise accurate results for both $\Delta M_K$ and $\epsilon_K$, including long distance effects.
  • We review the status of recent calculations by the RBC-UKQCD collaboration of the complex amplitude $A_2$, corresponding to the decay of a kaon to a two pion state with total isospin 2. In particular, we present preliminary results from two new ensembles: $48^3 \times 96$ with $a^{-1}=1.73$ GeV and $64^3 \times 128$ with $a^{-1}=2.3$ GeV, both at physical kinematics. Both ensembles were generated Iwasaki gauge action and domain wall fermion action with 2+1 flavours. These results, in comparison to our earlier ones on a $32^3$ DSDR lattice with $a^{-1}=1.36$ GeV, enable us to significantly reduce the discretization errors. The partial cancellation between the two dominant contractions contributing to Re($A_2$) has been confirmed and we believe that this cancellation is a major contribution to the $\Delta I=1/2$ rule.
  • There has been much speculation as to the origin of the \Delta I = 1/2 rule (Re A_0/Re A_2 \simeq 22.5). We find that the two dominant contributions to the \Delta I=3/2, K \to \pi \pi{} correlation functions have opposite signs leading to a significant cancellation. This partial cancellation occurs in our computation of Re A_2 with physical quark masses and kinematics (where we reproduce the experimental value of A_2) and also for heavier pions at threshold. For Re A_0, although we do not have results at physical kinematics, we do have results for pions at zero-momentum with m_\pi{} \simeq 420 MeV (Re A_0/Re A_2=9.1(2.1)) and m_\pi{} \simeq 330 MeV (Re A_0/Re A_2=12.0(1.7)). The contributions which partially cancel in Re A_2 are also the largest ones in Re A_0, but now they have the same sign and so enhance this amplitude. The emerging explanation of the \Delta I=1/2 rule is a combination of the perturbative running to scales of O(2 GeV), a relative suppression of Re A_2 through the cancellation of the two dominant contributions and the corresponding enhancement of Re A_0. QCD and EWP penguin operators make only very small contributions at such scales.
  • We present results on both the restoration of the spontaneously broken chiral symmetry and the effective restoration of the anomalously broken U(1)_A symmetry in finite temperature QCD at zero chemical potential using lattice QCD. We employ domain wall fermions on lattices with fixed temporal extent N_\tau = 8 and spatial extent N_\sigma = 16 in a temperature range of T = 139 - 195 MeV, corresponding to lattice spacings of a \approx 0.12 - 0.18 fm. In these calculations, we include two degenerate light quarks and a strange quark at fixed pion mass m_\pi = 200 MeV. The strange quark mass is set near its physical value. We also present results from a second set of finite temperature gauge configurations at the same volume and temporal extent with slightly heavier pion mass. To study chiral symmetry restoration, we calculate the chiral condensate, the disconnected chiral susceptibility, and susceptibilities in several meson channels of different quantum numbers. To study U(1)_A restoration, we calculate spatial correlators in the scalar and pseudo-scalar channels, as well as the corresponding susceptibilities. Furthermore, we also show results for the eigenvalue spectrum of the Dirac operator as a function of temperature, which can be connected to both U(1)_A and chiral symmetry restoration via Banks-Casher relations.
  • We report on the first realistic \emph{ab initio} calculation of a hadronic weak decay, that of the amplitude $A_2$ for a kaon to decay into two \pi-mesons with isospin 2. We find Re$A_2=(1.436\pm 0.063_{\textrm{stat}}\pm 0.258_{\textrm{syst}})\,10^{-8}\,\textrm{GeV}$ in good agreement with the experimental result and for the hitherto unknown imaginary part we find {Im}$\,A_2=-(6.83 \pm 0.51_{\textrm{stat}} \pm 1.30_{\textrm{syst}})\,10^{-13}\,{\rm GeV}$. Moreover combining our result for Im\,$A_2$ with experimental values of Re\,$A_2$, Re\,$A_0$ and $\epsilon^\prime/\epsilon$, we obtain the following value for the unknown ratio Im\,$A_0$/Re\,$A_0$ within the Standard Model: $\mathrm{Im}\,A_0/\mathrm{Re}\,A_0=-1.63(19)_{\mathrm{stat}}(20)_{\mathrm{syst}}\times10^{-4}$. One consequence of these results is that the contribution from Im\,$A_2$ to the direct CP violation parameter $\epsilon^{\prime}$ (the so-called Electroweak Penguin, EWP, contribution) is Re$(\epsilon^\prime/\epsilon)_{\mathrm{EWP}} = -(6.52 \pm 0.49_{\textrm{stat}} \pm 1.24_{\textrm{syst}}) \times 10^{-4}$. We explain why this calculation of $A_2$ represents a major milestone for lattice QCD and discuss the exciting prospects for a full quantitative understanding of CP-violation in kaon decays.
  • We study the region of the QCD phase transition using 2+1 flavors of domain wall fermions (DWF) and a $16^3 \times 8$ lattice volume with a fifth dimension of $L_s = 32$. The disconnected light quark chiral susceptibility, quark number susceptibility and the Polyakov loop suggest a chiral and deconfining crossover transition lying between 155 and 185 MeV for our choice of quark mass and lattice spacing. In this region the lattice scale deduced from the Sommer parameter $r_0$ is $a^{-1} \approx 1.3$ GeV, the pion mass is $\approx 300$ MeV and the kaon mass is approximately physical. The peak in the chiral susceptibility implies a pseudo critical temperature $T_c = 171(10)(17)$ MeV where the first error is associated with determining the peak location and the second with our unphysical light quark mass and non-zero lattice spacing. The effects of residual chiral symmetry breaking on the chiral condensate and disconnected chiral susceptibility are studied using several values of the valence $L_s$.
  • We present results on the equation of state in QCD with two light quark flavors and a heavier strange quark. Calculations with improved staggered fermions have been performed on lattices with temporal extent Nt =4 and 6 on a line of constant physics with almost physical quark mass values; the pion mass is about 220 MeV, and the strange quark mass is adjusted to its physical value. High statistics results on large lattices are obtained for bulk thermodynamic observables, i.e. pressure, energy and entropy density, at vanishing quark chemical potential for a wide range of temperatures, 140 MeV < T < 800 MeV. We present a detailed discussion of finite cut-off effects which become particularly significant for temperatures larger than about twice the transition temperature. At these high temperatures we also performed calculations of the trace anomaly on lattices with temporal extent Nt=8. Furthermore, we have performed an extensive analysis of zero temperature observables including the light and strange quark condensates and the static quark potential at zero temperature. These are used to set the temperature scale for thermodynamic observables and to calculate renormalized observables that are sensitive to deconfinement and chiral symmetry restoration and become order parameters in the infinite and zero quark mass limits, respectively.
  • We present the first results for neutral kaon mixing using 2+1 flavors of domain wall fermions. A new approach is used to extrapolate to the physical up and down quark masses from our numerical studies with pion masses in the range 240 -- 420 MeV; only $SU(2)_L \times SU(2)_R$ chiral symmetry is assumed and the kaon is not assumed to be light. Our main result is $B_K^{\bar{\rm MS}}(2 \mathrm{GeV}) = 0.524(10)(28)$ where the first error is statistical and the second incorporates estimates for all systematic errors.
  • We present a study of the flavor symmetry breaking in the pion spectrum for various improved staggered fermion actions. To study the effects of link fattening and tadpole improvement, we use three different variants of the p4 action - p4fat3, p4fat7, and p4fat7tad. These are compared to Asqtad and also to naive staggered. To study the pattern of symmetry breaking, we measure all 15 meson masses in the 4-flavor staggered theory. The measurements are done on a quenched gauge background, generated using a one-loop improved Symanzik action with $\beta=10/g^2 = 7.40, 7.75,$ and 8.00, corresponding to lattice spacings of approximately a = .31 fm., .21 fm., and .14 fm. We also study how the lattice scale set by the $\rho$ mass on each of these ensembles compares to one set by the static quark potential.
  • We present results for the static interquark potential, light meson and baryon masses, and light pseudoscalar meson decay constants obtained from simulations of domain wall QCD with one dynamical flavour approximating the $s$ quark, and two degenerate dynamical flavours with input bare masses ranging from $m_s$ to $m_s/4$ approximating the $u$ and $d$ quarks. We compare these quantities obtained using the Iwasaki and DBW2 improved gauge actions, and actions with larger rectangle coefficients, on $16^3\times32$ lattices. We seek parameter values at which both the chiral symmetry breaking residual mass due to the finite lattice extent in the fifth dimension and the Monte Carlo time history for topological charge are acceptable for this set of quark masses at lattice spacings above 0.1 fm. We find that the Iwasaki gauge action is best, demonstrating the feasibility of using QCDOC to generate ensembles which are good representations of the QCD path integral on lattices of up to 3 fm in spatial extent with lattice spacings in the range 0.09-0.13 fm. Despite large residual masses and a limited number of sea quark mass values with which to perform chiral extrapolations, our results for light hadronic physics scale and agree with experimental measurements within our statistical uncertainties.
  • We study the finite temperature transition in QCD with three flavors of equal masses using the R and RHMC algorithm on lattices with temporal extent N_{\tau}=4 and 6. For the transition temperature in the continuum limit we find r_0 T_c=0.429(8) for the light pseudo-scalar mass corresponding to the end point of the 1st order transition region. When comparing the results obtained with the R and RHMC algorithms for p4fat3 action we see no significant step-size errors down to a lightest pseudo-scalar mass of m_{ps} r_0=0.4.
  • We present a detailed calculation of the transition temperature in QCD with two light and one heavier (strange) quark mass on lattices with temporal extent N_t =4 and 6. Calculations with improved staggered fermions have been performed for various light to strange quark mass ratios in the range, 0.05 <= m_l/m_s <= 0.5, and with a strange quark mass fixed close to its physical value. From a combined extrapolation to the chiral (m_l -> 0) and continuum (aT = 1/N_t -> 0) limits we find for the transition temperature at the physical point T_c r_0 = 0.457(7) where the scale is set by the Sommer-scale parameter r_0 defined as the distance in the static quark potential at which the slope takes on the value, (dV_qq(r)/dr)_r=r_0 = 1.65/r_0^2. Using the currently best known value for r_0 this translates to a transition temperature T_c = 192(7)(4)MeV. The transition temperature in the chiral limit is about 3% smaller. We discuss current ambiguities in the determination of T_c in physical units and also comment on the universal scaling behavior of thermodynamic quantities in the chiral limit.
  • We present numerical results for the kaon B-parameter, B_K, determined in the quenched approximation of lattice QCD. Our simulations are performed using domain-wall fermions and the renormalization group improved, DBW2 gauge action which combine to give quarks with good chiral symmetry at finite lattice spacing. Operators are renormalized non-perturbatively using the RI/MOM scheme. We study scaling by performing the simulation on two different lattices with a^{-1} = 1.982(30) and 2.914(54) GeV. We combine this quenched scaling study with an earlier calculation of B_K using two flavors of dynamical, domain-wall quarks at a single lattice spacing to obtain B_K(MS,NDR,mu=2GeV)=0.563(21)(39)(30), were the first error is statistical, the second systematic (without quenching errors) and the third estimates the error due to quenching.
  • We discuss the implementation and properties of the quenched approximation in the calculation of the left-right, strong penguin contributions (i.e. Q_6) to epsilon'/epsilon. The coefficient of the new chiral logarithm, discovered by Golterman and Pallante, which appears at leading order in quenched chiral perturbation theory is evaluated using both the method proposed by those authors and by an improved approach which is free of power divergent corrections. The result implies a large quenching artifact in the contribution of Q_6 to epsilon'/epsilon. This failure of the quenched approximation affects only the strong penguin operators and so does not affect the Q_8 contribution to epsilon'/epsilon nor Re A_0, Re A_2 and thus the Delta I=1/2 rule at tree level in chiral perturbation theory.
  • The architecture of a new class of computers, optimized for lattice QCD calculations, is described. An individual node is based on a single integrated circuit containing a PowerPC 32-bit integer processor with a 1 Gflops 64-bit IEEE floating point unit, 4 Mbyte of memory, 8 Gbit/sec nearest-neighbor communications and additional control and diagnostic circuitry. The machine's name, QCDOC, derives from ``QCD On a Chip''.