Showing 1–50 of 56 results for author: Izubuchi, T

$ΔI = 3/2$ and $ΔI = 1/2$ channels of $K\toππ$ decay at the physical point with periodic boundary conditions

Authors: Thomas Blum, Peter A. Boyle, Daniel Hoying, Taku Izubuchi, Luchang Jin, Chulwoo Jung, Christopher Kelly, Christoph Lehner, Amarjit Soni, Masaaki Tomii

Abstract: 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 pr… ▽ More 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. △ Less

Submitted 3 June, 2024; v1 submitted 11 June, 2023; originally announced June 2023.

Journal ref: Phys,Rev,D.,108,094517 (2023)

An update of Euclidean windows of the hadronic vacuum polarization

Authors: T. Blum, P. A. Boyle, M. Bruno, D. Giusti, V. Gülpers, R. C. Hill, T. Izubuchi, Y. -C. Jang, L. Jin, C. Jung, A. Jüttner, C. Kelly, C. Lehner, N. Matsumoto, R. D. Mawhinney, A. S. Meyer, J. T. Tsang

Abstract: We compute the standard Euclidean window of the hadronic vacuum polarization using multiple independent blinded analyses. We improve the continuum and infinite-volume extrapolations of the dominant quark-connected light-quark isospin-symmetric contribution and address additional sub-leading systematic effects from sea-charm quarks and residual chiral-symmetry breaking from first principles. We fin… ▽ More We compute the standard Euclidean window of the hadronic vacuum polarization using multiple independent blinded analyses. We improve the continuum and infinite-volume extrapolations of the dominant quark-connected light-quark isospin-symmetric contribution and address additional sub-leading systematic effects from sea-charm quarks and residual chiral-symmetry breaking from first principles. We find $a_μ^{\rm W} = 235.56(65)(50) \times 10^{-10}$, which is in $3.8σ$ tension with the recently published dispersive result of Colangelo et al., $a_μ^{\rm W} = 229.4(1.4) \times 10^{-10}$, and in agreement with other recent lattice determinations. We also provide a result for the standard short-distance window. The results reported here are unchanged compared to our presentation at the Edinburgh workshop of the g-2 Theory Initiative in 2022. △ Less

Submitted 20 January, 2023; originally announced January 2023.

Comments: 24 pages, 15 figures

Transverse momentum-dependent parton distributions for longitudinally polarized nucleons from domain wall fermion calculations at the physical pion mass

Authors: M. Engelhardt, N. Hasan, T. Izubuchi, C. Kallidonis, S. Krieg, S. Meinel, J. Negele, A. Pochinsky, G. Silvi, S. Syritsyn

Abstract: Previous Lattice QCD calculations of nucleon transverse momentum-dependent parton distributions (TMDs) focused on the case of transversely polarized nucleons, and thus did not encompass two leading-twist TMDs associated with longitudinal polarization, namely, the helicity TMD and the worm-gear TMD corresponding to transversely polarized quarks in a longitudinally polarized nucleon. Based on a defi… ▽ More Previous Lattice QCD calculations of nucleon transverse momentum-dependent parton distributions (TMDs) focused on the case of transversely polarized nucleons, and thus did not encompass two leading-twist TMDs associated with longitudinal polarization, namely, the helicity TMD and the worm-gear TMD corresponding to transversely polarized quarks in a longitudinally polarized nucleon. Based on a definition of TMDs via hadronic matrix elements of quark bilocal operators containing staple-shaped gauge connections, TMD observables characterizing the aforementioned two TMDs are evaluated, utilizing an RBC/UKQCD domain wall fermion ensemble at the physical pion mass. △ Less

Submitted 15 January, 2023; originally announced January 2023.

Comments: 9 pages, 4 figures, to appear in the proceedings of the 39th International Symposium on Lattice Field Theory, LATTICE2022, PoS LATTICE2022, 103

Snowmass Theory Frontier Report

Authors: N. Craig, C. Csáki, A. X. El-Khadra, Z. Bern, R. Boughezal, S. Catterall, Z. Davoudi, A. de Gouvêa, P. Draper, P. J. Fox, D. Green, D. Harlow, R. Harnik, V. Hubeny, T. Izubuchi, S. Kachru, G. Kribs, H. Murayama, Z. Ligeti, J. Maldacena, F. Maltoni, I. Mocioiu, E. T. Neil, S. Pastore, D. Poland , et al. (16 additional authors not shown)

Abstract: This report summarizes the recent progress and promising future directions in theoretical high-energy physics (HEP) identified within the Theory Frontier of the 2021 Snowmass Process. This report summarizes the recent progress and promising future directions in theoretical high-energy physics (HEP) identified within the Theory Frontier of the 2021 Snowmass Process. △ Less

Submitted 12 December, 2022; v1 submitted 10 November, 2022; originally announced November 2022.

Comments: Contribution to the US Community Study on the Future of Particle Physics (Snowmass 2021), v2: fixed typos

Report of the Snowmass 2021 Topical Group on Lattice Gauge Theory

Authors: Zohreh Davoudi, Ethan T. Neil, Christian W. Bauer, Tanmoy Bhattacharya, Thomas Blum, Peter Boyle, Richard C. Brower, Simon Catterall, Norman H. Christ, Vincenzo Cirigliano, Gilberto Colangelo, Carleton DeTar, William Detmold, Robert G. Edwards, Aida X. El-Khadra, Steven Gottlieb, Rajan Gupta, Daniel C. Hackett, Anna Hasenfratz, Taku Izubuchi, William I. Jay, Luchang Jin, Christopher Kelly, Andreas S. Kronfeld, Christoph Lehner , et al. (13 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 21 September, 2022; originally announced September 2022.

Comments: 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

Report number: UMD-PP-022-08, LA-UR-22-29361, FERMILAB-CONF-22-703-T

Lattice QCD and Particle Physics

Authors: Andreas S. Kronfeld, Tanmoy Bhattacharya, Thomas Blum, Norman H. Christ, Carleton DeTar, William Detmold, Robert Edwards, Anna Hasenfratz, Huey-Wen Lin, Swagato Mukherjee, Konstantinos Orginos, Richard Brower, Vincenzo Cirigliano, Zohreh Davoudi, Bálint Jóo, Chulwoo Jung, Christoph Lehner, Stefan Meinel, Ethan T. Neil, Peter Petreczky, David G. Richards, Alexei Bazavov, Simon Catterall, Jozef J. Dudek, Aida X. El-Khadra , et al. (57 additional authors not shown)

Abstract: Contribution from the USQCD Collaboration to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021). Contribution from the USQCD Collaboration to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021). △ Less

Submitted 2 October, 2022; v1 submitted 15 July, 2022; originally announced July 2022.

Comments: 27 pp. main text, 4 pp. appendices, 29 pp. references, 1 p. index

Report number: FERMILAB-CONF-22-531-T

Prospects for precise predictions of $a_μ$ in the Standard Model

Authors: G. Colangelo, M. Davier, A. X. El-Khadra, M. Hoferichter, C. Lehner, L. Lellouch, T. Mibe, B. L. Roberts, T. Teubner, H. Wittig, B. Ananthanarayan, A. Bashir, J. Bijnens, T. Blum, P. Boyle, N. Bray-Ali, I. Caprini, C. M. Carloni Calame, O. Catà, M. Cè, J. Charles, N. H. Christ, F. Curciarello, I. Danilkin, D. Das , et al. (57 additional authors not shown)

Abstract: We discuss the prospects for improving the precision on the hadronic corrections to the anomalous magnetic moment of the muon, and the plans of the Muon $g-2$ Theory Initiative to update the Standard Model prediction. We discuss the prospects for improving the precision on the hadronic corrections to the anomalous magnetic moment of the muon, and the plans of the Muon $g-2$ Theory Initiative to update the Standard Model prediction. △ Less

Submitted 29 March, 2022; originally announced March 2022.

Comments: Contribution to the US Community Study on the Future of Particle Physics (Snowmass 2021)

Report number: FERMILAB-CONF-22-236-T, LTH 1303, MITP-22-030

Discovering new physics in rare kaon decays

Authors: Thomas Blum, Peter Boyle, Mattia Bruno, Norman Christ, Felix Erben, Xu Feng, Vera Guelpers, Ryan Hill, Raoul Hodgson, Danel Hoying, Taku Izubuchi, Yong-Chull Jang, Luchang Jin, Chulwoo Jung, Joe Karpie, Christopher Kelly, Christoph Lehner, Antonin Portelli, Christopher Sachrajda, Amarjit Soni, Masaaki Tomii, Bigeng Wang, Tianle Wang

Abstract: 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 latt… ▽ More 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. △ Less

Submitted 21 March, 2022; originally announced March 2022.

Comments: 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)

Electric dipole moments and the search for new physics

Authors: Ricardo Alarcon, Jim Alexander, Vassilis Anastassopoulos, Takatoshi Aoki, Rick Baartman, Stefan Baeßler, Larry Bartoszek, Douglas H. Beck, Franco Bedeschi, Robert Berger, Martin Berz, Hendrick L. Bethlem, Tanmoy Bhattacharya, Michael Blaskiewicz, Thomas Blum, Themis Bowcock, Anastasia Borschevsky, Kevin Brown, Dmitry Budker, Sergey Burdin, Brendan C. Casey, Gianluigi Casse, Giovanni Cantatore, Lan Cheng, Timothy Chupp , et al. (118 additional authors not shown)

Abstract: Static electric dipole moments of nondegenerate systems probe mass scales for physics beyond the Standard Model well beyond those reached directly at high energy colliders. Discrimination between different physics models, however, requires complementary searches in atomic-molecular-and-optical, nuclear and particle physics. In this report, we discuss the current status and prospects in the near fu… ▽ More Static electric dipole moments of nondegenerate systems probe mass scales for physics beyond the Standard Model well beyond those reached directly at high energy colliders. Discrimination between different physics models, however, requires complementary searches in atomic-molecular-and-optical, nuclear and particle physics. In this report, we discuss the current status and prospects in the near future for a compelling suite of such experiments, along with developments needed in the encompassing theoretical framework. △ Less

Submitted 4 April, 2022; v1 submitted 15 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021; updated with community edits and endorsements

Quark spin-orbit correlations in the proton

Authors: M. Engelhardt, J. R. Green, N. Hasan, T. Izubuchi, C. Kallidonis, S. Krieg, S. Liuti, S. Meinel, J. Negele, A. Pochinsky, A. Rajan, G. Silvi, S. Syritsyn

Abstract: Generalized transverse momentum-dependent parton distributions (GTMDs) provide a comprehensive framework for imaging the internal structure of the proton. In particular, by encoding the simultaneous distribution of quark transverse positions and momenta, they allow one to directly access longitudinal quark orbital angular momentum, and, moreover, to correlate it with the quark helicity. The releva… ▽ More Generalized transverse momentum-dependent parton distributions (GTMDs) provide a comprehensive framework for imaging the internal structure of the proton. In particular, by encoding the simultaneous distribution of quark transverse positions and momenta, they allow one to directly access longitudinal quark orbital angular momentum, and, moreover, to correlate it with the quark helicity. The relevant GTMD is evaluated through a lattice calculation of a proton matrix element of a quark bilocal operator (the separation in which is Fourier conjugate to the quark momentum) featuring a momentum transfer (which is Fourier conjugate to the quark position), as well as the Dirac structure appropriate for capturing the quark helicity. The weighting by quark transverse position requires a derivative with respect to momentum transfer, which is obtained in unbiased fashion using a direct derivative method. The lattice calculation is performed directly at the physical pion mass, using domain wall fermions to mitigate operator mixing effects. Both the Jaffe-Manohar as well as the Ji quark spin-orbit correlations are extracted, yielding evidence for a strong quark spin-orbit coupling in the proton. △ Less

Submitted 26 December, 2021; originally announced December 2021.

Comments: 8 pages, 2 figures, to appear in the proceedings of the 38th International Symposium on Lattice Field Theory, LATTICE2021, PoS LATTICE2021, 413

The anomalous magnetic moment of the muon in the Standard Model

Authors: T. Aoyama, N. Asmussen, M. Benayoun, J. Bijnens, T. Blum, M. Bruno, I. Caprini, C. M. Carloni Calame, M. Cè, G. Colangelo, F. Curciarello, H. Czyż, I. Danilkin, M. Davier, C. T. H. Davies, M. Della Morte, S. I. Eidelman, A. X. El-Khadra, A. Gérardin, D. Giusti, M. Golterman, Steven Gottlieb, V. Gülpers, F. Hagelstein, M. Hayakawa , et al. (107 additional authors not shown)

Abstract: We review the present status of the Standard Model calculation of the anomalous magnetic moment of the muon. This is performed in a perturbative expansion in the fine-structure constant $α$ and is broken down into pure QED, electroweak, and hadronic contributions. The pure QED contribution is by far the largest and has been evaluated up to and including $\mathcal{O}(α^5)$ with negligible numerical… ▽ More We review the present status of the Standard Model calculation of the anomalous magnetic moment of the muon. This is performed in a perturbative expansion in the fine-structure constant $α$ and is broken down into pure QED, electroweak, and hadronic contributions. The pure QED contribution is by far the largest and has been evaluated up to and including $\mathcal{O}(α^5)$ with negligible numerical uncertainty. The electroweak contribution is suppressed by $(m_μ/M_W)^2$ and only shows up at the level of the seventh significant digit. It has been evaluated up to two loops and is known to better than one percent. Hadronic contributions are the most difficult to calculate and are responsible for almost all of the theoretical uncertainty. The leading hadronic contribution appears at $\mathcal{O}(α^2)$ and is due to hadronic vacuum polarization, whereas at $\mathcal{O}(α^3)$ the hadronic light-by-light scattering contribution appears. Given the low characteristic scale of this observable, these contributions have to be calculated with nonperturbative methods, in particular, dispersion relations and the lattice approach to QCD. The largest part of this review is dedicated to a detailed account of recent efforts to improve the calculation of these two contributions with either a data-driven, dispersive approach, or a first-principle, lattice-QCD approach. The final result reads $a_μ^\text{SM}=116\,591\,810(43)\times 10^{-11}$ and is smaller than the Brookhaven measurement by 3.7$σ$. The experimental uncertainty will soon be reduced by up to a factor four by the new experiment currently running at Fermilab, and also by the future J-PARC experiment. This and the prospects to further reduce the theoretical uncertainty in the near future-which are also discussed here-make this quantity one of the most promising places to look for evidence of new physics. △ Less

Submitted 13 November, 2020; v1 submitted 8 June, 2020; originally announced June 2020.

Comments: 196 pages, 103 figures, version published in Phys. Rept., bib files for the citation references are available from: https://muon-gm2-theory.illinois.edu

Report number: FERMILAB-PUB-20-207-T, INT-PUB-20-021, KEK Preprint 2020-5, MITP/20-028, CERN-TH-2020-075, IFT-UAM/CSIC-20-74, LMU-ASC 18/20, LTH 1234, LU TP 20-20, MAN/HEP/2020/003, PSI-PR-20-06, UWThPh 2020-14, ZU-TH 18/20

Journal ref: Phys. Rept. 887 (2020) 1-166

Computing Nucleon Electric Dipole Moment from lattice QCD

Authors: Taku Izubuchi, Hiroshi Ohki, Sergey Syritsyn

Abstract: Electric dipole moments (EDMs) of nucleons and nuclei are actively considered as direct evidence of the CP violation. Calculations of nucleon EDMs on lattice are required to connect the quark- and hadron- level effective CP violating interactions within QCD or other CP violating sources in new physics beyond the standard model. Among them, the theta-induced nucleon EDM, that is the only such renor… ▽ More Electric dipole moments (EDMs) of nucleons and nuclei are actively considered as direct evidence of the CP violation. Calculations of nucleon EDMs on lattice are required to connect the quark- and hadron- level effective CP violating interactions within QCD or other CP violating sources in new physics beyond the standard model. Among them, the theta-induced nucleon EDM, that is the only such renormalizable interaction, has widely been investigated on a lattice. In the report, we review recent developments of the lattice calculations of nucleon EDM induced QCD theta term. △ Less

Submitted 22 April, 2020; originally announced April 2020.

Comments: Proceedings of the 37th International Symposium on Lattice Field Theory -- Lattice2019, 16-22 June 2019, Wuhan, China; 15 pages, 11 figures

Classically Emulated Digital Quantum Simulation of the Schwinger Model with Topological Term via Adiabatic State Preparation

Authors: Bipasha Chakraborty, Masazumi Honda, Taku Izubuchi, Yuta Kikuchi, Akio Tomiya

Abstract: 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 $θ$-term known as the Schwinger model. We construct the true vacuum state of a lattice Schwinger model using adiabatic state preparation which, in turn,… ▽ More 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 $θ$-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. △ Less

Submitted 22 April, 2022; v1 submitted 2 January, 2020; originally announced January 2020.

Comments: 13 pages, 6 figures; v3: accepted version in PRD

The hadronic light-by-light scattering contribution to the muon anomalous magnetic moment from lattice QCD

Authors: Thomas Blum, Norman Christ, Masashi Hayakawa, Taku Izubuchi, Luchang Jin, Chulwoo Jung, Christoph Lehner

Abstract: 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öbius 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… ▽ More 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öbius 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_μ^{\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. △ Less

Submitted 23 March, 2020; v1 submitted 19 November, 2019; originally announced November 2019.

Journal ref: Phys. Rev. Lett. 124, 132002 (2020)

Hadronic light-by-light contribution to the muon anomalous magnetic moment from lattice QCD

Authors: Thomas Blum, Norman Christ, Masashi Hayakawa, Taku Izubuchi, Luchang Jin, Chulwoo Jung, Christoph Lehner

Abstract: 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öbius 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_μ^{\rm HLbL} = (7.41\pm6.33)\times 10^{-10}$ 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öbius 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_μ^{\rm HLbL} = (7.41\pm6.33)\times 10^{-10}$ △ Less

Submitted 1 July, 2019; originally announced July 2019.

Comments: 8 pages, 19 figures, to appear in the proceedings of Rencontres de Moriond (EW)

Valence parton distribution function of pion from fine lattice

Authors: Taku Izubuchi, Luchang Jin, Christos Kallidonis, Nikhil Karthik, Swagato Mukherjee, Peter Petreczky, Charles Shugert, Sergey Syritsyn

Abstract: We present a lattice QCD study of valence parton distribution inside the pion within the framework of Large Momentum Effective Theory. We use a mixed action approach with 1-HYP smeared valence Wilson clover quarks on 2+1 flavor HISQ sea with the valence quark mass tuned to 300 MeV pion mass. We use $48^3 \times 64$ lattice at a fine lattice spacing $a=0.06$ fm for this computation. We renormalize… ▽ More We present a lattice QCD study of valence parton distribution inside the pion within the framework of Large Momentum Effective Theory. We use a mixed action approach with 1-HYP smeared valence Wilson clover quarks on 2+1 flavor HISQ sea with the valence quark mass tuned to 300 MeV pion mass. We use $48^3 \times 64$ lattice at a fine lattice spacing $a=0.06$ fm for this computation. We renormalize the quasi-PDF matrix element in the non-perturbative RI-MOM scheme. As a byproduct, we test the validity of 1-loop matching procedure by comparing the RI-MOM renormalized quasi-PDF matrix element with off-shell quark external states as computed in the continuum 1-loop perturbation theory with the lattice results at $a=0.04$ and 0.06 fm. By applying the RI-MOM to ${\bar{\rm MS}}$ one-loop matching, implemented through a fit to phenomenologically motivated PDFs, we obtain the valence PDF of pion. △ Less

Submitted 5 September, 2019; v1 submitted 15 May, 2019; originally announced May 2019.

Comments: 25 pages, 22 figures; version published in PRD

Journal ref: Phys. Rev. D 100, 034516 (2019)

The Role of Lattice QCD in Searches for Violations of Fundamental Symmetries and Signals for New Physics

Authors: Vincenzo Cirigliano, Zohreh Davoudi, Tanmoy Bhattacharya, Taku Izubuchi, Phiala E. Shanahan, Sergey Syritsyn, Michael L. Wagman

Abstract: This document is one of a series of whitepapers from the USQCD collaboration. Here, we discuss opportunities for Lattice Quantum Chromodynamics (LQCD) in the research frontier in fundamental symmetries and signals for new physics. LQCD, in synergy with effective field theories and nuclear many-body studies, provides theoretical support to ongoing and planned experimental programs in searches for e… ▽ More This document is one of a series of whitepapers from the USQCD collaboration. Here, we discuss opportunities for Lattice Quantum Chromodynamics (LQCD) in the research frontier in fundamental symmetries and signals for new physics. LQCD, in synergy with effective field theories and nuclear many-body studies, provides theoretical support to ongoing and planned experimental programs in searches for electric dipole moments of the nucleon, nuclei and atoms, decay of the proton, $n$-$\overline{n}$ oscillations, neutrinoless double-$β$ decay of a nucleus, conversion of muon to electron, precision measurements of weak decays of the nucleon and of nuclei, precision isotope-shift spectroscopy, as well as direct dark matter detection experiments using nuclear targets. This whitepaper details the objectives of the LQCD program in the area of Fundamental Symmetries within the USQCD collaboration, identifies priorities that can be addressed within the next five years, and elaborates on the areas that will likely demand a high degree of innovation in both numerical and analytical frontiers of the LQCD research. △ Less

Submitted 21 April, 2019; originally announced April 2019.

Comments: A whitepaper by the USQCD Collaboration, 30 pages, 9 figures

Journal ref: Eur. Phys. J. A (2019) 55: 197

Calculation of Nucleon Electric Dipole Moments Induced by Quark Chromo-Electric Dipole Moments and the QCD $θ$-term

Authors: S. Syritsyn, T. Izubuchi, H. Ohki

Abstract: Electric dipole moments (EDMs) of nucleons and nuclei, which are sought as evidence of CP violation, require lattice calculations to connect constraints from experiments to limits on the strong CP violation within QCD or CP violation introduced by new physics from beyond the standard model. Nucleon EDM calculations on a lattice are notoriously hard due to large statistical noise, chiral symmetry v… ▽ More Electric dipole moments (EDMs) of nucleons and nuclei, which are sought as evidence of CP violation, require lattice calculations to connect constraints from experiments to limits on the strong CP violation within QCD or CP violation introduced by new physics from beyond the standard model. Nucleon EDM calculations on a lattice are notoriously hard due to large statistical noise, chiral symmetry violating effects, and potential mixing of the EDM and the anomalous magnetic moment of the nucleon. In this report, details of ongoing lattice calculations of proton and neutron EDMs induced by the QCD $θ$-term and the quark chromo-EDM, the lowest-dimension effective CP-violating quark-gluon interaction are presented. Our calculation employs chiral-symmetric fermion discretization. An assessment of feasibility of nucleon EDM calculations at the physical point is discussed. △ Less

Submitted 16 January, 2019; originally announced January 2019.

Comments: proceedings XIIIth Quark Confinement and Hadron Spectrum conference (Maynooth University, Ireland, July 31-August 6, 2018); 11 pages, 13 figures. arXiv admin note: substantial text overlap with arXiv:1810.03721

Report number: RBRC-1304

Proton decay matrix element on the lattice with physical pion mass

Authors: Jun-Sik Yoo, Yasumichi Aoki, Taku Izubuchi, Sergey Syritsyn

Abstract: Proton decay is one of possible signatures of baryon number violation, which has to exist to explain the baryon asymmetry and the existence of nuclear matter. Proton decays must be mediated through effective low-energy baryon number violating operators made of three quarks and a lepton. We calculate matrix elements of these operators between the proton and various meson final states using the dire… ▽ More Proton decay is one of possible signatures of baryon number violation, which has to exist to explain the baryon asymmetry and the existence of nuclear matter. Proton decays must be mediated through effective low-energy baryon number violating operators made of three quarks and a lepton. We calculate matrix elements of these operators between the proton and various meson final states using the direct method. We report on preliminary results of matrix element calculation done with the 2+1 dynamical flavor domain wall fermions at the physical point for the first time. △ Less

Submitted 21 December, 2018; originally announced December 2018.

Comments: 7 pages, 5 figures, 1 table

Report number: RBRC-130, BNL-209806-2018-COPA

Progress in the Nucleon Electric Dipole Moment Calculations in Lattice QCD

Authors: S. Syritsyn, T. Izubuchi, H. Ohki

Abstract: Electric dipole moments (EDMs), which are sought as evidence of CP violation, require lattice calculations to connect constraints from experiments to limits on the strong CP violation within QCD or CP violation in new physics beyond the standard model. Nucleon EDM calculations on a lattice are notoriously hard due to large statistical noise, chiral symmetry violating effects, and potential mixing… ▽ More Electric dipole moments (EDMs), which are sought as evidence of CP violation, require lattice calculations to connect constraints from experiments to limits on the strong CP violation within QCD or CP violation in new physics beyond the standard model. Nucleon EDM calculations on a lattice are notoriously hard due to large statistical noise, chiral symmetry violating effects, and potential mixing of the EDM and the anomalous magnetic moment of the nucleon. In this report, details of ongoing lattice calculations of proton and neutron EDMs induced by the QCD $θ$-term and the quark chromo-EDM, the lowest-dimension effective CP-violating quark-gluon interaction are presented. Our calculation employs chiral-symmetric fermion discretization. An assessment of feasibility of nucleon EDM calculations at the physical point is discussed. △ Less

Submitted 8 October, 2018; originally announced October 2018.

Comments: proceedings for conference CIPANP 2018 (Palm Springs); 13 pages, 8 figures

Report number: RBRC-1292

Finite-volume correction on the hadronic vacuum polarization contribution to muon g-2 in lattice QCD

Authors: Taku Izubuchi, Yoshinobu Kuramashi, Christoph Lehner, Eigo Shintani

Abstract: We study the finite-volume correction on the hadronic vacuum polarization contribution to the muon g-2 ($a_μ^{\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 acti… ▽ More We study the finite-volume correction on the hadronic vacuum polarization contribution to the muon g-2 ($a_μ^{\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_μ^{\rm hvp}$ between $(5.4~{\rm fm})^4$ and $(8.1~{\rm fm})^4$ lattice in 146 MeV pion. △ Less

Submitted 2 October, 2018; v1 submitted 11 May, 2018; originally announced May 2018.

Journal ref: Phys. Rev. D 98, 054505 (2018)

Novel $|V_{us}|$ Determination Using Inclusive Strange $τ$ Decay and Lattice HVPs

Authors: Peter Boyle, Renwick James Hudspith, Taku Izubuchi, Andreas Jüttner, Christoph Lehner, Randy Lewis, Kim Maltman, Hiroshi Ohki, Antonin Portelli, Matthew Spraggs

Abstract: 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 error… ▽ More 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. △ Less

Submitted 19 November, 2018; v1 submitted 19 March, 2018; originally announced March 2018.

Comments: 10 pages, 11 figures. References added and minor modifications. Accepted for publication in Physical Review Letters

Journal ref: Phys. Rev. Lett. 121, 202003 (2018)

Calculation of the hadronic vacuum polarization contribution to the muon anomalous magnetic moment

Authors: T. Blum, P. A. Boyle, V. Gülpers, T. Izubuchi, L. Jin, C. Jung, A. Jüttner, C. Lehner, A. Portelli, J. T. Tsang

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

Submitted 22 January, 2018; originally announced January 2018.

Comments: 12 pages, 11 figures

Journal ref: Phys. Rev. Lett. 121, 022003 (2018)

Factorization Theorem Relating Euclidean and Light-Cone Parton Distributions

Authors: Taku Izubuchi, Xiangdong Ji, Luchang Jin, Iain W. Stewart, Yong Zhao

Abstract: 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 p… ▽ More 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. △ Less

Submitted 15 September, 2018; v1 submitted 11 January, 2018; originally announced January 2018.

Comments: 19 pages, 4 figures

Report number: MIT-CTP 4960

Journal ref: Phys. Rev. D 98, 056004 (2018)

Improved lattice computation of proton decay matrix elements

Authors: Yasumichi Aoki, Taku Izubuchi, Eigo Shintani, Amarjit Soni

Abstract: 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_π=0.34$--0.69 GeV on 2.7~fm… ▽ More 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_π=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πe^+$ or from 20--40\% for $p\rightarrow K \barν$. The accuracy of the low energy constants $α$ and $β$ in the leading-order baryon chiral perturbation theory (BChPT) of proton decay are also improved. The relevant form factors of $p\rightarrow π$ 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 $α$ and $β$. 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. △ Less

Submitted 3 May, 2017; originally announced May 2017.

Comments: 28 pages, 14 figures, 6 tables

Journal ref: Phys. Rev. D 96, 014506 (2017)

Connected and leading disconnected hadronic light-by-light contribution to the muon anomalous magnetic moment with physical pion mass

Authors: Thomas Blum, Norman Christ, Masashi Hayakawa, Taku Izubuchi, Luchang Jin, Chulwoo Jung, Christoph Lehner

Abstract: 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… ▽ More 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_μ^{\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. △ Less

Submitted 14 October, 2016; originally announced October 2016.

Comments: 6 pages

Journal ref: Phys. Rev. Lett. 118, 022005 (2017)

Erratum: Standard-model prediction for direct CP violation in $K\toππ$ decay

Authors: Z. Bai, T. Blum, P. A. Boyle, N. H. Christ, J. Frison, N. Garron, T. Izubuchi, C. Jung, C. Kelly, C. Lehner, R. D. Mawhinney, C. T. Sachrajda, A. Soni, D. Zhang

Abstract: In this document we address an error discovered in the ensemble generation for our calculation of the $I=0$ $K\toππ$ amplitude (Phys. Rev. Lett. 115, 212001 (2015), arXiv:1505.07863) whereby the same random numbers were used for the two independent quark flavors, resulting in small but measurable correlations between gauge observables separated by 12 units in the y-direction. We conclude that the… ▽ More In this document we address an error discovered in the ensemble generation for our calculation of the $I=0$ $K\toππ$ amplitude (Phys. Rev. Lett. 115, 212001 (2015), arXiv:1505.07863) whereby the same random numbers were used for the two independent quark flavors, resulting in small but measurable correlations between gauge observables separated by 12 units in the y-direction. We conclude that the effects of this error are negligible compared to the overall errors on our calculation. △ Less

Submitted 8 March, 2016; originally announced March 2016.

Comments: 2 pages, 1 figure

Lattice calculation of the leading strange quark-connected contribution to the muon $g-2$

Authors: T. Blum, P. A. Boyle, L. Del Debbio, R. J. Hudspith, T. Izubuchi, A. Jüttner, C. Lehner, R. Lewis, K. Maltman, M. Krstić Marinković, A. Portelli, M. Spraggs

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

Submitted 3 May, 2016; v1 submitted 4 February, 2016; originally announced February 2016.

Comments: 21 pages, 5 figures

Journal ref: JHEP 1604 (2016) 063

Calculation of the hadronic vacuum polarization disconnected contribution to the muon anomalous magnetic moment

Authors: T. Blum, P. A. Boyle, T. Izubuchi, L. Jin, A. Jüttner, C. Lehner, K. Maltman, M. Marinkovic, A. Portelli, M. Spraggs

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

Submitted 30 December, 2015; originally announced December 2015.

Comments: 8 pages, 13 figures

Journal ref: Phys. Rev. Lett. 116, 232002 (2016)

Neutron and proton electric dipole moments from $N_f=2+1$ domain-wall fermion lattice QCD

Authors: Eigo Shintani, Thomas Blum, Taku Izubuchi, Amarjit Soni

Abstract: 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 $θ$ 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 pi… ▽ More 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 $θ$ 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. △ Less

Submitted 1 December, 2015; originally announced December 2015.

Comments: 33 pages, 22 figures

Journal ref: Phys. Rev. D 93, 094503 (2016)

On calculating disconnected-type hadronic light-by-light scattering diagrams from lattice QCD

Authors: Masashi Hayakawa, Thomas Blum, Norman Christ, Taku Izubuchi, Luchang Jin, Christoph Lehner

Abstract: 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… ▽ More 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. △ Less

Submitted 18 November, 2015; v1 submitted 4 November, 2015; originally announced November 2015.

Comments: 7 pages, 12 figures. 2 tables, Proceedings of the 33rd International Symposium on Lattice Field Theory, July 14-18, 2015, Kobe, Japan

Lattice Calculation of Hadronic Light-by-Light Contribution to the Muon Anomalous Magnetic Moment

Authors: Thomas Blum, Norman Christ, Masashi Hayakawa, Taku Izubuchi, Luchang Jin, Christoph Lehner

Abstract: 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 thes… ▽ More 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. △ Less

Submitted 23 October, 2015; originally announced October 2015.

Comments: 52 pages

Journal ref: Phys. Rev. D 93, 014503 (2016)

Standard-model prediction for direct CP violation in $K\toππ$ decay

Authors: Z. Bai, T. Blum, P. A. Boyle, N. H. Christ, J. Frison, N. Garron, T. Izubuchi, C. Jung, C. Kelly, C. Lehner, R. D. Mawhinney, C. T. Sachrajda, A. Soni, D. Zhang

Abstract: We report the first lattice QCD calculation of the complex kaon decay amplitude $A_0$ with physical kinematics, using a $32^3\times 64$ lattice volume and a single lattice spacing $a$, with $1/a= 1.3784(68)$ GeV. We find Re$(A_0) = 4.66(1.00)(1.26) \times 10^{-7}$ GeV and Im$(A_0) = -1.90(1.23)(1.08) \times 10^{-11}$ GeV, where the first error is statistical and the second systematic. The first va… ▽ More We report the first lattice QCD calculation of the complex kaon decay amplitude $A_0$ with physical kinematics, using a $32^3\times 64$ lattice volume and a single lattice spacing $a$, with $1/a= 1.3784(68)$ GeV. We find Re$(A_0) = 4.66(1.00)(1.26) \times 10^{-7}$ GeV and Im$(A_0) = -1.90(1.23)(1.08) \times 10^{-11}$ GeV, where the first error is statistical and the second systematic. The first value is in approximate agreement with the experimental result: Re$(A_0) = 3.3201(18) \times 10^{-7}$ GeV while the second can be used to compute the direct CP violating ratio Re$(\varepsilon'/\varepsilon)=1.38(5.15)(4.59)\times 10^{-4}$, which is $2.1σ$ below the experimental value $16.6(2.3)\times 10^{-4}$. The real part of $A_0$ is CP conserving and serves as a test of our method while the result for Re$(\varepsilon'/\varepsilon)$ provides a new test of the standard-model theory of CP violation, one which can be made more accurate with increasing computer capability. △ Less

Submitted 20 January, 2016; v1 submitted 28 May, 2015; originally announced May 2015.

Comments: 9 pages, 3 figures. Updated to match published version

Report number: RBRC 1141

Journal ref: Phys. Rev. Lett. 115, 212001 (2015)

$B \to π\ell ν$ and $B_s \to K \ell ν$ form factors and $|V_{ub}|$ from 2+1-flavor lattice QCD with domain-wall light quarks and relativistic heavy quarks

Authors: J. M. Flynn, T. Izubuchi, T. Kawanai, C. Lehner, A. Soni, R. S. Van de Water, O. Witzel

Abstract: We calculate the $B \toπ\ellν$ and $B_s \to K \ellν$ form factors in dynamical lattice QCD. We use the (2+1)-flavor RBC-UKQCD gauge-field ensembles generated with the domain-wall fermion and Iwasaki gauge actions. For the $b$ quarks we use the anisotropic clover action with a relativistic heavy-quark interpretation. We analyze two lattice spacings $a \approx 0.11, 0.086$ fm and unitary pion masses… ▽ More We calculate the $B \toπ\ellν$ and $B_s \to K \ellν$ form factors in dynamical lattice QCD. We use the (2+1)-flavor RBC-UKQCD gauge-field ensembles generated with the domain-wall fermion and Iwasaki gauge actions. For the $b$ quarks we use the anisotropic clover action with a relativistic heavy-quark interpretation. We analyze two lattice spacings $a \approx 0.11, 0.086$ fm and unitary pion masses as light as $M_π\approx 290$ MeV. We simultaneously extrapolate our numerical results to the physical light-quark masses and to the continuum and interpolate in the pion/kaon energy using SU(2) "hard-pion" chiral perturbation theory. We provide complete error budgets for the form factors $f_+(q^2)$ and $f_0(q^2)$ at three momenta that span the $q^2$ range accessible in our numerical simulations. We extrapolate these results to $q^2 = 0$ using a model-independent $z$-parametrization and present our final form factors as the $z$-coefficients and the matrix of correlations between them. Our results agree with other lattice determinations using staggered light quarks and provide important independent cross-checks. Both $B \toπ\ellν$ and $B_s \to K \ellν$ decays enable a determination of the CKM matrix element $|V_{ub}|$. To illustrate this, we perform a combined $z$-fit of our numerical $B\toπ\ellν$ form-factor data with the experimental branching-fraction measurements leaving the relative normalization as a free parameter; we obtain $|V_{ub}| = 3.61(32) \times 10^{-3}$, where the error includes statistical and systematic uncertainties. This approach can be applied to $B_s\to K \ellν$ decay to determine $|V_{ub}|$ once the process has been measured experimentally. Finally, in anticipation of future measurements, we make predictions for $B \to π\ellν$ and $B_s\to K \ellν$ Standard-Model differential branching fractions and forward-backward asymmetries. △ Less

Submitted 26 April, 2015; v1 submitted 21 January, 2015; originally announced January 2015.

Comments: V3: References and discussions added. Large tables move to appendix. 39 pages, 22 tables, 23 figures. Version published in Physical Review D

Journal ref: Phys. Rev. D 91, 074510 (2015)

Hadronic light-by-light scattering contribution to the muon anomalous magnetic moment from lattice QCD

Authors: Thomas Blum, Saumitra Chowdhury, Masashi Hayakawa, Taku Izubuchi

Abstract: 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… ▽ More 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. △ Less

Submitted 10 July, 2014; originally announced July 2014.

Comments: 5 pages, 5 figures

Journal ref: Phys. Rev. Lett. 114, 012001 (2015)

Neutral $B$ meson mixings and $B$ meson decay constants with static heavy and domain-wall light quarks

Authors: Yasumichi Aoki, Tomomi Ishikawa, Taku Izubuchi, Christoph Lehner, Amarjit Soni

Abstract: 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… ▽ More 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 $ξ=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. △ Less

Submitted 29 June, 2015; v1 submitted 24 June, 2014; originally announced June 2014.

Comments: 33 pages, 18 figures, v3: published version

Report number: RBRC-1080

Journal ref: Phys. Rev. D 91, 114505 (2015)

B-meson decay constants from 2+1-flavor lattice QCD with domain-wall light quarks and relativistic heavy quarks

Authors: Norman H. Christ, Jonathan M. Flynn, Taku Izubuchi, Taichi Kawanai, Christoph Lehner, Amarjit Soni, Ruth S. Van de Water, Oliver Witzel

Abstract: 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 wit… ▽ More 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. △ Less

Submitted 2 February, 2015; v1 submitted 17 April, 2014; originally announced April 2014.

Comments: 25 pages, 13 tables, 12 figures. Version 2 revised for PRD. Sign in Eq. (27) corrected. Chiral-continuum extrapolations of fB and fBs/fB changed. Results for fB and fBs/fB, plots, and text updated accordingly. Updated references

Report number: FERMILAB-PUB-14-100-T

Journal ref: Phys. Rev. D 91, 054502 (2015)

Neutral B meson mixing with static heavy and domain-wall light quarks

Authors: Tomomi Ishikawa, Yasumichi Aoki, Taku Izubuchi, Christoph Lehner, Amarjit Soni

Abstract: 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 ar… ▽ More 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 ξ=1.222(60) in the static limit of b quark. △ Less

Submitted 3 December, 2013; originally announced December 2013.

Comments: 7 pages, 3 figures, talk presented at 31st International Symposium on Lattice Field Theory (LATTICE 2013), July 29 - August 3, 2013, Mainz, Germany

Report number: RBRC-1052

Hadronic corrections to the muon anomalous magnetic moment from lattice QCD

Authors: T. Blum, M. Hayakawa, T. Izubuchi

Abstract: After a brief self-contained introduction to the muon anomalous magnetic moment, (g-2), we review the status of lattice calculations of the hadronic vacuum polarization contribution and present first results from lattice QCD for the hadronic light-by-light scattering contribution. The signal for the latter is consistent with model calculations. While encouraging, the statistical error is large and… ▽ More After a brief self-contained introduction to the muon anomalous magnetic moment, (g-2), we review the status of lattice calculations of the hadronic vacuum polarization contribution and present first results from lattice QCD for the hadronic light-by-light scattering contribution. The signal for the latter is consistent with model calculations. While encouraging, the statistical error is large and systematic errors are mostly uncontrolled. The method is applied first to pure QED as a check. △ Less

Submitted 11 January, 2013; originally announced January 2013.

Comments: 15 pages, 11 figures. Proceedings of talk presented at Lattice 2012, Cairns, Australia, June 29, 2012

Lattice determination of the $K \to (ππ)_{I=2}$ Decay Amplitude $A_2$

Authors: T. Blum, P. A. Boyle, N. H. Christ, N. Garron, E. Goode, T. Izubuchi, C. Jung, C. Kelly, C. Lehner, M. Lightman, Q. Liu, A. T. Lytle, R. D. Mawhinney, C. T. Sachrajda, A. Soni, C. Sturm

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

Submitted 22 June, 2012; originally announced June 2012.

Comments: 59 pages, 11 figures

Report number: CU-TP-1202, Edinburgh 2012/10, MPP-2012-101, SHEP-1217

Full QED+QCD Low-Energy Constants through Reweighting

Authors: Tomomi Ishikawa, Thomas Blum, Masashi Hayakawa, Taku Izubuchi, Chulwoo Jung, Ran Zhou

Abstract: The effect of sea quark electromagnetic charge on meson masses is investigated, and first results for full QED+QCD low-energy constants are presented. The electromagnetic charge for sea quarks is incorporated in quenched QED+full QCD lattice simulations by a reweighting method. The reweighting factor, which connects quenched and unquenched QED, is estimated using a stochastic method on 2+1 flavor… ▽ More The effect of sea quark electromagnetic charge on meson masses is investigated, and first results for full QED+QCD low-energy constants are presented. The electromagnetic charge for sea quarks is incorporated in quenched QED+full QCD lattice simulations by a reweighting method. The reweighting factor, which connects quenched and unquenched QED, is estimated using a stochastic method on 2+1 flavor dynamical domain-wall quark ensembles. △ Less

Submitted 17 August, 2012; v1 submitted 27 February, 2012; originally announced February 2012.

Comments: 5 pages, 9 figures, REVTeX 4.1, v2: published version

Report number: RBRC-945

Journal ref: Phys. Rev. Lett. 109, 072002 (2012)

The $K\to(ππ)_{I=2}$ Decay Amplitude from Lattice QCD

Authors: T. Blum, P. A. Boyle, N. H. Christ, N. Garron, E. Goode, T. Izubuchi, C. Jung, C. Kelly, C. Lehner, M. Lightman, Q. Liu, A. T. Lytle, R. D. Mawhinney, C. T. Sachrajda, A. Soni, C. Sturm

Abstract: We report on the first realistic \emph{ab initio} calculation of a hadronic weak decay, that of the amplitude $A_2$ for a kaon to decay into two π-mesons with isospin 2. We find Re$A_2=(1.436\pm 0.063_{\textrm{stat}}\pm 0.258_{\textrm{syst}})\,10^{-8}\,\textrm{GeV}$ in good agreement with the experimental result and for the hitherto unknown imaginary part we find {Im}… ▽ More We report on the first realistic \emph{ab initio} calculation of a hadronic weak decay, that of the amplitude $A_2$ for a kaon to decay into two π-mesons with isospin 2. We find Re$A_2=(1.436\pm 0.063_{\textrm{stat}}\pm 0.258_{\textrm{syst}})\,10^{-8}\,\textrm{GeV}$ in good agreement with the experimental result and for the hitherto unknown imaginary part we find {Im}$\,A_2=-(6.83 \pm 0.51_{\textrm{stat}} \pm 1.30_{\textrm{syst}})\,10^{-13}\,{\rm GeV}$. Moreover combining our result for Im\,$A_2$ with experimental values of Re\,$A_2$, Re\,$A_0$ and $ε^\prime/ε$, we obtain the following value for the unknown ratio Im\,$A_0$/Re\,$A_0$ within the Standard Model: $\mathrm{Im}\,A_0/\mathrm{Re}\,A_0=-1.63(19)_{\mathrm{stat}}(20)_{\mathrm{syst}}\times10^{-4}$. One consequence of these results is that the contribution from Im\,$A_2$ to the direct CP violation parameter $ε^{\prime}$ (the so-called Electroweak Penguin, EWP, contribution) is Re$(ε^\prime/ε)_{\mathrm{EWP}} = -(6.52 \pm 0.49_{\textrm{stat}} \pm 1.24_{\textrm{syst}}) \times 10^{-4}$. We explain why this calculation of $A_2$ represents a major milestone for lattice QCD and discuss the exciting prospects for a full quantitative understanding of CP-violation in kaon decays. △ Less

Submitted 7 November, 2011; originally announced November 2011.

Comments: 5 pages, 1 figure

$K$ to $ππ$ Decay amplitudes from Lattice QCD

Authors: T. Blum, P. A. Boyle, N. H. Christ, N. Garron, E. Goode, T. Izubuchi, C. Lehner, Q. Liu, R. D. Mawhinney, C. T. Sachrajda, A. Soni, C. Sturm, H. Yin, R. Zhou

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

Submitted 14 June, 2011; originally announced June 2011.

Comments: 40 pages, 12 figures

Continuum Limit of $B_K$ from 2+1 Flavor Domain Wall QCD

Authors: Y. Aoki, R. Arthur, T. Blum, P. A. Boyle, D. Brömmel, N. H. Christ, C. Dawson, T. Izubuchi, C. Jung, C. Kelly, R. D. Kenway, M. Lightman, R. D. Mawhinney, Shigemi Ohta, C. T. Sachrajda, E. E. Scholz, A. Soni, C. Sturm, J. Wennekers, R. Zhou

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

Submitted 19 December, 2010; originally announced December 2010.

Comments: 65 pages

Journal ref: Phys.Rev.D84:014503,2011

Continuum Limit Physics from 2+1 Flavor Domain Wall QCD

Authors: Y. Aoki, R. Arthur, T. Blum, P. A. Boyle, D. Brommel, N. H. Christ, C. Dawson, J. M. Flynn, T. Izubuchi, X-Y. Jin, C. Jung, C. Kelly, M. Li, A. Lichtl, M. Lightman, M. F. Lin, R. D. Mawhinney, C. M. Maynard, S. Ohta, B. J. Pendleton, C. T. Sachrajda, E. E. Scholz, A. Soni, J. Wennekers, J. M. Zanotti , et al. (1 additional authors not shown)

Abstract: We present physical results obtained from simulations using 2+1 flavors of domain wall quarks and the Iwasaki gauge action at two values of the lattice spacing $a$, ($a^{-1}$=\,1.73\,(3)\,GeV and $a^{-1}$=\,2.28\,(3)\,GeV). On the coarser lattice, with $24^3\times 64\times 16$ points, the analysis of ref.[1] is extended to approximately twice the number of configurations. The ensembles on the fine… ▽ More We present physical results obtained from simulations using 2+1 flavors of domain wall quarks and the Iwasaki gauge action at two values of the lattice spacing $a$, ($a^{-1}$=\,1.73\,(3)\,GeV and $a^{-1}$=\,2.28\,(3)\,GeV). On the coarser lattice, with $24^3\times 64\times 16$ points, the analysis of ref.[1] is extended to approximately twice the number of configurations. The ensembles on the finer $32^3\times 64\times 16$ lattice are new. We explain how we use lattice data obtained at several values of the lattice spacing and for a range of quark masses in combined continuum-chiral fits in order to obtain results in the continuum limit and at physical quark masses. We implement this procedure at two lattice spacings, with unitary pion masses in the approximate range 290--420\,MeV (225--420\,MeV for partially quenched pions). We use the masses of the $π$ and $K$ mesons and the $Ω$ baryon to determine the physical quark masses and the values of the lattice spacing. While our data are consistent with the predictions of NLO SU(2) chiral perturbation theory, they are also consistent with a simple analytic ansatz leading to an inherent uncertainty in how best to perform the chiral extrapolation that we are reluctant to reduce with model-dependent assumptions about higher order corrections. Our main results include $f_π=124(2)_{\rm stat}(5)_{\rm syst}$\,MeV, $f_K/f_π=1.204(7)(25)$ where $f_K$ is the kaon decay constant, $m_s^{\bar{\textrm{MS}}}(2\,\textrm{GeV})=(96.2\pm 2.7)$\,MeV and $m_{ud}^{\bar{\textrm{MS}}}(2\,\textrm{GeV})=(3.59\pm 0.21)$\,MeV\, ($m_s/m_{ud}=26.8\pm 1.4$) where $m_s$ and $m_{ud}$ are the mass of the strange-quark and the average of the up and down quark masses respectively, $[Σ^{\msbar}(2 {\rm GeV})]^{1/3} = 256(6)\; {\rm MeV}$, where $Σ$ is the chiral condensate, the Sommer scale $r_0=0.487(9)$\,fm and $r_1=0.333(9)$\,fm. △ Less

Submitted 4 July, 2011; v1 submitted 3 November, 2010; originally announced November 2010.

Comments: 129 pages, 59 figures, Published version containing an extended discussion of reweighting and including a new appendix (Appendix C)

Journal ref: Phys.Rev.D83:074508,2011

Electromagnetic mass splittings of the low lying hadrons and quark masses from 2+1 flavor lattice QCD+QED

Authors: T. Blum, T. Doi, M. Hayakawa, T. Izubuchi, S. Uno, N. Yamada, R. Zhou

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

Submitted 2 September, 2010; v1 submitted 7 June, 2010; originally announced June 2010.

Comments: Significant additions: finite volume analysis of kaon mass (Sec. II and the Appendix) which is used to obtain new central values, extended discussion of the breaking of Dashen's theorem (VII), additional phenomenological results (VII), and a new section (VIII) on isospin breaking effects on the kaon decay constant. No conclusions or results from the first version changed significantly

Report number: KEK CP-235, RBRC-844

Journal ref: Phys.Rev.D82:094508,2010

The $η$ and $η^\prime$ mesons from Lattice QCD

Authors: N. H. Christ, C. Dawson, T. Izubuchi, C. Jung, Q. Liu, R. D. Mawhinney, C. T. Sachrajda, A. Soni, R. Zhou

Abstract: The large mass of the ninth pseudoscalar meson, the $η^\prime$, is believed to arise from the combined effects of the axial anomaly and the gauge field topology present in QCD. We report a realistic, 2+1 flavor, lattice QCD calculation of the $η$ and $η^\prime$ masses and mixing which confirms this picture. The physical eigenstates show small octet-singlet mixing with a mixing angle of… ▽ More The large mass of the ninth pseudoscalar meson, the $η^\prime$, is believed to arise from the combined effects of the axial anomaly and the gauge field topology present in QCD. We report a realistic, 2+1 flavor, lattice QCD calculation of the $η$ and $η^\prime$ masses and mixing which confirms this picture. The physical eigenstates show small octet-singlet mixing with a mixing angle of $θ= -14.1(2.8)^\circ$. Extrapolation to physical light quark mass gives, with statistical errors only, $m_η=573(6)$ MeV and $m_{η^\prime}=947(142)$ MeV, consistent with the experimental values of 548 MeV and 958 MeV. △ Less

Submitted 19 October, 2010; v1 submitted 15 February, 2010; originally announced February 2010.

Comments: 4 pages, 5 figures

Journal ref: Phys.Rev.Lett.105:241601,2010

Renormalization of quark bilinear operators in a momentum-subtraction scheme with a nonexceptional subtraction point

Authors: C. Sturm, Y. Aoki, N. H. Christ, T. Izubuchi, C. T. C. Sachrajda, A. Soni

Abstract: We extend the Rome-Southampton regularization independent momentum-subtraction renormalization scheme(RI/MOM) for bilinear operators to one with a nonexceptional, symmetric subtraction point. Two-point Green's functions with the insertion of quark bilinear operators are computed with scalar, pseudoscalar, vector, axial-vector and tensor operators at one-loop order in perturbative QCD. We call thi… ▽ More We extend the Rome-Southampton regularization independent momentum-subtraction renormalization scheme(RI/MOM) for bilinear operators to one with a nonexceptional, symmetric subtraction point. Two-point Green's functions with the insertion of quark bilinear operators are computed with scalar, pseudoscalar, vector, axial-vector and tensor operators at one-loop order in perturbative QCD. We call this new scheme RI/SMOM, where the S stands for "symmetric". Conversion factors are derived, which connect the RI/SMOM scheme and the MSbar scheme and can be used to convert results obtained in lattice calculations into the MSbar scheme. Such a symmetric subtraction point involves nonexceptional momenta implying a lattice calculation with substantially suppressed contamination from infrared effects. Further, we find that the size of the one-loop corrections for these infrared improved kinematics is substantially decreased in the case of the pseudoscalar and scalar operator, suggesting a much better behaved perturbative series. Therefore it should allow us to reduce the error in the determination of the quark mass appreciably. △ Less

Submitted 22 April, 2010; v1 submitted 19 January, 2009; originally announced January 2009.

Comments: 15 pages, 4 figures; final version accepted for publication in the journal

Report number: CU-TP-1186, KANAZAWA-09-01, RBRC-771, SHEP 09/02

Journal ref: Phys.Rev.D80:014501,2009

Determination of light quark masses from the electromagnetic splitting of pseudoscalar meson masses computed with two flavors of domain wall fermions

Authors: Thomas Blum, Takumi Doi, Masashi Hayakawa, Taku Izubuchi, Norikazu Yamada

Abstract: We determine the light quark masses from lattice QCD simulations incorporating the electromagnetic interaction of valence quarks. The meson masses are calculated on lattice QCD configurations generated by the RBC Collaboration for two flavors of dynamical domain wall fermions, which are combined with QED configurations generated via quenched non-compact lattice QED. The electromagnetic part of t… ▽ More We determine the light quark masses from lattice QCD simulations incorporating the electromagnetic interaction of valence quarks. The meson masses are calculated on lattice QCD configurations generated by the RBC Collaboration for two flavors of dynamical domain wall fermions, which are combined with QED configurations generated via quenched non-compact lattice QED. The electromagnetic part of the pion mass splitting is found to be $m_{π^+}-m_{π^0}=4.12(21)$ MeV, where only the statistical error is quoted, and similarly for the kaon, 1.443(55) MeV. Our results for the light quark masses are $m_u^{\rm\bar{MS}}$(2 GeV)=$3.02(27)(19)$ MeV, $m_d^{\rm\bar{MS}}$(2 GeV)=$5.49(20)(34)$ MeV, and $m_s^{\rm\bar{MS}}$(2 GeV)=$119.5(56)(74)$ MeV, where the first error is statistical and the second systematic. By averaging over $\pm e$ to cancel ${\cal O}(e)$ noise exactly on each combined gauge field configuration, we are able to work at physical $α=1/137$ and obtain very small statistical errors. In our calculation, several sources of systematic error remain, including finite volume, non-zero lattice spacing, chiral extrapolation, quenched QED, and quenched strange quark, which may be more significant than the errors quoted above. △ Less

Submitted 10 August, 2007; v1 submitted 3 August, 2007; originally announced August 2007.

Comments: 38 pages, 7 figures. Fixed typo in eq. 47, added footnote to Table I. Version submitted to Phys. Rev. D

Journal ref: Phys.Rev.D76:114508,2007

Neutral kaon mixing from 2+1 flavor Domain Wall QCD

Authors: D. J. Antonio, P. A. Boyle, T. Blum, N. H. Christ, S. D. Cohen, C. Dawson, T. Izubuchi, R. D. Kenway, C. Jung, S. Li, M. F. Lin, R. D. Mawhinney, J. Noaki, S. Ohta, B. J. Pendleton, E. E. Scholz, A. Soni, R. J. Tweedie, A. Yamaguchi

Abstract: We present the first results for neutral kaon mixing using 2+1 flavors of domain wall fermions. A new approach is used to extrapolate to the physical up and down quark masses from our numerical studies with pion masses in the range 240 -- 420 MeV; only $SU(2)_L \times SU(2)_R$ chiral symmetry is assumed and the kaon is not assumed to be light. Our main result is… ▽ More We present the first results for neutral kaon mixing using 2+1 flavors of domain wall fermions. A new approach is used to extrapolate to the physical up and down quark masses from our numerical studies with pion masses in the range 240 -- 420 MeV; only $SU(2)_L \times SU(2)_R$ chiral symmetry is assumed and the kaon is not assumed to be light. Our main result is $B_K^{\bar{\rm MS}}(2 \mathrm{GeV}) = 0.524(10)(28)$ where the first error is statistical and the second incorporates estimates for all systematic errors. △ Less

Submitted 22 October, 2007; v1 submitted 5 February, 2007; originally announced February 2007.

Comments: 4 pages, 4 figures. RBC and UKQCD collaborations

Report number: Edinburgh 2007/2, RBRC-638,CU-TP-1174,KEK-TH-1131

Journal ref: Phys.Rev.Lett.100:032001,2008