We assume that New Physics effects are parametrized within the Standard Model Effective Field Theory (SMEFT) written in a complete basis of gauge invariant operators up to dimension 6, commonly referred to as "Warsaw basis". We discuss all steps necessary to obtain a consistent transition to the spontaneously broken theory and several other important aspects, including the BRST-invariance of the SMEFT action for linear R ξ -gauges. The final theory is expressed in a basis characterized by SM-like propagators for all physical and unphysical fields. The effect of the non-renormalizable operators appears explicitly in triple or higher multiplicity vertices. In this mass basis we derive the complete set of Feynman rules, without resorting to any simplifying assumptions such as baryon-, lepton-number or CP conservation. As it turns out, for most SMEFT vertices the expressions are reasonably short, with a noticeable exception of those involving 4, 5 and 6 gluons. We have also supplemented our set of Feynman rules, given in an appendix here, with a publicly available Mathematica code working with the FeynRules package and producing output which can be integrated with other symbolic algebra or numerical codes for automatic SMEFT amplitude calculations.

The forthcoming precision data on lepton flavour violating (LFV) decays require precise and efficient calculations in New Physics models. In this article lepton flavour violating processes within the Minimal Supersymmetric Standard Model (MSSM) are calculated using the method based on the Flavour Expansion Theorem, a recently developed technique performing a purely algebraic mass-insertion expansion of the amplitudes. The expansion in both flavour-violating and flavour-conserving off-diagonal terms of sfermion and supersymmetric fermion mass matrices is considered. In this way the relevant processes are expressed directly in terms of the parameters of the MSSM Lagrangian. We also study the decoupling properties of the amplitudes. The results are compared to the corresponding calculations in the mass eigenbasis (i.e. using the exact diagonalization of the mass matrices). Using these methods, we consider the following processes: → γ, → 3 , → 2 , h → as well as µ → e conversion in nuclei. In the numerical analysis we update the bounds on the flavour changing parameters of the MSSM and examine the sensitivity to the forthcoming experimental results. We find that flavour violating muon decays provide the most stringent bounds on supersymmetric effects and will continue to do so in the future. Radiative → γ decays and leptonic three-body decays → 3 show an interesting complementarity in eliminating "blind spots" in the parameter space. In our analysis we also include the effects of non-holomorphic A-terms which are important for the study of LFV Higgs decays.