
We study the finitevolume correction on the hadronic vacuum polarization
contribution to the muon g2 in lattice QCD at physical pion mass using two
different volumes: $(5.4~{\rm fm})^4$ and $(8.1~{\rm fm})^4$. We use an
optimized AMA technique for noise reduction on $N_f=2+1$ PACS gauge
configurations with cloverWilson fermion action and stoutsmeared gluon action
at a single lattice cutoff $a^{1}=2.33$ GeV. The calculation is performed for
the quarkconnected lightquark contribution in the isospin symmetric limit. We
take into account the effects of backward propagating states by extending a
temporal boundary condition. In addition we study a quarkmass correction to
tune to the exactly same physical pion mass on different volume and compare its
correction with chiral perturbation.

We propose and apply a new approach to determining $V_{us}$ using
dispersion relations with weight functions having poles at Euclidean
(spacelike) momentum which relate strange hadronic $\tau$ decay distributions
to hadronic vacuum polarization functions (HVPs) obtained from lattice QCD. We
show examples where spectral integral contributions from the region where
experimental data have large errors or do not exist are strongly suppressed but
accurate determinations of the relevant lattice HVP combinations remain
possible. The resulting $V_{us}$ agrees well with determinations from $K$
physics and 3family CKM unitarity. Advantages of this new approach over the
conventional hadronic $\tau$ decay determination employing flavorbreaking sum
rules are also discussed.

In a large momentum nucleon state, the matrix element of a gaugeinvariant
Euclidean Wilson line operator accessible from lattice QCD can be related to
the standard lightcone parton distribution function through the largemomentum
effective theory (LaMET) expansion. This relation is given by a factorization
theorem with a nontrivial matching coefficient. Using the operator product
expansion we prove the largemomentum factorization of the quasiparton
distribution function in LaMET, and show that the more recently discussed
Ioffetime distribution approach also obeys an equivalent factorization
theorem. Explicit results for the coefficients are obtained and compared at
oneloop. Our proof clearly demonstrates that the matching coefficients in the
$\overline{\rm MS}$ scheme depend on the large partonic momentum rather than
the nucleon momentum.

We present an improved result of lattice computation of the proton decay
matrix elements in $N_f=2+1$ QCD. In this study, the significant improvement of
statistical accuracy by adopting the error reduction technique of
Allmodeaveraging, is achieved for relevant form factor to proton (and also
neutron) decay on the gauge ensemble of $N_f=2+1$ domainwall fermions in
$m_\pi=0.34$0.69 GeV on 2.7~fm$^3$ lattice as used in our previous work
\cite{Aoki:2013yxa}. We improve total accuracy of matrix elements to 1015\%
from 3040\% for $p\rightarrow\pi e^+$ or from 2040\% for $p\rightarrow K
\bar\nu$. The accuracy of the low energy constants $\alpha$ and $\beta$ in the
leadingorder baryon chiral perturbation theory (BChPT) of proton decay are
also improved. The relevant form factors of $p\rightarrow \pi$ estimated
through the "direct" lattice calculation from threepoint function appear to be
1.4 times smaller than those from the "indirect" method using BChPT with
$\alpha$ and $\beta$. It turns out that the utilization of our result will
provide a factor 23 larger proton partial lifetime than that obtained using
BChPT. We also discuss the use of these parameters in a dark matter model.

In our previous work, the connected and leading disconnected hadronic
lightbylight contributions to the muon anomalous magnetic moment (g  2) have
been computed using lattice QCD ensembles corresponding to physical pion mass
generated by the RBC/UKQCD collaboration. However, the calculation is expected
to suffer from a significant finite volume error that scales like $1/L^2$ where
$L$ is the spatial size of the lattice. In this paper, we demonstrate that this
problem is cured by treating the muon and photons in infinite volume, continuum
QED, resulting in a weighting function that is precomputed and saved with
affordable cost and sufficient accuracy. We present numerical results for the
case when the quark loop is replaced by a muon loop, finding the expected
exponential approach to the infinite volume limit and consistency with the
known analytic result. We have implemented an improved weighting function which
reduces both discretization and finite volume effects arising from the hadronic
part of the amplitude.

We present initial results of computing nucleon electric dipole moment
induced by quark chromoEDM, CPviolating quarkgluon coupling. Using
chirallysymmetric domain wall and M\"obius fermions with pion mass
$m_\pi=172\text{ MeV}$, we calculate the connected part of the electric dipole
form factor $F_3(Q^2)$. In addition, we perform an exploratory study of the
method to calculate EDM using uniform background electric field on a lattice
introduced without breaking the periodicity in the time direction.

Semileptonic $B$ decays provide promising channels to test the Standard
Model, search for signs of new physics, or determine fundamental parameters
like CKM matrix elements. We present an update on our calculation of short
distance contributions to GIM suppressed rare $B$ decays focusing in particular
on $B_s\to \phi \ell^+ \ell^$ decays. Furthermore we show first results for
our calculation of $B_{(s)}\to D_{(s)}^{(*)}\ell\nu$ semileptonic decays
involving $b\to c$ transitions. Our calculations are based on RBCUKQCD's 2+1
flavor domainwall fermion and Iwasaki gauge field configurations featuring
three lattice spacings in the range $1.73$ GeV $\le a^{1} \le 2.77$ GeV and
pion masses down to the physical value. We calculate the form factors by
simulating $b$quarks using the relativistic heavy quark action, create light
$u/d$ and $s$ quarks with standard domainwall kernel, and use optimised
M\"obius domainwall fermions for charm quarks.

We report our recent lattice calculation of hadronic lightbylight
contribution to muon $g2$ using our recently developed moment method. The
connected diagrams and the leading disconnected diagrams are included. The
calculation is performed on a $48^3 \times 96$ lattice with physical pion mass
and 5.5 fm box size. We expect sizable finite volume and finite lattice spacing
corrections to the results of these calculations which will be estimated in
calculations to be carried out over the next 12 years.

We report a lattice QCD calculation of the hadronic lightbylight
contribution to the muon anomalous magnetic moment at physical pion mass. The
calculation includes the connected diagrams and the leading,
quarklinedisconnected diagrams. We incorporate algorithmic improvements
developed in our previous work. The calculation was performed on the $48^3
\times 96$ ensemble generated with a physicalpionmass and a 5.5 fm spatial
extent by the RBC and UKQCD collaborations using the chiral, domain wall
fermion (DWF) formulation. We find $a_\mu^{\text{HLbL}} = 5.35 (1.35) \times
10^{ 10}$, where the error is statistical only. The finitevolume and finite
latticespacing errors could be quite large and are the subject of ongoing
research. The omitted disconnected graphs, while expected to give a correction
of order 10\%, also need to be computed.

The application of the Dirac operator on a spinor field, the Dslash
operation, is the most computationintensive part of the lattice QCD
simulations. It is often the key kernel to optimize to achieve maximum
performance on various platforms. Here we report on a project to optimize the
domain wall fermion Dirac operator in Columbia Physics System (CPS) using the
RStream sourcetosource compiler. Our initial target platform is the Intel PC
clusters. We discuss the optimization strategies involved before and after the
automatic code generation with RStream and present some preliminary benchmark
results.

We present a lattice calculation of the neutron and proton electric dipole
moments (EDM's) with $N_f=2+1$ flavors of domainwall fermions. The neutron and
proton EDM form factors are extracted from threepoint functions at the
nexttoleading order in the $\theta$ vacuum of QCD. In this computation, we
use pion masses 0.33 and 0.42 GeV and 2.7 fm$^3$ lattices with Iwasaki gauge
action and a 0.17 GeV pion and 4.6 fm$^3$ lattice with IDSDR gauge action, all
generated by the RBC and UKQCD collaborations. The allmodeaveraging technique
enables an efficient and high statistics calculation. Chiral behavior of
lattice EDM's is discussed in the context of baryon chiral perturbation theory.
In addition, we also show numerical evidence on relationship of three and
twopoint correlation function with local topological distribution.

For reliable comparison of the standard model prediction to the muon g2 with
its experimental value, the hadronic lightbylight scattering (HLbL)
contribution must be calculated by lattice QCD simulation. HLbL contribution
has many types of disconnectedtype diagrams. Here, we start with recalling the
point that must be taken care of in every method to calculate them by lattice
QCD, and present one concrete method called nonperturbative QED method.

The current measurement of muonic $g  2$ disagrees with the theoretical
calculation by about 3 standard deviations. Hadronic vacuum polarization (HVP)
and hadronic light by light (HLbL) are the two types of processes that
contribute most to the theoretical uncertainty. The current value for HLbL is
still given by models. I will describe results from a firstprinciples lattice
calculation with a 139 MeV pion in a box of 5.5 fm extent. Our current
numerical strategies, including noise reduction techniques, evaluating the HLbL
amplitude at zero external momentum transfer, and important remaining
challenges, in particular those associated with finite volume effects, will be
discussed.

Computing plays an essential role in all aspects of high energy physics. As
computational technology evolves rapidly in new directions, and data throughput
and volume continue to follow a steep trendline, it is important for the HEP
community to develop an effective response to a series of expected challenges.
In order to help shape the desired response, the HEP Forum for Computational
Excellence (HEPFCE) initiated a roadmap planning activity with two key
overlapping drivers  1) software effectiveness, and 2) infrastructure and
expertise advancement. The HEPFCE formed three working groups, 1) Applications
Software, 2) Software Libraries and Tools, and 3) Systems (including systems
software), to provide an overview of the current status of HEP computing and to
present findings and opportunities for the desired HEP computational roadmap.
The final versions of the reports are combined in this document, and are
presented along with introductory material.

The quarkconnected part of the hadronic lightbylight scattering
contribution to the muon's anomalous magnetic moment is computed using lattice
QCD with chiral fermions. We report several significant algorithmic
improvements and demonstrate their effectiveness through specific calculations
which show a reduction in statistical errors by more than an order of
magnitude. The most realistic of these calculations is performed with a
nearphysical, $171$ MeV pion mass on a $(4.6\;\mathrm{fm})^3$ spatial volume
using the $32^3\times 64$ Iwasaki+DSDR gauge ensemble of the RBC/UKQCD
Collaboration.

The anomalous magnetic moment of muon, $g2$, is a very precisely measured
quantity. However, the current measurement disagrees with standard model by
about 3 standard deviations. Hadronic vacuum polarization and hadronic light by
light are the two types of processes that contribute most to the theoretical
uncertainty. I will describe how lattice methods are wellsuited to provide a
firstprinciple's result for the hadronic light by light contribution, the
various numerical strategies that are presently being used to evaluate it, our
current results and the important remaining challenges which must be overcome.

We present a new class of statistical error reduction techniques for
MonteCarlo simulations. Using covariant symmetries, we show that correlation
functions can be constructed from inexpensive approximations without
introducing any systematic bias in the final result. We introduce a new class
of covariant approximation averaging techniques, known as allmode averaging
(AMA), in which the approximation takes account of contributions of all
eigenmodes through the inverse of the Dirac operator computed from the
conjugate gradient method with a relaxed stopping condition. In this paper we
compare the performance and computational cost of our new method with
traditional methods using correlation functions and masses of the pion,
nucleon, and vector meson in $N_f=2+1$ lattice QCD using domainwall fermions.
This comparison indicates that AMA significantly reduces statistical errors in
MonteCarlo calculations over conventional methods for the same cost.

Neutral $B$ meson mixing matrix elements and $B$ meson decay constants are
calculated. Static approximation is used for $b$ quark and domainwall fermion
formalism is employed for light quarks. The calculations are carried out on
$2+1$ flavor dynamical ensembles generated by RBC/UKQCD Collaborations with
lattice spacings $0.086$fm ($a^{1}\sim 2.3$GeV) and $0.11$fm ($1.7$GeV), and a
fixed physical spatial volume of about $(2.7{\rm fm})^3$. In the static quark
action, linksmearings are used to improve the signaltonoise ratio. We employ
two kinds of linksmearings, HYP1 and HYP2, and their results are combined in
taking the continuum limit. For the matching between the lattice and the
continuum theory, oneloop perturbative $O(a)$ improvements are made to reduce
discretization errors. As the most important quantity of this work, we obtain
SU(3) breaking ratio $\xi=1.208(60)$, where the error includes statistical and
systematic one. (Uncertainty from infinite $b$ quark mass is not included.) We
also find other neutral $B$ meson mixing quantities
$f_B\sqrt{\hat{B}_B}=240(22)$MeV, $f_{B_s}\sqrt{\hat{B}_{B_s}}=290(22)$MeV,
$\hat{B}_B=1.17(22)$, $\hat{B}_{B_s}=1.22(13)$ and $B_{B_s}/B_B=1.028(74)$, $B$
meson decay constants $f_B=219(17)$MeV, $f_{B_s}=264(19)$MeV and
$f_{B_s}/f_B=1.193(41)$, in the static limit of $b$ quark, which do not include
infinite $b$ quark mass uncertainty.

We present a method to couple finitevolume QCD to infinitevolume QED by an
appropriate twistaveraging procedure. We demonstrate the prescription
numerically for the leadingorder hadronic contribution to the anomalous
magnetic moment of the muon and the electromagnetic pion mass splitting.

We calculate the Bmeson decay constants f_B, f_Bs, and their ratio in
unquenched lattice QCD using domainwall light quarks and relativistic
bquarks. We use gaugefield ensembles generated by the RBC and UKQCD
collaborations using the domainwall fermion action and Iwasaki gauge action
with three flavors of light dynamical quarks. We analyze data at two lattice
spacings of a ~ 0.11, 0.086 fm with unitary pion masses as light as M_pi ~ 290
MeV; this enables us to control the extrapolation to the physical lightquark
masses and continuum. For the bquarks we use the anisotropic clover action
with the relativistic heavyquark interpretation, such that discretization
errors from the heavyquark action are of the same size as from the lightquark
sector. We renormalize the lattice heavylight axialvector current using a
mostly nonperturbative method in which we compute the bulk of the matching
factor nonperturbatively, with a small correction, that is close to unity, in
lattice perturbation theory. We also improve the lattice heavylight current
through O(alpha_s a). We extrapolate our results to the physical lightquark
masses and continuum using SU(2) heavymeson chiral perturbation theory, and
provide a complete systematic error budget. We obtain f_B0 = 199.5(12.6) MeV,
f_B+ = 195.6(14.9) MeV, f_Bs = 235.4(12.2) MeV, f_Bs/f_B0 = 1.197(50), and
f_Bs/f_B+ = 1.223(71), where the errors are statistical and total systematic
added in quadrature. These results are in good agreement with other published
results and provide an important independent cross check of other threeflavor
determinations of $B$meson decay constants using staggered light quarks.

The form factor that yields the lightbylight scattering contribution to the
muon anomalous magnetic moment is computed in lattice QCD+QED and QED. A
nonperturbative treatment of QED is used and is checked against perturbation
theory. The hadronic contribution is calculated for unphysical quark and muon
masses, and only the diagram with a single quark loop is computed.
Statistically significant signals are obtained. Initial results appear
promising, and the prospect for a complete calculation with physical masses and
controlled errors is discussed.

The real and imaginary parts of the $K_LK_S$ mixing matrix receive
contributions from all three charge2/3 quarks: up, charm and top. These give
both short and longdistance contributions which are accessible through a
combination of perturbative and lattice methods. We will discuss a strategy to
compute both the mass difference, $\Delta M_K$ and $\epsilon_K$ to subpercent
accuracy, looking in detail at the contributions from each of the three CKM
matrix element products $V_{id}^*V_{is}$ for $i=u, c$ and $t$ as described in
Ref. [1]

Neutral B meson mixing matrix elements and B meson decay constants are
calculated. Static approximation is used for b quark and domainwall fermion
formalism is employed for light quarks. The calculations are done on 2+1 flavor
dynamical ensembles, whose lattice spacings are 0.086 fm and 0.11 fm with a
fixed physical spatial volume of about (2.7 fm)^3. In the static quark action,
linksmearings are used to improve the signaltonoise ratio. We employ two
kinds of linksmearings and their results are combined in taking a continuum
limit. For the matching between the lattice and the continuum theory, oneloop
perturbative calculations are used including O(a) improvements to reduce
discretization errors. We obtain SU(3) braking ratio \xi=1.222(60) in the
static limit of b quark.

We report on a calculation of the effects of isospin breaking in Lattice
QCD+QED. This involves using Chiral Perturbation Theory with Electromagnetic
corrections to find the renormalized, nondegenerate, light quark masses. The
calculations are carried out on QCD ensembles generated by the RBC and UKQCD
collaborations using Domain Wall Fermions and the Iwasaki and Iwasaki+DSDR
Gauge Actions with unitary pion masses down to 170 MeV. Noncompact QED is
treated in the quenched approximation. The simulations use a $32^3$ lattice
size with $a^{1}=2.28(3)$ GeV (Iwasaki) and 1.37(1) (Iwasaki+DSDR). This
builds on previous work from the RBC/UKQCD collaboration with lattice spacing
$a^{1}=1.78(4)$ GeV.

We demonstrate the new class of variance reduction techniques for hadron
propagator and nucleon isovector form factor in the realistic lattice of
$N_f=2+1$ domainwall fermion. Allmode averaging (AMA) is one of the powerful
tools to reduce the statistical noise effectively for wider varieties of
observables compared to existing techniques such as lowmode averaging (LMA).
We adopt this technique to hadron twopoint functions and threepoint
functions, and compare with LMA and traditional sourceshift method in the same
ensembles. We observe AMA is much more cost effective in reducing statistical
error for these observables.