Kaon mixing beyond the standard model with physical masses

INSPIRE: 2773765
arXiv: 2404.02297
DOI: 10.1103/PhysRevD.110.034501

Authors: Boyle, Peter A., Erben, Felix, Flynn, Jonathan M., Garron, Nicolas, Kettle, Julia, Mukherjee, Rajnandini, Tsang, J. Tobias

Submitted: 2 April 2024

Subjects:

Journal reference: Phys.Rev.D 110 (3) 034501 (2024)

Citations: 20

Abstract

We present nonperturbative results for beyond the standard model kaon mixing matrix elements in the isospin symmetric limit (<math display="inline"><msub><mi>m</mi><mi>u</mi></msub><mo>=</mo><msub><mi>m</mi><mi>d</mi></msub></math>) of QCD, including a complete estimate of all dominant sources of systematic error. Our results are obtained from numerical simulations of lattice QCD with <math display="inline"><msub><mi>N</mi><mi>f</mi></msub><mo>=</mo><mn>2</mn><mo>+</mo><mn>1</mn></math> flavors of dynamical domain wall fermions. For the first time, these quantities are simulated directly at the physical pion mass <math display="inline"><msub><mi>m</mi><mi>π</mi></msub><mo>∼</mo><mn>139</mn><mtext> </mtext><mtext> </mtext><mi>MeV</mi></math> for two different lattice spacings. We include data at three lattice spacings in the range <math display="inline"><mrow><mi>a</mi><mo>=</mo><mn>0.11</mn><mi>–</mi><mn>0.07</mn><mtext> </mtext><mtext> </mtext><mi>fm</mi></mrow></math> and with pion masses ranging from the physical value up to 450 MeV. Compared to our earlier work, we have added both direct calculations at physical quark masses and a third lattice spacing making the removal of discretization effects significantly more precise and eliminating the need for any significant mass extrapolation beyond the range of simulated data. We renormalize the lattice operators nonperturbatively using RI-SMOM off-shell schemes. These schemes eliminate the need to model and subtract nonperturbative pion poles that arises in the RI-MOM scheme and, since the calculations are performed with domain wall fermions, the unphysical mixing between chirality sectors is suppressed. Our results for the bag parameters in the <math display="inline"><mover accent="true"><mi>MS</mi><mo stretchy="true">¯</mo></mover></math> scheme at 3 GeV are <math display="inline"><msub><mi>B</mi><mi>K</mi></msub><mo>≡</mo><msub><mi mathvariant="script">B</mi><mn>1</mn></msub><mo>=</mo><mn>0.5240</mn><mo stretchy="false">(</mo><mn>17</mn><mo stretchy="false">)</mo><mo stretchy="false">(</mo><mn>54</mn><mo stretchy="false">)</mo></math>, <math display="inline"><msub><mi mathvariant="script">B</mi><mn>2</mn></msub><mo>=</mo><mn>0.4794</mn><mo stretchy="false">(</mo><mn>25</mn><mo stretchy="false">)</mo><mo stretchy="false">(</mo><mn>35</mn><mo stretchy="false">)</mo></math>, <math display="inline"><msub><mi mathvariant="script">B</mi><mn>3</mn></msub><mo>=</mo><mn>0.746</mn><mo stretchy="false">(</mo><mn>13</mn><mo stretchy="false">)</mo><mo stretchy="false">(</mo><mn>17</mn><mo stretchy="false">)</mo></math>, <math display="inline"><msub><mi mathvariant="script">B</mi><mn>4</mn></msub><mo>=</mo><mn>0.897</mn><mo stretchy="false">(</mo><mn>02</mn><mo stretchy="false">)</mo><mo stretchy="false">(</mo><mn>10</mn><mo stretchy="false">)</mo></math> and <math display="inline"><msub><mi mathvariant="script">B</mi><mn>5</mn></msub><mo>=</mo><mn>0.6882</mn><mo stretchy="false">(</mo><mn>78</mn><mo stretchy="false">)</mo><mo stretchy="false">(</mo><mn>94</mn><mo stretchy="false">)</mo></math>, where the first error is from lattice uncertainties and the second is the uncertainty due to the perturbative matching to <math display="inline"><mover accent="true"><mi>MS</mi><mo stretchy="true">¯</mo></mover></math>.