• We study the finite-volume correction on the hadronic vacuum polarization contribution to the muon g-2 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 clover-Wilson fermion action and stout-smeared gluon action at a single lattice cut-off $a^{-1}=2.33$ GeV. The calculation is performed for the quark-connected light-quark 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 quark-mass 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 (space-like) 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 3-family CKM unitarity. Advantages of this new approach over the conventional hadronic $\tau$ decay determination employing flavor-breaking sum rules are also discussed.
  • In a large momentum nucleon state, the matrix element of a gauge-invariant Euclidean Wilson line operator accessible from lattice QCD can be related to the standard light-cone parton distribution function through the large-momentum effective theory (LaMET) expansion. This relation is given by a factorization theorem with a non-trivial matching coefficient. Using the operator product expansion we prove the large-momentum factorization of the quasi-parton distribution function in LaMET, and show that the more recently discussed Ioffe-time distribution approach also obeys an equivalent factorization theorem. Explicit results for the coefficients are obtained and compared at one-loop. 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 All-mode-averaging, is achieved for relevant form factor to proton (and also neutron) decay on the gauge ensemble of $N_f=2+1$ domain-wall 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 10--15\% from 30--40\% for $p\rightarrow\pi e^+$ or from 20--40\% for $p\rightarrow K \bar\nu$. The accuracy of the low energy constants $\alpha$ and $\beta$ in the leading-order 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 three-point 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 2--3 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 light-by-light 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 pre-computed 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 chromo-EDM, CP-violating quark-gluon coupling. Using chirally-symmetric 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.
  • Semi-leptonic $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$ semi-leptonic decays involving $b\to c$ transitions. Our calculations are based on RBC-UKQCD's 2+1 flavor domain-wall 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 domain-wall kernel, and use optimised M\"obius domain-wall fermions for charm quarks.
  • We report our recent lattice calculation of hadronic light-by-light contribution to muon $g-2$ 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 1-2 years.
  • We report a lattice QCD calculation of the hadronic light-by-light contribution to the muon anomalous magnetic moment at physical pion mass. The calculation includes the connected diagrams and the leading, quark-line-disconnected diagrams. We incorporate algorithmic improvements developed in our previous work. The calculation was performed on the $48^3 \times 96$ ensemble generated with a physical-pion-mass 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 finite-volume and finite lattice-spacing errors could be quite large and are the subject of on-going 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 computation-intensive 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 R-Stream source-to-source compiler. Our initial target platform is the Intel PC clusters. We discuss the optimization strategies involved before and after the automatic code generation with R-Stream 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 domain-wall fermions. The neutron and proton EDM form factors are extracted from three-point functions at the next-to-leading 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 I-DSDR gauge action, all generated by the RBC and UKQCD collaborations. The all-mode-averaging 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 two-point correlation function with local topological distribution.
  • For reliable comparison of the standard model prediction to the muon g-2 with its experimental value, the hadronic light-by-light scattering (HLbL) contribution must be calculated by lattice QCD simulation. HLbL contribution has many types of disconnected-type 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 first-principles 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 trend-line, 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 (HEP-FCE) initiated a roadmap planning activity with two key overlapping drivers -- 1) software effectiveness, and 2) infrastructure and expertise advancement. The HEP-FCE 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 quark-connected part of the hadronic light-by-light 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 near-physical, $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, $g-2$, 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 well-suited to provide a first-principle'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 Monte-Carlo 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 all-mode 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 domain-wall fermions. This comparison indicates that AMA significantly reduces statistical errors in Monte-Carlo 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 domain-wall 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, link-smearings are used to improve the signal-to-noise ratio. We employ two kinds of link-smearings, HYP1 and HYP2, and their results are combined in taking the continuum limit. For the matching between the lattice and the continuum theory, one-loop 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 finite-volume QCD to infinite-volume QED by an appropriate twist-averaging procedure. We demonstrate the prescription numerically for the leading-order hadronic contribution to the anomalous magnetic moment of the muon and the electro-magnetic pion mass splitting.
  • We calculate the B-meson decay constants f_B, f_Bs, and their ratio in unquenched lattice QCD using domain-wall light quarks and relativistic b-quarks. We use gauge-field ensembles generated by the RBC and UKQCD collaborations using the domain-wall 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 light-quark masses and continuum. For the b-quarks we use the anisotropic clover action with the relativistic heavy-quark interpretation, such that discretization errors from the heavy-quark action are of the same size as from the light-quark sector. We renormalize the lattice heavy-light axial-vector 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 heavy-light current through O(alpha_s a). We extrapolate our results to the physical light-quark masses and continuum using SU(2) heavy-meson 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 three-flavor determinations of $B$-meson decay constants using staggered light quarks.
  • The form factor that yields the light-by-light scattering contribution to the muon anomalous magnetic moment is computed in lattice QCD+QED and QED. A non-perturbative 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_L-K_S$ mixing matrix receive contributions from all three charge-2/3 quarks: up, charm and top. These give both short- and long-distance 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 sub-percent 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 domain-wall 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, link-smearings are used to improve the signal-to-noise ratio. We employ two kinds of link-smearings and their results are combined in taking a continuum limit. For the matching between the lattice and the continuum theory, one-loop 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, non-degenerate, 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. Non-compact 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$ domain-wall fermion. All-mode 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 low-mode averaging (LMA). We adopt this technique to hadron two-point functions and three-point functions, and compare with LMA and traditional source-shift method in the same ensembles. We observe AMA is much more cost effective in reducing statistical error for these observables.