We present a master formula describing the neutrinoless-double-beta decay (0νββ) rate induced by lepton-number-violating (LNV) operators up to dimension nine in the Standard Model Effective Field Theory. We provide an end-to-end framework connecting the possibly very high LNV scale to the nuclear scale, through a chain of effective field theories. Starting at the electroweak scale, we integrate out the heavy Standard Model degrees of freedom and we match to an SU (3) c ⊗ U (1) em effective theory. After evolving the resulting effective Lagrangian to the QCD scale, we use chiral perturbation theory to derive the lepton-number-violating chiral Lagrangian. The chiral Lagrangian is used to derive the twonucleon 0νββ transition operators to leading order in the chiral power counting. Based on renormalization arguments we show that in various cases short-range two-nucleon operators need to be enhanced to leading order. We show that all required nuclear matrix elements can be taken from existing calculations. Our final result is a master formula that describes the 0νββ rate in terms of phase-space factors, nuclear matrix elements, hadronic low-energy constants, QCD evolution factors, and high-energy LNV Wilson coefficients, including all the interference terms. Our master formula can be easily matched to any model where LNV originates at energy scales above the electroweak scale. As an explicit example, we match our formula to the minimal left-right-symmetric model in which contributions of operators of different dimension compete, and we discuss the resulting phenomenology.

Within the framework of chiral effective field theory, we discuss the leading contributions to the neutrinoless double-beta decay transition operator induced by light Majorana neutrinos. Based on renormalization arguments in both dimensional regularization with minimal subtraction and a coordinate-space cutoff scheme, we show the need to introduce a leading-order short-range operator, missing in all current calculations. We discuss strategies to determine the finite part of the short-range coupling by matching to lattice QCD or by relating it via chiral symmetry to isospin-breaking observables in the two-nucleon sector. Finally, we speculate on the impact of this new contribution on nuclear matrix elements of relevance to experiment.

We compute the one-loop matching between the Standard Model Effective Field Theory and the low-energy effective field theory below the electroweak scale, where the heavy gauge bosons, the Higgs particle, and the top quark are integrated out. The complete set of matching equations is derived including effects up to dimension six in the power counting of both theories. We present the results for general flavor structures and include both the CP -even and CP -odd sectors. The matching equations express the masses, gauge couplings, as well as the coefficients of dipole, three-gluon, and four-fermion operators in the low-energy theory in terms of the parameters of the Standard Model Effective Field Theory. Using momentum insertion, we also obtain the matching for the CP -violating theta angles. Our results provide an ingredient for a model-independent analysis of constraints on physics beyond the Standard Model. They can be used for fixed-order calculations at one-loop accuracy and represent a first step towards a systematic next-to-leading-log analysis. E SMEFT operator basis 56 F LEFT operator basis 58 G Conventions for the supplemental material 59 References 61

We analyze neutrinoless double beta decay (0νββ) within the framework of the Standard Model Effective Field Theory. Apart from the dimension-five Weinberg operator, the first contributions appear at dimension seven. We classify the operators and evolve them to the electroweak scale, where we match them to effective dimension-six, -seven, and -nine operators. In the next step, after renormalization group evolution to the QCD scale, we construct the chiral Lagrangian arising from these operators. We develop a power-counting scheme and derive the two-nucleon 0νββ currents up to leading order in the power counting for each lepton-number-violating operator. We argue that the leadingorder contribution to the decay rate depends on a relatively small number of nuclear matrix elements. We test our power counting by comparing nuclear matrix elements obtained by various methods and by different groups. We find that the power counting works well for nuclear matrix elements calculated from a specific method, while, as in the case of light Majorana neutrino exchange, the overall magnitude of the matrix elements can differ by factors of two to three between methods. We calculate the constraints that can be set on dimension-seven lepton-number-violating operators from 0νββ experiments and study the interplay between dimension-five and -seven operators, discussing how dimension-seven contributions affect the interpretation of 0νββ in terms of the effective Majorana mass m ββ .