Determining the leading-order contact term in neutrinoless double β decay
Vincenzo Cirigliano, Wouter Dekens, Jordy de Vries, Martin Hoferichter, Emanuele Mereghetti
Abstract
A bstract We present a method to determine the leading-order (LO) contact term contributing to the nn → ppe − e − amplitude through the exchange of light Majorana neutrinos. Our approach is based on the representation of the amplitude as the momentum integral of a known kernel (proportional to the neutrino propagator) times the generalized forward Compton scattering amplitude n ( p 1 ) n ( p 2 ) W + ( k ) → $$ p\left({p}_1^{\prime}\right)p\left({p}_2^{\prime}\right){W}^{-}(k) $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mi>p</mml:mi> <mml:mfenced> <mml:msubsup> <mml:mi>p</mml:mi> <mml:mn>1</mml:mn> <mml:mo>′</mml:mo> </mml:msubsup> </mml:mfenced> <mml:mi>p</mml:mi> <mml:mfenced> <mml:msubsup> <mml:mi>p</mml:mi> <mml:mn>2</mml:mn> <mml:mo>′</mml:mo> </mml:msubsup> </mml:mfenced> <mml:msup> <mml:mi>W</mml:mi> <mml:mo>−</mml:mo> </mml:msup> <mml:mfenced> <mml:mi>k</mml:mi> </mml:mfenced> </mml:math> , in analogy to the Cottingham formula for the electromagnetic contribution to hadron masses. We construct model-independent representations of the integrand in the low- and high-momentum regions, through chiral EFT and the operator product expansion, respectively. We then construct a model for the full amplitude by interpolating between these two regions, using appropriate nucleon factors for the weak currents and information on nucleon-nucleon ( NN ) scattering in the 1 S 0 channel away from threshold. By matching the amplitude obtained in this way to the LO chiral EFT amplitude we obtain the relevant LO contact term and discuss various sources of uncertainty. We validate the approach by computing the analog I = 2 NN contact term and by reproducing, within uncertainties, the charge-independence-breaking contribution to the 1 S 0 NN scattering lengths. While our analysis is performed in the $$ \overline{\mathrm{MS}} $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mover> <mml:mi>MS</mml:mi> <mml:mo>¯</mml:mo> </mml:mover> </mml:math> scheme, we express our final result in terms of the scheme-independent renormalized amplitude $$ {\mathcal{A}}_{\nu}\left(\left|\mathbf{p}\right|,\left|\mathbf{p}^{\prime}\right|\right) $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:msub> <mml:mi>A</mml:mi> <mml:mi>ν</mml:mi> </mml:msub> <mml:mfenced> <mml:mfenced> <mml:mi>p</mml:mi> </mml:mfenced> <mml:mfenced> <mml:mrow> <mml:mi>p</mml:mi> <mml:mo>′</mml:mo> </mml:mrow> </mml:mfenced> </mml:mfenced> </mml:math> at a set of kinematic points near threshold. We illustrate for two cutoff schemes how, using our synthetic data for $$ {\mathcal{A}}_{\nu } $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:msub> <mml:mi>A</mml:mi> <mml:mi>ν</mml:mi> </mml:msub> </mml:math> , one can determine the contact-term contribution in any regularization scheme, in particular the ones employed in nuclear-structure calculations for isotopes of experimental interest.