Lectures on SMEFT

This document contains lectures on SMEFT, which is an effective field theory of the degrees of freedom of the Standard Model. The material is at a basic, introductory level, without assuming any prior knowledge of effective field theory techniques. The main focus is on phenomenological applications of SMEFT in collider, flavor, and low-energy physics.


Goals, notation, conventions
This is a write-up of my lectures given on several occasions (Saclay'17, Orsay'20, Les Houches'21, Bhubaneswar'22, Orsay'22, Les Houches'23) about a specific effective field theory (EFT) called the Standard Model Effective Field Theory or SMEFT in short. My intention was to prepare these lectures at a very basic, introductory level, without assuming any prior knowledge of EFT techniques. On the other hand, I assume the reader is versed in Quantum Field Theory (QFT) roughly at the Peskin-Schroeder level [1]. The focus is on phenomenological applications of SMEFT, especially in Higgs, electroweak, flavor, and low-energy physics. This document is not meant to be a SMEFT review. It leaves out or barely touches upon many important topics (running, matching, Hilbert series, collider simulations, on-shell techniques, …), it does not try to summarize all recent developments in this field, and it does not attempt to provide references to all the important papers in the vast SMEFT literature. For a more general introduction to EFT, I recommend Refs. [2][3][4]. For a broader scope of SMEFT topics and more references see Refs. [5,6].
Here is the layout of these lectures. Section 1 describes the place of SMEFT in the ladder of effective theories, from extremely low energies to large scales that are not directly accessible by experiment. In Sect. 2 I will discuss in some details the assumptions under which SMEFT is the relevant formalism to describe physics above the electroweak scale. Next, in Sect. 3 I explain how to systematically construct the SMEFT Lagrangian. Various equivalent representations of the SMEFT Lagrangian, the so-called bases, are discussed in Sect. 4. In Sect. 5 I will discuss how the observables, at the Large Hadron Collider (LHC) and in other experiments, depend on the Wilson coefficients of higher-dimensional SMEFT operators. Finally, Sect. 6 is devoted to the phenomenological importance of new SMEFT sources of C P violation.
I work with the mostly minus Minkowski metric η μν = (1, −1, −1, −1). The metric is used to raise and lower indices, e.g. A μ = η μν A ν . As usual, repeated Lorentz and other indices are implicitly summed over, unless otherwise noted. Since the Lorentz contractions are unambiguous, sometimes I may write contracted Lorentz indices on the same level (e.g. A μ A μ instead of A μ A μ ) if this improves the aesthetics (usually when there are many other indices). The sign convention for the totally anti-symmetric Levi-Civita tensor μνρα is 0123 = 1, which implies 0123 = −1. When I refer to a vector I always mean a Lorentz 4-vector. For 3-vectors I use the bold notation rather than an arrow top: x ≡ x.
I always use the natural units,h = c = 1. Energy, momentum, area, distance, time, etc. are expressed in appropriate powers of electronvolts (eV).
The SU (3) C × SU (2) L × U (1) Y gauge fields of the Standard Model (SM) are denoted by G a μ , W k μ , B μ , where a = 1 . . . 8, k = 1 . . . 3. The corresponding gauge couplings are called g s , g L , g Y , and the corresponding field strengths are defined as G a I use the plus sign covariant derivative convention: D μ X = ∂ μ X + i G a μ T a X + i W k μ σ k 2 X + i B μ Y X X . The (sine of the) weak mixing angle is related to the electroweak coupling as sin θ W = g Y g 2 L +g 2 Y , and the electromagnetic coupling is . After electroweak symmetry breaking, the photon field is denoted by A μ , and its field strength by F μν . The massive electroweak vector bosons are denoted W ± μ and Z μ , and in this case I define V μν ≡ ∂ μ V ν − ∂ ν V μ , without any non-abelian piece. The vector boson eigenstates are related to the SU (2) L × U (1) Y gauge fields by I use the 2-component spinor formalism, following the conventions of Ref. [7]. A Dirac fermion is described by a pair anti-commuting fields f α ,f ċ α transforming respectively under the first and the second component of the SU (2) ⊗ SU (2) Lorentz algebra. The spinor index can be raised and lowered by the anti-symmetric tensor, f α = αβ f β , 12 If you're not familiar with this notation…that's very bad, you should learn this as soon as possible, it's an essential part of modern education of a particle physicist. But if you don't want to learn, you can always quickly translate to the 4-component Dirac fermion using the map For example,fσ μ ∂ μ f =Fγ μ ∂ μ P L F, f c σ μ ∂ μf c = Fγ μ ∂ μ P R F, f c f =F P L F,ff c =F P R F, where P L ,R = 1∓γ 5 2 are the Dirac chirality projectors. The 1σ uncertainty on theoretical or experimental quantities is often expressed either using the bracket notation, e.g. x = 1.234 (56) is the same as x = 1.234 ± 0.056. The former notation is especially useful when precision reaches many digits.
All abbreviations are defined the first time they are introduced, but in case you forget they are all collected in Appendix A.

EFT ladder
In the previous century the Holy Grail of theoretical particle physics was the Theory of Everything. Physicist would imagine that internal consistency of quantum theories of matter and gravity selects essentially a unique theory, with very few or none at all free parameters. The Theory of Everything would be valid at all energies, up to the Planck scale and beyond, and it would lead to the SM as its low-energy approximation, possibly with some intermediate supersymmetric or grand-unified theories emerging between the electroweak and the Planck scales. Alas, this top-down approach has not quite delivered, and the quest for the Theory of Everything is now largely abandoned. Around the turn of century the focus shifted to the less ambitious but more practical Theories of Something. These theories are meant to be valid only in a restricted energy range, and often the degrees of freedom they describe are emergent rather than fundamental. For these reasons they are commonly referred as effective fields theories, or EFTs in short.
The central idea behind EFT is that things may appear simpler when viewed from a distance. For countless physical systems complexity is dramatically reduced by focusing on the large-scale behavior. Take for example a system of many static electric charges confined to a region of space of size R, see Fig. 1. A near observer positioned at a distance L ∼ R must trace the position of each charge to accurately determine the electric field in her vicinity. However, for a Fig. 1 Simple illustration of the EFT idea using the example of multipole expansion of the potential produced by electric charges far observer at r R the details of the charge distribution are not essential. Instead, the electric field at large r can be described using the multipole moments of the charge distribution: the total charge, the dipole moment, the quadrupole moment, etc. The error of this approximation is controlled by the ratio (R/r ) n , where n is the number of multipoles taken into account. For large enough r only a first few multipoles need to be included to adequately describe the electric field, and this way the description of a possibly complex system with many degrees of freedom is reduced to a small number of discrete parameters.
This reduction of degrees of freedom with increasing distance is so pervasive in physics that is often taken for granted. Indeed, the far observer could be a humble engineer tuning his antenna, who would be shocked when told that his actions have anything to do with EFT. It is easy to evoke other familiar examples: a gas in equilibrium, where the enormous mess of gazillions of atoms bouncing against each other can be summarized by a small number of thermodynamic quantities like temperature, pressure, entropy; planets rotating around the star, all of which may be complicated objects, with a finite radius, not-exactly-spherical shapes, and non-trivial density profiles, but for the purpose of planetary dynamics they can be perfectly approximated as point particles, while tidal corrections due to their sizes are calculable in a quickly converging expansion; and then the whole universe at the scales larger than that of galaxy clusters is described by the simple Friedmann equations depending only on the density of matter, radiation, and dark energy. We could go on with similar classical examples for hours. Instead, we will now head straight to particle physics and quantum field theory (QFT).
In relativistic QFT, instead of distance scales L, it is more convenient to refer to energy scales E. The two are simply connected via the uncertainty relation, E ∼ 1/L (in natural unitsh = c = 1), thus long distance translates low energy, or infrared (IR) in our jargon, while short distance translates to high energy, or ultraviolet (UV). As understood long ago by Wilson, Weinberg, and other giants of the past century, changes in complexity of QFT as we move towards lower energies can be nicely formalized in the language of "integrating out the UV degrees of freedom". The concept is perhaps most succinctly summarized using the path integral formulation of QFT. Consider a QFT with low-energy degrees of freedom denoted collectively as φ, and with high-energy degrees of freedom denoted as H . Here, φ and H can refer to, respectively, light and heavy particles in the theory, or to low-and high-frequency modes of the same particle. Quite generally, the full UV theory of φ and H can be defined by the Lagrangian L UV from which the partition function Z UV is calculated: All correlation functions of φ's and H 's (and thus all Smatrix elements) can be obtained by differentiating Z UV with respect to the auxiliary sources J . At low energies, such that the H modes cannot be excited, we only need the correlators of φ to calculate observables, hence we can set J H = 0. From this IR perspective we can define (1. 2) The first equality is a definition of the EFT partition function Z EFT . It is adequate for our purpose since Z EFT , trivially, leads to the same correlation functions of φ as Z UV . However, this is a tad formal as in QFT we rarely know the full partition function. The second equality is far more useful for practitioners. It defines the EFT Lagrangian L EFT , which captures the effective interactions of the light degrees of freedom. 1 These interactions should reproduce the correlation functions and scattering amplitudes in the full theory. An important point is that L EFT can be determined algorithmically, order by order in perturbation theory, if L UV is known and weakly coupled. The effective Lagrangian is sufficient to calculate all low-energy observables involving φ without ever referring to H . This often leads to conceptual and calculational simplifications compared to working with the full UV theory.
In general, L EFT defined by Eq. (1.2) is a very complicated non-local object. 2 Indeed, even though it depends only on φ, it must somehow contain information about the scattering amplitudes in the full UV theory, including the effects of propagating virtual H degrees of freedom at all loop levels. However, a dramatic simplification occurs in the presence of scale separation, when the mass M (or the characteristic frequency) of H is parametrically larger than the relevant energy scale, E M. In such a case, L EFT can be approximated by a local Lagrangian, that is by polynomial of φ and its derivatives. Mathematically speaking, the full L EFT will contain non-local expressions such as e.g. ( + M 2 ) −1 , which for E M can be Taylor-expanded into local ones, ( + M 2 ) −1 ≈ M −2 − M −4 + · · · . This is where the power of EFT shows up: the possibly very complicated shortdistance physics mediated by H is summarized by a discrete set of local interactions of φ suppressed by increasing powers of M, in close analogy with the multipole expansion in electrodynamics. The expansion can be truncated at some fixed order n in 1/M, depending on the precision required.
The philosophy sketched above has been applied over and over again in many areas of particle physics. One important example is the SM below the electroweak scale. At energies E m W 80 GeV (giga-electronvolt, that is 10 9 eV), one can integrate out the W and Z bosons, together with the Higgs boson and the top quark. The resulting EFT, which I will call Weak Effective Field Theory (WEFT, also known as WET or LEFT in the literature), has the photon, the gluon octet, the 3 generations of SM leptons, and the 5 lightest flavors of SM quarks as the degrees of freedom, and is valid in the range 2 GeV E m W . In this case the UV theory is known and weakly coupled, therefore the effective Lagrangian defined by Eq. (1.2) can be calculated and all Wilson coefficients can be determined as functions of the SM parameters. Because the massive electroweak gauge bosons are absent, the WEFT Lagrangian is invariant only under the SU (3) C ×U (1) em subgroup of the SM SU (3) C × ×SU (2) W × U (1) Y gauge symmetry. The memory of the larger gauge symmetry survives only in the specific pattern of the WEFT interactions and corresponding Wilson coefficients. In particular, the weak interactions, which in the SM are mediated by the exchange of W and Z , emerge in WEFT as 4-fermion effective interactions between the quarks and leptons. Many of these effective operators violate the (approximate) flavor symmetry and thus can mediate transitions between different quark flavors. WEFT is thus the theory underlying the vast phenomenology of flavor physics, that is the studies of transitions between mesons and baryons made of different quarks. Another well-known example is the so-called Chiral Perturbation Theory (ChPT) describing the physics of light mesons (pions, kaons, eta) at energies below the ρ resonance mass, E m ρ 775 MeV (mega-electronvolt = 10 6 eV). In this case the UV theory is known, but because the SM SU (3) C interactions become strongly coupled at E ∼ 2 GeV, the effective Lagrangian cannot be calculated analytically. Nevertheless, the approximate chiral symmetry of Quantum Chromodynamics (QCD) with light quarks, which resurfaces in another form in ChPT, allows one to systematically construct L EFT as a derivative expansion in ∂/(4π F), where F ∼ 100 MeV is called the pion decay constant. A less familiar example is the so-called General Relativity EFT (GREFT), which is an EFT extension of the Einstein theory of general relativity. Here the quantum field encoding the gravitational degrees of freedom is the spacetime metric g μν , which describes, in the limit of the flat Minkowski background, a massless spin-2 particle called the graviton. The Lagrangian is invariant under general coordinate transformations, which is necessary to decouple the unphysical degrees of freedom in the metric. The lowest order term is the Einstein-Hilbert Lagrangian, L GREFT ⊃ 1 2 M 2 Pl R, corresponding in the classical limit to the usual Einsteinian general relativity. Higher order corrections are constructed from powers of the Riemann tensor R μναβ (more precisely, from its Weyl tensor part), with the EFT expansion organized in powers of R μν / , where may or may not be equal to the Planck scale. The validity regime of this EFT is 0 E min(M Pl , ). 3 Unlike in the two previous examples, we only have vague speculations about the UV completion of GREFT: it may be some form of string theory, or something completely different. For this reason, we do not know the coupling constants multiplying the higherderivative interactions terms in the GREFT Lagrangian; they have to be treated as free parameters to be determined one day from experiment.
Yet another important example is SMEFT, which is the main topic of these lectures. The SMEFT philosophy has been employed in high-energy physics since more than 40 years [8], but only quite recently, around the year 2010, the theory gained large prominence. SMEFT is an EFT of the SM degrees of freedom: the photon, the gluon octet, the W and Z bosons, the Higgs boson, and the 3 generations of quarks and leptons. Much as in the SM, the action is exactly invariant under the local (gauge) SU (3) × SU (2) × U (1) symmetry. The SMEFT Lagrangian contains the SM one, but also an infinite set of higher-dimensional gauge-invariant interaction Fig. 2 The ladder of EFTs describing nature, assuming SMEFT is a valid description of nature in the regime 100 GeV E , with 100 GeV terms, 4 The latter interactions, which are non-renormalizable in the old parlance, describe the effects of heavy particles from beyond the SM. Under very broad assumptions, which will be spelled out in Sect. 2, SMEFT is the theory of fundamental interactions in the energy range 100 GeV E , where m W is the scale at which non-SM particles appear. The Lagrangian is organized in a systematic expansion based on the canonical dimensions of the interaction terms, with the operators of canonical dimension D suppressed by D−4 . The operators with D = 5 and D = 6 are expected to provide the leading deformations of the SM Lagrangian. Most often, the expansion is truncated at D = 6, with the D > 6 operators deemed as irrelevant at the currently available energies.
Let us discuss the place of SMEFT in the larger scheme of things. The ladder of effective theories is sketched in Fig. 2. On top of the ladder, for E , we have a hypothetical theory that UV completes SMEFT (which itself may be an EFT of another, more fundamental theory). At the time of writing, we have no clue what it is, what its degrees of freedom are, or what its mass scale is. Nevertheless, given the lack of discovery of non-SM particles at the LHC, it is a reasonable assumption that 1 TeV (tera-electronvolt = 10 12 eV). The validity domain of SMEFT is 100 GeV E , and our assumption about implies in particular that SMEFT is 4 In the EFT jargon, these higher-dimensional interactions terms are often referred to as operators for no good reason. The coupling constants multiplying these operators are often referred to as the Wilson coefficients. I will use this jargon in the following. the relevant theory to describe processes at the LHC collider. Since the UV completion of SMEFT is unknown, much as for GREFT, the Wilson coefficients in the SMEFT Lagrangian should be treated as free parameters to be determined from experiment. Below the electroweak scale, for 1 GeV E 100 GeV, SMEFT reduces to WEFT mentioned earlier, but now due to more general assumptions regarding the UV completion the WEFT Wilson coefficients should be treated as free parameters. Moving to lower energies, something dramatic happens around 1 GeV. The number of degrees of freedom and complexity explodes due to the onset of strong QCD coupling and emergence of baryons and hadrons as bound states of quarks. When the smoke clears, for 100 MeV E m ρ , we are left with ChPT describing the lightest mesons coupled to electrons, muons, photons, and neutrinos. The local symmetry is reduced to U (1) em , while the color SU (3) C and the associate gluons are no longer relevant degrees of freedom at these energies. Below 100 MeV we have a series of EFTs with less and less degrees of freedom and complexity. First we have an EFT extension of QED coupled to neutrinos. Then electrons can be decoupled, and we are left with massless photons, and almost massless neutrinos, interacting with the former via highly suppressed dipole interactions. Finally, below the mass of the lightest neutrinos (which is also unknown, but I'm assuming here it is nonzero) we have the EFT of pure light. As far as we know, photons are exactly massless, therefore the theory preserves the U (1) em local symmetry, and its validity extends down to the scales of the order of the inverse size of the universe. In this ultimate EFT, photons still interact with each other, albeit very weakly, via dimension-8 and higher operators in the so-called Euler-Heisenberg (EH) Lagrangian.
Each rung of this ladder deserves a series of lectures on its own, but in the following I will focus almost exclusively on SMEFT.

Assumptions behind SMEFT
In theory, SMEFT is a perfectly consistent EFT of the SM degrees of freedom. However, it is not guaranteed that there is any energy range where SMEFT is the relevant EFT to describe physical processes. For this to happen, several broad assumptions have to be satisfied. Let me first list these assumption, and then we will discuss them in some detail. #1 QFT. Physics above the electroweak scale is described by a manifestly Poincaré-invariant local quantum theory. #2 Mass Gap. The mass scale of the non-SM particles is much larger than the electroweak scale, m W . #3 Gauge Symmetry. The Lagrangian describing interactions above the electroweak scale is invariant under the SM gauge symmetry SU In Assumption #1, Poincaré-invariant entails Lorentz and translational invariance in four spacetime dimensions. This is often taken for granted. QFT is a surprisingly rigid structure, and there are very few ways to modify it without wrecking some fundamental principles, such as unitarity and causality for example. Within QFT, one consistent departure from Assumption #1 would be to introduce extra space-time dimensions, but these are very unlikely to be relevant anywhere near the electroweak scale. Consistent non-QFT quantum frameworks are rare and far between; one such example is string theory, but, again, it is unlikely to be relevant anywhere near the electroweak scale. At the same time, there is no single hint from experiment that the standard QFT techniques may break down at the energy scales available in the foreseeable future. All in all, the first assumption seems a very safe one. Even though we do not know until which energy scale QFT can be used, we will assume it works at all the scales relevant for these lectures.
Assumption #2 is more tricky because, strictly speaking, it is false. Indeed, the degrees of freedom at the electroweak scale include not only the SM spectrum, but also a massless spin-2 particle called the graviton, which mediates the gravitational interactions. Thus, in order to describe all known physics at that scale we should also include the graviton in our EFT, which leads to the construction called GRSMEFT [9]. Nevertheless, gravity is expected to be very weak around the electroweak scale. Consistency of the theory requires the leading order coupling of matter to gravitons to be universal and controlled by the scale M Pl 10 18 GeV, leading to the suppression factor of TeV/M Pl ∼ 10 −15 at the LHC energies. Subleading graviton couplings are controlled by the GRSMEFT expansion scale , which is unknown, but the (rather safe) assumption here is that m W , perhaps even ∼ M Pl . If that is satisfied, graviton emission is totally irrelevant at the LHC and in other experiments that focus on non-gravitational interactions. For those experiments, SMEFT provides an adequate description. On the other hand, for observables where gravity plays a central role, for example for gravitational wave emission and detection, GREFT or GRSMEFT should be used.
Are there any other light non-SM degrees of freedom except for the graviton? This is an open question at present. Theorist have hypothesized countless light particles, some of which are even well motivated, and sometimes even hinted at by some experiments. As examples one could mention the sterile neutrinos, the axion, and a light dark matter particle. An affirmative answer to our question will be provided if we are very lucky and such a particle is discovered in some ongoing or future experiment. However a negative answer may never be established, because in many scenarios the coupling of the new particle to the SM matter is a free parameter that can be adjusted to arbitrary small values. From our point of view, a more immediate question is whether the non-SM degrees of freedom are relevant at the LHC energies. Again, this is an open question that may be difficult to settle in the near future. For all we know, a new light particle could for example couple to the Higgs boson, and could lead to an invisible Higgs branching fraction up to O(10)%. Using the SMEFT framework one misses such a possibility. All in all, it is reasonable to assume that the graviton is the only non-SM light degree of freedom, however it certainly requires a certain leap of faith. SMEFT practitioners should always keep their eyes and minds open and follow experimental developments in collider physics and elsewhere. In case the existence of a new light particle is established, the SMEFT approach may have to be abandoned.
We arrive at Assumption #3, which is the most mysterious one. In the SM, the action is exactly invariant under the SU (3) C × SU (2) W × U (1) Y local symmetry, which in the global limit acts as a linear transformation on the fields in the Lagrangian. At the level of the spectrum this symmetry is not visible, because it is spontaneously broken by a vacuum expectation value (VEV) of the Higgs field. With some experimental input about the quantum numbers of SM matter, the gauge principle has led to highly non-trivial and successful predictions. For example, the interactions strength of all left-handed fermions with the W boson are predicted to be universal (in the tree-level approximation) and controlled by the SU (2) L gauge coupling g L , while the interactions with the Z boson are predicted non-universal but controlled only by the fermion's quantum numbers and one universal parameter called the weak mixing angle sin θ W .
All in all, gauge symmetry has proved to be one of the deepest foundational ideas in QFT, and the SM gauge symmetry has time and again proved to be extremely successful phenomenologically. That's all very impressive, but why should SMEFT respect the same gauge symmetry as the SM? In the end, the goal of SMEFT is to provide a model independent description of heavy new physics beyond the SM. The discussion is further complicated by the fact that, in the modern view, gauge symmetry is not a real symmetry of the physical system, but merely a redundancy of ifs description. Why do we insist on imposing that particular redundancy on SMEFT?
First, let us recall what is the true purpose of gauge symmetry, or gauge redundancy [10]. The point is that a consistent, unitary QFT that is manifestly Lorentz invariant and contains massless spin-1 particles must be equipped with gauge redundancy, one generator for each massless spin-1 particle. Heuristically, this is because a spin-1 particle is described in QFT by a 4-component vector field A μ , μ = 0 . . . 3, or equivalently by the associated polarization wave function μ ( p). Since, an on-shell massless spin-1 particle has 2 degrees of freedom, corresponding to the two helicities, two of the four components must be somehow projected from μ ( p). One can be taken care of in a Lorentz invariant way by the transversality condition p μ μ ( p) = 0. It turns out that the only Lorentz invariant way to project out the other spurious degree of freedom is to identify the states described by the polarization wave functions μ ( p) and μ ( p) + p μ , that is by imposing gauge redundancy on the theory.
In the SMEFT we have two kinds of massless spin-1 particles: a photon and a gluon octet. Accordingly, we need 9 generators of local symmetry to have a consistent and manifestly Lorentz-invariant theory. An input from phenomenology is needed to identify that SU (3) C × U (1) em provides a correct description of these degrees of freedom, because the gluons all self-interact with each other, thus they are described by the non-abelian SU (3) factor, while the photons do not have self-interactions, thus they are described by the abelian U (1) factor. But this raises another question: why do we insist on the larger SU is enough to satisfy the consistency principles of QFT?
In fact, an EFT for the SM degrees of freedom, where only the SU (3) C × U (1) em gauge symmetry is realized linearly, does exist and is most often referred to as HEFT (as in Higgs EFT). In HEFT, the generators of the larger are realized as a non-linear transformation of the scalar Goldstone bosons eaten by W and Z , akin to the realization of the SU (2) L × SU (2) R /SU (2) V in ChPT. While the formal difference between HEFT and SMEFT is clear, the physical difference between the two EFTs is more subtle and was elucidated only recently [11,12]. The long story short: HEFT is an effective theory for nondecoupling UV physics, that is for theories where the masses of non-SM particles are dominated by contributions from electroweak symmetry breaking. A simple toy model for such a UV completion is a real scalar field S without a mass term but with the quartic interaction with the Higgs field: L ⊃ −λ|H | 2 S 2 . After electroweak symmetry breaking S acquires mass m 2 S = 2λ|H | 2 , which can be large if the quartic coupling λ is O(1) or larger. Integrating out S will lead to an EFT described by the HEFT framework rather than SMEFT. Another less artificial example is the SM with 4 generations of chiral fermions, in which case all fermions are massless in the limit of the Higgs VEV going to zero. Integrating out the 4th generation will again lead to HEFT rather than SMEFT. On the other hand, integrating out the 4th generation of vector-like fermions, where the masses of the non-SM fermions are dominated by a vector-like mass term M v, will lead to SMEFT rather than HEFT. In the end, the gauge symmetry Assumption #3 turns out to be closely related to the mass gap Assumption #2. Indeed, in non-decoupling theories masses of non-SM particles are of the form m i ∼ g i v, where g i is some gauge or Yukawa coupling. Since couplings are restricted by perturbativity to be |g i | 4π , the masses are m i 4πv. This means the new particles in non-decoupling theories are within the reach of the LHC or just around the corner. Conversely, if new physics enters at the scale 4πv ∼ 3 TeV, then the physics below is necessarily described by SMEFT and not HEFT. By imposing Assumption #3 we make an implicit decision to neglect the possibility of non-decoupling UV completions. Note that large swathes of non-decoupling theories have already been experimentally excluded; for example, the chiral 4th generation was definitely excluded by the Higgs production rate measurements at the LHC. Even though, at present, one cannot formally exclude the existence of nondecoupling new physics, and some wiggle room remains for certain constructions, it is a very unlikely possibility in my opinion. Focusing on decoupling new physics, and thus restricting our scope to SMEFT, seems a very reasonable assumption.
Note that assumptions #1-#3 do not restrict the SMEFT Lagrangian to be renormalizable. There was a time in the history of particle physics when renormalizability was hailed as a sacred principle that every successful quantum theory should obey. Now the pendulum has swung in the opposite direction, and we think that every QFT description of a reallife physical system corresponds to a non-renormalizable EFT. Now, in some case that EFT may be well approximated by a renormalizable QFT, as is the case for physics at the electroweak scale. We think of this as an accident due to a large separation between the electroweak scale and the scale suppressing the non-renormalizable interactions. However we expect that these non-renormalizable interactions are present in the Lagrangian, and that they will become apparent when enough experimental precision is achieved.

Constructing SMEFT
This section reviews a systematic prescription to construct the SMEFT Lagrangian. The fields corresponding to the SM particles and their representations under the gauge symmetry are summarized in Table 1. Using these fields as building blocks, we will write down the most general Lagrangian consistent with the assumptions spelled out in Sect. 2.

Power counting
Because the SMEFT Lagrangian is non-renormalizable, it contains an infinite number of interaction terms. Even if we wanted to arbitrarily restrict to a finite number of interactions, loop corrections would force us to introduce an infinite number of counterterms to cancel the UV divergences. In order to make the theory usable in practice we need power counting, which is the EFT jargon for an organizing principle that allows us to establish a relative importance of different interaction terms. In SMEFT, a natural power counting is based on the canonical dimension of an interaction. We organize the SMEFT Lagrangian as where each term L D in this series contains operators O i,D of canonical dimension D: where c i,D are dimensionless, and is a common mass scale entering all Wilson coefficients. At this point Eq. Then the dimensionless coefficients c i,D are functions of the couplings and mass ratios in the UV completion of SMEFT, as well as of the SM couplings. Now, the standard SMEFT power counting relies on the assumption that |c i,D | ∼ 1, that is to say which is basically dimensional analysis. In such a case we have a simple estimate of the relative relevance of different Wilson coefficients. Matching the dimensions in tree-level scattering amplitudes (which are dimensionless) one finds that, for the relevant scattering energy E much larger than the particles' mass, a Wilson coefficient at a given D will enter as For example, the effects of dimension-4 operators are unsuppressed, the effects of dimension-5 operators are suppressed by E/ , the effects of dimension-6 operators are suppressed by (E/ ) 2 , and so on. The higher the dimension of the operator, the larger is the suppression. Thus, operators with lower dimensions will have a larger impact on phenomenology, assuming E , that is when SMEFT is used at the energy scale well below the mass scale of the UV completion. We can thus truncate the SMEFT Lagrangian at some particular D, ignoring the contributions of all but a finite number of operators. Conversely, for E ∼ the suppression of higher-dimensional operators is no more, and one should take into account the whole infinite series of operators in the Lagrangian to correctly evaluate the amplitude. Obviously, in this regime SMEFT in unusable, and thus is the cutoff scale of SMEFT, beyond which it should be replaced by a more fundamental theory.
One important consequence of the standard power counting is that it allows one to define SMEFT at the quantum level. Recall that SMEFT is non-renormalizable, thus in principle an infinite number of unknown counterterms has to be introduced to properly define loop corrections to amplitudes of physical processes. However, working at E , we can declare that we drop from the amplitudes all the contributions that are O( 4−D max −1 ) or smaller. By dimensional analysis it is easy to see that the counterterms corresponding to operators of dimension D max + 1 are moot and we can neglect them in our analysis. This leaves a finite number of operators of dimension D ≤ D max , together with the associated counterterms. Thus, SMEFT with the standard power counting and truncated at a finite D max is as renormalizable as the renormalizable theories in the standard sense (D max = 4). From the SM it differs only by a larger number of counterterms (if D max > 4), thus a larger number of free parameters that have to be fixed by experiments. The standard power counting sketched above has the advantage of being simple and self-consistent. One should remember however that it is not the only option, and it may not be the most sound one from the physics point of view. A run-of-the-mill UV completion will not generate all Wilson coefficients universally; typically it will generate a handful of operators at tree level, while others will be suppressed by loop factors, leading to hierarchies not captured by Eq. (3.5). Moreover, certain types of operators can never be generated at tree level, independently of the UV completion. Next, flavor or other symmetries in the UV completion may lead to special patterns in SMEFT, leading to additional suppression of Wilson coefficients. For example, Eq. (3.5) suggests that Wilson coefficients corresponding to analogous operators involving say, up and top quarks scale in the same way, however if the UV completion incorporates something akin to SM flavor hierarchies (which is very likely) one expects the former will be suppressed compared to the latter by a small factor (m u /m t ) n . Finally, Eq. (3.5) ignores the dependence of the Wilson coefficients on the coupling strength in the UV theory. Consider a UV theory with a single coupling g * . Very often, Wilson coefficients of dimension-6 and -8 operators will scale C i,6 ∼ g 2 * 2 and C i,8 ∼ g 2 * 4 . In the standard power counting, C 2 i,6 is always of the same order as C i, 8 , which is indeed the case for g * ∼ 1. But for g * 1 we have C i, 8 C 2 i,6 , whereas for 1 g * 4π we have 6 , in both case the parametric hierarchy being missed in the standard power counting.
Nevertheless, let us brush aside these caveats for the time being and proceed under the assumption that the canonical dimension of an operator is the central determinant of its relevance for the low-energy phenomenology at E . Consequently, we will build the SMEFT Lagrangian starting from the operators of lowest dimensions, and working up towards higher D.

D ≤ 4
The sum in Eq. (3.1) starts at D = 2 because there is nothing at lower dimensions: D = 0 would be a field-independent constant, which has no physical consequences in non-gravitational theories, while there is no gauge invariant D = 1 operators because there are no singlet scalars in the spectrum in Table 1. At D = 2 there is a single gauge invariant operator, the Higgs mass squared: The Wilson coefficient in this case has mass dimension 2 and is denoted as μ 2 H . According to our power counting in Eq. (3.4), we should have μ H ∼ v. In reality we expect μ H v because the Higgs mass term triggers electroweak symmetry breaking by the Higgs VEV. In the SM, where there are no free unknown parameters anymore, we know precisely the tree level value μ H 88 GeV. In SMEFT I cannot give you a number for μ H because unknown higher dimensional operators also affect the Higgs VEV. Nevertheless, μ H v would be unnatural as it would require large cancelations between μ H and higher-dimensional operators to arrive at the correct value of v. We thus have a puzzle. On one hand, power counting predicts μ H ∼ v. On the other hand, phenomenological and naturalness arguments imply μ H v. This clash is nothing else but the hierarchy problem. 5 Not so long ago, the hierarchy problem was considered an almost certain indication that there are new degrees of freedom at the electroweak scale, for example the supersymmetric partners or the Kaluza-Klein modes of the SM particles. If that were the case, SMEFT would not be a useful theory in any energy range. However, the results from the LHC strongly suggest that the SM degrees of freedom are all there is near the electroweak scale, and that SMEFT is the correct description of physics, at least in the energy range from 100 GeV up to a few TeV. That's good for SMEFT and fortunate for my lectures, however the hierarchy problem remains puzzling. Have we somehow missed the degrees of freedom responsible for stabilizing the electroweak scale? Can the hierarchy problem be addressed with no new degrees of freedom at the electroweak scale? Do we misunderstand something about how QFT works? Is the SM more fundamental than we think? It is fair to say that no one has presented a convincing solution so far.
So, we start with high ideals: everything is EFT, physics is basically dimensional analysis, etc., but at the first opportunity reality slaps us in the face…Nevertheless let us press on and apply the standard power counting to SMEFT operators of dimensions higher than two. At D = 3 again there are no gauge invariant operators because there are no sin-glet fermions in the spectrum in Table 1. 6 At D = 4 there are multiple gauge-invariant operators. Here is the complete list 7 , and Y f are 3 × 3 matrices in the generation space. Dimensional analysis dictates that all the couplings in the dimension-4 Lagrangian: the gauge couplings g X , the Yukawa couplings Y f , and the quartic coupling λ, are dimensionless. The standard power counting in Eq. (3.4) treats them all as O( 0 ) couplings. In reality, this is reasonably well borne out for the gauge and quartic couplings, but not for most of the elements of Y f . Clearly Eq. (3.4) does not know about flavor hierarchies. Some of the D = 4 Wilson coefficients are extremely suppressed, e.g. [Y e ] 11 3×10 −6 (in a convenient basis). It is conceivable that contributions of some D > 4 operators to certain scattering amplitudes will be larger than the effects proportional to the electron Yukawa coupling, which would represent another break down of the standard power counting. But, overall, the standard power counting is a very successful principle at D = 4: all but the last term in Eq. (3.7) have been experimentally shown to exist (again assuming that hey are not somehow mimicked by higher-dimensional operators). Of course, L D=2 + L D=4 is nothing else than the SM Lagrangian, so the success of SMEFT with the standard power counting is to reproduce the SM as the leading terms in its EFT expansion. Concerning the last term in Eq. (3.7), the current constraints are |θ| 10 −12 . The lack of experimental evidence for the θ term, which is referred to as the strong CP problem, is as much puzzling from the EFT perspective as it is within the SM. Fortunately, unlike for the hierarchy problem, we have some reasonable ideas about the solution. The smallness ofθ likely means that there is a new particle called the QCD axion, which effec- 6 Dimension-3 operators are present for example in the extension of SMEFT featuring singlet (right-handed) neutrinos. 7 We only wrote the so-called θ-term G a μνG a μν for the SU (3) C gauge bosons because analogous terms for other group factors have no physical effect. The θ-term is moot for U (1) gauge groups, while for SU (2) W the term W μν : kW k μν can be redefined away via a chiral transformation.
tively makesθ a dynamical quantity settled in a minimum whereθ ∼ 0. But at this point one cannot completely exclude the possibility that we misunderstand something fundamental about QCD, and in reality theθ term has no physical effects. 8 Or that the parameterθ is very small by pure accident.

D = 5
We move to D = 5, that is beyond the SM. At this order in the SMEFT expansion we have the following gauge-invariant interactions [8]: The Wilson coefficients C 5 form a 3 × 3 matrix in the generation space. Here and in most of the following, the generation indices are implicitly contracted, so that one should read Incidentally, neutrinos are known to be massive particles. 9 While we do not know the absolute values of the masses, we know the mass differences (squared) with a good accuracy, see e.g. [15]. Given this, one can estimate C 5 v 2 ∼ 10 −1 eV, that is to say C 5 ∼ 1 10 15 GeV . One cannot emphasize enough what an enormous success of the SMEFT paradigm this is. In SMEFT, the most relevant phenomenological effects at E are expected from the D = 2 and D = 4 operators, which are those of the 8 Such claims appear on arXiv in regular intervals. However, these are at odds with the preliminary evidence from lattice calculations [13,14], which observe nucleon electric dipole moments (EDMs) proportional toθ at large pion masses. 9 There is no doubt that at least two neutrinos have masses, but their precise nature is experimentally an open question. There are two different mechanisms to implement the neutrino masses in the Lagrangian. The simplest option is to write down the so-called Majorana mass term for a left-handed neutrino ν: L = − 1 2 m M νν + h.c.. Another option is to add a new degree of freedom to the SM -the right-handed neutrino ν c -together with the Dirac mass term L = −m D ν c ν + h.c.. The two options lead to different contributions to the neutrinoless double beta decay. The jury is still out whether the SM neutrino masses are of the Majorana, or the Dirac, or the mixed type. In the following I will be assuming without any further comment that the masses are purely Majorana. SM, and which are indeed seen in nature. Furthermore, the standard power counting predicts that the most relevant deviations from the SM should be due to D = 5 operators. This prediction was spectacularly confirmed by the discovery of neutrino masses in the Super-Kamiokande detector in 1998 [16], almost 20 years after Weinberg's paper [8].
At the same time, this very success carries a premonition of doom. The neutrino masses turn out to be quite small, leading to the appearance of a very large scale in the denominator of C 5 . Since in the standard power counting C 5 ∼ −1 , it would be most natural to conclude that the SMEFT expansion parameter = 10 15 GeV. This would not be a problem for SMEFT -on the contrary, it would mean that the expansion is very quickly convergent, and thus the operators up to D = 5, maybe plus a handful of operators at D = 6 are enough to describe all physics at available energy scales. But this would be a problem for you and for me. It would mean that the gap between the electroweak scale and the new physics scale is enormous, which would make the options for fundamental research very limited. The directions worth pursuing would be neutrino physics, and perhaps proton decay. Otherwise one could switch to astrophysics, cosmology, quantum computing, nuclear fusion, climate science, or banking. Not much point for future colliders, flavor physics, charged lepton flavor violation, which would only serve to confirm ad nauseam the SM predictions.
This may be the future, but it does not have to be. Even within the SMEFT paradigm (no new light degrees of freedom), it is quite possible that the expansion parameter is much smaller than 10 15 GeV. New physics responsible for the operators in Eq. (3.8) may be much lighter, perhaps even near the TeV scale, but coupled very weakly to the SM fermions. A sharper argument can be formulated by noticing that the operators is Eq. (3.8) are very special, as they violate the lepton number symmetry acting as L → e iα L, E c → e −iα E c . This is an accidental symmetry at the D ≤ 4 level, as one simply cannot construct a gauge invariant operators with D ≤ 4 that violates it, and thus D = 5 is the lowest dimension where lepton-number-violating operators can appear. One can modify the standard power counting by assuming that there are two scales governing the SMEFT expansion. One, call it L , corresponds to the mass scale of B − L-violating new physics, and it happens to be very high, L ∼ 10 15 GeV. Another, let's keep calling it without a sub-index, corresponds to the mass scale of B−L-conserving new physics. It is then perfectly natural to have a huge gap between these two scales, L . Symmetry consideration forbid new physics at the scale to generate D = 5 operators, and the lowest dimension it can show up is D = 6. This assumption of the two-scale expansion gives us a rationale for exploring the SMEFT Lagrangian at D = 6 and higher, and we will tacitly adopt this point of view in all of the following.

D = 6
We have arrived at dimension-6 operators, which is the nexus of the SMEFT research. At D = 2 there is a single operator; the D = 4 Lagrangian can fit a t-shirt; at D = 5 there is basically a single operator but, taking into account the generation structure, it counts as 12 operators. 10 At D = 6, all hell breaks loose: we have…wait for it…3045 independent operators. They contribute to phenomenology in virtually all areas of particle physics, such as Higgs physics, electroweak precision observables, flavor physics, nuclear physics, electric dipole moments, and much more. Below I will present a quick survey of dimension-6 operators using the set proposed in Ref. [17] and known under the name of the Warsaw basis. To organize the presentation, let me divide them into several classes: . (3.10) The bosonic operators, as the name suggest, are constructed out of the SM gauge and Higgs fields, without involving any fermionic fields. In the Warsaw basis there are 15 bosonic operators: where ≡ ∂ μ ∂ μ and σ k are the three Pauli matrices. Already this relatively small subset of dimension-6 operators contains rich phenomenology. C H changes the shape of the Higgs potential, in particular it affects the cubic Higgs boson self-coupling -perhaps the last landmark measurement to be delivered by the LHC. C H contributes to the Higgs boson kinetic term and thus, indirectly, affects universally all Higgs boson production and decay rates. The following two operators contribute to electroweak precision observables measured long ago by the LEP collider. C H D contributes to the Z boson mass, while C H W B contributes to the kinetic mixing between the photon and the Z boson.
Through these intermediaries, they affect the whole lot of electroweak precision observables. In fact, these two are just the famous oblique S and T parameters of Peskin and Takeuchi [18] in another (more modern) guise. Furthermore, C H W B as well as the Wilson coefficients C H G , C H W , C H B in the second line contribute to the ever important Higgs boson interaction strengths with gluons, W , Z , and photons, which are measured at the LHC. In the third line, C W and C G induce 3-derivative anomalous cubic interactions of electroweak gauge bosons and gluons, respectively. The final two lines contain CP violating interactions. They can be searched for in colliders, but more easily discernible effects appear via their loop contributions to electric dipole moments of the electron or the neutron.
The next class of dimension-6 operators we discuss are Yukawa-like interactions: Each C f H is a 3 × 3 complex matrix in the generation space, thus each comes with 18 free parameters, which makes 54 parameters overall. These operators contribute to the fermion masses, but that is unobservable because it merely renormalizes the unknown Yukawa matrices in Eq. (3.7). The observables effect is the modification of the Higgs boson Yukawa couplings to the fermions. In the SM, the Yukawa coupling is not a free parameter but it is uniquely fixed by the fermion's mass. In the presence of the operator is Eq. (3.12) that relation no longer holds, and the Higgs boson couplings to fermions become free parameters independent of fermion masses. Moreover, a qualitatively new effect of flavor violation in Higgs interactions may appear. That is to say, the Higgs boson can couple to two fermions from different generations, e.g. L SMEFT ⊃ hēμ c , which does not occur in the SM. Next, we have what I call the current operators: where The Wilson coefficient C H f are matrices in the generation space, but now only C Hud is a general complex matrix, while the remaining ones are Hermitian matrices (thus with 9 free parameters each). This adds up to 81 free parameters in Eq. (3.13). These operators contribute to the W and Z bosons interactions with fermions, which have been precisely measured in the LEP, Tevatron, and LHC colliders. Several qualitatively new effects are introduced by Eq. (3.13). One is the W boson couplings to right-handed quarks, e.g. L SMEFT ⊃ W μ (t c σ μbc ), whereas in the SM W couples only to left-handed quarks. Another is tree-level flavor-changing neutral currents, that is Z boson couplings to quarks or leptons of different generations, e.g. L SMEFT ⊃ Z μ (bσ μ s).
Next, we have the dipole operators (3.14) Given that C f V are 3 × 3 complex matrices in the generation space, the above introduces 144 free parameters. An important effect of the operators in Eq. (3.13) is their contribution to the anomalous magnetic dipole moments of fundamental particles. In particular, the Wilson coefficients [C eW ] 22 and [C eB ] 22 contribute to the muon g − 2 which, at the time of writing, may or may not deviate from the SM prediction. The imaginary parts of these Wilson coefficients contribute to electric dipole moments. Moreover, the operators in Eq. (3.13) can mediate certain processes that are forbidden in the SM, e.g. the μ → eγ decay. The dimension-6 operators introduced so far come with 15 + 54 + 81 + 144 = 294 free parameters. It follows that a large majority of dimension-6 operators are hiding in the last term in Eq. (3.10), which contains 4-fermion operators. For the sake of this discussion let me split them further into four sub-classes: defined by the number of lepton and of quark fields. The first sub-class in Eq. (3.15) is the 4-lepton operators: Thus, for example, Eq. (3.16) contains Therefore the components [C ll ] J K K J and [C ee ] J K K J can be declared real, as their imaginary parts do not enter the Lagrangian. 11 Four-lepton operators containing electron fields are relevant for physics at LEP-2, where e + e − pair were collided with the center-of-mass energy above the Z pole. A subset of interactions in Eq. (3.16) mediate tree-level charge-lepton-flavor violating processed, where the overall lepton number is conserved, but the separate electron, muon, or tau numbers are not. Such processes are forbidden in the SM, while they are mediated at loop level via the D = 5 intermediaries but with very suppressed rates due to the smallness of the neutrino masses. For example, [C ee ] 1112 mediates the μ − → e − e − e + decay, which is subject to current experimental searches [20]. Finally, it is worth mentioning that [C ll ] 1221 contributes to the usual muon decay μ − → e −ν e ν μ , which in the SM is a standard candle to determine the Fermi constant. By disrupting this standard candle, [C ll ] 1221 indirectly affects SM predictions for countless precision measurements.
The next sub-class in Eq. (3.15) are semi-leptonic operators, that is 4-fermion operators containing two quark and two lepton fields: 11 Another convention existing in the literature, see e.g. [19], is to set [C ll ] J K K J and [C ee ] J K K J to zero for J > K . This leads to a factor of two difference in the dependence of observables on these Wilson coefficients with J = K , as compared to the convention used in these lectures.
The semi-leptonic operators affect myriads of important precision observables: hadronic cross sections at LEP-2, Drell-Yan production of leptons in hadron colliders, electric dipole moments, beta decays, and so on. They also play a major role in flavor physics, where they contribute to semileptonic flavor transitions. These are often under reasonable theoretical control, such that reliable SM predictions can be established, and thus stringent constraints on the dimension-6 operators can de derived. Moreover, the importance of the operators contributing to flavor-changing neutral currents is amplified by the suppression of these processes in the SM. One of many relevant examples of this kind is the B s → μ + μ − decay, whose branching fraction is currently measured with 10% precision, and the SM prediction is know with a similar accuracy. That decay rate is affected, among others, by the Wilson coefficient [C (1) lq ] 2232 . The third sub-class singled out in Eq. (3.15) are four-quark operators: These play arguably a lesser role in phenomenology. The reason is that their effects have to compete with QCD processes, which are typically abundant and poorly controlled theoretically, especially at hadron colliders. Nevertheless, some of the operators in Eq. (3.19) will appear later in our story in the context of precision observables. The final sub-class in Eq. (3.15) is perhaps the most exciting one, as it consists of operators violating the baryon and lepton numbers: Above, the quark color indices are implicitly contracted by epsilon tensors, e.g. qqq ≡ abc q a q b q c . These operators violate the baryon number B defined as the global  Fig. 3.

Beyond
The sum in Eq. (3.1) extends to D = ∞, and one could press on, but, more often than not, the discussion in SMEFT stops at dimension six. What lies beyond? At D = 7 we have 1542 independent operators, at D = 8 the number is 44807, at D = 9 it grows to 90456, and at D = 10 we have a whopping 2092441 operators (two million!) [21]. There is a good chance, however, that in your research you won't ever deal with this cornucopia of higher-dimension operators. The exponential growth of the number of operators with increasing D is one reason, as it quickly makes any systematic analysis difficult. But that practical difficulty is not the only reason. The philosophy of SMEFT with the standard power counting is that, the higher the dimension of the operator, the more suppressed its effect is. Since at present we do not have any observational evidence of dimension-6 operators, it is hard to believe that D > 6 operators might show up in any experiment in a foreseeable future. This general conclusion should hold for other reasonable power counting beyond the standard one. There are a few exceptions, however, that one should be aware of. First of all, if new physics is close to the electroweak scale, the effect of higher-dimensional operators may be non-negligible, especially for observables probing the high-energy tail of differential distributions at the LHC, such as for example the Drell-Yan production of leptons, pp → + − . But that is of course also the situation where the SMEFT itself is least useful, since the central assumption of the mass gap is not quite satisfied. A more relevant case is when a qualitatively new phenomenon, which cannot be induced by D ≤ 6 operators, appears at D > 6. For example, tree-level contributions to light-by-light scattering enter at D = 8 from operators such as (B μν B μν ) 2 ; another D = 8 operator (l 1 Hσ μH l 1 )(u c σ μd c ) might play an important role in C P violation in nuclear beta decay [22]; neutronantineutron oscillations arise at D = 9 from operators such as (ū cdc ) 3 [23]. When analyzing this kind of observables, one should however pay attention whether new physics generating these higher-dimensional operators might not show up in other observables via D ≤ 6 operators.
All in all, while every particle physicist should have at least a rudimentary knowledge of the dimension-6 SMEFT operators, excursions beyond that are necessary and justified only in rare circumstances. You should be aware of the existence of higher-dimensional operators, and you should be able to use power counting to estimate their effects relative to the ones coming from lower-dimensional operators. But it is quite unlikely that the higher-dimensional operators can play a useful role in discovering or constraining new physics. The situation will of course change if dimension-6 operators are observed in nature. In such a case, going to higher dimensions will sharpen our low-energy description of new physics, and may help to pinpoint the scale of new physics [24].

Summary of the SMEFT Lagrangian
To summarize this section, SMEFT is a QFT with the spectrum given in Table 1 and the Lagrangian of the form where the consecutive terms are given in Eqs. (3.6) to (3.8) and (3.10) and the dots stand for higher dimensional interactions. The interactions in L D=5 and L D=6 describe various effects of the virtual exchange of exotic particles not included in Table 1, and the scale inherent in their Wilson coefficients is related to the mass scale M of these exotic particles. Almost any BSM model you will encounter will induce some D = 5 and D = 6 interactions. 13 The Wilson coefficients of the D ≥ 7 interactions are suppressed by higher powers 13 One exception is when the BSM degrees of freedom are massive particles with spin two and higher. This would be the situation in the Randall-Sundrum model [25] with all SM matter confined to one of the branes.  Table 1 emerge. At the energy scale E ∼ m Z it is convenient to integrate out the electroweak bosons, top quark, and Higgs, to derive another effective theory called WEFT with a smaller spectrum and a smaller gauge symmetry.
Since we do not know what is the UV completion of SMEFT, the usual philosophy is to include all possible independent interactions in Eq. (3.21) with completely arbitrary Wilson coefficients. This is a sane approach, up to small caveats. First, we should pay attention that the validity range is not null. This could happen if the Wilson coefficients are too large, as then the SMEFT amplitudes may violate perturbative unitarity already at E ∼ v. Second, the existing experimental results already put important bounds on some combinations of Wilson coefficients; more about it in Sect. 5. Third, several interaction terms in Eq. (3.21), such as e.g. klm W k μν W l νρ W m ρμ orL Hσ μνĒ c B μν arise from integrating out new particles at one loop level. Therefore we should not expect their Wilson coefficients to be very small, e.g. C W ∼ 1 v 2 would imply BSM particles below the electroweak scale, in contradiction with the assumptions behind SMEFT. Finally, there exist more subtle constraints on EFT Wilson coefficient arising from the mere assumption that the UV completion is local, unitary, and Poincaré invariant [26]. These typically imply a definite sign of a Wilson coefficient or a combination thereof, therefore they are referred to as the positivity constraints. With completely general assumptions about the high-energy behavior of the amplitudes in the UV completion, the positivity bounds are relevant for dimension-8 operators and higher. Therefore their application to SMEFT is limited to very specific cases, for example to the study of quartic gauge boson interactions [27,28]. In order to apply the positivity constraints to dimension-6 operators, additional assumptions about the UV completion are necessary [29].

Bases of SMEFT
A basis is a complete, non-redundant set of operators in the Lagrangian at a given order in the EFT expansion. Complete means that the full parameter space of the EFT is accounted for at that order. Non-redundant means that none of the operators in that set can be expressed as a linear combination of the remaining ones. In the previous section we introduced the dimension-6 SMEFT operators using one particular representation called the Warsaw basis. In this section we take a more general view, and discuss the technical side of establishing equivalence classes of higher-dimensional operator. The issue is non-trivial because there is a huge redundancy at the core of the Lagrangian formulation of QFT. We are so used to working with Lagrangians that we unconsciously assign a physical meaning to its interactions terms. However Lagrangians are merely tools -very useful, but also very crude ones. Seemingly distinct sets of interaction terms can lead to the same on-shell scattering amplitudes, and thus to exactly the same observable consequences.
One can distinguish three main ways we can reshuffle the interaction terms in the Lagrangian: 1. Integration by parts, 2. Fierz transformations, 3. Field redefinitions.
These operations do not change on-shell scattering amplitudes, which define the physical content of the theory, 14 Below I discuss them in more detail.
Starting with integration by parts, a simple example is This transformation, which consist in reshuffling how the derivatives act on the fields, changes the apparent form of the interaction term. However, both interaction terms above lead to the same Feynman rules once momentum conservation is taken into account. At a more mathematical level, integration by parts does not change physics because it shifts the Lagrangian by a total derivative ∂ μ F, which shifts the action S = d 4 xL by a boundary term, which does not affect scattering amplitudes. 15 Moving to the second point, Fierz transformations are rearrangements of the way the spinor or gauge indices are contracted. Let us illustrate them in a simple example, which is relevant for constructing a set of independent dimension-6 operators. Naively, it seems one could supplement the fourlepton operators in Eq. (3.16) by another one: where the σ k matrices act on the SU (2) W indices of the lepton fields (unlikeσ μ , which act on the Lorentz group spinor indices). This is a legitimate interaction term, in particular it is gauge invariant. However, one can show that it is redundant, in the sense it is already contained in the term proportional to C ll in Eq. (3.16). In order to show this, we first use the Fierz identity for the Pauli matrices, This allows us to rewrite I only display the SU (2) W indices in the first term in the bracket, where they contracted in an unobvious way. The second term in the bracket has already the same form as the L 4 term in Eq. (3.16), but we still need to work on the first term. To this end, we need another Fierz identity to rearrange the spinor indices. For anti-commuting spinors ψ i one can . Using lence classes, the so-called Green's bases where operators related by field redefinitions are not considered equivalent. See Refs. [30,31] for a construction of a Green's basis for dimension six and eight operators. 15 The D = 4 SMEFT Lagrangian in Eq. (3.7) contains a total derivative -the famous θ term -which does not affect perturbative amplitudes but does affect physics via non-perturbative instanton contributions to the path integral. To my knowledge, this is the only total derivative in the SMEFT Lagrangian that has physical consequences.
this we rewrite This completes the proof that the operator in Eq. (4.2) is equivalent to the L 4 term in Eq. (3.16), and we also found the mapping between the two tensor Wilson coefficients: Finally, we discuss the last point, which is the least trivial one. In QFT, the so-called equivalence theorem [32,33] stipulates that physics is invariant under arbitrary non-linear redefinitions of bosonic and fermionic fields of the form where F (i) X ( ) are function of all fields in the theory with the same Lorentz and gauge transformation properties as X . The expansion in some scale is to ensure that one can apply this transformation perturbatively, such that kinetic terms in the Lagrangian do not change, but only the interaction terms do. In SMEFT, can be naturally identified with the cutoff scale controlling the dimensional expansion of the Lagrangian. If we once again evoke the path integral formulation, invariance under field redefinitions seems obvious, as Eq. (4.5) amounts to a change of variables under the integral. In the context of effective theories the invariance of on-shell Green's functions under field redefinitions, also at the loop level, was formally proven in Ref. [34]. However, at the Lagrangian level this invariance seems far from obvious and often counter-intuitive. Indeed, field redefinitions can produce new interaction terms that were not present in the original Lagrangian. Let as consider an extreme example of a free scalar theory, and make a field redefinition φ → φ + 1 φ 2 . The shifted Lagrangian reads This naively looks like a Lagrangian of an interacting theory, with derivative cubic and quartic scalar self-interactions. However we know it is a free theory in disguise, as it was constructed via a field redefinition from the free Lagrangian in Eq. (4.6). And indeed, if you derive Feynman rules from Eq. (4.7) and calculate scattering amplitudes, you will always obtain zero, whatever the number of loops or external legs.
It would be however difficult to search whether there exists a field redefinition to a free theory, each time somebody hands you a new Lagrangian, Fortunately, there is a simple hack to quickly see that Eq. (4.7) is a free Lagrangian in disguise.
Using integration by parts we can rewrite the cubic interaction term as . This exercise demonstrates the truth that, at least at the leading order in the 1/ expansion, applying field redefinitions is equivalent to using equations of motions at a lower order in 1/ . In fact, in the literature this kind of equivalent Lagrangians are commonly referred to as related via equations of motions, rather than via field redefinitions. The two formulations are largely equivalent 16 but the former may be puzzling for some non-practitioners: why should quantum fields be constrained to obey the equations of motions, also in the context of loop calculations when they are off-shell? The formulation via field redefinitions appears much more intuitive. All in all, QFT Lagrangians have this huge redundancy due to field redefinitions. The redundancy exists also in the SM but it is rarely discussed in this context, as there exists a canonical choice of field variables where the SM Lagrangian appears renormalizable, that is without higher-dimensional operators. However in SMEFT, as in any garden-variety EFT, the redundancy is very relevant, because there is no canonical choice of field variables. As a consequence, Lagrangians with apparently distinct interaction terms may belong to the same equivalence classes, that is to say, they lead to the same scattering amplitudes. The equivalence classes can be confusing at times; for example a purely bosonic operator can be equivalent to a linear combination of operators containing fermions! Let us see in a concrete example how field redefinitions relate different dimension-6 SMEFT operators. The bosonic operator in the Warsaw basis, cf. Eq. (3.11), contains two quartic two-derivative self-interaction terms: (H  † H ) (H  † H ), and |H † D μ H | 2 . There is another interaction term of this type: 16 See Ref. [35] for the discussion of the differences between the use of field redefinitions and equations of motion beyond the D = 6 level of SMEFT.
which is perfectly legal and gauge invariant. However, using integration by parts and field redefinitions one can show its effects are already accounted for by other dimension-6 operators in the Warsaw basis. First we rewrite the operator multiplying C H D as In the first two lines we integrated by parts in two different fashions, and in the third line we summed the first two lines with weight 1/2 each. We also used a generalized version of the Leibnitz rule: The first term above can be absorbed into the C H Wilson coefficient in Eq. (3.11), but the second term requires more work. To get rid of it we can redefine the Higgs field as 17 where I put the prime on the right-hand side to mark that it is a redefined field, but in the following I will omit it to reduce clutter. Then the last term in Eq. (4.10) is canceled by the shift of the Higgs kinetic term However, the field redefinition in Eq. (4.11) also shifts other D = 2 and D = 4 terms (it also shift D = 5 and higher 17 Equivalently, we can replace D 2 H in Eq. (4.10) using the Higgs operators, but that does not affect dimension-6 operators, which we are focused on here). One finds Fortunately, all the terms that appear above are already in the Warsaw basis. At the end of the day, the effects of the dimension-6 operator in Eq. (4.9) can be simulated by concerted efforts of several dimension-4 and -6 operators in the Warsaw basis with the Wilson coefficients (4.14) Note that multiple bosonic and fermionic operators are needed to simulate a single bosonic operator. This may be counter-intuitive, but again can be confirmed by an explicit calculation of scattering amplitudes. By the sort of exercises as the ones above, one can show that every dimension-6 operator can be expressed by the ones already present in the Warsaw basis. A proof seems like a non-trivial task, and indeed it took a long time to complete. While dimension-6 operators have been frequently used in the literature since the influential work of Buchmüller and Wyler in 1986 [36], a complete basis of dimension-6 operators appeared only in 2010 [17]. A more systematic road to constructing a basis is offered by the Hilbert series techniques [21,37].
The Warsaw basis is by no means a unique representation of the parameter space of dimension-6 SMEFT operators. For example, one could drop from the Warsaw basis the operator proportional to C ll in Eq. (3.16) and replace it by the one proportional to C ll defined in Eq. (4.2). Or one could drop the operator proportional to C H in Eq. (3.11) and replace it with the one proportional to C H D defined in Eq. (4.9). Or keep an arbitrary linear combination of the two. And so on. This way, there are infinite ways to write a basis of dimension-6 SMEFT operators. Any such basis will have 3045 operators, and will lead to completely equivalent predictions of physical observables.
The Warsaw basis is widely used in the literature, and it is nowadays an indispensable reference point for all phenomenological studies of dimension-6 operators. Nevertheless, different bases may be more convenient for specific applications. Another popular choice of operators is the socalled SILH basis [38]. 18 In this basis, the bosonic operators are  16), and the four-quark operator multiplying [C uu ] 3333 , while keeping all other two-and four-fermions operators the same as in the Warsaw basis. From this discussion you can gather that the SILH basis may be a bit awkward for the sake of describing new physics with a non-trivial flavor structure. On the other hand, it is convenient to describe the so-called oblique new physics, that is the one where the leading non-SM effects enter via corrections to gauge boson propagators. Indeed, the SILH basis connects in a simple manner to the oblique S, T , W , and Y parameters defined in Refs. [18,41]: the S parameter corresponds to a linear combination of the S W and S B Wilson coefficients, while T , W , and Y translated directly into S T , S 2W , and S 2B , respectively.
A somewhat different approach to parametrizing the space of dimension-6 SMEFT operators was proposed in Ref. [42]. The connection between the Wilson coefficients in the Warsaw or SILH basis and observables can be quite complicated (see later in Sect. 5). For phenomenological applications it is often easier to work with the SMEFT Lagrangian expressed in terms of the fields that are mass eigenstates after electroweak symmetry breaking: the photon, the W and Z bosons, the Higgs boson. The idea is to pick a set of independent couplings in the Lagrangian -the so-called BSM primaries, which, at least at tree level, are simply related to experimentally observable quantities. Typically, such a set would include the Higgs boson couplings to matter (related to the Higgs signal strength measured at the LHC) and the W an Z boson couplings to fermions (related to cross sections, branching fractions, and forward-backward asymmetries precisely measured at LEP). With a judicious choice of this set (see for example the Higgs basis [43]), there exists an invertible map between the BSM primaries and the Wilson coefficients of the Warsaw or SILH basis, Then one can translate the theoretical predictions from the Wilson coefficient language to the BSM primaries language, and vice-versa. While the two are completely equivalent, the latter can be more convenient for certain applications. In global SMEFT fits to electroweak precision observables, the BSM primaries allow one to easily separate strongly constrained directions in the parameter space (typically, coming from the LEP measurements at the Z pole) from the weakly constrained direction (for example, coming from the Higgs measurements at the LHC), which improves stability of the fit.
Bases of SMEFT operators are also known beyond D = 6. For D = 7 operators (which are B − L-violating, much as those at D = 5 and at any other odd dimension) a basis was constructed in Ref. [44]. More recently, bases at D = 8, 9, 10, 11, 12 were constructed in Refs. [45][46][47]. A code to generate a SMEFT basis at any dimension is now available [48]. As I discussed in the previous section, operators of dimension higher than six currently have limited phenomenological applications, and no one even dreams of using a complete set of independent operators beyond D = 6. Nevertheless, these constructions serve as proof-of-principle demonstrations and useful look-up tables for specific applications.
Let me close this section by a couple of loosely related comments. The first concerns the choice of basis in the flavor space. Even when a specific operator basis is chosen, for example the Warsaw basis as in Sect. 3, there still remains freedom of making U (3) rotations separately for each fermionic 3-vector in the generation space: Q, U c , D c , L, E c . The power of these U (3) 5 rotations can be used for good. The most convenient choice is to make the fermion mass terms as diagonal as possible. This last clause is due to the fact that we cannot simultaneously make the up and down quark masses diagonal by these rotations without destroying the manifest SU (2) L symmetry of the SMEFT Lagrangian, that is without rotating differently the left-handed U and D fields sitting in Q. One convenient choice, which is adopted in these lectures, is the so-called down-type basis. In this basis, the fields U c , D c , L, E c are rotated such that their components are mass eigenstates, for example U c = (u c , c c , t c ) contains the right-handed up, charm, and top quark fields that do not have any flavor-off-diagonal mass terms. Concerning the quark doublet field Q, in the down-type basis it takes the form where U = (u, c, t) and D = (d, s, b) contain the lefthanded quark mass eigenstates, and V is the unitary 3 × 3 CKM matrix. The latter must appear here so that, after electroweak symmetry breaking, the kinetic term i Qσ μ D μ Q yields the usual quark-mixing interactions with the W boson: Another comment concerns on-shell amplitude methods. As discussed earlier in this section, the freedom of field reparametrization adds a lot of complications to constructing an EFT basis. It turns out that this excess baggage can be dropped: one can practice QFT without referring to fields, Lagrangians, gauge invariance, or other off-shell concepts. Instead, the theory can be formulated at the level of scattering amplitudes with all involved particles being on the mass shell (that is to say, with their respective momenta p constrained by p 2 = m 2 ). See Ref. [49] for a nice introduction. Unsurprisingly, the on-shell approach greatly simplifies the task of constructing an EFT basis, see e.g. [50,51]. It has been applied specifically to SMEFT [52], allowing one to swap a basis of higher-dimensional operators for a more tractable basis of independent amplitudes. This has lead to a more transparent derivation of some known results, as well as to several novel results. The advantages appear particularly clear in relation to renormalization group running, see e.g. [53][54][55], facilitating extension of the program beyond one loop.

From operators to observables
Currently, the SM is the reference point for any physical theory at the electroweak scale. It enjoys huge success, correctly accounting for a host of phenomena measured in colliders and low-energy precision experiments. SMEFT differs from the SM by the presence of interactions originating from operators with dimension D = 5, 6, and higher. In this section we discuss, in more precise and quantitative terms, the observable effects of these operators.
There are many ways in which higher-dimensional operators can affect observables. To organize the following discussion, it is convenient to divide them into three broad classes:

New vertices: interaction vertices in the SMEFT
Lagrangian that do not occur in the SM Lagrangian, due to symmetries or accidental reasons. In this following I will discuss each of these classes in turn.

New vertices
The most spectacular effects of SMEFT occur when higherdimensional operators violate an exact global symmetry of the SM. One important example of this kind is baryon and lepton number violation 19 by the dimension-6 operators in Eq. (3.20). Let us take one of these operators at random, say where I abbreviated C duu ≡ [C duu ] 1111 and explicitly displayed both the operator and its Hermitian conjugate. The latter mediates the quark-level process uu →de + . In a bit hand-waving but intuitive way, one can think of this process as transforming two up quarks from the proton (uud) into a down antiquark, leading to a dd meson state. Consequently, the operator can mediate proton decay into meson states. One important example is p → e + π 0 -the decay to a positron and a neutral pion (who is a combination of dd and uū quark states, the latter pair you can think of as being pulled from the vacuum sea during the decay). To calculate the rate for this process, one needs to take the on-shell matrix elements of 19 Strictly speaking, baryon or lepton number are an exact symmetry only at the perturbative level in the SM, but they are both violated by non-perturbative effects. Only one linear combination of the two, B − L, is conserved at the non-perturbative level. This subtlety is however irrelevant for the discussion in this section. the operator in Eq. (5.1) between the initial and final states: where T is the non-trivial part of the S matrix, S = 1 + iT , p 1 is the incoming momentum of the proton, k 2 , k 3 are the outgoing momenta of the pion and positron, and we also define q = p 1 −k 2 . The electron field acts on the annihilation operator of the positron final state, leaving the spinor wave functionx 3 corresponding to the momentum k 3 . We get The remaining matrix element between the proton and pion states is non-perturbative, and we cannot calculate it using the familiar textbook methods. Nevertheless, the Poincaré and little group covariance of the S-matrix tell us that it has to be proportional to a linear combination of the spinor wave functions of the incoming proton (pion is a scalar particle): We can thus parametrize where W 0 and W 1 encode the information about the nonperturbative brown muck. To know its value you have to ask your lattice friends, and they may reply W 0 ≈ 0.15 GeV 2 , W 1 ≈ −0.13 GeV 2 , with roughly a 20% error [56]. Thus or trading the matrix element for the more familiar amplitude via T = (2π) 4 δ 4 ( i p i − j k j )M: The rest is standard QFT manipulations. Taking the square of the amplitude, summing/averaging over the positron/proton spins, and plugging the result into the formula for the decay width one gets (5.7) Above I approximated the electron mass by zero, which is perfectly legitimate given the uncertainty on W 0 (actually, even approximating the much larger m π 0 as zero would be fine). In this limit the contribution of the W 1 form factor drops out. I presented this calculation here so that you can get familiar with hadronic matrix element, but as long as you are not interested in precision calculation you could easily obtain the order of magnitude of the result via dimensional analysis: ( p → e + π 0 ) ∼ m 5 p |C duu | 2 /16π . Now, a sneak peak into the Particle Data Group booklet [57] reveals that the limit on this proton decay channel is ( p → e + π 0 ) ≤ 1.3 × 10 −66 GeV at 90% confidence level (CL), which translates into the limit on the Wilson coefficient (5.8)

This limit is valid assuming only a single baryon-number violating operator is present in the Lagrangian; otherwise proton decay constrains a linear combination of various baryonnumber violating Wilson coefficients.
It is mind-blowing that low-energy experiments searching for proton decay allow us to probe new physics up to scales not so far from the Planck scale! To my knowledge, this is the highest scale we can indirectly access via low-energy experiments within a sane theoretical framework (unitary, causal, local, Lorentz-symmetric). The reason for this extreme sensitivity is that it is feasible to amass astronomical number of protons for a long period of time in a controlled setting, for example in a tank filled with water and surrounded by photo-detectors [58]. Moreover, the detection capabilities are impressive and would allow us to see the signal even if a handful of the protons in the tank decayed. Thanks to this combination of favorable circumstances, the limit on the proton lifetime can be orders of magnitude larger than the age of the universe! The final fact making the limit so strong is that we search for an effect that is predicted to be zero in the SM, so we do not have to face the uncertainty due to imprecise theory predictions.
Baryon and lepton number violation is certainly the most spectacular prediction of SMEFT. Nevertheless, higherdimensional operators can also break other exact or approximate global symmetries of the SM, with quite interesting consequences. A nice example is the decay μ → eγ . In the SM, not only the overall lepton number L is conserved, but also the individual lepton numbers L α for each generation, α = e, μ, τ . The process μ → eγ preserves L, but breaks L μ and L e , that is to say, these quantum numbers are dif-ferent for the initial and final states. Therefore in the SM μ → eγ is forbidden, and the predicted branching ratio is exactly zero. In the SMEFT, dimension-5 operators break L and consequently each L α , but the smallness of the neutrino masses (translating to the large scale suppressing the dimension-5 Wilson coefficients) suppresses their contribution to μ → eγ to an unobservable level. On the other hand, there are many dimension-6 operators that break L e and L μ , and their contributions may be more significant. In particular, one of the dipole operators in Eq. (3.14) reads L D=6 ⊃ [C eB ] 12 (l 1 Hσ αβμc )B αβ + h.c., (5.9) where l 1 = (ν e , e) is the doublet of left-handed first generation leptons. After electroweak symmetry breaking, this operator leads to the interaction term mediating μ → eγ at tree level, with cos θ W ≈ 0.89 being the cosine of the Weinberg angle. In the presence of this interaction, the branching ratio for the μ → eγ can be calculated to be where μ ≈ 3×10 −19 GeV is the total muon decay width. 20 We again peek into Particle Data Group, and find the 90% CL experimental constraint Br(μ → eγ ) ≤ 4.2 × 10 −13 coming from the MEG experiment [59]. This translates into the constraint on the Wilson coefficient The scale probed by μ → eγ is less impressive than that probed by proton decay, but it is nevertheless several orders of magnitude above the direct reach of the LHC. Again, it helps that we consider a process forbidden in the SM, so we avoid dealing with theoretical errors on the SM prediction. Moreover, muons are relatively long-lived ( μ in the denominator of Eq. (5.11) is small), and that we can easily produce and handle large amounts of them. There is one caveat concerning the scale probed by μ → eγ . In typical BSM models, the mass of the new particles that can be excluded by this constraint will be much smaller. One reason is that in perturbative models the operator in Eq. (5.9) cannot be generated at tree level, thus it will appear with at least one loop suppression factor. Furthermore, the BSM model is likely to have some form of the chiral symmetry, with some small parameters suppressing the transitions between left-and right-handed fermions. If that symmetry is akin to the one in the SM (for example, if chirality is violated only by the SM Yukawa interactions), the operator in Eq. (5.9) will be generated with the y μ ∼ 10 −3 suppression factor. But even if [C eB ] 12 = y μ e 16π 2 2 , with identified as the scale of BSM particles, μ → eγ still probes ∼ 10 5 GeV, comfortably above the LHC reach.
New SMEFT vertices violating established SM symmetries may be our best path to new physics at high energy scales. Nevertheless, not all new SMEFT vertices are of this type. There are many examples of SMEFT interactions that preserve all SM symmetries but do not appear in the SM Lagrangian, usually due to renormalizability of the latter. As an example, consider this dimension-6 operator: Its effect is to induce the Higgs boson couplings to gluons: which permits the Higgs boson to mutate into two gluons or vice-versa. Such a contact interaction term between the Higgs and gluons is absent in the SM. However, the process where two gluons collide to produce a Higgs boson does appear in the SM at the one-loop level, see Fig. 4. This is in fact the most common way the Higgs is created at the LHC. The SM loops and the contact interaction in Eq. (5.14) are in principle distinguishable experimentally, in particular they lead to a different p T distribution of the Higgs production in hadron colliders. In practice, we can best distinguish them indirectly via global fits to Higgs data, as the SM process and the Eq. (5.14) are differently correlated with the associated tth production rate. The resulting bounds are in the ballpark of |C H G | 1 (10 TeV) 2 , see e.g. Ref. [60]. This is visibly less spectacular than the bounds discussed previously in this subsection. The sensitivity to O(10) TeV scale is rather typical for new SMEFT vertices not violating any exact or approximate symmetries of the SM.

New Lorentz structures
We turn to another class of effects of higher-dimensional operators, which are related to interaction terms with different Lorentz structures compared to those in the SM. Perhaps the most iconic example in this class are the interactions contributing to the anomalous magnetic and electric moments of elementary particles. The dimension-6 SMEFT Lagrangian contains In the presence of this operator, the Lagrangian after electroweak symmetry breaking contains The first line in Eq. (5.16) contains the usual kinetic and mass terms for the electron. In the second line we have electron's interaction with the electromagnetic field. The minimal coupling proportional to q e is, of course, already present in the SM. The effect of the operator in Eq. (5.15) is to introduce another vertex with two electrons and one photon, but with a different Lorentz structure, in particular with the photon entering via the field strength F μν . As the name suggests, d e in Eq. (5.16) is the electric dipole moment of the electron, while μ e is its anomalous magnetic dipole moment. It is worth pausing here to justify this identification. The notion of dipole moments stems from non-relativistic physics, therefore it is convenient to match Eq. (5.16) to the non-relativistic EFT for the electron. This can be achieved by this change of variables: Here, ψ and ψ c are non-relativistic fermionic fields describing a spin-1/2 particle (electron) and a spin 1/2 antiparticle In other words, ψ and ψ c do not mix at the quadratic level, and they both satisfy the Schrödinger equation in the nonrelativistic limit ∇/m e 1.  (5.20) where V ≡ A 0 is the Coulomb potential, E ≡ −∇V − ∂ t A is the electric field, and B ≡ ∇ × A is the magnetic field. The first two terms describe the interaction of electron's charge with the Coulomb and vector potentials. The remaining two terms describes the interaction of its dipole moments with the external electric and magnetic fields. We can compare them with the textbook expressions for the dipole Hamiltonian: We can thus identify 21 the magnetic and electric moments of the electron at tree level: μ e = q e e m e + μ e s, d e = d e s, (5.21) where the spin vector is defined as s k = ψ † σ k 2 ψ. Furthermore, defining the g-factor for a charged particle via the 21 Recall that the sign of the potential terms is flipped between the Hamiltonian, H = T + V , and the Lagrangian, L = T − V . relation μ = g qe 2 m s, we have + μ e m e q e e , (5.22) where g SM e − 2 encodes the SM loop effects, and the contributions from higher-dimensional operators are contained in μ e . In the following let us assume that [C eB ] 11 is real (the imaginary part will be discussed later in Sect. 6), so that d e = 0 while μ e is allowed to linger. We want to determine the constraint on [C eB ] 11 from the measurement of the anomalous magnetic moment of the electron. These days, g e is predicted and measured with the incredible 10 −13 accuracy, which is often hailed as a triumph of both experimental and theoretical prowess and a powerful demonstration of the robustness of the QFT framework. The most recent experimental result comes from Ref. [61]: = 0.00115965218059 (13). To constrain new physics we also need the SM prediction g SM e . This can be calculated perturbatively, mainly in function of the fine structure constant α. Unfortunately, at this point in time there is some confusion about α. Measurements using rubidium atoms in Ref. [62] lead to 1/α(0) = 137.035999206(11), while those using cesium atoms in Ref. [63] find 1/α(0) = 137.035999046 (27), the two disagreeing at more than 5 sigma. Clearly, one or both experiments underestimated their systematic errors. To deal with this kind of situations, Particle Data Group developed a completely ad-hoc but nevertheless very useful procedure. The idea is to punish both groups indiscriminately by inflating the error bars to the point where the two measurements become consistent with each other at 1 sigma. Using this procedure, I find the combined value 1/α(0) = 137.035999183 (56), where the errors are inflated by S = 5.5. This is a loss of precision by a factor of 5(2) compared to the more (less) precise input, but, undeniably, this combination better reflects our current knowledge of α than the two individual results with the smaller errors. With this value of α one can obtain the prediction g SM e /2 = 1.00115965218045 (48). Note that the theoretical error is now almost 4 times larger than the experimental one. 22 At this point we have all the ingredients to constrain the Wilson coefficient [C eB ] 11 . Using Eq. (5.22) and replacing μ e using Eq. (5.17) we get [C eB ] 11 = 1.4(5.0) × 10 −13 GeV −2 or 11 1 (940 TeV) 2 (5.23) at 95% CL. We can see that the anomalous magnetic moment of the electron probes very high scales, although not as high as, say, μ → eγ . Moreover, similar caveat as the one discussed below Eq. (5.12) applies: in natural BSM models the chiral symmetry will typically be implemented, leading to [C eB ] 11 ∼ y e 16π 2 2 . The same scaling is true when [C eB ] 11 is induced by other dimension-6 operators via renormalization group running. If that is the case, measurements of g e currently probe the very unimpressive new physics scale ∼ 100 GeV, such that the validity range of SMEFT would be null in this scenario. The most accurate precision experiment in physics may not be accurate enough to reach new physics above a TeV! Let me also discuss a less spectacular but more typical example of new Lorentz structures in the SMEFT Lagrangian. Consider one of the bosonic dimension-6 operators in Eq. (3.11): It contributes to self-interactions of the electroweak gauge bosons, in particular, it affects the W W Z vertex. Cubic selfinteractions between electroweak gauge bosons are predicted as a consequence of the non-abelian structure of the electroweak gauge group, and the existence of the W W Z vertex was spectacularly confirmed at the turn of century by the LEP-2 experiment [64]. However, the SM predicts W W Z interactions with a single derivative acting on one of the fields, whereas Eq. (5.24) yields one derivative acting on each field. Let us put this in a wider context. The CP-conserving triple interactions of electroweak gauge bosons are traditionally parametrized as [65] where g 1,Z , κ V and λ V are real and referred to as the triple gauge couplings. It is worth knowing this parametrization because experimental collaboration more often than not present their results in this language. One physical interpretation of κ γ and λ γ is that they parametrize the dipole and quadrupole moments of the W boson: μ = e 2m W (κ γ − λ γ )s, Q i j = − 3e 2m 2 W (κ γ + λ γ )s i s j . In the SM, g Z 1 = κ γ = κ Z = 1, and λ γ = λ Z = 0, which corresponds to the W boson gfactor being equal to two at tree level, and a particular nonzero value of the quadrupole moment. In SMEFT, dimension-6 operators may lead to deviations from these predictions, while respecting the pattern κ γ = κ Z and λ γ = λ Z . For this reason, let us define g Z 1 = 1 + δg Z 1 , κ γ = κ Z = 1 + δκ, λ γ = λ Z = λ, and we will treat δg Z 1 , δκ and λ as O( −2 ) parameters in the SMEFT power counting. In particular, in the presence of just the single dimension-6 operator in Eq. (5.24) one has δg Z 1 = δκ = 0 and How can one constrain the parameters δg 1,Z , δκ and λ in Eq. (5.25)? Unlike for the electron or muon, measuring electromagnetic moments for such a short-lived particle as the W boson is not feasible. Instead, one good handle is the energy dependence of the ff → W + W − and ff → W ± Z amplitudes. In SMEFT with dimension-6 operators for E m W one expects (1) and M SM asymptoting to a constant at large energies. Therefore, the effects of the non-SM triple gauge couplings should increase with the collision energy, which provides a way to distinguish them from the SM pieces. It is fortunate that the LEP-2 phase of the CERN's e + e − collider was operating at several different center-of-mass energies. The measured e − e + → W − W + production cross section is tabulated in 10 different energy bins in the 161-207 GeV energy range [64]. Moreover, the triple gauge couplings affect differently the polarization of the produced W bosons, which is reflected in the distribution of the W production angle. All in all, there is enough information to disentangle all the three triple gauge couplings. It is a straightforward if tedious to calculate the differential e − e + → W − W + cross sections using the Lagrangian in Eq. (5.25) and then to compare them with the data in Ref. [64]. Once the dust settles one finds the result ⎛ (220 GeV) 2 at 95% CL. This is rather disappointing, and also surprising, given that the typical precision of the LEP-2 experiment is O(1%), suggesting the reach up to ∼ 3 TeV. The reason for this loose constraint is that LEP-2 alone, for certain accidental reasons, cannot efficiently disentangle all the three triple gauge couplings. The final likelihood has an approximately flat direction, as witnessed by the large off-diagonal entries in the correlation matrix (more precisely, the determinant of the correlation matrix is very close to zero).
We can nevertheless obtain stronger constraints on C 3W assuming that this is the only Wilson coefficient affecting W W production. All the necessary information for this purpose is already encoded in Eq. (5.27). First, from that equation we can reconstruct the complete likelihood: is the error matrix built from the 1σ uncertainties x and the correlation matrix ρ in Eq. (5.27), and σ −2 is the inverse of the error matrix. In this 3-dimensional likelihood function we set to zero δg Z 1 and δκ, to obtain 1-dimensional likelihood for λ. Minimizing that likelihood one finds λ = −0.049 (27), a much stronger constraint than in Eq. (5.27), although still weaker than naively expected, which has to do with the suppressed interference terms between the SM and λ contributions [66]. This translates to |C W | 1 (760 GeV) 2 (5.28) at 95% CL. The LHC also probes this Wilson coefficient. There the relevant processes at the parton level are qq → W + W − and qq → W ± Z , while the scanning of the centerof-mass energy is automatically delivered via to the parton distribution functions of the colliding protons. The recent analysis by ATLAS quotes |C W | 1 (1.6 TeV) 2 at 95% CL, a factor of four better than the LEP-2 bound in Eq. (5.28).

Modified couplings
We turn to yet another important class of effects of higherdimensional operators. The SM has merely 18 free parameters (not counting the theta term), and in terms of those it predicts countless interaction strengths between particles in the SM Lagrangian. For example, all interactions of the Higgs boson are uniquely predicted in terms of the Higgs VEV v and the SM particles masses: All of the parameters above are well known. The experimental precision varies (e.g. m Z is measured with a relative 10 −4 error, while for the up and down quark the accuracy is closer to 10%) but is invariably better than what is needed to adequately predict the LHC rates of Higgs production and decay. Staying within the SM paradigm, the LHC measurements of Higgs cross sections and branching ratios teach us next to nothing about fundamental interactions. Things are completely different in SMEFT, where literally every interaction strength in Eq. (5.29) can be altered by higher-dimensional operators. As an example, consider the dimension-6 operator in Eq. (3.11) that modifies the Higgs boson coupling to tau leptons: After electroweak symmetry breaking this becomes The h 2 and h 3 are new vertices in the nomenclature of this section. They have currently very limited phenomenological relevance, so let us leave them aside. The first two terms shift the τ mass term and the Higgs Yukawa coupling to τ :  33 [C eH ] 33 τ c τ +ττ c = m τ 1 + δy τ h τ c τ +ττ c ,  33 destroys the correlation between the Higgs boson Yukawa coupling to the τ lepton and the τ lepton mass. In other words, that coupling is modified, deviating from the SM prediction. This is in fact the best way to constrain [C eH ] 33 . Particle Data Group performs an average of the ATLAS and CMS bounds on the h → τ τ signal strength, finding (h → τ τ )/ (h → τ τ ) SM = 1.15 ± 0.15 [57]. This translates to δy τ = 0.08 ± 0.08, or 33 1 (5 TeV) 2 (5.34) at 95% CL. The reach up to a few TeV is representative to what one can currently squeeze out of Higgs physics. There are two main lessons from this simple example. One is the importance of precision measurements. C eH is just one of many dimension-6 operators that shift interaction strengths away from the SM value. Searching for such effects relies not only on improving experimental accuracy, but also on a good control of the theoretical predictions. This is often challenging, but the payoff is important: increasing precision of the measurements directly translates into increased scale of higher-dimensional operators (thus, increased scale of new physics) that we can probe. The second lesson concerns the importance of properly identifying the input parameters in SMEFT. In the above example, one had to take into account that [C eH ] 33 contributes not only to the hτ τ Yukawa, but also to the τ mass term. Had we forgotten about it, and just naively looked at the Yukawa term in Eq. (5.32), we would have obtained a wrong answer for δy τ . In this case the error committed would be of order one, but it can be much more dramatic. Consider another example, where we switch on the four-fermion operator where [C ll ] 1221 is real in our conventions (see the comment below Eq. (3.16)). This operator contains the interaction terms L SMEFT ⊃ [C ll ] 1221 (ēσ ρ μ)(ν μσρ ν e ) + h.c., which contributes to muon decay, (μ − → e −ν e ν μ )/ (μ − → e −ν e ν μ ) SM = 1 − v 2 [C ll ] 1221 . Now, muon decay is measured with impressive precision, with the relative error of order 10 −6 . One might naively jump to the conclusion that [C ll ] 1221 is stringently constrained, v 2 |[C ll ] 1221 | 10 −6 , that is to say |[C ll ] 1221 | 1 (300 TeV) 2 . This would be terribly wrong. The reason is that, in the usual approach to SM precision tests, muon decay is the standard candle that determines one of the unknown parameters of the SM -the Higgs VEV v. Indeed, the tree-level formula 23 SM μ→eνν = m 5 μ /(384π 3 v 4 ) allows one to precisely fix v, given that m μ is known with an even better accuracy. On the other hand, in SMEFT in the presence of [C ll ] 1221 one has Since both v and [C ll ] 1221 are a-priori unknown parameters, muon decay does not fix either, but just one combination of the two. Let us repeat it loudly and clearly: muon decay alone leads to no constraint at all on C ll ! Nevertheless, the effects of [C ll ] 1221 do not jus disappear, when regarded from a more global perspective. To understand how [C ll ] 1221 re-emerges we need to do a small detour first, and discuss the input parameters for electroweak precision tests. The latter can be defined as a set of observables that, in the SM at tree level, depend on the parameters g L , g Y and v in the electroweak sector. The numerical values of these parameters are traditionally fixed by three precisely known input observables: 1. The Fermi constant G F , extracted from the measured muon lifetime using the formula SM μ→eνν = G 2 F m 5 μ /(192π 3 ) plus radiative corrections [67]. 2. The electromagnetic structure constant α, currently best extracted from the spectroscopy of rubidium [62] and cesium [63] atoms. 3. The Z boson mass m Z , extracted from the position of the corresponding resonance in e + e − scattering in the LEP-1 collider [68].
In SMEFT, in the presence of [C ll ] 1221 , working at O( −2 ) and at tree 24 level, these input observables are connected to the SMEFT parameters as where I re-labeled the parameters from the SM Lagrangian using the subscript zero to distinguish them from the g L , g Y and v parameters in the following, which will differ by O( −2 ), and which will be assigned definite numerical values. Note at, at O( −2 ), it does not matter whether I write v 2 C ll or v 2 0 C ll . If [C ll ] 1221 = 0 then, as in the SM, Eq. (5.36) relates 3 parameters to 3 observables and we can readily solve for v 0 , g L0 , and g Y 0 . However when [C ll ] 1221 = 0 we have 23 Of course, at this level of precision, one should also take into account the radiative and O(m 2 e /m 3 μ ) corrections when relating v to the observable decay width. This does not interfere with the following discussion. 24 Once again, radiative corrections from D ≤ 4 operators must be taken into account in the matching of the input parameters to observables in order to meet the required precision level. The procedure is in fact very similar to our treatment of higher-dimensional effects.
3 equations for 4 parameters. In this case it is convenient to use a trick: we can get rid of [C ll ] 1221 from Eq. (5.36) by absorbing it into the other parameters. These can be achieved by defining which means that g L , g Y , and v are related to observables in exactly the same way as the corresponding SM parameters, and therefore they can be assigned exactly the same numerical values. Plugging in the numbers into Eq.

One more thing
In this section, for the sake of simplicity, the discussion of phenomenological effects of higher-dimensional operators is divided into "new vertices", "new Lorentz structures", and "modified couplings" parts. It is however important to mention that sometimes there is no invariant way to make this distinction. Consider the following example of self-interactions of the Higgs boson h: where the dots denote terms with 4 and more Higgs bosons. Two possible effects of higher-dimensional operators appear above. The one proportional to δ 1 changes the magnitude of the triple Higgs self-coupling, which is already present in the SM Lagrangian but with the magnitude strictly fixed by the Higgs boson mass. This is a modified coupling in our nomenclature. The other effect proportional to δ 2 is a two-derivative Higgs self-interaction term which does not appear in the SM Lagrangian in its canonical form. This is a new Lorentz structure in our nomenclature. Both δ i can be generated by dimension-6 operators, therefore we will treat δ i as O( −2 ). For example, switching on the Wilson coefficients C H and C H in Eq. (3.11), one gets terms resulting from the redefinition. Seemingly, the Lagrangians in Eqs. (5.43) and (5.45) are different, as they contain different interaction terms. However, the equivalence theorem ensures that field redefinitions cannot change the physical content of the theory, as was discussed in Sect. 4 in the context of SMEFT bases. Therefore, the two Lagrangians give exactly the same predictions for physical observables, at any order in the perturbation theory, which can be verified by explicit calculations. This demonstrates that we can get rid of some of the new Lorentz structures generated by higher-dimensional operators and absorb them into modified couplings. Conversely, a shift h → h − δ 1 2v h 2 applied to Eq. (5.43) would erase the modified coupling in favor of the new Lorentz structure. Once the fermions are included in the discussion, the above shift of the Higgs boson field produces Yukawa-like interactions between two powers of the Higgs boson and the SM fermions. Such interactions are absent in the SM Lagrangian, thus they are "new vertices" in our nomenclature. Therefore the boundaries between new vertices and modified couplings or new Lorentz structures are also blurred by field redefinitions. The sharp boundary persists only for the new vertices violating the SM symmetries, as those can never be redefined away into modified couplings or new Lorentz structures.

On the importance of global fits
In this section I have presented several examples of experimental constraints on Wilson coefficients of dimension-6 operators. Most of the time, in the derivation I was assuming that just a single Wilson coefficient was present at a time. This was for the sake of simplicity and clarity. My main goal was to sketch the path between the SMEFT Lagrangian and observable quantities. One should be aware, however, that such one-at-a-time results are rarely useful in practice. The higherdimensional SMEFT operators are meant to describe effects of heavy BSM particles. But BSM models will never produce a single operator. Typical ones will produce dozens already at tree level. Even if not, renormalization group running and other loop effects will make sure to generate additional operators. Note also that one-operator-at-a-time is inherently a basis-dependent statement: what is a single operator in one basis may map to several operators with correlated Wilson coefficients in another basis.
For the practical purposes what you need is a global likelihood function. This should be derived without assuming that any particular Wilson coefficient in the SMEFT Lagrangian vanishes. A single observable typically depends on several Wilson coefficients but, at least at tree level and truncating the analysis at O( −2 ), the relevant linear combination of Wilson coefficient can be determined by the similar techniques as discussed in this section. Beyond tree level it may be more challenging to proceed in full generality, but the more and more widespread availability of automatization tools should help to overcome the difficulties. The real power of this approach is highlighted when multiple observables are combined. More often than not, different observables will depend on different linear combinations of Wilson coefficients, therefore they are complementary. If a large enough number of observables is included in the global fit, multiple Wilson coefficients can be independently constrained without any simplifying assumptions. The taste of this procedure was given around Eq. (5.27), where constraints from the LEP-2 W W cross section measured in different angular and energy bins were combined to simultaneously determine three triple gauge couplings of electroweak gauge bosons.
The likelihood function obtained in such global fits can be easily shared in a numerical form. If the likelihood is approximately Gaussian, the full information can be succinctly encoded by giving the central values, uncertainties, and the correlation matrix for the constrained Wilson coefficients. This likelihood can then be used by other researchers to constrain their preferred BSM models. There are many groups extracting such likelihoods from different classes of observables: Higgs physics, electroweak precision tests, flavor physics, etc. Each such work can be regarded as a lego block to create an even larger likelihood function covering more aspects of physics beyond the SM. Combining different likelihoods often brings forward important synergies. For example, within the SMEFT framework there is an important complementarity between measuring the triple gauge couplings and the Higgs coupling. Combining the two sets of observables allows one to lift approximate flat directions that occur in each standalone set. An advanced example of such brinkmanship is Ref. [69], where a wide selection of electroweak precisions tests at low and high energies is used to simultaneously constrain 73 linear combinations of dimension-6 Wilson coefficients. This is still a bit less than the complete set of 3045 dimension-6 Wilson coefficients in the SMEFT Lagrangian, but we're getting there:)

CP violation in SMEFT
Charge conjugation (C) and parity (P) are badly violated in the SM. Their combination, C P, is not a symmetry of the Lagrangian either, however, due to certain accidents of nature, the observable effects of the SM sources of C P violation are feeble. This leaves ample room for C P-violating higher-dimensional operators to leave an imprint.

Formalism of CP violation
Let us begin by discussing how C P acts on the quantum fields in SMEFT. The Higgs field H is scalar 26 and the transformation is quite simple: The action of P is responsible for x → −x, and the action of C is responsible for the complex conjugation on the right-hand side. This in particular implies that the Higgs boson, which lives in the real part of H , is C P even: For the gauge fields, P and C act as The action of parity is the usual one for vectors. For electromagnetic fields it implies the well known fact that the electric field E is P-odd, while the magnetic field B is P-even. The minus in the action of C is intuitive: changing a particle into its anti-particle flips the electric charge and thus the sign of the potential. It follows that C P acts as It is easy to show that the kinetic term V μν V μν is C P even. On the other hand V μνṼ μν is C P odd: In the next-to-last step we changed the integration variables x → −x. This shows that the θ term in the D = 4 Lagrangian Eq. (3.7) violates C P. Likewise, the bosonic interactions proportional to the Wilson coefficients C H G , C H W , C H B , and C H W B in Eq. (3.11) are C P odd. By an analogous calculation one can show that also the ones proportional to C W and C G are C P violating. Note that in these cases C P violation is not associate with complex couplings in the Lagrangian. Let us move to the fermionic fields. One can show that parity and charge conjugation act on two-component fermion fields as , x), x). (6.6) Parity, intuitively, exchanges left-and right-handed fermions.
The action of C provides the rationale for the · c notation. It follows that C P acts as Now we are ready to investigate how C P acts on a fermionic mass term in the action: x). (6.8) In the last step we changed the integration variables x → −x, and used that for anti-commuting fermion fields we have ψ α ψ c α = −ψ c α ψ α = ψ c α ψ α . It follows that complex mass terms violate C P. This is however not the whole story. The chiral transformation: with an arbitrary phase η, removes the complex phase of the mass term and restores C P conservation at the level of the kinetics terms. Thus C P violation can only arise due to unremovable relative phases between fermionic mass and interaction terms. For example, by exactly the same calculation as in Eq. (6.8), C P transformation of a Yukawa interaction of the fermion with a real scalar field h is given by , x). (6.10) If Arg[m] = Arg[y] then the relative phase cannot be removed and C P violation is physical. This simple example illustrates the fact that C P is a collective endeavour. Phases in a single interaction are not physical, as they can always be transferred elsewhere by a phase redefinition of fermionic fields. For this reason, the mathematically proper language to discuss C P involves invariants under phase and other unitary transformations. In the SM there is only one such object: the famous Jarlskog invariant. SMEFT proliferates the number of physical phases and invariants; at the level of dimension-6 operators they were recently classified in Ref. [70].

CP violation at D = 4
To set the stage, I will now briefly discuss C P violation at the level of D = 4 operators in SMEFT, which is the same as C P violation in the SM. This will be the springboard to describe observable C P-violating effects of selected higherdimensional SMEFT operators in the following subsections. The D = 4 SMEFT Lagrangian in Eq. (3.7) contains the so-called theta term: μν . (6.11) This is a total derivative, which does not affect any perturbative amplitudes. However, at the non-perturbative level this term does contribute to C P-violating hadronic observables, in particular to the neutron EDM. The most recent lattice evaluation in Ref. [14] estimates d n = 0.158(36)θ e fm. (6.12) The current experimental limit is d n = (0.0 ± 1.1) × 10 −13 e fm [71], which leads to the 95% CL limit |θ | 1.4 × 10 −12 . (6.13) The smallness of C P violation due the theta term is one of the biggest mysteries of contemporary particle physics.
The D = 4 SMEFT Lagrangian also contains the Yukawa interactions: L D=4 ⊃QHY uŪ c +Q HY dD c +L HY eĒ c + h.c. (6.14) As discussed around Eq. (6.10), Yukawa terms transform under C P as Y f → Y * f , therefore complex Yukawa phases may lead to C P violation. However, phases in standalone Yukawa terms are removable, therefore they have to be discussed in conjunction with other interactions in the SMEFT Lagrangian. The leptonic part will be discussed later in the context of D = 5 operators, and we focus on the quark Yukawas for now. The usual way to isolate the physical phases is to first write down the quark mass term after electroweak symmetry breaking:  27 This transformation leaves invariant most of the terms originating from the D = 4 Lagrangian, in particular R f cancel out completely. The only place where L f emerge is in the charged current interactions between the W boson and left-handed quarks: (6.15) where the CKM matrix V is related to the left rotations as V = L † u L d and therefore it is unitary, V † V = 1. By an analogous calculation as the one in Eq. (6.8) one can show that the C P transformation effectively takes V → V * , thus phases of the CKM matrix may lead to C P violation. A general unitary matrix has 9 free parameters: 3 angles and 6 phases, however most of the phases are removable. Indeed, we can make the transformations u J → e iφ J u J , for I = J = K . All in all, for 3 generations there remains a single order parameter for C P violation via the CKM matrix, the so-called Jarlskog invariant: (6.16) while the other components of J are given by the formula J I J K L = J M,N I K M J L N . In the SM, this single invariant controls C P violation observed in multiple processes with mesons containing strange, charm, and bottom quarks. From these observations, assuming negligible contributions from higher-dimensional operators, one can extract J = 3.08(14) × 10 −5 [57]. This is a small number but, actually, the phase itself in the CKM matrix is large, order one. In particular, in the widely used Wolfenstein parametrization of the CKM matrix, where the phase is located mostly in the V ub and V td elements, one gets Arg −0.4. However the Jarlskog invariant is suppressed by the smallness of the absolute values of the offdiagonal CKM elements: |V us | ∼ λ, |V cb | ∼ λ 2 , |V ub | ∼ λ 3 , where λ ∼ 0.2 is the Cabibbo angle. Due to this hierarchical structure, J ∼ λ 6 , leading to a severe suppression of observable signals of C P violation. Although C P violation was first discovered some 60 years ago [72], to this day experimental detections of C P-violating quantities are few and far between, and each new detection is announced with fanfare.
The partly mysterious and partly accidental suppression of C P violation at the level of the D = 4 Lagrangian opens the playground for higher-dimensional C P-violating operators. D = 5 operators introduce three more C P invariants, including the so-called PMNS phase, which is a subject of an intense ongoing experimental program in neutrino physics. At the level of D = 6 operators all hell breaks loose, and there are literally hundreds of new C P invariants [70]. Dimension-6 operators can provide tree-level contributions to basically any C P-violating observable, including the EDMs of elementary (electron) and composite (neutron, nuclei, molecules) particles, and the neutral meson mixing (KK , DD, BB). In the rest of this section we discuss selected C P-violating observables and their dependence on higher-dimensional SMEFT operators.

CP violation at D = 5
The dimension-5 operator in Eq. (3.8) after electroweak symmetry breaking yields the neutrino mass terms. Together with the charged lepton masses from the D = 4 Yukawa interactions, the lepton masses are given by An analogous calculation as in Eq. (6.8) shows that C P takes Y e → Y * e , C 5 → C * 5 , therefore the phases in C 5 and Y e are potentially C P-violating. As in the quark case earlier, in order to identify the physical phases it is convenient to make the masses diagonal and real. One can decompose Y e = L e M e R † e , C 5 = L T ν M ν L ν (recall that C 5 is symmetric), where L e , R e and L ν are unitary, while M e and M ν are diagonal, real, and positive. Rotating the fields as ν → L ν ν, E → L e E, E c → E c R † e , the phases then re-emerge in the W boson interactions with left-handed leptons: (6.18) where the PMNS matrix U = L † e L ν is unitary. In general, a unitary matrix has 6 phases. We can remove 3 phases from U by redefining e J → e iφ J e J , e c J → e iφ J e c J . Contrary to the CKM case, we cannot remove further two phases, because the transformations ν n → e iη n ν n would feed the complex phases back into the neutrino mass terms. Thus, there are three physical phases in U . One corresponds to the usual Jarlskog invariant, The phase controllingĴ is called the Dirac phase, and it remains physical in the limit of vanishing neutrino masses. It is measurable in oscillation experiments as it affects the oscillation probabilities. Ignoring contributions from higherdimensional operators, these can be written as where 2 kl ≡ m 2 ν k − m 2 ν l are the differences of the squared masses of neutrino eigenstates. The C P-violating observable is the difference between oscillation probabilities of neutrinos and anti-neutrinos of a given flavor: This explicitly shows that C P violation in neutrino oscillations at the D = 5 level is controlled by the single Jarlskog invariantĴ defined in Eq. (6.19). Determining whetherĴ is non-zero is subject to an intensive experimental program, with several ongoing experiments (T2K, NOνA, Super-K) placing non-trivial limits. Some of these even provide preliminary hints thatĴ = 0, and that the phase in the PMNS matrix, much like in the CKM matrix, is large [73]. The matter should be ultimately resolved in the coming years, or at latest in the next decade, when the future DUNE and Hyper-K experiments collect enough data. The importance of these measurements go beyond measuring yet another parameter in the PMNS matrix. IfĴ = 0, it will be another (after the CKM phase) demonstration that phases are present in the fundamental Lagrangians. Philosophically, it will be a strong hint that phases are generic, that is to say, they arise whenever they are allowed to, and there is no physical principles suppressing them. Adopting this point of view would have tremendous consequences for BSM model building. 28 The other two physical phases in the PMNS matrix can be isolated by rewriting U J n =Û J n e iγ n , whereÛ J n contains only the Dirac phase. The γ n phases cancel out inĴ J K mn but the two relative phases (the differences γ n − γ m ) are physical. They are called the Majorana phases and they enter into predictions for physical observables, such as neutrino oscillations into anti-neutrinos, or neutrinoless double beta decay. However they become unphysical in the limit of vanishing neutrino masses, as in this limit one is free to use phase redefinitions of the neutrino fields to get rid of these phases. Therefore, their effects on physical observables must always be proportional to the neutrino masses, and therefore extremely suppressed. For the moment, there is little hope of ever measuring them in low-energy experiments.

Electric dipole moments
Electric dipole moments, abbreviated as EDMs, are mundane in classical electrodynamics. Take two opposite charges ±Q, put them very close together at a distance R, and what you get from a sufficiently large distance r R looks like a particle without an electric charge but with the electric dipole moment |d| = Q R. On the other hand, EDMs of elementary particles are dramatic because they signal CP violation. To date, no such EDMs have been detected, in spite of the existence of CP violation in the SM -the contributions proportional to the CKM phase are still many orders of magnitude below the experimental sensitivity. The same remains 28 The usual sales pitch for advertising neutrino C P violation is the question of the origin of matter and its dominance over anti-matter. This is however very controversial: there is no model-independent or natural way to connect the phase measured in neutrino oscillations to the phases responsible for the asymmetry between matter and anti-matter. true concerning the EDMs of sufficiently simple composite particles: protons, neutrons, atoms in ground state in vacu-um…To observe permanent EDMs one has to go up to the molecule level, the prime example being the water molecule with d H 2 O 2 × 10 −5 eV −1 .
Let us begin the discussion with the simple case of the electron. The SMEFT Lagrangian contains the dimension-6 operators 11 (l 1 Hσ μνēc )B μν + h.c. (6.22) where l 1 = (ν e , e) is the first generation left-handed lepton doublet. Once the Higgs gets a VEV this leads to where (6.24) As the notation suggests, μ e is the anomalous magnetic dipole moment of the electron, and d e is its electric dipole moment. Therefore, d e = 0 in the presence of imaginary Wilson coefficients in Eq. (6.22). See the derivation of Eq. (5.20) to understand the precise connection between the coefficients of the F μν f c σ μν f interaction, and the dipole moments of spin-1/2 particles. The current experimental constraint on the electric dipole moment of the electron is d e = (4.3±4.0)×10 −30 e cm [74], which is d e = (6.6 ± 6.2) × 10 −17 GeV −1 in natural units. Assuming C eW = 0, using Eq. (6.24) we can translate the experimental constraint into a bound on the SMEFT Wilson coefficient C eB : Similar constraints can be obtained for C H W , and some components of C lequ are also strongly constrained by the electron's EDM. The power of the electron EDM is such that even operators contributing via two-loop effects are subject to relevant constraints [76]. We move to discussing another important EDM. This time it is about a composite particle -the neutron. Because of the strong interactions involved, connecting the neutron's EDM to SMEFT Wilson coefficients is a complicated business, involving difficult calculations in Chiral Perturbation Theory, lattice determination of nuclear matrix elements, phenomenological extraction of nuclear parameters, etc. I am not daring enough to review the necessary technology here. Instead, I will just quote the result from Ref. [77]: d n = (43 ± 27)Im C us us 1 L R + (210 ± 130)Im C us us where the Wilson coefficients above are defined by the following Lagrangian below the electroweak scale, Thus, in sheer numbers, the scale of the SMEFT operators probed by the neutron EDM is some two orders of magnitude below the scale probed by the electron EDM, cf. Eq. (6.25). Let us also study another problem: what are the neutron EDM constraints on the C Hud Wilson coefficient in Eq. (3.13) ? This may be a surprising question, since that operator multiplied by C Hud , which isH † D μ H (u c σ μdc ), couples the Higgs field to right-handed quarks, and at first sight it has little to do with EDMs. It turns out that the effect on the neutron EDM originates from the fact that this operator induces a coupling of the W boson to the right-handed up and down quarks:  11 . Applying again Eq. (6.28), ignoring the theoretical error for the sake of sim-plicity, and comparing with the experimental constraint we find |Im [C Hud ] 11 | 1 (100 TeV) 2 (6.32) at 95% CL. This is again an impressive sensitivity. As remarked above, the result may be somewhat unexpected, because the connection between C Hud and EDMs is not obvious. Finally, we also discuss another class of observables, namely EDM of diamagnetic atoms (such as 199 Hg) and paramagnetic molecules (such as ThO). These are in fact the same observables from which the electron EDM is extracted. However, they are sensitive to other interactions beside the electron-photon ones. In particular, they are instrumental to probing the following SMEFT operators: Above q and l are the 1st generation quark and lepton doublet fields. These operators simultaneously generate charged current interactions (relevant e.g. for nuclear beta decay, and for pp → ν signals at colliders) and the neutral current ones: If the Wilson coefficients of these interactions have imaginary parts, at even lower energies they match to the nucleonelectron contact interactions i(N N)(ēγ 5 e), which in turn contributes to the EDM of diamagnetic atoms and paramagnetic molecules. Using the theoretical expressions from Ref. [78] and the EDM measurement using the ThO molecule [74] we find the 95% CL constraints Once again we are sensitive to scales 4 orders of magnitude above the LHC reach. Assuming the operators in Eq. (6.33) have generic phases, the CP constraints exclude the possibility that they have any visible effects of in nuclear beta decay or at the LHC,

Meson mixing
Mesons are bound states of a quark and an anti-quark. Of particular interest in this section are neutral kaons, who are spinzero mesons made of down (d) and strange (s) (anti-)quarks. One can construct two neutral kaon states: |K 0 ≡sd and its antiparticle |K 0 ≡ sd. From these two one can construct the CP eigenstates |K 0 ± = 1 √ 2 |K 0 ± |K 0 , such that C P|K 0 ± = ±|K 0 ± . Now, due to CP violating interactions, the observed mass eigenstates K S and K L are not C P eigenstates, but instead The "short" eigenstate K S decays mostly to two pions with a lifetime of about 10 −10 second, whereas the "long" eigenstate K L decays mostly to three pions and lives about 1000 times longer due to the phase space suppression. However, due to CP violation parametrized by K , one in a thousand of the long eigenstates decays to two pions as well, as observed in the Cronin and Fitch experiment in 1964. From the existing K L → ππ data one can extract | K | = 2.228(11) × 10 −3 [57]. The decay K S → 3π has not been observed to date. In the EFT below the electroweak scale, obtained from the SM by integrating out the W bosons at one loop, one generates the four-fermion operator with the numerical coefficient in front c ∼ 2.3. This operator violates the strangeness quantum number by two units, therefore it has good quantum numbers to contribute to K -K mixing. Furthermore, the Wilson coefficient in Eq. (6.37) inherits the complex phase from the CKM matrix, therefore it violates CP. Consequently, K is predicted to be non-zero in the SM: | SM K | = 2.027(195) × 10 −3 [79]. 29 While there is quite some lattice hocus-pocus involved to translate Eq. (6.37) into K , the error of the SM prediction is largely due to the uncertainty on the CKM elements. Note that the theoretical error dominates the experimental one by an order of magnitude. 29 In the Wolfenstein parametrization of the CKM matrix, (V * ts V td ) 2 ∼ λ 10 A 4ρη V 4 cb λ 2ρη . As a result, the prediction for K depends significantly on the value of V cb we insert as the input. There are competing determinations of V cb : the so-called inclusive and exclusive ones, yielding results in tension with each other, which has a large (magnified by the 4th power) impact on SM K . The quoted prediction uses the inclusive determination of V cb , but note that the most recent exclusive determination [80] is in good agreement with the inclusive value. All in all | K |/| SM K | = 1.099(97), that is the ratio is known with a 10% error.
The suppression scale of the effective operator in Eq. (6.37) is by the whopping O(10) PeV, rather than by 100 GeV characteristic of the particles being integrated out. The additional suppression, apart from the one-loop factor, is due to the socalled GIM mechanism [81]. This forces the main contribution to involve a top quark, and the effect is suppressed by the small CKM mixing between the first two generations and the third. The suppression of the SM contribution opens the playing field for new physics. In SMEFT there are several dimension-6 four-fermion operators changing strangeness by two units [17]: where the quark doublets above are q 1 = (u k V * kd , d), q 2 = (u k V * ks , s) (that is, we use the down-type basis). The operators in the first line directly correct the Wilson coefficient of the WEFT operator Eq. (6.37). The remaining operators yield different chirality and/or color structures, but they also contribute to K at tree level. Reference [82] estimates the contributions of the SMEFT Wilson coefficients in Eq. (6.38) to K as Assuming only one Wilson coefficient is switched on at a time, one finds the following 95% CL bounds Im [C (8) qd ] 2121 1 (240 PeV) 2 , (6.40) where the Wilson coefficients are defined at the electroweak scale. These quantitative results again rely on some lattice inputs, but their order of magnitude could have been guessed without calculations. C (1) qq leads to the same form of the interaction below the electroweak scale as the SM loop, and the bound is similar as the magnitude of the operator in Eq. (6.37) given the SM prediction is known only at the 10% level. Idem for C (3) qq . For C dd , the low-energy operator differs from the one in Eq. (6.37) by trading left-handed quarks for righthanded, and the bound is again the same as for C (1) qq because QCD is vector-like and cannot tell left from right. It is a bit less trivial to understand why the bounds on C (1) qd and C (8) qd are an order of magnitude stronger. This is because these operators marry left-and right-handed quarks, and gain in the course of low-energy running thanks to not having to rely on quark masses for chirality violation.
All in all, CP violation in kaon mixing impressively probes the scales up to almost EeV, five orders of magnitude above the direct reach of the LHC. New physics at the TeV scale typically requires some elaborate mechanism to suppress the operators in Eq. (6.38), causing much headache for pre-LHC model building. In supersymmetry it has motivated for example the baroque edifices around gauge mediated supersymmetry breaking. The situation was not much better in the Randall-Sundrum family of models, even in spite of the natural mechanism of suppressing flavor-changing neutral currents that these models offered.
In a similar fashion, one can constrain the SMEFT operators affecting the mixing of neutral B d and B s mesons. The analogous loop diagram to the one generating Eq. (6.37) yields the 4-fermion WEFT operators violating the b number: From this one can guess that CP violation in the neutral Bmeson mixing probes dimension-6 SMEFT operators suppressed by 1 PeV, only a tad lower than the kaon mixing. This is indeed borne out by more careful calculations [83]. Furthermore, one also generates the flavor-changing neutral current operator (cσ μ u) 2 , with the suppression scale similar to the one Eq. (6.37). Thus, CP violation in DD mixing also probes O(10) PeV scale [82]. The sensitivity of neutral meson mixing to new physics is visualized in Fig. 5.

Nuclear beta decay
A less frequented path to CP violation leads through nuclear beta decay. We will see it is a less promising place to find CP violating new physics. Nevertheless, the discussion is useful to highlight correlations and interdependencies between different probes of CP violation in the context of SMEFT. Consider the process N → N e −ν e (β − decay) or N → N e + ν e (β + decay). Here N and N are called the parent and daughter nuclei respectively, the electron and the positron are jointly referred to as the beta particle, and both neutrino and anti-neutrino are referred to simply as neutrino in the following. The simplest example is the neutron decay, where N is the neutron and N is the proton, but many more processes where N and N are heavier nuclei are studied in laboratories. We are dealing with a 3-body decay process, thus there is a lot of kinematics and differential distributions to play with. For the so-called allowed beta decays, at the leading (zero-th) order in expansion in 1/m N , after summing over beta particle and daughter nucleus polarizations, the differential distribution of the decay products takes the most general form [84] where m e is the electron mass, J is the polarization vector of the parent nucleus and J is its spin, j is the unit vector in the polarization direction, and k e , k ν , E e , E ν are the 3-momenta and the energies of the beta particle and of the neutrino.
The correlation coefficients in Eq. (6.42) are denoted somewhat inaptly as b, a, A, B,ĉ, and D. All of these coefficients except forĉ have been measured experimentally for many transitions, sometimes with an accuracy better than per-mille. Historically, these measurements have played a crucial role in the development of the theory of fundamental interactions. Since momenta are parity odd, whereas angular momenta (hence J) are parity even, non-zero values of A and B signal parity violation. The famous Wu experiment [85], which demonstrated parity violation in weak interaction, was in fact a measurement of the A coefficient in the beta decay of 60 Co.
For the sake of this section we will focus on the D coefficient. It quantifies the triple correlation between the polarization of the parent nucleus and the momenta of the beta particle and the neutrino. Since both momenta and angular momenta flip signs under T , a non-zero D signals violation of time reversal invariance. Due to C PT conservation in QFT, C P violation implies T violation, therefore C P violating interactions can contribute to the D parameter. Within the SM, D is very suppressed and has not been observed in any transition yet. 30 This is just as good, as it leaves a lot of room to spot the possible contributions to D from non-SM CP violation, in particular from complex phases of dimension-6 operators in SMEFT.
There are several SMEFT operators that contribute to the D parameter at tree level. Below we only discuss a couple of them, but the lesson we will draw is more general. Consider first the dimension-6 operator † D μ H (u c σ μdc ) + h.c. (6.43) This operator induces a coupling of the W boson to the righthanded up and down quarks: [C Hud ] 11 u c σ μdc + h.c. (6.44) in addition to the couplings to left-handed leptons and quarks present in the SM. To understand the consequences for beta decay we need to descend a few steps in the EFT ladder.
Integrating out the W boson at tree level, in the EFT below the electroweak scale one obtains the quark-level chargedcurrent interactions mediating beta decay: L WEFT ⊃ − 2V ud v 2 (ēσ μ ν)(ūσ μ d) + R (ēσ μ ν e )(u c σ μdc ) +h.c., Beta decays occurs at the level of nucleons or nuclei, thus we need to descend further. The nucleon-level Lagrangian relevant for beta decay matching to Eq. (6.45) is 30 It is not completely zero because CP conserving interactions also can contribute to D at one loop via the so-called final state interactions [86], that is to say, electromagnetic interactions between the daughter nucleus and the beta particle. This process leads to |D SM | 10 −4 , strongly depending on the transition. Contribution to D SM from the CP violation in the SM is completely negligible. L nucleon ⊃ − V ud v 2 g V (ēσ μ ν e )(pγ μ n) −g A R (ēσ μ ν e )(pγ μ γ 5 n) + h.c., (6.47) where for nucleon we use the Dirac fermion notation. Our ignorance of the non-perturbative aspects of the matching has been swiped under the carpet of the g V and g A parameters. Now, one can show that g V = 1 in the limit of unbroken isospin symmetry, and moreover g V ≈ 1 up to O(10 −4 ) corrections quadratic in isospin breaking [87]. On the other hand, g A is not protected by isospin symmetry, and we have to resort to lattice determinations, g A = 1.246 (28) [88]. From Eq. (6.47), using standard techniques to calculate amplitudes and distributions, a tedious but not impossible calculation yields (6.48) The current limit on the D parameter in neutron decay is D n = −1.2(2.0) × 10 −4 [57], while the SM contribution predicts D SM n ∼ 10 −5 . It appears that neutron decay is sensitive to non-standard CP violation, probing the scales of Im [C Hud ] 11 on the order of 10 TeV. There is an issue however: the same operator is better probed by EDMs [22]. We discussed this in the context of the neutron EDM, finding |v 2 Im [C Hud ] 11 | 6 × 10 −6 at 95% CL, see Eq. (6.32). We can rewrite Eq. (6.48) in the form |D n | ≈ 2 × 10 −6 v 2 |Im [C Hud ] 11 | 6 × 10 −6 , (6.49) which demonstrates that |D n | 2 ×10 −6 . This is two orders of magnitude below the current experimental sensitivity, and unlikely to be reached in the near future. One should not expect a discovery of CP violation in nuclear beta decay, at least in the scenario sketched above. Can we avoid this negative conclusion? One way is to ask for help our good old friend the fine-tuning. A larger D parameter can be achieved if one allows for some fine-tuning between different contributions to the neutron EDM. For example, one can arrange for a partial cancellation between the contributions proportional to Im [C Hud ] 11 and those proportional to the QCD θ parameter [89]. Another way is to resort to a more baroque model building.

Closing comments
In this section we have studied a couple of examples illustrating sensitivity of C P-violating observables to new physics. These observables probe D = 6 operators suppressed by scales much higher than the reach of current colliders -sometimes as large as EeV! Indeed, C P violation is the most sensitive probe of generic new physics that does not violate the SM symmetries. 31 The non-discovery of non-SM C P violation foreshadowed the non-discovery of new physics at the LHC. It has been known for decades that new particles with TeV scale masses and order one couplings to the SM can easily be at odds with the EDM constraints. This is especially true for frameworks addressing the hierarchy problem (supersymmetry, composite Higgs, extra dimensions, etc.) as they invariably come with many new interactions, which potentially contain many new C P-violating phases. The past expectations of dozens new particles popping out at the LHC were always relying on the belief that new physics, in stark contrast to the SM, is non-generic, that is to say, that the new C P invariants associated to the new particle interactions vanish. Well…now we are older and wiser. While new physics can still be around the corner, and be discovered tomorrow at the LHC, in my opinion a more probable scenario is that the new particles are very heavy, perhaps 100 TeV, perhaps more. At the same time, my expectation is that new physics, much like the SM, comes with generic C P-violating phases. If these expectations are correct, it is very likely that new physics will first be discovered via C P-violating observables, in particular via the electron or neutron EDM.
One more comment is in order. We have seen that distinct C P-violating observables may probe the same operator, but the sensitivity may be different. For example, both the neutron EDM and neutron beta decay constrain the C Hud Wilson coefficient, with the former having O(100) better sensitivity. The gap is even larger, O(10 5 ), concerning the sensitivity to C H B from the electron EDM and from Higgs physics. The logical conclusion is that some experiments are more destined to discover new physics. Indeed, if C P violation was detected in, say, Higgs decays to Z bosons, my theorist bias would strongly suspect an unfortunate fluctuation or an experimental error. One should however not jump to the conclusion that only EDMs are worth pursuing. Less sensitive experiments also play an important role. The SMEFT parameter space is multi-dimensional, and one needs many different probes to disentangle the different relevant parameters. For example, the neutron EDM is sensitive (among others) to a linear combination of C Hud and the theta term θ , while the D parameter measured in nuclear beta decay is sensitive only to the former. We cannot exclude the possibility of large cancellations (accidental or implied by the structure of the BSM theory) between C Hud andθ , and the flat direction in this 2-parameter fit is only lifted by the D parameter measurement. The less sensitive experiments will become even more vital if non-SM C P violation is finally observed in an EDM experiment, as they will play the crucial role in the task of identifying the UV completion of SMEFT.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 . SCOAP 3 supports the goals of the International Year of Basic Sciences for Sustainable Development.

Appendix: Abbreviations
• BSM: Beyond the Standard Model. In these lectures I define it as any theory with new particles in addition to the ones present in the SM. • ChPT: Chiral Perturbation Theory. The effective theory of pions and sometimes kaons and eta mesons, valid below the QCD resonance scale m ρ 770 MeV. • CL: Confidence level. This refers to the size of the confidence intervals in the standard Gaussian statistics. In particular, 68% CL is roughly one sigma, while 95% CL translates roughly to 1.96 sigma. • eV: Electronvolt. It cannot be logically explained why this is the basic unit in particle physics. EFT of a single massless spin-2 particle. • GRSMEFT: General Relativity Standard Model Effective Field Theory. The EFT of a single massless spin-2 particle coupled to the SM particles.
• IR: Infrared. In the particle physics literature this is a posh way to say low-energy.