Search

The separation of the hadronic vacuum polarization (HVP) contribution to the muon anomalous magnetic moment into Euclidean windows allows for a tailored approach to address the different dominant challenges at short, intermediate, and long distances. We present a novel approach to compute the short-distance window without the need for using perturbative QCD. We combine a quenched continuum extrapolation using 18 lattice spacings ($1.6 \,\text{GeV} \lesssim a^{-1} \lesssim 6.1\,\text{GeV}$) with a separate continuum extrapolation of the sea-quark effects. This method allows for the computationally expensive sea-quark effects to be estimated using only a smaller number of ensembles at coarser lattice spacings, while largely confining the logarithmic dependency of the continuum extrapolation to the quenched component., Comment: 19 pages, 20 figures

We present a lattice calculation of the $K\to\pi\pi$ matrix elements and amplitudes with both the $\Delta I = 3/2$ and 1/2 channels and $\varepsilon'$, the measure of direct $CP$ violation. We use periodic boundary conditions (PBC), where the correct kinematics of $K\to\pi\pi$ can be achieved via an excited two-pion final state. To overcome the difficulty associated with the extraction of excited states, our previous work \cite{Bai:2015nea,RBC:2020kdj} successfully employed G-parity boundary conditions, where pions are forced to have non-zero momentum enabling the $I=0$ two-pion ground state to express the on-shell kinematics of the $K\to\pi\pi$ decay. Here instead we overcome the problem using the variational method which allows us to resolve the two-pion spectrum and matrix elements up to the relevant energy where the decay amplitude is on-shell. In this paper we report an exploratory calculation of $K\to\pi\pi$ decay amplitudes and $\varepsilon'$ using PBC on a coarser lattice size of $24^3\times64$ with inverse lattice spacing $a^{-1}=1.023$ GeV and the physical pion and kaon masses. The results are promising enough to motivate us to continue our measurements on finer lattice ensembles in order to improve the precision in the near future.

We demonstrate that gauge-equivariant pooling and unpooling layers can perform as well as traditional restriction and prolongation layers in multigrid preconditioner models for lattice QCD. These layers introduce a gauge degree of freedom on the coarse grid, allowing for the use of explicitly gauge-equivariant layers on the coarse grid. We investigate the construction of coarse-grid gauge fields and study their efficiency in the preconditioner model. We show that a combined multigrid neural network using a Galerkin construction for the coarse-grid gauge field eliminates critical slowing down., Comment: 8 pages, 5 figures

The hadronic light-by-light scattering contribution to the muon anomalous magnetic moment, $(g-2$)/2, is computed in the infinite volume QED framework with lattice QCD. We report $a_\mu^\text{HLbL}=12.47(1.15)(0.99) \times 10^{-10}$ where the first error is statistical and the second systematic. The result is mainly based on the 2+1 flavor M\"obius domain wall fermion ensemble with inverse lattice spacing $a^{-1} = 1.73~\mathrm{GeV}$, lattice size $L=5.5~\mathrm{fm}$, and $m_\pi = 139~\mathrm{MeV}$, generated by the RBC-UKQCD collaborations. The leading systematic error of this result comes from the lattice discretization. This result is consistent with previous determinations., Comment: 29 pages

The two pion channel in Lattice QCD has long been a primary testing ground for studying multiparticle scattering in finite volume QCD. With the development of sophisticated techniques such as distillation, it is possible to carefully study two-pion scattering in order to constrain associated low-energy constants. In this work, correlation functions with multiparticle interpolating operators are constructed to compute pion scattering phase shifts and scattering lengths in the isospin 0 and 2 channels with both sea and valence quarks at physical mass. Contamination from vacuum and thermal contributions are explicitly quantified with dedicated calculations and the results obtained after subtracting these nuisance terms are compared with the traditional correlator time series subtraction method. Two physical point ensembles with different lattice actions are used, and our finest ensemble gives results for scattering lengths and phase shifts consistent with phenomenology to within the reported statistical uncertainty.

We demonstrate that a state-of-the art multi-grid preconditioner can be learned efficiently by gauge-equivariant neural networks. We show that the models require minimal re-training on different gauge configurations of the same gauge ensemble and to a large extent remain efficient under modest modifications of ensemble parameters. We also demonstrate that important paradigms such as communication avoidance are straightforward to implement in this framework., Comment: 12 pages, 12 figures

We present a study of lattice-QCD methods to determine the relevant hadronic form factors for radiative leptonic decays of pseudoscalar mesons. We provide numerical results for $D_s^+ \to \ell^+ \nu \gamma$. Our calculation is performed using a domain-wall action for all quark flavors and on a single RBC/UKQCD lattice gauge-field ensemble. The first part of the study is how to best control two sources of systematic error inherent in the calculation, specifically the unwanted excited states created by the meson interpolating field, and unwanted exponentials in the sum over intermediate states. Using a 3d sequential propagator allows for better control over unwanted exponentials from intermediate states, while using a 4d sequential propagator allows for better control over excited states. We perform individual analyses of the 3d and 4d methods as well as a combined analysis using both methods, and find that the 3d sequential propagator offers good control over both sources of systematic uncertainties for the smallest number of propagator solves. From there, we further improve the use of a 3d sequential propagator by employing an infinite-volume approximation method, which allows us to calculate the relevant form factors over the entire allowed range of photon energies. We then study improvements gained by performing the calculation using a different three-point function, using ratios of three-point functions, averaging over positive and negative photon momentum, and using an improved method for extracting the structure-dependent part of the axial form factor. The optimal combination of methods yields results for the $D_s^+ \to \ell^+ \nu \gamma$ structure-dependent vector and axial form factors in the entire kinematic range with statistical plus fitting uncertainties of order 5%, using 25 gauge configurations with 64 samples per configuration., Comment: 31 pages, 14 figures

A study of two-pion scattering for the isospin channels, $I=0$ and $I=2$, using lattice QCD is presented. M\"obius domain wall fermions on top of the Iwasaki-DSDR gauge action for gluons with periodic boundary conditions are used for the lattice computations which are carried out on two ensembles of gauge field configurations generated by the RBC and UKQCD collaborations with physical masses, inverse lattice spacings of 1.023 and 1.378 GeV, and spatial extents of $L=4.63$ and 4.58 fm, respectively. The all-to-all propagator method is employed to compute a matrix of correlation functions of two-pion operators. The generalized eigenvalue problem (GEVP) is solved for a matrix of correlation functions to extract phase shifts with multiple states, two pions with a non-zero relative momentum as well as two pions at rest. Our results for phase shifts for both $I=0$ and $I=2$ channels are consistent with and the Roy Equation and chiral perturbation theory, though at this preliminary stage our errors for $I=0$ are large. An important outcome of this work is that we are successful in extracting two-pion excited states, which are useful for studying $K\to\pi\pi$ decay, on physical-mass ensembles using GEVP.

We construct an approximate trivializing map by using a Schwinger-Dyson equation. The advantage of this method is that: (1) The basis for the flow kernel can be chosen arbitrarily by hand. (2) It can be applied to the general action of interest. (3) The coefficients in the kernel are determined by lattice estimates of the observables, which does not require analytic calculations beforehand. We perform the HMC with the effective action obtained by the Schwinger-Dyson method, and show that we can have better control of the effective action than the known $t$-expansion construction. However, the algorithmic overhead is still large and overwhelming the gain though faster decorrelation is observed for long-range observables in some cases. This contribution reports the preliminary results of this attempt., Comment: 19 pages, 12 figures, talk presented at the 39th International Symposium on Lattice Field Theory (Lattice 2022)

Lattice gauge theory continues to be a powerful theoretical and computational approach to simulating strongly interacting quantum field theories, whose applications permeate almost all disciplines of modern-day research in High-Energy Physics. Whether it is to enable precision quark- and lepton-flavor physics, to uncover signals of new physics in nucleons and nuclei, to elucidate hadron structure and spectrum, to serve as a numerical laboratory to reach beyond the Standard Model, or to invent and improve state-of-the-art computational paradigms, the lattice-gauge-theory program is in a prime position to impact the course of developments and enhance discovery potential of a vibrant experimental program in High-Energy Physics over the coming decade. This projection is based on abundant successful results that have emerged using lattice gauge theory over the years: on continued improvement in theoretical frameworks and algorithmic suits; on the forthcoming transition into the exascale era of high-performance computing; and on a skillful, dedicated, and organized community of lattice gauge theorists in the U.S. and worldwide. The prospects of this effort in pushing the frontiers of research in High-Energy Physics have recently been studied within the U.S. decadal Particle Physics Planning Exercise (Snowmass 2021), and the conclusions are summarized in this Topical Report., Comment: 57 pages, 1 figure. Submitted to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021). Topical Group Report for TF05 - Lattice Gauge Theory

Contribution from the USQCD Collaboration to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021)., Comment: 27 pp. main text, 4 pp. appendices, 29 pp. references, 1 p. index

Lattice quantum chromodynamics has proven to be an indispensable method to determine nonperturbative strong contributions to weak decay processes. In this white paper for the Snowmass community planning process we highlight achievements and future avenues of research for lattice calculations of weak $b$ and $c$ quark decays, and point out how these calculations will help to address the anomalies currently in the spotlight of the particle physics community. With future increases in computational resources and algorithmic improvements, percent level (and below) lattice determinations will play a central role in constraining the standard model or identifying new physics., Comment: contribution to Snowmass 2021; 19 pages; v2 corrected typo and added references

The decays and mixing of $K$ mesons are remarkably sensitive to the weak interactions of quarks and leptons at high energies. They provide important tests of the standard model at both first and second order in the Fermi constant $G_F$ and offer a window into possible new phenomena at energies as high as 1,000 TeV. These possibilities become even more compelling as the growing capabilities of lattice QCD make high-precision standard model predictions possible. Here we discuss and attempt to forecast some of these capabilities., Comment: submitted to the Rare Processes and Precision, Theory and Computational Frontiers for the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021)

We discuss a substantial update to the Grid software library for Lattice QCD, enabling it to port to multiple GPU architectures while retaining CPU vectorisation and SIMD execution within OpenMP threads. The GPU environments supported include vendor specific Nvidia CUDA and AMD HIP environments and a (mostly) standards based SYCL implementation. This is performed by an internal abstraction interface giving single source cross-platform performance portability across all number of planned Exascale architectures, and all those planned by the US Department of Energy., Comment: 10 pages, 6 figures

We introduce a new non-perturbative method to tune the parameters of the Columbia formulation of an anisotropic, clover-improved relativistic heavy-quark (RHQ) action. By making use of suitable observables which can be computed at a sequence of heavy-quark mass values, employing an $O(a)$-improved discretized action with domain-wall chiral fermion, and safely interpolated between the accessible heavy-quark mass region and the static point predicted by heavy-quark effective theory, we are able to precisely determine the unknown coefficients of the RHQ action. In this proof-of-principle study we benefit from the RBC/UKQCD Iwasaki gauge configurations with $2+1$ flavors of dynamical quarks, at three values of the lattice spacing varying from $0.11$ to $0.062$ fm. Preliminary results and applications to bottom spectroscopy are also presented., Comment: 11 pages, 5 figures, 1 table. Poster presented at The 38th International Symposium on Lattice Field Theory, LATTICE2021 26-30 July, 2021, Zoom/Gather@Massachusetts Institute of Technology

Two important sources of systematic errors in lattice QCD calculations of radiative leptonic decays are unwanted exponentials in the sum over intermediate states and unwanted excited states created by the meson interpolating field. Performing the calculation using a 3d sequential propagator allows for better control over the systematic uncertainties from intermediate states, while using a 4d sequential propagator allows for better control over the systematic uncertainties from excited states. We calculate form factors using both methods and compare how reliably each controls these systematic errors. We also employ a hybrid approach involving global fits to data from both methods., Comment: 10 pages, 4 figures, talk given at the 38th International Symposium on Lattice Field Theory (LATTICE2021), 26th-30th July 2021, Zoom/Gather@Massachusetts Institute of Technology

We present a lattice QCD calculation of the $\Delta I=1/2$, $K\to\pi\pi$ decay amplitude $A_0$ and $\varepsilon'$, the measure of direct CP-violation in $K\to\pi\pi$ decay, improving our 2015 calculation of these quantities. Both calculations were performed with physical kinematics on a $32^3\times 64$ lattice with an inverse lattice spacing of $a^{-1}=1.3784(68)$ GeV. However, the current calculation includes nearly four times the statistics and numerous technical improvements allowing us to more reliably isolate the $\pi\pi$ ground-state and more accurately relate the lattice operators to those defined in the Standard Model. We find ${\rm Re}(A_0)=2.99(0.32)(0.59)\times 10^{-7}$ GeV and ${\rm Im}(A_0)=-6.98(0.62)(1.44)\times 10^{-11}$ GeV, where the errors are statistical and systematic, respectively. The former agrees well with the experimental result ${\rm Re}(A_0)=3.3201(18)\times 10^{-7}$ GeV. These results for $A_0$ can be combined with our earlier lattice calculation of $A_2$ to obtain ${\rm Re}(\varepsilon'/\varepsilon)=21.7(2.6)(6.2)(5.0) \times 10^{-4}$, where the third error represents omitted isospin breaking effects, and Re$(A_0)$/Re$(A_2) = 19.9(2.3)(4.4)$. The first agrees well with the experimental result of ${\rm Re}(\varepsilon'/\varepsilon)=16.6(2.3)\times 10^{-4}$. A comparison of the second with the observed ratio Re$(A_0)/$Re$(A_2) = 22.45(6)$, demonstrates the Standard Model origin of this "$\Delta I = 1/2$ rule" enhancement., Comment: Updated to published version. 95 pages, 12 figures

There are emerging tensions for theory results of the hadronic vacuum polarization contribution to the muon anomalous magnetic moment both within recent lattice QCD calculations and between some lattice QCD calculations and R-ratio results. In this paper we work towards scrutinizing critical aspects of these calculations. We focus in particular on a precise calculation of Euclidean position-space windows defined by RBC/UKQCD that are ideal quantities for cross-checks within the lattice community and with R-ratio results. We perform a lattice QCD calculation using physical up, down, strange, and charm sea quark gauge ensembles generated in the staggered formalism by the MILC collaboration. We study the continuum limit using inverse lattice spacings from $a^{-1}\approx 1.6$ GeV to $3.5$ GeV, identical to recent studies by FNAL/HPQCD/MILC and Aubin et al. and similar to the recent study of BMW. Our calculation exhibits a tension for the particularly interesting window result of $a_\mu^{\rm ud, conn.,isospin, W}$ from $0.4$ fm to $1.0$ fm with previous results obtained with a different discretization of the vector current on the same gauge configurations. Our results may indicate a difficulty related to estimating uncertainties of the continuum extrapolation that deserves further attention. In this work we also provide results for $a_\mu^{\rm ud,conn.,isospin}$, $a_\mu^{\rm s,conn.,isospin}$, $a_\mu^{\rm SIB,conn.}$ for the total contribution and a large set of windows. For the total contribution, we find $a_\mu^{\rm HVP~LO}=714(27)(13) 10^{-10}$, $a_\mu^{\rm ud,conn.,isospin}=657(26)(12) 10^{-10}$, $a_\mu^{\rm s,conn.,isospin}=52.83(22)(65) 10^{-10}$, and $a_\mu^{\rm SIB,conn.}=9.0(0.8)(1.2) 10^{-10}$, where the first uncertainty is statistical and the second systematic. We also comment on finite-volume corrections for the strong-isospin-breaking corrections., Comment: 27 pages, 15 figures

We report the first result for the hadronic light-by-light scattering contribution to the muon anomalous magnetic moment with all errors systematically controlled. Several ensembles using 2+1 flavors of physical mass M\"obius domain-wall fermions, generated by the RBC/UKQCD collaborations, are employed to take the continuum and infinite volume limits of finite volume lattice QED+QCD. We find $a_\mu^{\rm HLbL} = 7.87(3.06)_\text{stat}(1.77)_\text{sys}\times 10^{-10}$. Our value is consistent with previous model results and leaves little room for this notoriously difficult hadronic contribution to explain the difference between the Standard Model and the BNL experiment.

The Hadronic Vacuum Polarization (HVP) is a dominant contribution to the theoretical uncertainty of the muon anomalous magnetic moment. The uncertainty in a lattice QCD calculation of the connected light-quark contribution to the HVP is dominated by the long-distance region of the vector correlation function. Explicit studies of the exclusive channels of the HVP diagram make it possible to reconstruct the long-distance behavior of the correlation function. This removes most of the statistical uncertainty of the correlation function. In these proceedings, preliminary results of an exclusive study of the isospin symmetric connected-only vector-vector correlation function using a hybrid of distillation and A2A techniques are presented. The computation is performed on 2+1 flavor M\"obius Domain Wall Fermion ensembles with physical pion mass. Reconstruction of the long-distance correlation function will enable lattice-only calculations of the HVP to achieve precision similar to estimates of the HVP from the R-ratio method on the timescale of the new experimental measurements of the muon anomalous magnetic moment., Comment: 7 pages, 4 figures; Presented at the 37th International Symposium on Lattice Field Theory, 16-22 June 2019, Wuhan, China

We report preliminary results for the hadronic light-by-light scattering contribution to the muon anomalous magnetic moment. Several ensembles using 2+1 flavors of M\"obius domain-wall fermions, generated by the RBC/UKQCD collaborations, are employed to take the continuum and infinite volume limits of finite volume lattice QED+QCD. We find $a_\mu^{\rm HLbL} = (7.41\pm6.33)\times 10^{-10}$, Comment: 8 pages, 19 figures, to appear in the proceedings of Rencontres de Moriond (EW)

Adding a hard photon to the final state of a leptonic pseudoscalar-meson decay lifts the helicity suppression and can provide sensitivity to a larger set of operators in the weak effective Hamiltonian. Furthermore, radiative leptonic $B$ decays at high photon energy are well suited to constrain the first inverse moment of the $B$-meson light-cone distribution amplitude, an important parameter in the theory of nonleptonic $B$ decays. We demonstrate that the calculation of radiative leptonic decays is possible using Euclidean lattice QCD, and present preliminary numerical results for $D_s^+ \to \ell^+ \nu\gamma$ and $K^- \to \ell^-\bar{\nu}\gamma$., Comment: 7 pages, 5 figures, Proceedings of the 37th Annual International Symposium on Lattice Field Theory (Lattice 2019), June 16-22, Wuhan, China

This document is one of a series of whitepapers from the USQCD collaboration. Here, we discuss opportunities for lattice QCD in quark and lepton flavor physics. New data generated at Belle II, LHCb, BES III, NA62, KOTO, and Fermilab E989, combined with precise calculations of the relevant hadronic physics, may reveal what lies beyond the Standard Model. We outline a path toward improvements of the precision of existing lattice-QCD calculations and discuss groundbreaking new methods that allow lattice QCD to access new observables., Comment: USQCD whitepaper

A determination of the hadronic vacuum polarization contribution to the anomalous magnetic moment of the muon from lattice QCD aiming at a precision of $1\%$ requires to include isospin breaking corrections in the computation. We present a lattice calculation of the QED and strong isospin breaking corrections to the hadronic vacuum polarization with Domain Wall fermions. The results are obtained using quark masses which are tuned such that pion and kaon masses agree with their physical values including isospin breaking corrections., Comment: 7 pages, 4 figures, presented at the 36th Annual International Symposium on Lattice Field Theory - LATTICE2018 - 22-28 July, 2018, Michigan State University, East Lansing, Michigan, USA

Hadronic spectral functions of $\tau$ decays have been used in the past to provide an alternative determination of the LO Hadronic Vacuum Polarization relevant for the (g-2) of the muon. Following recent developments and results in Lattice QCD+QED calculations, we explore the possibility of studying the isospin breaking corrections of $\tau$ spectral functions for this prediction. We present preliminary results at physical pion mass based on Domain Wall Fermion ensembles generated by the RBC/UKQCD collaboration., Comment: 7 pages, 6 figures. Talk presented at The 36th Annual International Symposium on Lattice Field Theory

We study the finite-volume correction on the hadronic vacuum polarization contribution to the muon g-2 ($a_\mu^{\rm hvp}$) in lattice QCD at (near) 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 stout-smeared clover-Wilson fermion action and Iwasaki gauge 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 state propagation 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 those correction with chiral perturbation. We find $10(26)\times10^{-10}$ difference for light quark $a_\mu^{\rm hvp}$ between $(5.4~{\rm fm})^4$ and $(8.1~{\rm fm})^4$ lattice in 146 MeV pion.

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., Comment: 10 pages, 11 figures. References added and minor modifications. Accepted for publication in Physical Review Letters

We propose a method to compute the Wilson coefficients of the weak effective Hamiltonian to all orders in the strong coupling constant using Lattice QCD simulations. We perform our calculations adopting an unphysically light weak boson mass of around $2~\mathrm{GeV}$. We demonstrate that systematic errors for the Wilson coefficients $C_1$ and $C_2$, related to the current-current four-quark operators, can be controlled and present a path towards precise determinations in subsequent works., Comment: 16 pages, 9 figures

We discuss practical methods to ensure near wirespeed performance from clusters with either one or two Intel(R) Omni-Path host fabric interfaces (HFI) per node, and Intel(R) Xeon Phi(TM) 72xx (Knight's Landing) processors, and using the Linux operating system. The study evaluates the performance improvements achievable and the required programming approaches in two distinct example problems: firstly in Cartesian communicator halo exchange problems, appropriate for structured grid PDE solvers that arise in quantum chromodynamics simulations of particle physics, and secondly in gradient reduction appropriate to synchronous stochastic gradient descent for machine learning. As an example, we accelerate a published Baidu Research reduction code and obtain a factor of ten speedup over the original code using the techniques discussed in this paper. This displays how a factor of ten speedup in strongly scaled distributed machine learning could be achieved when synchronous stochastic gradient descent is massively parallelised with a fixed mini-batch size. We find a significant improvement in performance robustness when memory is obtained using carefully allocated 2MB "huge" virtual memory pages, implying that either non-standard allocation routines should be used for communication buffers. These can be accessed via a LD\_PRELOAD override in the manner suggested by libhugetlbfs. We make use of a the Intel(R) MPI 2019 library "Technology Preview" and underlying software to enable thread concurrency throughout the communication software stake via multiple PSM2 endpoints per process and use of multiple independent MPI communicators. When using a single MPI process per node, we find that this greatly accelerates delivered bandwidth in many core Intel(R) Xeon Phi processors., Comment: 17 pages, 5 figures

We present results for the QED and strong isospin breaking corrections to the hadronic vacuum polarization using $N_f=2+1$ Domain Wall fermions. QED is included in an electro-quenched setup using two different methods, a stochastic and a perturbative approach. Results and statistical errors from both methods are directly compared with each other., Comment: 8 pages, 6 figures, presented at the 35th International Symposium on Lattice Field Theory (Lattice 2017), Granada, Spain, June 18-24, 2017

In this talk I present the current status of a precise first-principles calculation of the quark connected, quark disconnected, and leading QED and strong isospin-breaking contributions to the leading-order hadronic vacuum polarization by the RBC and UKQCD collaborations. The lattice data is also combined with experimental $e^+ e^-$ scattering data, consistency between the two datasets is checked, and a combined result with smaller error than the lattice data and $e^+ e^-$ scattering data individually is presented., Comment: 15 pages, 17 figures; Talk given at the 35th International Symposium on Lattice Field Theory, 18-24 June 2017, Granada, Spain

We present a Lanczos algorithm utilizing multiple grids that reduces the memory requirements both on disk and in working memory by one order of magnitude for RBC/UKQCD's 48I and 64I ensembles at the physical pion mass. The precision of the resulting eigenvectors is on par with exact deflation., Comment: 6 pages, 7 figures; Talk given at the 35th International Symposium on Lattice Field Theory, 18-24 June 2017, Granada, Spain

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., Comment: 12 pages

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., Comment: 13 pages, 2 tables, 6 figures, presented at the 34th annual International Symposium on Lattice Field Theory (Lattice2016), 24-30 July 2016, University of Southampton, UK

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., Comment: 34th annual International Symposium on Lattice Field Theory, 24-30 July 2016, University of Southampton, UK. 7 pages

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., Comment: 6 pages

We present preliminary results for the strange leading-order hadronic contribution to the anomalous magnetic moment of the muon using RBC/UKQCD physical point domain wall fermions ensembles. We discuss various analysis strategies in order to constrain the systematic uncertainty in the final result.

We report on progress to renormalize non-pertubatively the static heavy quark theory on the lattice. In particular, we present first results for position-space renormalization scheme for heavy-light bilinears. We test our approach on RBC's 16^3 x 32 lattice ensemble with m_pi = 420 MeV, Iwasaki gauge action and domain wall light fermions., Comment: 7 pages, 6 figures; contribution at the 33rd International Symposium on Lattice Field Theory, 14 -18 July 2015, Kobe, Japan

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., Comment: The 33rd International Symposium on Lattice Field Theory, 14 -18 July 2015, Kobe International Conference Center, Kobe, Japan

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., Comment: 7 pages, 12 figures. 2 tables, Proceedings of the 33rd International Symposium on Lattice Field Theory, July 14-18, 2015, Kobe, Japan

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., Comment: 52 pages

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., Comment: Twelfth Conference on the Intersections of Particle and Nuclear Physics. Vail Colorado at the Vail Marriott from May 19-24, 2015

It has been clear for past many years that in low energy observables new physics can only appear as a perturbation. Therefore precise theoretical predictions and precise experimental measurements have become mandatory. Here we draw attention to the significant advances that have been made on the lattice in recent years in $K\to \pi \pi$, $\Delta M_K$, the long-distance part of $\varepsilon$ and rare K-decays. Thus, in conjunction with experiments, the construction of a unitarity triangle purely from Kaon physics should soon become feasible. Along with the B-unitarity triangle, this should allow for more stringent tests of the Standard Model and tighter constraints on new physics., Comment: 9 pages, 6 figures, Version to appear on PLB

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., Comment: Presented at the 32nd International Symposium on Lattice Field Theory (Lattice 2014), 23-28 June 2014, New York, NY, USA

Measurements of the $\Lambda_b \to p \ell^- \bar{\nu}_\ell$ and $\Lambda_b \to \Lambda_c \ell^- \bar{\nu}_\ell$ decay rates can be used to determine the magnitudes of the CKM matrix elements $V_{ub}$ and $V_{cb}$, provided that the relevant hadronic form factors are known. Here we present a precise calculation of these form factors using lattice QCD with 2+1 flavors of dynamical domain-wall fermions. The $b$ and $c$ quarks are implemented with relativistic heavy-quark actions, allowing us to work directly at the physical heavy-quark masses. The lattice computation is performed for six different pion masses and two different lattice spacings, using gauge-field configurations generated by the RBC and UKQCD collaborations. The $b \to u$ and $b \to c$ currents are renormalized with a mostly nonperturbative method. We extrapolate the form factor results to the physical pion mass and the continuum limit, parametrizing the $q^2$-dependence using $z$-expansions. The form factors are presented in such a way as to enable the correlated propagation of both statistical and systematic uncertainties into derived quantities such as differential decay rates and asymmetries. Using these form factors, we present predictions for the $\Lambda_b \to p \ell^- \bar{\nu}_\ell$ and $\Lambda_b \to \Lambda_c \ell^- \bar{\nu}_\ell$ differential and integrated decay rates. Combined with experimental data, our results enable determinations of $|V_{ub}|$, $|V_{cb}|$, and $|V_{ub}/V_{cb}|$ with theory uncertainties of 4.4%, 2.2%, and 4.9%, respectively., Comment: 40 pages, 15 figures, form factor parameters included as ancillary files. v3: new definitions of z-expansion parameters t_+ and other minor improvements