Equations of motion, symmetry currents and EFT below the electroweak scale
Andreas Helset, Michael Trott

TL;DR
This paper derives the equations of motion and symmetry currents for the low-energy effective field theory below the electroweak scale, incorporating higher dimensional operators and analyzing their impact on conserved currents and gauge interactions.
Contribution
It provides a detailed construction of equations of motion and symmetry currents in the EFT below the electroweak scale, including effects of operators up to dimension six.
Findings
Higher dimensional operators modify equations of motion and conserved currents.
Surface terms relate to multipole charge distributions sourcing gauge fields.
Electromagnetic current corrections reflect multipole contributions.
Abstract
The low-energy effective field theory is constructed by integrating out Standard Model states with masses proximate to the electroweak scale. We report the equations of motion for this theory, including corrections due to higher dimensional operators up to mass dimension six. We construct the corresponding symmetry currents, and discuss how the symmetry, and global symmetries, are manifested when Standard Model states are integrated out. Including contributions from higher dimensional operators to the equations of motion modifies the interpretation of conserved currents. We discuss the corrections to the electromagnetic current as an example, showing how modifications to the equation of motion, and corresponding surface terms, have a direct interpretation in terms of multipole charge distributions that act to source gauge fields.
Click any figure to enlarge with its caption.
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6
Figure 7
Figure 8
Figure 9
Figure 10
Figure 11
Figure 12
Figure 13Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Equations of motion, symmetry currents and EFT below the electroweak scale
Andreas Helset and Michael Trott
Niels Bohr International Academy and Discovery Center, Niels Bohr Institute, University of Copenhagen, Blegdamsvej 17, DK-2100 Copenhagen, Denmark
Abstract
The low-energy effective field theory is constructed by integrating out Standard Model states with masses proximate to the electroweak scale. We report the equations of motion for this theory, including corrections due to higher dimensional operators up to mass dimension six. We construct the corresponding symmetry currents, and discuss how the symmetry, and global symmetries, are manifested when Standard Model states are integrated out. Including contributions from higher dimensional operators to the equations of motion modifies the interpretation of conserved currents. We discuss the corrections to the electromagnetic current as an example, showing how modifications to the equation of motion, and corresponding surface terms, have a direct interpretation in terms of multipole charge distributions that act to source gauge fields.
pacs:
Introduction. Assuming physics beyond the Standard Model (SM) at scales , the embedding of the discovered “Higgs-like” scalar into an scalar doublet (), and the absence of hidden states with couplings to the SM and masses , the SM can be extended into the Standard Model Effective Field Theory (SMEFT). Current LHC results are consistent with interpreting data in this framework, where an infinite tower of higher dimensional operators is added to the SM. The lack of any direct discovery of new physics resonances indicating beyond the SM states with masses also supports the assumption that . As a result, the SMEFT expansion in terms of local contact operators is a useful and predictive formalism to employ studying measurements with characteristic scales .
The SMEFT has the same field content as the SM, and reduces to the later by taking . As the SM is falsified due to the evidence of neutrino masses from neutrino oscillations, we assume that neutrino masses are generated by the dimension five SMEFT operator.
The LHC is providing large amounts of data measured around the scale to search indirectly for physics beyond the SM. These efforts are important to combine with experimental measurements at scales , where the Low-Energy Effective Field Theory (LEFT) is the appropriate EFT description.111The notation indicates that this expectation value includes the effects of possible higher dimensional operators. The LEFT is built out of the field content of the SM, but as the Higgs, , , and top have masses , these states are integrated out in sequence. The gauged and linearly realized symmetries of the LEFT are and . To perform EFT studies that combine data sets at scales and , one matches the SMEFT onto the LEFT, and uses renormalization group evolution to run between the different scales. For recent results to this end, see Refs. Jenkins et al. (2018a, b).
When considering matching onto the LEFT at sub-leading order, it is usually necessary to take into account corrections to the equations of motion (EOM) that occur due to the local contact operators present in this theory. In Ref. Barzinji et al. (2018), such corrections for the SMEFT were determined. In this paper, we determine these corrections for the LEFT up to operators of mass dimension six.
The pattern of local operator corrections to the EOM encodes a (non-manifest) symmetry, when this symmetry is assumed to be present in the UV completion of the LEFT. In this paper, we also construct the corresponding symmetry currents and explain the way that the SM gauge symmetries, and global symmetries such as lepton number, are encoded in the LEFT.
Modifying the equations of motion of SM fields by higher dimensional operators challenges the standard interpretation of conserved currents which is appropriate for, and limited to, renormalizable theories. The generalized currents encode symmetry constraints that still constrain an EFT. We also discuss how higher dimensional operator corrections to the equation of motion have a direct interpretation in terms of multipole charge distributions that act to source the corresponding gauge fields. We use the electromagnetic current as an example of this phenomena, and redefine the source in Gauss’s law.
Effective field theory taxonomy. This paper is concerned with the connection between three effective theories: the Standard Model, the SMEFT and the LEFT. Our SM notation is defined in Ref. Barzinji et al. (2018). The SMEFT extends the SM with higher dimensional operators of mass dimension ,
[TABLE]
The operators are suppressed by powers of the cut-off scale and the are the Wilson coefficients. The are constructed out of all of the SM fields and the mass dimension label on the operators is suppressed. We use the non-redundant Warsaw basis Grzadkowski et al. (2010) for , which removed some redundancies in the result reported in Ref. Buchmüller and Wyler (1986). (See also Refs. Aguilar-Saavedra (2011); Alonso et al. (2014).)
The LEFT is given by
[TABLE]
where
[TABLE]
The dual fields are defined with the convention with . The dimension four mass terms are
[TABLE]
is the field strength of . Here labels the fermion fields. In the chiral basis for the we use, charge conjugation is given by . This is not to be confused with a Wilson coefficient . As chiral projection and charge conjugation do not commute, we fix notation . has been rescaled by and has mass dimension zero.
The are constructed out of the SM fields except the Higgs, , and the chiral top fields . The dimensionfull cut off scale of the operators has been chosen to be in the LEFT. The relative couplings required to transform this scale into the mass of a particle integrated out (or a numerical factor in the case of ) are absorbed by the Wilson coefficients.
Equations of Motion. The SM, the SMEFT and the LEFT are all consistent field theories defined by actions
[TABLE]
Each theory contains field variables, here generically denoted . The meaning of the field variables, even those with the same notational label, differs in these theories. A field is redefined order by order in an EFT power counting expansion to remove redundancies of description out of the Lagrangian. As a result, the extremum of the action under variations of field configurations,
[TABLE]
is also redefined order by order. The descendent EOM for then depend on the local contact operators that are present in the EFT expansion. Asymptotic states can be considered to be free field solutions to the modified EOM. The corrections to the EOM modify matching to sub-leading order onto an EFT Elgaard-Clausen and Trott (2017); Jenkins et al. (2018b); Barzinji et al. (2018), and modify the sources of gauge fields. Obviously, one must be careful to include all effects when dealing with higher orders in the power counting expansion.
For the LEFT the gauge fields have the expanded EOM
[TABLE]
Here we have used the adjoint derivative with definition
[TABLE]
For the fermions, the EOM take the form
[TABLE]
Each up to is given in the Appendix.
Symmetry Currents. A continuous transformation of a field,
[TABLE]
under a deformation , with an associated infinitesimal parameter , is a symmetry of if is invariant under this transformation, up to the possible generation of a surface term. The EOM defined by the variations of field configurations in the action –– is unchanged by this transformation. The EOM are defined with surface terms neglected, and the surface terms themselves are defined to be those of the form
[TABLE]
generated by . The Lagrangian is then invariant under , up to a possible total derivative
[TABLE]
for some . Associated with each symmetry defined in this manner is a conserved current Noether (1918). The definition of the current is
[TABLE]
The conservation of the current corresponds to
[TABLE]
Due to the presence of an EFT power counting expansion, it is interesting to examine how symmetry currents are defined when non-renormalizable operators are included, and how these currents encode symmetry constraints.
Basis dependence. The symmetry currents are basis dependent in an EFT, but still meaningful. They receive corrections due to the local contact operators in a particular basis through the modification of the EOM. The basis dependence of the symmetry currents can be made clear by considering a space-time symmetry. For an infinitesimal translation of this form
[TABLE]
up to . Comparing to Eqn. (15) identifies . Four separately conserved currents result, identified as the stress-energy tensor, given by
[TABLE]
The become basis dependent when redundant operators are removed from the EFT, leading to the chosen basis of operators for . The constructed from is also basis dependent as a result at the same order in the power counting. This should be unsurprising, as the currents are auxiliary operators, and sources and the related Green’s functions are not invariant under field redefinitions. For more detailed discussion on this point, see Refs. Politzer (1980); Brivio and Trott (2017). This basis dependence is similar to scheme dependence. It vanishes in relationships between a set of physical measured quantities (i.e. -matrix elements constructed with an LSZ procedure) defined via the same stress-energy tensor. Symmetry constraints between -matrix elements are basis independent, even though the symmetry current itself carries basis dependence.
Non-linear global symmetries. The effect of non-linear representations of the symmetries of the LEFT is straightforward in some cases. As a simple example, consider transforming the charged lepton fields as
[TABLE]
by some global phase . By inspection of the LEFT operator basis, the operators all respect this transformation, except . The charged lepton current is
[TABLE]
The kinetic terms are taken to a flavour diagonal form
[TABLE]
using the flavour space rotation matrix . In the remainder of the paper, the prime superscript is suppressed. descends from the kinetic terms and is also flavour diagonal after these rotations. can receive contributions from higher dimensional operators in a basis, as indicated by the ellipsis in the above expression. The LEFT basis of Refs. Jenkins et al. (2018a, b) removes derivative operators systematically so there are no contributions of this form due to the defined in these works. The divergence of the current including the EOM corrections is
[TABLE]
and similarly for . The mass terms are invariant under Eqn. (20) and cancel when the expressions are summed. We split the EOM correction and into lepton number conserving and violating parts, and . First, consider the lepton number conserving part of Eqn. (23). A significant degree of cancellation occurs in the resulting expression. The only Wilson coefficient remaining corresponds to , an operator which is not individually invariant under the charged lepton field transformation. The explicit expression is
[TABLE]
Similarly, we can define a neutrino current
[TABLE]
The lepton number conserving contributions to the divergence of the neutrino current are such that
[TABLE]
This is as expected, and provides a cross check of the EOM corrections in the Appendix. The total lepton field current is conserved by the subset of operators leading to
[TABLE]
where is the doublet field. Considering the transformation of only part of the lepton multiplet under a phase change also illustrates how a symmetry can be present in a Lagrangian, but non-linearly realized. The symmetry constraint is only made manifest when all terms corresponding to the linear symmetry multiplet are simultaneously included in the constructed symmetry current. This re-emphasizes the requirement to use a consistent LEFT with all operators retained when studying the data. Doing so ensures that the LEFT represents a consistent IR limit. Conversely, dropping operators can forbid non-linear realizations in the LEFT of UV symmetries, which can block a consistent IR limit of some UV completions being defined. For this reason (see also Ref. Jiang and Trott (2017)), experimental studies of constraints on higher dimensional operators done “one at a time” can result in misleading conclusions.
Linear representations of global symmetries. Operator dimension in the SMEFT is even (odd) if is even (odd) de Gouvea et al. (2014); Kobach (2016). Here and are respectively the baryon and lepton number violation of the operator considered. In , is an accidentally conserved quantity consistent with this constraint.
In the LEFT, incomplete SM multiplets are used to construct operators, and operators are not constructed to respect hypercharge. The relationship between operator dimension and global lepton and baryon number in the LEFT is different than in the SMEFT as a result. When considering arbitrary Wilson coefficients in the LEFT, the classes of , , and defined in Refs. Jenkins et al. (2018a, b) are present. These operators are not present in in the SMEFT, and these operators violate .
The SMEFT relationship between operator dimension and these global symmetries is projected onto the LEFT operator basis when the matching result of Ref. Jenkins et al. (2018a) is imposed. The corresponding - matchings that violate vanish exactly.
Hypercharge. The fermion hypercharge current of the SM is
[TABLE]
where . This current is manifestly not conserved in the LEFT
[TABLE]
In the LEFT, a hypercharge current can be defined as
[TABLE]
Here and the hypercharges are assigned as in the SM. Part of the non-conservation of the current stems from the fermion mass terms. In addition, the corrections also lead to the current not being conserved when the Wilson coefficients in the LEFT take arbitrary values. When the matching conditions on the Wilson coefficients to the SMEFT are imposed Jenkins et al. (2018a), many of the EOM corrections generating a non-vanishing are removed. The terms that remain are
[TABLE]
Here we have used the fact that in whole or in part, composite operators forms with have a corresponding vanishing contribution to the current. This condition being fulfilled also provides a cross check of the EOM corrections in the Appendix.
Enforcing matching constraints to the SMEFT is insufficient to make the hypercharge current manifest. The reason is that SM states are integrated out in constructing the LEFT, that carry this quantum number. Consider the definition of the full hypercharge current
[TABLE]
where for the Higgs field. Here, and later, we are using the Hermitian derivative defined by
[TABLE]
for a field . To make hypercharge conservation manifest, we include the transformation properties of the masses associated with states integrated out that depended on . This can be done in a spurion analysis. Rescaled Wilson coefficients and mass terms are promoted to spurion fields with tilde superscripts
[TABLE]
These spurion fields have the hypercharge assignments
[TABLE]
As the spurions are charged under hypercharge, we need to include them in the current in the LEFT
[TABLE]
where
[TABLE]
Here the flavour indices are suppressed. When promoting the Wilson coefficients to fields, we need to include kinetic terms,
[TABLE]
The EOM for the spurion fields are . Including these contributions, the hypercharge current is conserved: .
This provides a cross check of the EOM corrections in the Appendix and the results in Ref. Barzinji et al. (2018); Jenkins et al. (2018a).
** current.** The current in the SMEFT is defined as
[TABLE]
This definition of the current fixes the embedding of the LEFT states into doublets. Here are the generators (Pauli matrices) with normalization for . The fields and are left-handed quark and lepton doublets, which are absent in the LEFT as linear multiplets. To examine the current we need to combine terms in the LEFT into reconstructed multiplets and also introduce spurions to account for the transformation properties of . We illustrate the constraints of the current with an operator from the class as an example,
[TABLE]
where
[TABLE]
We have promoted the Wilson coefficient to a doublet field, and collected the left-handed leptons into a doublet. Analogous promotions can be made for all the operators in this class. The relevant terms in the equations of motion are
[TABLE]
The covariant derivative of gives
[TABLE]
To recover a conserved current, we perform a spurion analysis, similar to the one for hypercharge. We have the EOM for the spurion , in Eqns. (43) and (44). The spurion current is
[TABLE]
with flavour indices suppressed. The covariant divergence of the spurion current is
[TABLE]
Combining Eqns. (Equations of motion, symmetry currents and EFT below the electroweak scale) and (Equations of motion, symmetry currents and EFT below the electroweak scale), the new current is covariantly conserved for the chosen operator from the class ,
[TABLE]
The generalization to include quarks is straightforward.
For operators a similar spurion analysis that also includes the promotion of all of the fermion fields into the corresponding fermion multiplet of the SM is done. The procedure is straightforward. When imposing the - matching and performing this spurion analysis, the current is conserved.
Contraints due to non-manifest currents. The and currents are not conserved in the LEFT when the Wilson coefficients of this theory are treated as free parameters. Furthermore, the implication of these currents in the LEFT is distinct than in the SM or the SMEFT, as there is no manifest field corresponding to these currents when they are conserved. There is no direct construction of a Ward identity using a propagating gauge field as a result.
The conserved currents do constrain the LEFT by fixing relationships between otherwise free parameters of the theory. Matrix elements of the currents can be directly constructed, as they are composed of the fields of the LEFT. Constructing such a matrix element from the generalized Heisenberg current field, with a set of initial and final states denoted , and taking a total derivative gives
[TABLE]
A series of relationships between the Wilson coefficients then follows
[TABLE]
Formally, the measured -matrix elements must be constructed using an LSZ reduction formula. The constraints that follow for the Wilson coefficients are trivially satisfied only if the Wilson coefficients are already fixed by a UV matching preserving the corresponding symmetry.
** and the LEFT multipole expansion.** The classical limit of reproduces the well known physics of Maxwell’s equations, and in particular Gauss’s law Lagrange (1773) (see also Ref. Gauss (1845)). Gauss’s law relates the time component of the electromagnetic current to
[TABLE]
Here , is the electron charge in the usual SI units Mohr et al. (2016). In the LEFT, the electromagnetic current is also expected to be conserved
[TABLE]
without any of the subtleties of the previous sections as the are constructed to manifestly preserve .
The current is subject to its own set of subtleties. First, the naive understanding that being conserved directly leads to its non-renormalization requires some refinement. This issue was comprehensively addressed for QED in Ref. Collins et al. (2006), neglecting higher dimensional operators and considering a one electron state and the corresponding electron number current. Here we review the result of Refs. Collins et al. (2006); Manohar et al. (2017) and then directly extend this result into the LEFT.
The definition of the electromagnetic current is affected by the presence of a surface term Lurie (1968); Collins et al. (2006) introducing a renormalization of this current. We define the CP conserving QED Lagrangian as
[TABLE]
where all superscripted quantities are bare parameters. is introduced so that the renormalized coupling is dimensionless and is the charge of . We restrict our attention to for simplicity (even in loops) in the discussion below. Renormalized quantities are introduced above with a suppressed superscript, and we use as a subtraction scheme so that
[TABLE]
Here . The renormalization constants in QED are given by
[TABLE]
and at one loop. Here , following the notation of Ref. Collins et al. (2006). Hereon we define our subtractions in and suppress the corresponding constant terms, setting .
Standard arguments advanced to establish the non-renormalization of are concerned with Fig. 1 a)-b). Fig. 1 a) represents wavefunction renormalization, while the insertion of the current is represented with a circled cross in Fig. 1 b)-c). The divergence and finite terms of diagrams a)-b) cancel at zero momentum transfer for an on-shell state. For a one electron state, the Noether current corresponds to the electron number current, which we label as consistent with Ref. Collins et al. (2006). The usual textbook argument then concludes
[TABLE]
consistent with the current being conserved. However, the penguin diagram in Fig. 1 c) is divergent. This divergence is cancelled by a counter-term of the form shown in Fig. 1 d). This operator has a four divergence that identically vanishes (i.e. corresponds to a surface term). The EOM of the field is given by
[TABLE]
The EOM relates terms in a non-intuitive fashion when an extremum of the action is taken. receives a multiplicative renormalization generated from the nonzero anomalous dimension of the second term as a result. The current can be subsequently redefined to remove this effect and cancel the running, as shown in Ref. Collins et al. (2006).
Fig. 2 shows the need to further refine this argument in the presence of higher dimensional operators. These diagrams are the direct analogy to the arguments of Ref. Collins et al. (2006) leading to a redefinition of the current due to the mixing of the dipole operator with the counter-term multiplying . Inserting the dipole operator (indicated with a black box) with the electromagnetic current, indicated with a circled cross in Fig. 2a), gives mixing proportional to . Including the effect of the dipole operator in the current insertion is indicated by a “circled cross box” in Fig. 2b). Calculating the diagrams directly for an electron in the loop gives a contribution to the photon two point function of the form
[TABLE]
This divergence is cancelled by a counter-term Jenkins et al. (2018b) which leads to a modification of of the form . (The generalization to other charged leptons in the loop is trivial.) This is as expected as a corresponding divergence is present in the LEFT in Fig. 3 a)-b) and the external photon does not play a role that distinguishes the divergence obtained once the current is redefined.
We have calculated the diagrams in Fig. 3 and agree with the corresponding dipole operator results in Ref. Jenkins et al. (2018b).
The interpretation of this mixing effect is subtle in the LEFT. Varying with respect to gives
[TABLE]
The tree level contributions to the electron number operator of terms vanish at infinity by Stokes’ theorem.222We thank Mark Wise for discussions on this point. We define a -renormalized current
[TABLE]
The -renormalized current expressed in terms of bare quantities is
[TABLE]
The renormalization group flow of the current is
[TABLE]
The -renormalized current depends on the renormalization scale as in the SM case. The LEFT dipole corrections to the current fall off at infinity when considering the electron number operator. They also vanish from Eqn. (59) as separate terms, which is consistent with this fact. The dipole operators mix into proportional to , a correction with a natural interpretation of an electron dipole charge distribution in the LEFT. In order to extract a conserved electron number which is independent of the renormalization scale, we redefine the current, including the effect of dipole operators in direct analogy to Ref. Collins et al. (2006). We define
[TABLE]
where is the electron vacuum polarization in the LEFT, including the effects of operators of mass dimension greater than four. The electron vacuum polarization is still defined in the standard manner, and the current is modified by a redefinition at . It follows that
[TABLE]
In the scheme
[TABLE]
From these results one directly defines the time component of the physical current as
[TABLE]
which is the appropriate generalization of the source in Gauss’s law into the LEFT. This is a numerically small effect, as the electromagentic dipole operator is constrained Dekens et al. (2018).
To summarize, higher dimensional operators in the LEFT act to change the relationship between the Lagrangian parameter and experimental measurements in a manner that corresponds to dipole operators being present in the LEFT. This occurs through a modified source term in Gauss’s law that reflects the presence of a multipole expansion in the EFT. The tree level dipole contributions to the electron number operator vanish at infinity by Stokes’ theorem, but quantum effects necessitates a redefinition of the current.
I Conclusions
We have reported the equations of motion for the LEFT including corrections due to dimension six operators. These results are listed in the Appendix. These corrections lead directly to questions on the meaning of conserved currents in the LEFT. We have examined how the conserved currents of the LEFT encode symmetry constraints that are manifest or non-linearly realized. We have also generalized and embedded the source in Gauss’s law into the LEFT, incorporating the effects of electrically charged particles having dipole operator sources.
II acknowledgments
The authors acknowledge support from the Villum Fonden and the Danish National Research Foundation (DNRF91) through the Discovery center. MT thanks Caltech, San Diego and Perimeter Institute for hospitality when part of this work was completed. We thank Ilaria Brivio, Andy Jackson, Aneesh Manohar, Peter Stoffer, Anagha Vasudevan and Mark Wise for helpful discussions.
Appendix A Appendix
Our operator label notation for the LEFT is largely consistent with Refs. Jenkins et al. (2018a, b). We use a different sign convention on the charge conjugation operator, here , where as in Refs. Jenkins et al. (2018a, b) is defined with opposite sign. We further introduce the current notation
[TABLE]
where etc. We also define the currents where one of the fields is charge conjugated
[TABLE]
and similarly for etc.
Using these notational conventions, the EOM for the gauge fields from are
[TABLE]
The , contributions to the EOM from are as follows
[TABLE]
[TABLE]
The , contributions to the EOM from are
[TABLE]
, contributions to the EOM from are
[TABLE]
Finally, the dimension 3 and 5 LEFT operators contributing to the neutrino EOM give
[TABLE]
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Jenkins et al. (2018 a) E. E. Jenkins, A. V. Manohar, and P. Stoffer, JHEP 03 , 016 (2018 a) , ar Xiv:1709.04486 [hep-ph] . · doi ↗
- 2Jenkins et al. (2018 b) E. E. Jenkins, A. V. Manohar, and P. Stoffer, JHEP 01 , 084 (2018 b) , ar Xiv:1711.05270 [hep-ph] . · doi ↗
- 3Barzinji et al. (2018) A. Barzinji, M. Trott, and A. Vasudevan, (2018), ar Xiv:1806.06354 [hep-ph] .
- 4Grzadkowski et al. (2010) B. Grzadkowski, M. Iskrzynski, M. Misiak, and J. Rosiek, JHEP 1010 , 085 (2010) , ar Xiv:1008.4884 [hep-ph] . · doi ↗
- 5Buchmüller and Wyler (1986) W. Buchmüller and D. Wyler, Nucl.Phys. B 268 , 621 (1986) . · doi ↗
- 6Aguilar-Saavedra (2011) J. A. Aguilar-Saavedra, Nucl. Phys. B 843 , 638 (2011) , [Erratum: Nucl. Phys. B 851 443 (2011)], ar Xiv:1008.3562 [hep-ph] . · doi ↗
- 7Alonso et al. (2014) R. Alonso, H.-M. Chang, E. E. Jenkins, A. V. Manohar, and B. Shotwell, Phys. Lett. B 734 , 302 (2014) , ar Xiv:1405.0486 [hep-ph] . · doi ↗
- 8Elgaard-Clausen and Trott (2017) G. Elgaard-Clausen and M. Trott, JHEP 11 , 088 (2017) , ar Xiv:1703.04415 [hep-ph] . · doi ↗
