Version 32 (modified by richardphysics, 5 weeks ago) (diff)

Added new references

SM + Heavy N at NLO in QCD

Contact Author

Richard Ruiz

  • University of Durham (IPPP)
  • richard.ruiz AT

In collaboration with Daniel Alva and Tao Han [ 1 ], and Celine Degrande, Olivier Mattelear, and Jessica Turner [ 2 ].

Model Description

This effective model extends the Standard Model (SM) field content by introducing three right-handed (RH) neutrinos, which are singlets under the SM gauge symmetry (no color, weak isospin, or weak hypercharge charges). Each RH neutrino possesses one RH Majorana mass. After electroweak symmetry breaking, the Lagrangian with three heavy Majorana neutrinos Ni (for i=1,2,3) is given by [ 3 ]

\mathcal{L} = \mathcal{L}_{\rm SM} + \mathcal{L}_{N} + \mathcal{L}_{N~\text{Int.}}

The first term is the Standard Model Lagrangian. In the mass basis, i.e., after mixing with active neutrinos, the heavy Majorana neutrinos' kinetic and mass terms are

\mathcal{L}_{N} = \frac{1}{2}\overline{N_k} i\!\not\!\partial N_k - \frac{1}{2}m_{N_k} \overline{N_k}N_k, \quad k=1,\dots,3,

and its interactions with the Weak gauge and Higgs bosons are given by

\mathcal{L}_{N~\text{Int}} = 
&-&\frac{g}{\sqrt{2}} W_{\mu}^{+}\sum_{k=1}^{3}\sum_{\ell=e}^{\tau} \overline{N_k}V_{\ell k}^{*}\gamma^{\mu}P_{L}\ell^{-}
+{\rm H.c.}
&-&\frac{g}{2\cos\theta_W}Z_{\mu}\sum_{k=1}^{3}\sum_{\ell=e}^{\tau} \overline{N_k}V_{\ell k}^{*}\gamma^{\mu}P_{L}\nu_\ell
+{\rm H.c.}
&-&\frac{g m_N}{2 M_W}         h \sum_{k=1}^{3}\sum_{\ell=e}^{\tau} \overline{N_k}V_{\ell k}^{*}P_{L}\nu_\ell
+{\rm H.c.}

Neutrino masses (mNk) and mixing parameters (Vlk) between heavy mass eigenstate and (active) flavor eigenstates are taken to be independent, phenomenological parameters. This allows for maximum flexibility and model independence when calculating rates. Therefore, some care is required by the user. The lepton number- and flavor-violating interactions of the Lagrangian allow for modeling of the Type I, Inverse, and Linear seesaw mechanisms at both lepton, hadron, and lepto-hadron colliders.

QCD Corrections

The above Lagrangian with Goldstone boson couplings and in the Feynman Gauge was implemented into FeynRules 2.3.10. QCD renormalization and R2 rational counter terms were determined using NLOCT 1.02 and FeynArts? 3.8. Feynman rules were collected into a single UFO, available below. In the UFO file, five massless quarks are assumed as are zero off-diagonal CKM matrix entries. For additional details, see [ 2 ] and references therein. These additions permit tree-level calculations at LO and NLO in QCD and loop-induced calculations at LO in QCD using MadGraph_aMC@NLO.

Model Files

  • Main model file. Relies on (default FR model file) being declared elsewhere.
  • heavyN_NLO.nb: Mathematica notebook file that generates UFO file from FeynRules model files. Allows user to also run quick sanity checks (optional) on model.
  • SM_HeavyN_Files.tgz: Standalone package containing, heavyN_NLO.nb, massless.rst (default FR file), diagonalCKM.rst (default FR file), and (default FR file).


  • To download any of the packages and unpack via the terminal, use the commands:

~/Path $ wget

~/Path $ tar -zxvf SM_HeavyN_NLO_UFO.tgz

~/Path $ wget

~/Path $ tar -zxvf SM_HeavyN_FilesWithUFO.tgz

  • This model contains 15 free parameters:
    • Three masses: mN1, mN2, mN3. Defaults are 300 GeV, 500 GeV, and 1 TeV, respectively.
    • Three widths: WN1, WN2, WN3. Defaults are 0.303 GeV, 1.50 GeV, and 12.3 GeV, respectively.
    • Nine real (no CP violation) mixing parameters: Vlk for l = e, mu, tau and k = 1,2,3. Default values are Vlk = Identity(3x3), i.e., Ve1 = Vmu2 = Vta3 = 1 and all others zero.
    • Note: VlN are restricted to be real in the model file.
    • Note: Default parameters are set so "out-of-the-box" checks can be made with [ 1 ] and [ 2 ].


  • The model file was validated at LO in [ 1 ]; see Tables 2 and 5 for further validation checks.
  • The model file was validated at NLO in [ 2 ]; see Table 1 for further validation checks.
  • For very large masses, this model has been constructed to satisfy the Goldstone Equivalence Theorem: \begin{eqnarray}
\text{BR}(N_i \rightarrow W^+ \ell^-) &=& \text{BR}(N_i \rightarrow W^- \ell^+) =
\text{BR}(N_i \rightarrow Z \nu_\ell + Z \overline{\nu_\ell} ) &=&
\text{BR}(N_i \rightarrow h \nu_\ell + h \overline{\nu_\ell} ) = 25\%.

List of studies that have used the model file

  • Please email to update this space.
  • Golling, et al., Physics at a 100 TeV pp collider: beyond the Standard Model phenomena, arXiv:1606.00947
  • CMS, Search for heavy neutral leptons in events with three charged leptons in proton-proton collisions at 13 TeV, arXiv:1802.02965 [hep-ex]


  • Please cite [ 3 ] for the model and [ 1-2 ] for the FR/UFO files.

[1] D. Alva, T. Han, R. Ruiz, Heavy Majorana neutrinos from $W\gamma$ fusion at hadron colliders, JHEP 1502, 072 (2015), arXiv:1411.7305 [hep-ph]

[2] C. Degrande, O. Mattelear, R. Ruiz, J. Turner, Fully-Automated Precision Predictions for Heavy Neutrino Production Mechanisms at Hadron Colliders, PRD 94, 053002 (2016), arXiv:1602.06957 [hep-ph]

[3] A. Atre, T. Han, S. Pascoli and B. Zhang, The Search for Heavy Majorana Neutrinos, JHEP 0905, 030 (2009), arXiv:0901.3589 [hep-ph]

[4] S. Pascoli, R. Ruiz and C. Weiland, Safe Jet Vetoes, PLB 786, 106 (2018), arXiv:1805.09335 [hep-ph]

Attachments (13)

Download all attachments as: .zip