Neutral Kaon Mixing Beyond the Standard Model with $n_f=2+1$ Chiral Fermions Part 1: Bare Matrix Elements and Physical Results

We compute the hadronic matrix elements of the four-quark operators relevant for $K^0-{\bar K^0}$ mixing beyond the Standard Model. Our results are from lattice QCD simulations with $n_f=2+1$ flavours of domain-wall fermion, which exhibit continuum-like chiral-flavour symmetry. The simulations are performed at two different values of the lattice spacing ($a\sim0.08$ and $a\sim 0.11 \, \fm $) and with lightest unitary pion mass $\sim 300\, \MeV$. For the first time, the full set of relevant four-quark operators is renormalised non-perturbatively through RI-SMOM schemes; a detailed description of the renormalisation procedure is presented in a companion paper. We argue that the intermediate renormalisation scheme is responsible for the discrepancies found by different collaborations. We also study different normalisations and determine the matrix elements of the relevant four-quark operators with a precision of $\sim 5\%$ or better.