R.D. Kenway

The International Exascale Software Project roadmap

Over the last 20 years, the open-source community has provided more and more software on which the world’s high-performance computing systems depend for performance and productivity. The community has invested millions of dollars and years of effort to build key components. However, although the investments in these separate software elements have been tremendously valuable, a great deal of productivity has also been lost because of the lack of planning, coordination, and key integration of technologies necessary to make them work together smoothly and efficiently, both within individual petascale systems and between different systems. It seems clear that this completely uncoordinated development model will not provide the software needed to support the unprecedented parallelism required for peta/ exascale computation on millions of cores, or the flexibility required to exploit new hardware models and features, such as transactional memory, speculative execution, and graphics processing units. This report describes the work of the community to prepare for the challenges of exascale computing, ultimately combing their efforts in a coordinated International Exascale Software Project.

Physical results from 2+1 flavor domain wall QCD and SU(2) chiral perturbation theory

We have simulated QCD using 2+1 flavors of domain wall quarks on a (2.74fm)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 BK 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 ? baryon, and the ? and K mesons to set the lattice scale and determine the quark masses. We then find f?=124.1(3.6)stat(6.9)systMeV, fK=149.6(3.6)stat(6.3)systMeV and fK/f?=1.205(0.018)stat(0.062)syst. Using non-perturbative renormalization to relate lattice regularized quark masses to RI-MOM masses, and perturbation theory to relate these to MS¯ we find mMS¯ud(2GeV)=3.72(0.16)stat(0.33)ren(0.18)systMeV and mMS¯s(2GeV)=107.3(4.4)stat(9.7)ren(4.9)systMeV.

Effects of nonperturbatively improved dynamical fermions in QCD at fixed lattice spacing

We present results for the static inter-quark potential, lightest glueballs, light hadron spectrum and topological susceptibility using a non-perturbatively improved action on a

Domain wall QCD with physical quark masses

16^3\times 32

Monte Carlo renormalisation group studies of SU(3) lattice gauge theory

lattice at a set of values of the bare gauge coupling and bare dynamical quark mass chosen to keep the lattice size fixed in physical units (

Gauge-invariant smearing and matrix correlators using Wilson fermions at beta =6.2.

\sim 1.7

2 + 1 flavor domain wall QCD on a (2 fm)3 lattice : Light meson spectroscopy with Ls = 16

fm). By comparing these measurements with a matched quenched ensemble, we study the effects due to two degenerate flavours of dynamical quarks. With the greater control over residual lattice spacing effects which these methods afford, we find some evidence of charge screening and some minor effects on the light hadron spectrum over the range of quark masses studied (

Hadron Mass Calculations With Susskind and Wilson Fermions in the Fundamental Adjoint Plane

M_{PS}/M_{V}\ge0.58

Kl3 semileptonic form factor from 2 + 1 flavor lattice QCD

). More substantial differences between quenched and unquenched simulations are observed in measurements of topological quantities.

Neutral-Kaon Mixing from (2 + 1)-Flavor Domain-Wall QCD

We present results for several light hadronic quantities ( f π , f K , B K , m ud , m s , t 0 ½, w 0 ) obtained from simulations of 2+1 flavor domain wall lattice QCD with large physical volumes and nearly physical pion masses at two lattice spacings. We perform a short, O (3) %, extrapolation in pion mass to the physical values by combining our new data in a simultaneous chiral/continuum “global fit” with a number of other ensembles with heavier pion masses. We use the physical values of m π , m K and m Ω to determine the two quark masses and the scale—all other quantities are outputs from our simulations. We obtain results with subpercent statistical errors and negligible chiral and finite-volume systematics for these light hadronic quantities, including f π = 130.2(9) MeV; f K = 155.5(8) MeV; the average up/down quark mass and strange quark mass in the ‾MS scheme at 3 GeV, 2.997(49) and 81.64(1.17) MeV respectively; and the neutral kaon mixing parameter, B K , in the renormalization group invariant scheme, 0.750(15) and the ‾MS scheme at 3 GeV, 0.530(11).