SmeftFR -- Feynman rules generator for the Standard Model Effective Field Theory

TL;DR

SmeftFR is a Mathematica package that automates the generation of Feynman rules for the Standard Model Effective Field Theory, including all gauge-invariant operators up to dimension 6, facilitating theoretical and phenomenological analyses.

Contribution

It provides a comprehensive, user-friendly tool for deriving Feynman rules in SMEFT with flexible options for gauge fixing, neutrino treatment, and output formats, enhancing research efficiency.

Findings

Generates complete Feynman rules for SMEFT up to dimension 6.

Supports various gauges and neutrino treatments.

Exports in multiple formats and includes a LaTeX generator.

Abstract

We present SmeftFR, a Mathematica package designed to generate the Feynman rules for the Standard Model Effective Field Theory (SMEFT) including the complete set of gauge invariant operators up to dimension~6. Feynman rules are generated with the use of FeynRules package, directly in the physical (mass eigenstates) basis for all fields. The complete set of interaction vertices can be derived including all or any chosen subset of SMEFT operators. As an option, the user can also choose preferred gauge fixing, generating Feynman rules in unitary or $R_\xi$-gauges (the latter include generation of ghost vertices). Further options allow to treat neutrino fields as massless Weyl or massive Majorana fermions. The derived Lagrangian in the mass basis can be exported in various formats supported by FeynRules, such as UFO, FeynArts, etc. Initialisation of numerical values of $d=6$ Wilson coefficients used by SmeftFR is interfaced to WCxf format. The package also includes dedicated Latex generator allowing to print the result in clear human-readable form. SmeftFR can be downloaded from www.fuw.edu.pl/smeft.

Tables

Table 1. Table 1: Dimension-6 operators other than the four-fermion ones (from ref.  [ 9 ] ). For brevity we suppress fermion chiral indices L , R 𝐿 𝑅 L,R .

$X^3$		${\phi }^6$  and  ${\phi }^4{D}^2$		${\psi }^2{\phi }^3$
$Q_G$	$f^{ABC}{G}_{\mu }^{A\nu }{G}_{\nu }^{B\rho }{G}_{\rho }^{C\mu }$	$Q_{\phi }$	${\left({\phi }^{†}\phi \right)}^3$	$Q_{e\phi }$	$\left({\phi }^{†}\phi \right)\left({\overline{l}}_p^{\prime }{e}_r^{\prime }\phi \right)$
$Q_{\tilde{G}}$	$f^{ABC}{\tilde{G}}_{\mu }^{A\nu }{G}_{\nu }^{B\rho }{G}_{\rho }^{C\mu }$	$Q_{\phi \square }$	$\left({\phi }^{†}\phi \right)\square \left({\phi }^{†}\phi \right)$	$Q_{u\phi }$	$\left({\phi }^{†}\phi \right)\left({\overline{q}}_p^{\prime }{u}_r^{\prime }\tilde{\phi }\right)$
$Q_W$	${\epsilon }^{IJK}{W}_{\mu }^{I\nu }{W}_{\nu }^{J\rho }{W}_{\rho }^{K\mu }$	$Q_{\phi D}$	${\left({\phi }^{†}{D}^{\mu }\phi \right)}^{\ast }\left({\phi }^{†}{D}_{\mu }\phi \right)$	$Q_{d\phi }$	$\left({\phi }^{†}\phi \right)\left({\overline{q}}_p^{\prime }{d}_r^{\prime }\phi \right)$
$Q_{\tilde{W}}$	${\epsilon }^{IJK}{\tilde{W}}_{\mu }^{I\nu }{W}_{\nu }^{J\rho }{W}_{\rho }^{K\mu }$
$X^2{\phi }^2$		${\psi }^{2}X\phi$		${\psi }^2{\phi }^{2}D$
$Q_{\phi G}$	${\phi }^{†}\phi G_{\mu \nu }^A{G}^{A\mu \nu }$	$Q_{eW}$	$\left({\overline{l}}_p^{\prime }{\sigma }^{\mu \nu }{e}_r^{\prime }\right){\tau }^I\phi W_{\mu \nu }^I$	$Q_{\phi l}^{\left(1\right)}$	$\left({\phi }^{†}i{}_{}{}^{\leftrightarrow }D_{\mu }^{}\phi \right)\left({\overline{l}}_p^{\prime }{\gamma }^{\mu }{l}_r^{\prime }\right)$
$Q_{\phi \tilde{G}}$	${\phi }^{†}\phi {\tilde{G}}_{\mu \nu }^A{G}^{A\mu \nu }$	$Q_{eB}$	$\left({\overline{l}}_p^{\prime }{\sigma }^{\mu \nu }{e}_r^{\prime }\right)\phi B_{\mu \nu }$	$Q_{\phi l}^{\left(3\right)}$	$\left({\phi }^{†}i{}_{}{}^{\leftrightarrow }D_{\mu }^I\phi \right)\left({\overline{l}}_p^{\prime }{\tau }^I{\gamma }^{\mu }{l}_r^{\prime }\right)$
$Q_{\phi W}$	${\phi }^{†}\phi W_{\mu \nu }^I{W}^{I\mu \nu }$	$Q_{uG}$	$\left({\overline{q}}_p^{\prime }{\sigma }^{\mu \nu }{𝒯}^A{u}_r^{\prime }\right)\tilde{\phi }{G}_{\mu \nu }^A$	$Q_{\phi e}$	$\left({\phi }^{†}i{}_{}{}^{\leftrightarrow }D_{\mu }^{}\phi \right)\left({\overline{e}}_p^{\prime }{\gamma }^{\mu }{e}_r^{\prime }\right)$
$Q_{\phi \tilde{W}}$	${\phi }^{†}\phi {\tilde{W}}_{\mu \nu }^I{W}^{I\mu \nu }$	$Q_{uW}$	$\left({\overline{q}}_p^{\prime }{\sigma }^{\mu \nu }{u}_r^{\prime }\right){\tau }^I\tilde{\phi }{W}_{\mu \nu }^I$	$Q_{\phi q}^{\left(1\right)}$	$\left({\phi }^{†}i{}_{}{}^{\leftrightarrow }D_{\mu }^{}\phi \right)\left({\overline{q}}_p^{\prime }{\gamma }^{\mu }{q}_r^{\prime }\right)$
$Q_{\phi B}$	${\phi }^{†}\phi B_{\mu \nu }{B}^{\mu \nu }$	$Q_{uB}$	$\left({\overline{q}}_p^{\prime }{\sigma }^{\mu \nu }{u}_r^{\prime }\right)\tilde{\phi }{B}_{\mu \nu }$	$Q_{\phi q}^{\left(3\right)}$	$\left({\phi }^{†}i{}_{}{}^{\leftrightarrow }D_{\mu }^I\phi \right)\left({\overline{q}}_p^{\prime }{\tau }^I{\gamma }^{\mu }{q}_r^{\prime }\right)$
$Q_{\phi \tilde{B}}$	${\phi }^{†}\phi {\tilde{B}}_{\mu \nu }{B}^{\mu \nu }$	$Q_{dG}$	$\left({\overline{q}}_p^{\prime }{\sigma }^{\mu \nu }{𝒯}^A{d}_r^{\prime }\right)\phi G_{\mu \nu }^A$	$Q_{\phi u}$	$\left({\phi }^{†}i{}_{}{}^{\leftrightarrow }D_{\mu }^{}\phi \right)\left({\overline{u}}_p^{\prime }{\gamma }^{\mu }{u}_r^{\prime }\right)$
$Q_{\phi WB}$	${\phi }^{†}{\tau }^I\phi W_{\mu \nu }^I{B}^{\mu \nu }$	$Q_{dW}$	$\left({\overline{q}}_p^{\prime }{\sigma }^{\mu \nu }{d}_r^{\prime }\right){\tau }^I\phi W_{\mu \nu }^I$	$Q_{\phi d}$	$\left({\phi }^{†}i{}_{}{}^{\leftrightarrow }D_{\mu }^{}\phi \right)\left({\overline{d}}_p^{\prime }{\gamma }^{\mu }{d}_r^{\prime }\right)$
$Q_{\phi \tilde{W}B}$	${\phi }^{†}{\tau }^I\phi {\tilde{W}}_{\mu \nu }^I{B}^{\mu \nu }$	$Q_{dB}$	$\left({\overline{q}}_p^{\prime }{\sigma }^{\mu \nu }{d}_r^{\prime }\right)\phi B_{\mu \nu }$	$Q_{\phi ud}$	$i\left({\tilde{\phi }}^{†}{D}_{\mu }\phi \right)\left({\overline{u}}_p^{\prime }{\gamma }^{\mu }{d}_r^{\prime }\right)$

Table 2. Table 2: Four-fermion operators (from ref.  [ 9 ] ). Fermion chiral ( L , R ) 𝐿 𝑅 (L,R) indices are suppressed.

$\left(\overline{L}L\right)\left(\overline{L}L\right)$		$\left(\overline{R}R\right)\left(\overline{R}R\right)$		$\left(\overline{L}L\right)\left(\overline{R}R\right)$
$Q_{ll}$	$\left({\overline{l}}_p^{\prime }{\gamma }_{\mu }{l}_r^{\prime }\right)\left({\overline{l}}_s^{\prime }{\gamma }^{\mu }{l}_t^{\prime }\right)$	$Q_{ee}$	$\left({\overline{e}}_p^{\prime }{\gamma }_{\mu }{e}_r^{\prime }\right)\left({\overline{e}}_s^{\prime }{\gamma }^{\mu }{e}_t^{\prime }\right)$	$Q_{le}$	$\left({\overline{l}}_p^{\prime }{\gamma }_{\mu }{l}_r^{\prime }\right)\left({\overline{e}}_s^{\prime }{\gamma }^{\mu }{e}_t^{\prime }\right)$
$Q_{qq}^{\left(1\right)}$	$\left({\overline{q}}_p^{\prime }{\gamma }_{\mu }{q}_r^{\prime }\right)\left({\overline{q}}_s^{\prime }{\gamma }^{\mu }{q}_t^{\prime }\right)$	$Q_{uu}$	$\left({\overline{u}}_p^{\prime }{\gamma }_{\mu }{u}_r^{\prime }\right)\left({\overline{u}}_s^{\prime }{\gamma }^{\mu }{u}_t^{\prime }\right)$	$Q_{lu}$	$\left({\overline{l}}_p^{\prime }{\gamma }_{\mu }{l}_r^{\prime }\right)\left({\overline{u}}_s^{\prime }{\gamma }^{\mu }{u}_t^{\prime }\right)$
$Q_{qq}^{\left(3\right)}$	$\left({\overline{q}}_p^{\prime }{\gamma }_{\mu }{\tau }^I{q}_r^{\prime }\right)\left({\overline{q}}_s^{\prime }{\gamma }^{\mu }{\tau }^I{q}_t^{\prime }\right)$	$Q_{dd}$	$\left({\overline{d}}_p^{\prime }{\gamma }_{\mu }{d}_r^{\prime }\right)\left({\overline{d}}_s^{\prime }{\gamma }^{\mu }{d}_t^{\prime }\right)$	$Q_{ld}$	$\left({\overline{l}}_p^{\prime }{\gamma }_{\mu }{l}_r^{\prime }\right)\left({\overline{d}}_s^{\prime }{\gamma }^{\mu }{d}_t^{\prime }\right)$
$Q_{lq}^{\left(1\right)}$	$\left({\overline{l}}_p^{\prime }{\gamma }_{\mu }{l}_r^{\prime }\right)\left({\overline{q}}_s^{\prime }{\gamma }^{\mu }{q}_t^{\prime }\right)$	$Q_{eu}$	$\left({\overline{e}}_p^{\prime }{\gamma }_{\mu }{e}_r^{\prime }\right)\left({\overline{u}}_s^{\prime }{\gamma }^{\mu }{u}_t^{\prime }\right)$	$Q_{qe}$	$\left({\overline{q}}_p^{\prime }{\gamma }_{\mu }{q}_r^{\prime }\right)\left({\overline{e}}_s^{\prime }{\gamma }^{\mu }{e}_t^{\prime }\right)$
$Q_{lq}^{\left(3\right)}$	$\left({\overline{l}}_p^{\prime }{\gamma }_{\mu }{\tau }^I{l}_r^{\prime }\right)\left({\overline{q}}_s^{\prime }{\gamma }^{\mu }{\tau }^I{q}_t^{\prime }\right)$	$Q_{ed}$	$\left({\overline{e}}_p^{\prime }{\gamma }_{\mu }{e}_r^{\prime }\right)\left({\overline{d}}_s^{\prime }{\gamma }^{\mu }{d}_t^{\prime }\right)$	$Q_{qu}^{\left(1\right)}$	$\left({\overline{q}}_p^{\prime }{\gamma }_{\mu }{q}_r^{\prime }\right)\left({\overline{u}}_s^{\prime }{\gamma }^{\mu }{u}_t^{\prime }\right)$
		$Q_{ud}^{\left(1\right)}$	$\left({\overline{u}}_p^{\prime }{\gamma }_{\mu }{u}_r^{\prime }\right)\left({\overline{d}}_s^{\prime }{\gamma }^{\mu }{d}_t^{\prime }\right)$	$Q_{qu}^{\left(8\right)}$	$\left({\overline{q}}_p^{\prime }{\gamma }_{\mu }{𝒯}^A{q}_r^{\prime }\right)\left({\overline{u}}_s^{\prime }{\gamma }^{\mu }{𝒯}^A{u}_t^{\prime }\right)$
		$Q_{ud}^{\left(8\right)}$	$\left({\overline{u}}_p^{\prime }{\gamma }_{\mu }{𝒯}^A{u}_r^{\prime }\right)\left({\overline{d}}_s^{\prime }{\gamma }^{\mu }{𝒯}^A{d}_t^{\prime }\right)$	$Q_{qd}^{\left(1\right)}$	$\left({\overline{q}}_p^{\prime }{\gamma }_{\mu }{q}_r^{\prime }\right)\left({\overline{d}}_s^{\prime }{\gamma }^{\mu }{d}_t^{\prime }\right)$
				$Q_{qd}^{\left(8\right)}$	$\left({\overline{q}}_p^{\prime }{\gamma }_{\mu }{𝒯}^A{q}_r^{\prime }\right)\left({\overline{d}}_s^{\prime }{\gamma }^{\mu }{𝒯}^A{d}_t^{\prime }\right)$
$\left(\overline{L}R\right)\left(\overline{R}L\right)$ and $\left(\overline{L}R\right)\left(\overline{L}R\right)$		$B$-violating
$Q_{ledq}$	$\left({\overline{l}}_p^{{}_{}{}^{\prime }j}{e}_r^{\prime }\right)\left({\overline{d}}_s^{\prime }{q}_t^{{}_{}{}^{\prime }j}\right)$	$Q_{duq}$	${\epsilon }^{\alpha \beta \gamma }{\epsilon }_{jk}\left[{\left(d_p^{{}_{}{}^{\prime }\alpha }\right)}^Tℂu_r^{{}_{}{}^{\prime }\beta }\right]\left[{\left(q_s^{{}_{}{}^{\prime }\gamma j}\right)}^Tℂl_t^{{}_{}{}^{\prime }k}\right]$
$Q_{quqd}^{\left(1\right)}$	$\left({\overline{q}}_p^{{}_{}{}^{\prime }j}{u}_r^{\prime }\right){\epsilon }_{jk}\left({\overline{q}}_s^{{}_{}{}^{\prime }k}{d}_t^{\prime }\right)$	$Q_{qqu}$	${\epsilon }^{\alpha \beta \gamma }{\epsilon }_{jk}\left[{\left(q_p^{{}_{}{}^{\prime }\alpha j}\right)}^Tℂq_r^{{}_{}{}^{\prime }\beta k}\right]\left[{\left(u_s^{{}_{}{}^{\prime }\gamma }\right)}^Tℂe_t^{\prime }\right]$
$Q_{quqd}^{\left(8\right)}$	$\left({\overline{q}}_p^{{}_{}{}^{\prime }j}{𝒯}^A{u}_r^{\prime }\right){\epsilon }_{jk}\left({\overline{q}}_s^{{}_{}{}^{\prime }k}{𝒯}^A{d}_t^{\prime }\right)$	$Q_{qqq}$	${\epsilon }^{\alpha \beta \gamma }{\epsilon }_{jn}{\epsilon }_{km}\left[{\left(q_p^{{}_{}{}^{\prime }\alpha j}\right)}^Tℂq_r^{{}_{}{}^{\prime }\beta k}\right]\left[{\left(q_s^{{}_{}{}^{\prime }\gamma m}\right)}^Tℂl_t^{{}_{}{}^{\prime }n}\right]$
$Q_{lequ}^{\left(1\right)}$	$\left({\overline{l}}_p^{\prime j}{e}_r^{\prime }\right){\epsilon }_{jk}\left({\overline{q}}_s^{\prime k}{u}_t^{\prime }\right)$	$Q_{duu}$	${\epsilon }^{\alpha \beta \gamma }\left[{\left(d_p^{{}_{}{}^{\prime }\alpha }\right)}^Tℂu_r^{{}_{}{}^{\prime }\beta }\right]\left[{\left(u_s^{{}_{}{}^{\prime }\gamma }\right)}^Tℂe_t^{\prime }\right]$
$Q_{lequ}^{\left(3\right)}$	$\left({\overline{l}}_p^{{}_{}{}^{\prime }j}{\sigma }_{\mu \nu }{e}_r^{\prime }\right){\epsilon }_{jk}\left({\overline{q}}_s^{{}_{}{}^{\prime }k}{\sigma }^{\mu \nu }{u}_t^{\prime }\right)$

Equations

$$Q_{\nu \nu }={\epsilon }_{jk}{\epsilon }_{mn}{\phi }^j{\phi }^m(l_{Lp}^{\prime k}{)}^T\mathbb{C}{l}_{Lr}^{\prime n}\equiv ({\tilde{\phi }}^{†}{l}_{Lp}^{\prime }{)}^T\mathbb{C}({\tilde{\phi }}^{†}{l}_{Lr}^{\prime }),$$

$$\mathcal{L}={\mathcal{L}}_{\mathrm{S}M}^{(4)}+\frac{1}{\Lambda }{C}^{\nu \nu }{Q}_{\nu \nu }^{(5)}+\frac{1}{{\Lambda }^2}\sum _X{C}^X{Q}_X^{(6)}+\frac{1}{{\Lambda }^2}\sum _f{C}^{\prime f}{Q}_f^{(6)}.$$

$$\displaystyle\left(\begin{array}[]{c}\varphi^{+}\\ \varphi^{0}\end{array}\right)$$

$$\displaystyle\left(\begin{array}[]{c}{B}_{\mu}\\ {W}^{3}_{\mu}\end{array}\right)$$

$$W_{\mu }^1$$

$$W_{\mu }^2$$

$$G_{\mu }^A$$

$$\displaystyle\begin{array}[]{l}g\,=Z_{g}\bar{g}\qquad g^{\prime}=Z_{g^{\prime}}\bar{g}^{\prime}\qquad g_{s}=Z_{g_{s}}\bar{g}_{s}\,.\end{array}$$

$${\psi }_X^{\prime }=U_{{\psi }_X}{\psi }_X,$$

$$\begin{array}[]{cc}M^{\prime}_{\nu}=-v^{2}C^{\prime\nu\nu}\;,&M^{\prime}_{e}=\frac{v}{\sqrt{2}}\,\left(\Gamma_{e}-\frac{v^{2}}{2}C^{\prime e\varphi}\right),\\[5.69054pt] M^{\prime}_{u}=\frac{v}{\sqrt{2}}\,\left(\Gamma_{u}-\frac{v^{2}}{2}C^{\prime u\varphi}\right),&M^{\prime}_{d}=\frac{v}{\sqrt{2}}\,\left(\Gamma_{d}-\frac{v^{2}}{2}C^{\prime d\varphi}\right).\end{array}$$

$$K=U_{u_L}^{†}{U}_{d_L},U=U_{e_L}^{†}{U}_{{\nu }_L}.$$

Full text

††journal: Computer Physics Communications

PROGRAM SUMMARY

Manuscript Title:

SmeftFR – Feynman rules generator for the Standard Model Effective Field Theory

Authors: Athanasios Dedes, Michael Paraskevas, Janusz Rosiek, Kristaq Suxho, Lampros Trifyllis

Program Title: SmeftFR v2.0

Journal Reference:

Catalogue identifier:

Licensing provisions: None

Programming language: Mathematica 11.3 (earlier versions were reported to have problems running this code)

Computer: any running Mathematica

Operating system: any running Mathematica

RAM: allocated dynamically by Mathematica, at least 4GB total RAM suggested

Number of processors used: allocated dynamically by Mathematica

Supplementary material: None

Keywords: Standard Model Effective Field Theory, Feynman rules, unitary and $R_{\xi}$ gauges

*Classification:

11.1 General, High Energy Physics and Computing,

4.4 Feynman diagrams,

5 Computer Algebra.

External routines/libraries:* Wolfram Mathematica program

Subprograms used: FeynRules v2.3 package

*Nature of problem:

*Automatised generation of Feynman rules in physical (mass) basis for the Standard Model Effective Field Theory with user defined operator subset and gauge fixing terms selection.

*Solution method:

*Implementation (as the Mathematica package) of the results of Ref. [1] in the FeynRules package [2], including dynamic “model files” generation.

Restrictions: None

Unusual features: None

Additional comments: None

Running time: depending on control variable settings, from few minutes for a selected subset of few SMEFT operators and Feynman rules generation up to several hours for generating UFO output for maximum 60-operator set (using Mathematica 11.3 running on a personal computer)

1 Introduction

Despite lack of direct discoveries of new particles at LHC, indirect experiments searching for cosmological dark matter and dark energy, neutrino masses and mixing, heavy meson decays, anomalous magnetic moment of muon, point to existence of new physics phenomena that go beyond the Standard Model (SM) [1, 2, 3] of particle physics. It has become customary in recent years to parameterise the potential effects of such New Physics (NP) in terms of the so-called SM Effective Field Theory (SMEFT) [4, 5, 6]111For a recent review see ref. [7].. In this theory, SM is extended by a complete set of gauge invariant operators constructed only by the SM fields which are considered to be light in mass with respect to, yet unknown, particles responsible for NP phenomena. Such (non-renormalisable) operators can be classified according to their growing mass dimension, with couplings (usually called Wilson coefficients) suppressed by respectively inverse powers of a typical mass scale $\Lambda$ of the NP extension.

As observable effects of new operators are typically suppressed by powers of $v/\Lambda$ (where $v$ is the true vacuum expectation value of the Higgs field), in many analyses it is sufficient to restrict the maximal dimension of new operators to $d\leq 6$. Classification of the complete independent set of the gauge invariant operators up to $d=6$ was conducted in an older study by Buchmüller and Wyler [8] and more recently put in a non-redundant form in ref. [9], where the basis of such operators, usually referred as the “Warsaw basis”, has been given. Suppressing the flavour indices of the fields and not counting hermitian conjugated operators, Warsaw basis contains 59+1 baryon-number conserving and 4 baryon-number violating operators.

In its most general version, SMEFT is a hugely complicated model. Including all possible CP- and flavour-violating interactions, it contains 2499 free parameters. Due to large number and complicated structure of new terms in the Lagrangian, theoretical calculations of physical processes within the SMEFT can be very challenging — it is enough to notice that the number of primary vertices when SMEFT is quantised in $R_{\xi}$-gauges, printed for the first time in ref. [10], is almost 400 without counting the hermitian conjugates. Thus, it is important to develop technical methods and tools facilitating such calculations, starting from developing the universal set of the Feynman rules for propagators and vertices for physical fields, after spontaneous symmetry breaking (SSB) of the full effective theory. The initial version of relevant package, SmeftFR v1, was announced and briefly described for the first time in Appendix B of ref. [10]. In this paper we present SmeftFR v2.0, a Mathematica symbolic language package generating Feynman rules in several formats, based on the formulae developed in ref. [10]. Comparing to version 1, the new code has been extended with a number of additional options and capabilities. It has also been made and tested to be compatible with many other publicly available high-energy physics related computer codes accepting standardised input and output data formats. Features of the presented package contain:

SmeftFR is able to generate interactions in the most general form of the SMEFT Lagrangian, without any restrictions on the structure of flavour violating terms and on CP-, lepton- or baryon-number conservation.222However, we do restrict ourselves to linear realisations of the SSB. Feynman rules are expressed in terms of physical SM fields and canonically normalised Goldstone and ghost fields. Expressions for interaction vertices are analytically expanded in powers of inverse New Physics scale $1/\Lambda$, with all terms of dimension higher than $d=6$ consistently truncated.

• 2.

SmeftFR is written as an overlay to FeynRules package [11], used as the engine to generate Feynman rules.

• 3.

Including the full set of SMEFT parameters in model files for FeynRules may lead to very slow computations. SmeftFR can generate FeynRules model files dynamically, including only the user defined subset of higher dimension operators. It significantly speeds up the calculations and produces simpler final result, containing only the Wilson coefficients relevant for a process chosen to analyse.

• 4.

Feynman rules can be generated in the unitary or in linear $R_{\xi}$-gauges by exploiting four different gauge-fixing parameters $\xi_{\gamma},\xi_{Z},\xi_{W},\xi_{G}$ for thorough amplitude checks. In the latter case also all relevant ghost vertices are obtained.

• 5.

Feynman rules are calculated first in Mathematica/FeynRules format. They can be further exported in other formats: UFO [12] (importable to Monte Carlo generators like MadGraph5_aMC@NLO 5  [13], Sherpa [14], CalcHEP [15], Whizard [16, 17]), FeynArts [18] which generates inputs for loop amplitude calculators like FeynCalc [19], or FormCalc [20], and others output types supported by FeynRules .

• 6.

SmeftFR provides a dedicated Latex generator, allowing to display vertices and analytical expressions for Feynman rules in clear human readable form, best suited for hand-made calculations.

• 7.

SmeftFR is interfaced to the WCxf format [21] of Wilson coefficients. Numerical values of SMEFT parameters in model files can be read from WCxf JSON-type input produced by other computer packages written for SMEFT. Alternatively, SmeftFR can translate FeynRules model files to the WCxf format.

• 8.

Further package options allow to treat neutrino fields as massless Weyl or (in the case of non-vanishing dimension-5 operator) massive Majorana fermions, to correct signs in 4-fermion interactions not yet fully supported by FeynRules and to perform some additional operations as described later in this manual.

Feynman rules derived in ref. [10] using the SmeftFR package have been used successfully in many articles including refs. [22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32] and have passed certain non-trivial tests, such as gauge-fixing parameter independence of the $S$-matrix elements, validity of Ward identities, cancellation of infinities in loop calculations, etc.

We note here in passing, that there is a growing number of publicly available codes performing computations related to SMEFT. These include, Wilson [33], DSixTools [34], MatchingTools [35], which are codes for running and matching Wilson coefficients, SMEFTsim [36], a package for calculating tree level observables, CoDEx [37] or a version of SARAH code [38], that calculate Wilson Coefficients after the decoupling of a more fundamental theory, and finally, DirectDM [39], a code for dark matter EFT. To a degree, these codes (especially the ones supporting WCxf format) can be used in conjunction with SmeftFR. For example, some of them can provide the numerical input for Wilson coefficients of higher dimensional operators at scale $\Lambda$, while others, the running of these coefficients from that scale down to the EW one. Alternatively, Feynman rules evaluated by SmeftFR can be used with Monte Carlo generators to test the predictions of other packages.

The paper is organised as follows. After this Introduction, in Sec. 2 we define the notation and conventions, listing for reference the operator set in Warsaw basis and the formulae for transition to the mass basis. In Sec. 3, we present the structure of the code, installation procedure and available functions. Section 4, contains examples of programs generating the Feynman rules in various formats. We conclude in Sec. 5.

2 SMEFT Lagrangian in Warsaw and mass basis

The classification of higher order operators in SMEFT is done in terms of fields in electroweak basis, before the Spontaneous Symmetry Breaking. SmeftFR uses the so-called “Warsaw basis” [9] as a starting point to calculate physical states in SMEFT and their interactions. For easier reference we copy here from ref. [9] Tables 1 and 2 containing the full collection of $d=6$ operators in Warsaw basis.333We do not list here all details of conventions used — they are identical to these listed in refs. [9, 10]. In addition, we include the single lepton flavour violating dimension-5 operator:

[TABLE]

where $\mathbb{C}$ is the charge conjugation matrix.

The SMEFT Lagrangian is the sum of the dimension-4 terms and operators defined in Tables 1 and 2:

[TABLE]

Physical fields in SMEFT are obtained after the SSB. We follow here the systematic presentation (and notation) of ref. [10]. In the gauge and Higgs sectors physical and Goldstone fields ($h,G^{0},G^{\pm},W^{\pm}_{\mu},Z^{0}_{\mu},A_{\mu}$) are related to initial (Warsaw basis) fields ($\varphi,W_{\mu}^{i},B_{\mu},G_{\mu}^{A}$) by the normalisation constants:

[TABLE]

In addition, Feynman rules for physical fields are expressed in terms of effective gauge couplings, chosen to preserve the natural form of covariant derivative:

[TABLE]

In $d=6$ SMEFT $SU(2)$ and $SU(3)$ gauge field and gauge normalisation constants are equal, $Z_{g}=Z_{W}$, $Z_{g_{s}}=Z_{G}$. In addition $Z_{g^{\prime}}=1-v^{2}C^{\varphi B}$ and $Z_{G^{+}}=1$. Complete expressions for the field normalisation constants, $Z_{X}$, for the corrected Higgs field VEV, $v$, and for the gauge and Higgs boson masses, $M_{Z}$, $M_{W}$ and $M_{h}$, including corrections from the full set of dimension-5 and -6 operators, are given in ref. [10]. SmeftFR recalculates these corrections by taking into account only the subset of non-vanishing SMEFT Wilson coefficients chosen by the user, as described in Sec. 3.

The basis in the fermion sector is not fixed by the structure of gauge interactions and allows for unitary rotation freedom in the flavour space:

[TABLE]

with $\psi=\nu,e,u,d$ and $X=L,R$. We choose the rotations such that $\psi_{X}$ eigenstates correspond to real and non-negative eigenvalues of $3\times 3$ fermion mass matrices:

[TABLE]

The fermion flavour rotations can be adsorbed in redefinitions of Wilson coefficients, as a trace leaving CKM and PMNS matrices (denoted in SmeftFR as $K$ and $U$ respectively) multiplying them. It is important to note a departure from ref. [10] in the definition of these matrices, assumed here to be unitary products of the form444In ref. [10], these matrices are defined as non-unitary ones, $K=U_{u_{L}}^{\dagger}U_{d_{L}}(1+C^{\varphi q(3)})$ and $U=(1+C^{\varphi l(3)})U_{e_{L}}^{\dagger}U_{\nu_{L}}$. It is however more natural and economical to use the unitary definitions of eq. (2.18) instead. The conventions for $K$ and $U$ of ref. [10] can be still enforced in SmeftFR by editing the file code/smeft_variables.m and changing there the value of control variable SMEFT$Unitary=0 to SMEFT$Unitary=1.:

[TABLE]

In practice, comparing to results of ref. [10] the altered definition of eq. (2.18) affects only four vertices, namely couplings of $W$ boson and charged Goldstone to fermions: $W^{+}(G^{+})$-$\bar{u}^{I}$-$d^{J}$, $W^{+}(G^{+})$-$\bar{\nu}^{I}$-$e^{J}$ and their Hermitian conjugates.

The complete list of redefinitions of flavour-dependent Wilson coefficients is given in Table 4 of ref. [10]. After rotations, they are defined in so called “Warsaw mass” basis (as also described in WCxf standard [21]). SmeftFR assumes that the numerical values of Wilson coefficients are given in this particular basis.

In summary, Feynman rules generated by the SmeftFR package describe interactions of SMEFT physical (mass eigenstates) fields, with numerical values of Wilson coefficients defined in the same (“Warsaw mass”) basis.

It is also important to stress that in the general case of lepton number flavour violation, with non-vanishing Weinberg operator of eq. (2.1), neutrinos are massive Majorana spinors, whereas under the assumption of $L$-conservation they can be regarded as massless Weyl spinors. As described in the next Section, SmeftFR is capable to generate Feynman rules for neutrino interactions in both cases, depending on the choice of initial options. One should remember that treating neutrinos as Majorana particles requires special set of rules for propagators, vertices and diagram combinatorics. We follow here the treatment of refs. [40, 41, 10, 42].

3 Deriving SMEFT Feynman rules with SmeftFR package

3.1 Installation

SmeftFR package works using the FeynRules system, so both need to be properly installed first. A recent version and installation instructions for the FeynRules package can be downloaded from the address:

https://feynrules.irmp.ucl.ac.be

SmeftFR has been tested with FeynRules version 2.3.

Standard FeynRules installation assumes that the new models description is put into Model subdirectory of its main tree. We follow this convention, so that SmeftFR archive should be unpacked into

Models/SMEFT_N_NN

catalogue, where N_NN denotes the package version (currently version 2.00). After installation, Models/SMEFT_N_NN contains the following files and subdirectories listed in Table 3.

Before running the package, one needs to set properly the main FeynRules installation directory, defining the $FeynRulesPath variable at the beginning of smeft_init.m and smeft_outputs.m files. For non-standard installations (not advised!), also the variable SMEFT$Path has to be updated accordingly.

3.2 Code structure

The most general version of SMEFT, including all possible flavour violating couplings, is very complicated. Symbolic operations on the full SMEFT Lagrangian, including complete set of dimension 5 and 6 operators and with numerical values of all Wilson coefficients assigned are time consuming and can take hours or even days on a standard personal computer. For most of the physical applications it is sufficient to derive interactions only for a subset of operators.555Eventually, operators must be selected with care as in general they may mix under renormalisation [43, 44, 45].

To speed up the calculations, SmeftFR can evaluate Feynman rules for a chosen subset of operators only, generating dynamically the proper FeynRules “model files”. The calculations are divided in two stages, as illustrated in flowchart of Fig. 1. First, the SMEFT Lagrangian is initialised in Warsaw basis and transformed to mass eigenstates basis analytically, truncating all terms of the order ${\cal O}(\frac{1}{\Lambda^{3}})$ and higher. To speed up the program, at this stage all flavour parameters are considered to be tensors with indices without assigned numerical values (they are “Internal” parameters in FeynRules notation). The resulting mass basis Lagrangian and Feynman rules written in Mathematica format are stored on disk. In the second stage, the previously generated output can be used together with new “model file”, this time containing numerical values of (“External”) parameters, to export mass basis SMEFT interactions in various commonly used external formats such as Latex, WCxf and standard FeynRules supported interfaces – UFO, FeynArts and others.

3.3 Model initialisation

In the first step, the relevant FeynRules model files must be generated. This is done by calling the function:

SMEFTInitializeModel[Option1 $\to$ Value1, Option2 $\to$ Value2, $\ldots$]

with the allowed options listed in Table 4.

Names of operators used in SmeftFR are derived from the subscript indices of operators listed in Tables 1 and 2, with obvious transcriptions of “tilde” symbol and Greek letters to Latin alphabet. By default, all possible 59+1+4 SMEFT ($d=6$) operator classes are included in calculations, which is equivalent to setting:

Operators $\to$ { "G", "Gtilde", "W", "Wtilde", "phi", "phiBox", "phiD", "phiW", "phiB", "phiWB", "phiWtilde", "phiBtilde", "phiWtildeB", "phiGtilde", "phiG", "ephi", "dphi", "uphi", "eW", "eB", "uG", "uW", "uB", "dG", "dW", "dB", "phil1", "phil3", "phie", "phiq1", "phiq3", "phiu", "phid", "phiud", "ll", "qq1", "qq3", "lq1", "lq3", "ee", "uu", "dd", "eu", "ed", "ud1", "ud8", "le", "lu", "ld", "qe", "qu1", "qu8", "qd1", "qd8", "ledq", "quqd1", "quqd8", "lequ1", "lequ3", "vv", "duq", "qqu", "qqq", "duu" }

To speed up the derivation of Feynman rules and to get more compact expressions, the user can restrict the list above to any preferred subset of operators.

SmeftFR is fully integrated with the WCxf standard. Apart from numerically editing Wilson coefficients in FeynRules model files, reading them from the WCxf input is the only way of automatic initialisation of their numerical values. Such an input format is exchangeable between a larger set of SMEFT-related public packages [21] and may help to compare their results.

An additional advantage of using WCxf input format comes in the flavour sector of the theory. Here, Wilson coefficients are in general tensors with flavour indices, in many cases symmetric under various permutations. WCxf input requires initialisation of only the minimal set of flavour dependent Wilson coefficients, those which could be derived by permutations are also automatically properly set.666We would like to thank D. Straub for supplying us with a code for symmetrisation of flavour-dependent Wilson coefficients.

Further comments concern MajoranaNeutrino and Correct4Fermion options. They are used to modify the analytical expressions only for the Feynman rules, not at the level of the mass basis Lagrangian from which the rules are derived. This is because some FeynRules interfaces, like UFO, intentionally leave the relative sign of 4-fermion interactions uncorrected777B. Fuks, private communication., as it is later changed by Monte Carlo generators like MadGraph5. Correcting the sign before generating UFO output would therefore lead to wrong final result. Similarly, treatment of neutrinos as Majorana fields could not be compatible with hard coded quantum number definitions in various packages. On the other hand, in the manual or symbolic computations it is convenient to have from the start the correct form of Feynman rules, as done by SmeftFR when both options are set to their default values.

SMEFTInitializeModel routine does not require prior loading of FeynRules package. After execution, it creates in the output subdirectory three model files listed in Table 5. Parameter files generated by SMEFTInitializeModel contain also definitions of SM parameters, copied from templates definitions/smeft_par_head_WB.fr and definitions/smeft_par_head_MB.fr. The values of SM parameters can be best updated directly by editing the template files and the header of the code/smeft_variables.m file, otherwise they will be overwritten in each rerun of SmeftFR initialisation routines.

As mentioned above, in all analytical calculations performed by SmeftFR , terms suppressed by ${\cal O}(1/\Lambda^{3})$ or higher $1/\Lambda$ powers are always neglected. Therefore, the resulting Feynman rules can be consistently used to calculate physical observables, symbolically or numerically by Monte Carlo generators, up to the linear order in dimension-6 operators. This information is encoded in FeynRules SMEFT model files by assigning the “interaction order” parameter NP=1 to each Wilson coefficients and setting in smeft_field_WB.fr and smeft_field_MB.fr the limits:

[TABLE]

If one wishes to include (e.g. for testing purposes) terms quadratic in Wilson coefficients in the automatic cross section calculations performed by matrix element and Monte Carlo generators, one should edit both model files and increase the allowed order of NP contributions by changing the definition of M$InteractionOrderLimit, setting there {NP,2}.

An additional remark concerns the value of neutrino masses. In mass basis, the neutrino masses are equal to $-v^{2}C_{\nu\nu}^{II}$  [see eq. (2.17)]. Thus, the numerical values of $C_{\nu\nu}^{II}$ coefficients should be real and negative. If positive or complex values of $C_{\nu\nu}^{II}$ are given in the WCxf input file, then the SMEFTInitializeModel routine evaluates neutrino masses as $M_{\nu^{I}}=v^{2}|C_{\nu\nu}^{II}|$.

3.4 Calculation of mass basis Lagrangian and Feynman rules

By loading the FeynRules model files the derivation of SMEFT Lagrangian in mass basis is performed by calling the following sequence of routines:

[TABLE]

The calculation time may vary considerably depending on the choice of operator (sub-)set and gauge fixing conditions chosen. For the full list of SMEFT $d=5$ and $d=6$ operators and in $R_{\xi}$-gauges, one can expect CPU time necessary to evaluate all Feynman rules, from about an hour to many hours on a typical personal computer, depending on its speed capabilities.

One should note that when neutrinos are treated as Majorana particles, (as necessary in case of non-vanishing Wilson coefficient of $d=5$ Weinberg operator), their interactions involve lepton number non-conservation. When FeynRules is dealing with them it produces warnings of the form:

*QN::NonConserv: Warning: non quantum number conserving vertex encountered!

Quantum number LeptonNumber not conserved in vertex $\ldots$ *

Obviously such warnings should be ignored.

Evaluation of Feynman rules for vertices involving more than two fermions is not fully implemented yet in FeynRules. To our experience, apart from the issue of relative sign of four fermion diagrams mentioned earlier, particularly problematic was the correct automatic derivation of quartic interactions with four Majorana neutrinos and similar vertices which violate $B$- and $L$-quantum numbers. For these special cases, SmeftFR overwrites the FeynRules result with manually calculated formulae encoded in Mathematica format.

Another remark concerns the hermicity property of the SMEFT Lagrangian. For some types of interactions, e.g. four-fermion vertices involving two-quarks and two-leptons, the function CheckHermicity provided by FeynRules reports non-Hermitian terms in the Lagrangian. However, such terms are actually Hermitian if permutation symmetries of indices of relevant Wilson coefficients are taken into account. Such symmetries are automatically imposed if numerical values of Wilson coefficients are initialized with the use of SMEFTInitializeMB or SMEFTToWCXF routines (see Sections 3.5 and 3.5.1).

Results of the calculations are collected in file output/smeft_feynman_rules.m. The Feynman rules and pieces of the mass basis Lagrangian for various classes of interactions are stored in the variables with self-explanatory names listed in Table 6.

File output/smeft_feynman_rules.m contains also expressions for the normalisation factors relating Higgs and gauge fields and couplings in the Warsaw and mass basis. Namely, variables Hnorm, G0norm, GPnorm, AZnorm[i,j], Wnorm, Gnorm, correspond to, respectively, $Z_{h}^{-1}$, $Z_{G^{0}}^{-1}$, $Z_{G^{+}}^{-1}$ $\hat{Z}_{AZ}^{-1}$, $Z_{W}^{-1}$ and $Z_{G}^{-1}$ in (2.13). In addition, formulae for tree level corrections to SM mass parameters and Yukawa couplings are stored in variables SMEFT$vev, SMEFT$MH2, SMEFT$MW2, SMEFT$MZ2, SMEFT$YL[i,j], SMEFT$YD[i,j] and SMEFT$YU[i,j].

It is important to note that although at this point the Feynman rules for the mass basis Lagrangian are already calculated, definitions for fields and parameters used to initialise the SMEFT model in FeynRules are still given in Warsaw basis. To avoid inconsistencies, it is strongly advised to quit the current Mathematica kernel and start new one reloading the mass basis Lagrangian together with the compatible model files with fields defined also in mass basis, as described next in Sec. 3.5. All further calculations should be performed within this new kernel.

3.5 Interfaces

SmeftFR output in some of portable formats must be generated from the SMEFT Lagrangian transformed to mass basis, with all numerical values of parameters initialised. As FeynRules does not allow for two different model files loaded within a single Mathematica session, one needs to quit the kernel used to run routines necessary to obtain Feynman rules and, as described in previous Section, start a new Mathematica kernel. Within it, the user must reload FeynRules and SmeftFR packages and call the following routine:

SMEFTInitializeMB[ Options ]

Allowed options are given in Table 7. After call to SMEFTInitializeMB, mass basis model files are read and the mass basis Lagrangian is stored in a global variable named SMEFTMBLagrangian for further use by interface routines.

3.5.1 WCxf input and output

Translation between FeynRules model files and WCxf format is done by the functions SMEFTToWCXF and WCXFToSMEFT. They can be used standalone and do not require loading FeynRules and calling first SMEFTInitializeMB routine to work properly.

Exporting numerical values of Wilson coefficients of operators in the WCxf format is done by the function:

SMEFTToWCXF[ SMEFT_Parameter_File, WCXF_File ]

where the arguments SMEFT_Parameter_File, WCXF_File define the input model parameter file in the FeynRules format and the output file in the WCxf JSON format, respectively. The created JSON file can be used to transfer numerical values of Wilson coefficients to other codes supporting WCxf format. Note that in general, the FeynRules model files may contain different classes of parameters, according to the Value property defined to be a number (real or complex), a formula or even not defined at all. Only the Wilson coefficients with Value defined to be a number are transferred to the output file in WCxf format.

Conversely, files in WCxf format can be translated to FeynRules parameter files using:

WCXFToSMEFT[ WCXF_File, SMEFT_Parameter_File Options]

with the allowed options defined in Table 8.

3.5.2 Latex output

SmeftFR provides a dedicated Latex generator (not using the generic FeynRules Latex export routine). Its output has the following structure:

For each interaction vertex, the diagram is drawn, using the axodraw style [46]. Expressions for Feynman rules are displayed next to corresponding diagrams.

• 2.

In analytical expressions, all terms multiplying a given Wilson coefficient are collected together and simplified.

• 3.

Long analytical expressions are automatically broken into many lines using breakn style (this does not always work perfectly but the printout is sufficiently readable).

Latex output is generated by the function:

SMEFTToLatex[ Options ]

with the allowed options listed in Table 9. The function SMEFTToLatex assumes that the variables listed in Table 6 are initialised. It can be called either after executing relevant commands, described in Sec. 3.4, or after reloading the mass basis Lagrangian with the SMEFTInitializeMB routine, see Sec. 3.5.

Latex output is stored in output/latex subdirectory, split into smaller files each containing one primary vertex. The main file is named smeft_feynman_rules.tex. The style files necessary to compile Latex output are supplied with the SmeftFR distribution.

Note that the correct compilation of documents using “axodraw.sty” style requires creating intermediate Postscript file. Programs like pdflatex producing directly PDF output will not work properly. One should instead use e.g.:

latex smeft_feynman_rules.tex

dvips smeft_feynman_rules.dvi

ps2pdf smeft_feynman_rules.ps

The smeft_feynman_rules.tex does not contain analytical expressions for five and six gluon vertices. Such formulae are very long (multiple pages, hard to even compile properly) and not useful for hand-made calculations. If such vertices are needed, they should be rather directly exported in some other formats as described in the next subsection.

Other details not printed in the Latex output, such as, the form of field propagators, conventions for parameters and momenta flow in vertices (always incoming), manipulation of four-fermion vertices with Majorana fermions etc, are explained thoroughly in the Appendices A1–A3 of ref. [10].

3.5.3 Standard FeynRules interfaces

After calling the initialisation routine SMEFTInitializeMB, the output to UFO, FeynArts and other formats supported by FeynRules interfaces, can be generated using FeynRules commands and options from the mass basis Lagrangian stored in the SMEFTMBLagrangian variable. For instance, one could call:

WriteUFO[ SMEFTMBLagrangian, Output $\to$ ”output/UFO”, AddDecays $\to$ False, …]

WriteFeynArtsOutput[ SMEFTMBLagrangian, Output $\to$ ”output/FeynArts”, …]

and similarly for other formats.

It is important to note that FeynRules interfaces like UFO or FeynArts generate their output starting from the level of SMEFT mass basis Lagrangian. Thus, options of SMEFTInitializeModel function like MajoranaNeutrino and Correct4Fermion (see Table 4) have no effect on output generated by the interface routines. As explained in Section 3.3 they affect only the expressions for Feynman rules.

If four-fermion vertices are included in SMEFT Lagrangian, UFO produces warning messages of the form:

Warning: Multi-Fermion operators are not yet fully supported!

Therefore, the output for four-fermion interactions in UFO or other formats must be treated with care and limited trust — performing appropriate checks are left to users’ responsibility. To our experience, implementation in FeynRules of baryon and lepton number violating four-fermion interactions, with charge conjugation matrix appearing explicitly in vertices, is even more problematic. Thus, for safety in current SmeftFR version (2.00) such terms are never included in SMEFTMBLagrangian variable, eventually they can be passed to interface routines separately via the BLViolatingLagrangian variable.

Exporting to UFO or other formats can take a long time, even several hours for $R_{\xi}$-gauges and complete SMEFT Lagrangian with fully general flavour structure and all numerical values of parameters initialised.

Finally, it is important to stress here that our Feynman rules communicate properly with MadGraph5 and FeynArts. In particular, we ran without errors test simulations in MadGraph5 using UFO model files produced by SmeftFR v2. Similar tests were performed with amplitude generation for sample processes using SmeftFR v2 FeynArts output.

4 Sample programs

After setting the variable $FeynRulesPath to correct value, in order to evaluate mass basis SMEFT Lagrangian and analytical form of Feynman rules one can use the following sequence of commands:

SMEFT$MajorVersion = "2";

SMEFT$MinorVersion = "00";

SMEFTPath = FileNameJoin[{FeynRulesPath, "Models", "SMEFT_" <>

                                    SMEFT$MajorVersion <> "_" <> SMEFT$MinorVersion}];

Get[ FileNameJoin[$FeynRulesPath,"FeynRules.m"] ];

Get[ FileNameJoin[ SMEFT$Path, "code", "smeft_package.m"] ];

OpList = { "G", "Gtilde", "W", "Wtilde", "phi", "phiBox", "phiD", "phiW", "phiB", "phiWB", "phiWtilde", "phiBtilde", "phiWtildeB", "phiGtilde", "phiG", "ephi", "dphi", "uphi", "eW", "eB", "uG", "uW", "uB", "dG", "dW", "dB", "phil1", "phil3", "phie", "phiq1", "phiq3", "phiu", "phid", "phiud", "ll", "qq1", "qq3", "lq1", "lq3", "ee", "uu", "dd", "eu", "ed", "ud1", "ud8", "le", "lu", "ld", "qe", "qu1", "qu8", "qd1", "qd8", "ledq", "quqd1", "quqd8", "lequ1", "lequ3", "vv", "duq", "qqu", "qqq", "duu" };

SMEFTInitializeModel[ Operators -> OpList, Gauge -> Rxi, WCXFInitFile -> "wcxf_input_file_with_path.json" ];

SMEFTLoadModel[ ];

SMEFTFindMassBasis[ ];

SMEFTFeynmanRules[ ];

or alternatively rerun the supplied programs: the notebook SmeftFR-init.nb or the text script smeft_fr_init.m.

As described before, Latex, WCxf, UFO and FeynArts formats can be exported after rerunning first SmeftFR-init.nb or equivalent set of commands generating file smeft_feynman_rules.m containing the expressions for the mass basis Lagrangian. Then, the user needs to start a new Mathematica kernel and rerun the notebook file SmeftFR-interfaces.nb or the script smeft_fr_interfaces.m. Alternatively, one can manually type the commands, if necessary changing some of their options as described in previous sections:

Get[ FileNameJoin[{$FeynRulesPath,"FeynRules.m"}] ];

Get[ FileNameJoin[{SMEFT$Path, "code", "smeft_package.m"}] ];

SMEFTInitializeMB[ ];

SMEFTToWCXF[ FileNameJoin[{SMEFT$Path, "output", "smeft_par_MB.fr"}],

                        FileNameJoin[{SMEFT$Path, "output", "smeft_wcxf_MB.json"}] ];

SMEFToLatex[ ];

WriteUFO[ SMEFTMBLagrangian, $\ldots$ ];

WriteFeynArtsOutput[ SMEFTMBLagrangian, $\ldots$ ];

Feynman rules calculated for the full SMEFT operator set and two possible gauge choices (unitary and $R_{\xi}$), including also outputs in Latex, UFO and FeynArts formats, are stored in the subdirectories full_rxi_results and full_unitary_results. They can be directly used, without rerunning the SmeftFR package. However, such a general output with all operators and numerical values for all 2499 SMEFT parameters initialised, is huge and importing it directly to other SMEFT codes may cause them to work very slowly.

5 Summary

The proliferation of the primitive vertices in SMEFT, even in non-redundant basis of effective field operators, suggests a certain kind of computational aid. During the last few years, such an effort has been intensive between high energy physicists. Aiming at this direction, we present here a code, named SmeftFR, that generates Feynman rules in SMEFT with $\text{dimension}\leq 6$ operators given in “Warsaw basis” [9] without restricting to specific flavour structure, $CP$-, $B$- nor $L$-number conservation. We have exploited the quantisation steps of SMEFT in unitary and $R_{\xi}$-gauges following a procedure described in ref. [10]. SmeftFR has been written on top of the Mathematica package FeynRules.

In this article, we describe how to use SmeftFR package in order to produce Feynman rules for a selection of operators relevant to observable (or observables) under study, with further options to handle massive Majorana or massless Weyl neutrinos. The output of the package can be printed in Latex or exported in various formats supported by FeynRules, such as UFO, FeynArts, etc. Input parameters for Wilson coefficients used in SmeftFR can communicate with WCxf format for further numerical handling. Feynman rules are given in SMEFT mass basis and in both unitary or linear $R_{\xi}$-gauges for further computational checks.

The current version of SmeftFR code and its manual can be downloaded from the address

www.fuw.edu.pl/smeft

We believe that SmeftFR is bridging a gap between the effective SM Lagrangian all the way down to amplitude calculations required by the experimental analyses.

Acknowledgements

The work of JR was supported in part by Polish National Science Centre under research grants DEC-2015/19/B/ST2/02848, and DEC-2015/18/M/ST2/00054. MP acknowledges financial support by the Polish National Science Center under the Beethoven series grant number DEC-2016/23/G/ST2/04301. KS would like to thank the Greek State Scholarships Foundation (IKY) for full financial support through the Operational Programme “Human Resources Development, Education and Lifelong Learning, 2014-2020”. LT is supported by the Onassis Foundation — Scholarship ID: G ZO 029-1/2018-2019. JR would like to thank University of Ioannina and CERN for the hospitality during his stays there.