Search

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.

A study of two-pion scattering for the isospin channels, I=0 and I=2, using lattice QCD is presented. Möbius 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 nonzero relative momentum, as well as two pions at rest. Our results for phase shifts for both the I=0 and I=2 channels are consistent with 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→ππ decay, on physical-mass ensembles using the GEVP.

We present a lattice calculation of the $K\toππ$ matrix elements and amplitudes with both the $Δ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ππ$ 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ππ$ 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ππ$ 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.

The search for new physics requires a joint experimental and theoretical effort. Lattice QCD is already an essential tool for obtaining precise model-free theoretical predictions of the hadronic processes underlying many key experimental searches, such as those involving heavy flavor physics, the anomalous magnetic moment of the muon, nucleon-neutrino scattering, and rare, second-order electroweak processes. As experimental measurements become more precise over the next decade, lattice QCD will play an increasing role in providing the needed matching theoretical precision. Achieving the needed precision requires simulations with lattices with substantially increased resolution. As we push to finer lattice spacing we encounter an array of new challenges. They include algorithmic and software-engineering challenges, challenges in computer technology and design, and challenges in maintaining the necessary human resources. In this white paper we describe those challenges and discuss ways they are being dealt with. Overcoming them is key to supporting the community effort required to deliver the needed theoretical support for experiments in the coming decade., Contribution to Snowmass 2021. 22 pages

We report a first principle lattice calculation of the pion electric polarizability $\alpha_\pi$ at the physical pion mass. First, we derive the master formula, which relates the pion polarizabilities with the position space hadronic Compton tensor, $\langle \pi | J_\mu(x) J_\nu(0) | \pi \rangle$. The finite volume error of the master formula is exponentially suppressed by the spatial extent of the lattice. Then, the hadronic tensor is calculated using domain wall fermions (DWF) directly at physical pion mass. The gauge ensembles are generated by the RBC-UKQCD collaborations., Comment: 9 pages, 5 figures. Talk presented at the 38th International Symposium on Lattice Field Theory (Lattice 2021), 26-30 July, 2021, Zoom/Gather@Massachusetts Institute of Technology

We calculate $K\to\pi\pi$ matrix elements using periodic boundary conditions as an independent calculation from our previous study with G-parity boundary conditions. We present our preliminary results for $K\to\pi\pi$ three-point functions and matrix elements on a $24^3, a^{-1} = 1$~GeV, $2+1$-flavor M\"obius DWF ensemble at physical pion and kaon masses generated by the RBC and UKQCD collaborations and discuss the prospect for high-precision computation of $\varepsilon'$ with periodic boundary conditions., Proceedings of the 38th International Symposium on Lattice Field Theory (Lattice 2021), July 26-30, 2021, Zoom/Gather@MIT

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

We perform a digital quantum simulation of a gauge theory with a topological term in Minkowski spacetime, which is practically inaccessible by standard lattice Monte Carlo simulations. We focus on $1+1$ dimensional quantum electrodynamics with the $\theta$-term known as the Schwinger model. We construct the true vacuum state of a lattice Schwinger model using adiabatic state preparation which, in turn, allows us to compute an expectation value of the fermion mass operator with respect to the vacuum. Upon taking a continuum limit we find that our result in massless case agrees with the known exact result. In massive case, we find an agreement with mass perturbation theory in small mass regime and deviations in large mass regime. We estimate computational costs required to take a reasonable continuum limit. Our results imply that digital quantum simulation is already useful tool to explore non-perturbative aspects of gauge theories with real time and topological terms., Comment: 13 pages, 6 figures; v3: accepted version in PRD

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)

We report nucleon mass, isovector vector and axial-vector charges, and tensor and scalar couplings, calculated using two recent 2+1-flavor dynamical domain-wall fermions lattice-QCD ensembles generated jointly by the RIKEN-BNL-Columbia and UKQCD collaborations. These ensembles were generated with Iwasaki $\times$ dislocation-suppressing-determinant-ratio gauge action at inverse lattice spacing of 1.378(7) GeV and pion mass values of 249.4(3) and 172.3(3) MeV. The nucleon mass extrapolates to a value $m_N = 0.950(5)$ GeV at physical point. The isovector vector charge renormalizes to unity in the chiral limit, narrowly constraining excited-state contamination in the calculation. The ratio of the isovector axial-vector to vector charges shows a deficit of about ten percent. The tensor coupling no longer depends on mass and extrapolates to 1.04(5) in $\overline {\rm MS}$ 2-GeV renormalization at physical point, in a good agreement with the value obtained at the lightest mass in our previous calculations and other calculations that followed. The scalar charge, though noisier, does not show mass dependence and is in agreement with other calculations., Comment: Accepted for publication in Physical Review D, 11 pages, 13 figures and 7 tables. arXiv admin note: text overlap with arXiv:1710.06656

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 $��$ 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 $��$ decay determination employing flavor-breaking sum rules are also discussed., 10 pages, 11 figures. References added and minor modifications. Accepted for publication in Physical Review Letters

Hadronic spectral functions of $��$ 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 $��$ 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., 7 pages, 6 figures. Talk presented at The 36th Annual International Symposium on Lattice Field Theory

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., Comment: 33 pages, 22 figures

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

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., Comment: 33 pages, 18 figures, v3: published version

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, $��M_K$ and $��_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], Presented at the 31st International Symposium on Lattice Field Theory (Lattice 2013), 29 July - 3 August 2013, Mainz, Germany, 7 pages and 5 figures

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., Comment: Presented at the 31st International Symposium on Lattice Field Theory (Lattice 2013), 29 July - 3 August 2013, Mainz, Germany

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., Comment: Proceedings of the 30th International Symposium on Lattice Field Theory, June 24-29, 2012, Cairns, Australia

We discuss the theoretical input into the current RBC-UKQCD calculation of $f_{B_{d, s}}$ and $B_{B_{d, s}}$ using a smeared static heavy quark propagator, light domain-wall quarks and the Iwasaki gauge action. We present the complete one-loop, mean-field improved matching of heavy-light current and four-fermion lattice operators onto the static continuum theory renormalized in $\bar{\text{MS}}$(NDR). The large degeneracies present in a static calculation are addressed, and a method for extracting $f_B$ and $B_B$ using only box sources is described; implications for future calculations are discussed., Presented at the XXV International Symposium on Lattice Field Theory, Regensburg, Germany (30 July - 4 August 2007); 7 pages, 2 figures

We study the isospin breaking effect on octet baryons. Using the two-flavor dynamical domain-wall QCD configurations combined with the quenched non-compact QED configurations, the electromagnetic mass splittings between isomultiplets (p, n), (Sigma^+, Sigma^0, Sigma^-), (Xi^0, Xi^-) are investigated. We evaluate the main source of statistical fluctuations in the two-point correlation function, and find that the elimination of O(e) fluctuation (e: the QED charge) is essential to extract the signal. Preliminary results for $m_p - m_n$ as well as other mass splittings are presented. Possible origin of systematic uncertainty is also discussed., Comment: Talk given at 24th International Symposium on Lattice Field Theory (LATTICE 2006), Tucson, Arizona, USA, 23-28 Jul. 2006 (Hadron Spectroscopy), 7pages, 4figures

We calculate the static quark and anti-quark potential both in quenched and two-flavor dynamical quark lattice QCD using DBW2 gauge and domain-wall quark actions. Lattice spacings from Sommer scale are determined. We find (i) mixing of excited states is different in between quenched and dynamical, (ii) lattice spacing $a_{r_0}\sim a_{m_{\rho}}$ in dynamical and (iii) coefficient of Coulomb term being $\alpha_{N_f=0}, Comment: Talk presented at Lattice2004(spectrum), Fermilab, June 21-26, 2004

We study lattice QCD with two flavors of dynamical domain wall quarks. With renormalization group motivated actions, we find chiral symmetry can be preserved to a high degree at lattice cut off of a^{-1} \sim 2 GeV even for fifth dimension size as small as L_s = 12. In addition two new steps are introduced to improve the performance of the hybrid Monte Carlo simulation., Lattice2002(chiral), 2 figures

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)

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

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