
We present a firstprinciples lattice QCD+QED calculation at physical pion
mass of the leadingorder hadronic vacuum polarization contribution to the muon
anomalous magnetic moment. The total contribution of up, down, strange, and
charm quarks including QED and strong isospin breaking effects is found to be
$a_\mu^{\rm HVP~LO}=715.4(16.3)(9.2) \times 10^{10}$, where the first error is
statistical and the second is systematic. By supplementing lattice data for
very short and long distances with experimental Rratio data using the
compilation of Ref. [1], we significantly improve the precision of our
calculation and find $a_\mu^{\rm HVP~LO} = 692.5(1.4)(0.5)(0.7)(2.1) \times
10^{10}$ with lattice statistical, lattice systematic, Rratio statistical,
and Rratio systematic errors given separately. This is the currently most
precise determination of the leadingorder hadronic vacuum polarization
contribution to the muon anomalous magnetic moment. In addition, we present the
first lattice calculation of the lightquark QED correction at physical pion
mass.

We analyze commonly used expressions for computing the nucleon electric
dipole form factors (EDFF) $F_3$ and moments (EDM) on a lattice and find that
they lead to spurious contributions from the Pauli form factor $F_2$ due to
inadequate definition of these form factors when parity mixing of lattice
nucleon fields is involved. Using chirally symmetric domain wall fermions, we
calculate the proton and the neutron EDFF induced by the CPviolating quark
chromoEDM interaction using the corrected expression. In addition, we
calculate the electric dipole moment of the neutron using background electric
field that respects time translation invariance and boundary conditions, and
find that it decidedly agrees with the new formula but not the old formula for
$F_3$. Finally, we analyze some selected lattice results for the nucleon EDM
and observe that after the correction is applied, they either agree with zero
or are substantially reduced in magnitude, thus reconciling their difference
from phenomenological estimates of the nucleon EDM.

We review lattice results related to pion, kaon, D and Bmeson physics with
the aim of making them easily accessible to the particle physics community.
More specifically, we report on the determination of the lightquark masses,
the form factor f+(0), arising in the semileptonic K > pi transition at zero
momentum transfer, as well as the decay constant ratio fK/fpi and its
consequences for the CKM matrix elements Vus and Vud. Furthermore, we describe
the results obtained on the lattice for some of the lowenergy constants of
SU(2)LxSU(2)R and SU(3)LxSU(3)R Chiral Perturbation Theory. We review the
determination of the BK parameter of neutral kaon mixing as well as the
additional four B parameters that arise in theories of physics beyond the
Standard Model. The latter quantities are an addition compared to the previous
review. For the heavyquark sector, we provide results for mc and mb (also new
compared to the previous review), as well as those for D and Bmeson decay
constants, form factors, and mixing parameters. These are the heavyquark
quantities most relevant for the determination of CKM matrix elements and the
global CKM unitaritytriangle fit. Finally, we review the status of lattice
determinations of the strong coupling constant alpha_s.

We present results for the leading hadronic contribution to the muon
anomalous magnetic moment due to strange quarkconnected vacuum polarisation
effects. Simulations were performed using RBCUKQCD's $N_f=2+1$ domain wall
fermion ensembles with physical light sea quark masses at two lattice spacings.
We consider a large number of analysis scenarios in order to obtain solid
estimates for residual systematic effects. Our final result in the continuum
limit is $a_\mu^{(2)\,{\rm
had},\,s}=53.1(9)\left(^{+1}_{3}\right)\times10^{10}$.

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
ydirection. 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
standardmodel theory of CP violation, one which can be made more accurate with
increasing computer capability.

We report the first lattice QCD calculation of the hadronic vacuum
polarization disconnected contribution to the muon anomalous magnetic moment at
physical pion mass. The calculation uses a refined noisereduction technique
which enabled the control of statistical uncertainties at the desired level
with modest computational effort. Measurements were performed on the $48^3
\times 96$ physicalpionmass lattice generated by the RBC and UKQCD
collaborations. We find $a_\mu^{\rm HVP~(LO)~DISC} = 9.6(3.3)(2.3)\times
10^{10}$, where the first error is statistical and the second systematic.

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 twopion 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 report initial nucleon structure results computed on lattices with 2+1
dynamical M\"obius domain wall fermions at the physical point generated by the
RBC and UKQCD collaborations. At this stage, we evaluate only connected quark
contributions. In particular, we discuss the nucleon vector and axialvector
form factors, nucleon axial charge and the isovector quark momentum fraction.
From currently available statistics, we estimate the stochastic accuracy of the
determination of $g_A$ and $<x>_{ud}$ to be around 10%, and we expect to
reduce that to 5% within the next year. To reduce the computational cost of our
calculations, we extensively use acceleration techniques such as loweigenmode
deflation and allmodeaveraging (AMA). We present a method for choosing
optimal AMA parameters.

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
partiallyquenched 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 strangequark mass in the MSbarscheme 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
MSbarscheme 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}.

These reports present the results of the 2013 Community Summer Study of the
APS Division of Particles and Fields ("Snowmass 2013") on the future program of
particle physics in the U.S. Chapter 9, on Computing, discusses the computing
challenges for future experiments in the Energy, Intensity, and Cosmic
Frontiers, for accelerator science, and for particle theory, as well as
structural issues in supporting the intense uses of computing required in all
areas of particle physics.

This is the report of the Intensity Frontier Charged Lepton Working Group of
the 2013 Community Summer Study "Snowmass on the Mississippi", summarizing the
current status and future experimental opportunities in muon and tau lepton
studies and their sensitivity to new physics. These include searches for
charged lepton flavor violation, measurements of magnetic and electric dipole
moments, and precision measurements of the decay spectrum and parityviolating
asymmetries.

This is the report of the Computing Frontier working group on Lattice Field
Theory prepared for the proceedings of the 2013 Community Summer Study
("Snowmass"). We present the future computing needs and plans of the U.S.
lattice gauge theory community and argue that continued support of the U.S.
(and worldwide) latticeQCD effort is essential to fully capitalize on the
enormous investment in the highenergy physics experimental program. We first
summarize the dramatic progress of numerical latticeQCD simulations in the
past decade, with some emphasis on calculations carried out under the auspices
of the U.S. LatticeQCD Collaboration, and describe a broad program of
latticeQCD calculations that will be relevant for future experiments at the
intensity and energy frontiers. We then present details of the computational
hardware and software resources needed to undertake these calculations.

After a brief selfcontained introduction to the muon anomalous magnetic
moment, (g2), we review the status of lattice calculations of the hadronic
vacuum polarization contribution and present first results from lattice QCD for
the hadronic lightbylight scattering contribution. The signal for the latter
is consistent with model calculations. While encouraging, the statistical error
is large and systematic errors are mostly uncontrolled. The method is applied
first to pure QED as a check.

We describe the computation of the amplitude A_2 for a kaon to decay into two
pions with isospin I=2. The results presented in the letter Phys.Rev.Lett. 108
(2012) 141601 from an analysis of 63 gluon configurations are updated to 146
configurations giving Re$A_2=1.381(46)_{\textrm{stat}}(258)_{\textrm{syst}}
10^{8}$ GeV and Im$A_2=6.54(46)_{\textrm{stat}}(120)_{\textrm{syst}}10^{13}$
GeV. Re$A_2$ is in good agreement with the experimental result, whereas the
value of Im$A_2$ was hitherto unknown. We are also working towards a direct
computation of the $K\to(\pi\pi)_{I=0}$ amplitude $A_0$ but, within the
standard model, our result for Im$A_2$ can be combined with the experimental
results for Re$A_0$, Re$A_2$ and $\epsilon^\prime/\epsilon$ to give
Im$A_0/$Re$A_0= 1.61(28)\times 10^{4}$ . Our result for Im\,$A_2$ implies
that the electroweak penguin (EWP) contribution to $\epsilon^\prime/\epsilon$
is Re$(\epsilon^\prime/\epsilon)_{\mathrm{EWP}} = (6.25 \pm
0.44_{\textrm{stat}} \pm 1.19_{\textrm{syst}}) \times 10^{4}$.

The Proceedings of the 2011 workshop on Fundamental Physics at the Intensity
Frontier. Science opportunities at the intensity frontier are identified and
described in the areas of heavy quarks, charged leptons, neutrinos, proton
decay, new light weaklycoupled particles, and nucleons, nuclei, and atoms.

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 \pimesons
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 socalled 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
CPviolation in kaon decays.

We report a direct lattice calculation of the $K$ to $\pi\pi$ decay matrix
elements for both the $\Delta I=1/2$ and 3/2 amplitudes $A_0$ and $A_2$ on 2+1
flavor, domain wall fermion, $16^3\times32\times16$ lattices. This is a
complete calculation in which all contractions for the required ten, fourquark
operators are evaluated, including the disconnected graphs in which no quark
line connects the initial kaon and final twopion states. These lattice
operators are nonperturbatively renormalized using the RomeSouthampton method
and the quadratic divergences are studied and removed. This is an important but
notoriously difficult calculation, requiring high statistics on a large volume.
In this paper we take a major step towards the computation of the physical
$K\to\pi\pi$ amplitudes by performing a complete calculation at unphysical
kinematics with pions of mass 422\,MeV at rest in the kaon rest frame. With
this simplification we are able to resolve Re$(A_0)$ from zero for the first
time, with a 25% statistical error and can develop and evaluate methods for
computing the complete, complex amplitude $A_0$, a calculation central to
understanding the $\Delta =1/2$ rule and testing the standard model of CP
violation in the kaon system.

We determine the neutral kaon mixing matrix element $B_K$ in the continuum
limit with 2+1 flavors of domain wall fermions, using the Iwasaki gauge action
at two different lattice spacings. These lattice fermions have near exact
chiral symmetry and therefore avoid artificial lattice operator mixing.
We introduce a significant improvement to the conventional NPR method in
which the bare matrix elements are renormalized nonperturbatively in the
RIMOM scheme and are then converted into the MSbar scheme using continuum
perturbation theory. In addition to RIMOM, we introduce and implement four
nonexceptional intermediate momentum schemes that suppress infrared
nonperturbative uncertainties in the renormalization procedure. We compute the
conversion factors relating the matrix elements in this family of RISMOM
schemes and MSbar at oneloop order. Comparison of the results obtained using
these different intermediate schemes allows for a more reliable estimate of the
unknown higherorder contributions and hence for a correspondingly more robust
estimate of the systematic error. We also apply a recently proposed approach in
which twisted boundary conditions are used to control the Symanzik expansion
for offshell vertex functions leading to a better control of the
renormalization in the continuum limit.
We control chiral extrapolation errors by considering both the NLO SU(2)
chiral effective theory, and an analytic mass expansion. We obtain
$B_K^{\msbar}(3 GeV) = 0.529(5)_{stat}(15)_\chi(2)_{FV}(11)_{NPR}$. This
corresponds to $\hat{B}_K = 0.749(7)_{stat}(21)_\chi(3)_{FV}(15)_{NPR}$. Adding
all sources of error in quadrature we obtain $\hat{B}_K =
0.749(27)_{combined}$, with an overall combined error of 3.6%.

We report on an investigation of loading characteristics of deep microscopic
dipole traps. The dipole trap is loaded from a low density magneto optical trap
(MOT) containing $\approx 5\times 10^{6}$ atoms. We determine the loading
parameters that maximize the trapped atom number for a trap of waist 5 {\mu}m
with trap depths ranging from 3.5 mK to 10 mK. We determine the optimal trap
loading conditions and the loading rates, loss coefficients and temperature of
the trapped atoms under these conditions. We show that it is possible to load a
few hundred to thousand atoms in dipole traps of depth 3.5 mK to 8.5 mK under
the optimal loading conditions.

Results are presented for the electromagnetic mass splittings of the low
lying hadrons.
These are used to determine the nondegenerate light quark masses.
It is found that m_u=2.24(10)(34), m_d=4.65(15)(32), and $m_s=97.6(2.9)(5.5)$
MeV (MSbar scheme, 2 GeV scale). The first error is statistical and the second
systematic. We find the lowest order electromagnetic splitting
(m_pi+m_pi0)_QED=3.38(23) MeV, the splittings including nexttoleading order,
(m_pi+m_pi0)_QED=4.50(23) MeV, (m_K+m_K0)_QED=1.87(10) MeV, and the m_u !=
m_d contribution to the kaon mass difference, (m_K+m_K0)_(m_um_d)=5.840(96)
MeV.
All errors are statistical only, and the nexttoleading order pion splitting
is only approximate; it does not contain all nexttoleading order
contributions. We also computed the protonneutron mass difference, including
for the first time, QED interactions in a realistic 2+1 flavor calculation. We
find $(m_pm_n)_{\rm QED}=0.383(68)$ MeV, (m_pm_n)_(m_um_d)=2.51(14) MeV,
and the total m_pm_n=2.13(16)(70) MeV, where the first error is statistical,
and the second, part of the systematic error. We use domain wall fermions and
the Iwasaki gauge action (gauge coupling beta=2.13). We use two lattice sizes,
16^3 and 24^3, to address finite volume effects. Noncompact QED is treated in
the quenched approximation.
We present new results for the electromagnetic low energy constants in SU(3)
and SU(2) partiallyquenched chiral perturbation theory to the nexttoleading
order, obtained from fits to our data.
Detailed analysis of systematic errors in our results and methods for
improving them are discussed. Finally, new analytic results for SU(2)_L x
SU(2)_Rpluskaon chiral perturbation theory, including the oneloop logs
proportional to alpha_em*m, are given.

The exciting possibility of direct observation of QCD instantons in heavyion
collisions has recently been proposed by Kharzeev. The underlying phenomenon,
known as the chiral magnetic effect, may have been observed recently at RHIC,
and a first principles calculation is needed to confirm and understand the
results. The chiral magnetic effect is thought to be visible in the symmetric
phase, at temperatures above the QCD critical temperature, and in the presence
of an external magnetic field. We report on first 2+1 flavor, domain wall
fermion, QCD+QED dynamical simulations above the critical temperature, in a
fixed topological sector(s), which are used to study the electric charge
separation produced by the effect.

We discuss Quantum Chromodynamics calculations using the lattice regulator.
The theory of the strong force is a cornerstone of the Standard Model of
particle physics. We present USQCD collaboration results obtained on Argonne
National Lab's Intrepid supercomputer that deepen our understanding of these
fundamental theories of Nature and provide critical support to frontier
particle physics experiments and phenomenology.

We present results for the nucleon axial charge g_A at a fixed lattice
spacing of 1/a=1.73(3) GeV using 2+1 flavors of domain wall fermions on size
16^3x32 and 24^3x64lattices (L=1.8 and 2.7 fm) with length 16 in the fifth
dimension. The length of the Monte Carlo trajectory at the lightest m_\pi is
7360 units, including 900 for thermalization. We find finite volume effects are
larger than the pion mass dependence at m_\pi= 330 MeV. We also find that g_A
exhibits a scaling with the single variable m_\pi L which can also be seen in
previous twoflavor domain wall and Wilson fermion calculati ons. Using this
scaling to eliminate the finitevolume effect, we obtain g_A = 1.20(6)(4) at
the physical pion mass, m_\pi = 135 MeV, where the first and second errors are
statistical and systematic. The observed finitevolume scaling also appears in
similar quenched simulations, but disappear when V\ge (2.4 fm)^3. We argue this
is a dynamical quark effect.

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 nexttoleading 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
nonperturbative renormalization to relate lattice regularized quark masses to
RIMOM 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}$.