On expansions in neutrino effective field theory

We match the seesaw model for generating neutrino masses onto the Standard Model Effective Field Theory (SMEFT). We perform this matching at tree level up to dimension seven in the operator expansion. We explain how some of the perturbations of the neutrino mass matrix due to operators of mass dimension greater than five are tied to integrating out the heavy Majorana mass eigenstates in sequence. We demonstrate that the low energy limit of seesaw models are well described by the SMEFT, particularly when constructed using a flavour space expansion. Flavour space expansions of seesaw models are of interest as the coupling of the heavy states to the Standard Model, that are integrated out to generate neutrino masses, are through flavour space vectors ∈ ℂ3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\mathbb{C}}^3 $$\end{document}. We point out that neutrino phenomenology can be systematically developed as a perturbation around the unknown eigenvectors diagonalizing the charged lepton mass matrix using the fact that these eigenvectors also form a basis of ℂ3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\mathbb{C}}^3 $$\end{document}. This point holds in seesaw models and can also be applied to other models of neutrino mass generation to develop systematic expansions. We develop the algebra for this flavour space and discuss some phenomenology to illustrate this approach.


Introduction
Recently, there has been an escalation of theoretical efforts treating the Standard Model (SM) as a consistent low energy limit of a more fundamental theory. This is a natural result of the discovery of a dominantly J P = 0 + Higgs like boson at the Large Hadron Collider (LHC), and increased experimental indications that there is a mass gap between the electroweak scale (v ∼ 246 GeV) and any scale of new physics. It is reasonable to assume that the SM Lagrangian terms are the leading terms in the Standard Model Effective Field Theory (SMEFT) operator expansion [1][2][3][4][5][6][7][8][9]. Despite the power of the SMEFT formalism and recent systematic developments, there is little direct experimental evidence that higher dimensional operators supplementing the SM have non-vanishing Wilson coefficients. One exception is arguably supplied by the Wilson coefficient of the dimension five operator given by [2,3], (1.1)

JHEP11(2017)088
This operator 1 violates global Lepton number, U(1) L , which is accidentally preserved in the SMEFT operators of mass dimension less than or equal to four [11,12]. For this reason it is generally neglected in LHC studies of the SMEFT. A nonzero value of this operator's Wilson coefficient leads to Majorana neutrino mass terms, which are not present in the minimal SM Lagrangian. Masses for neutrinos are now strongly experimentally supported [13] which makes it appealing to obtain a nonzero Wilson coefficient for this operator. Arguably the simplest way to generate this Wilson coefficient is to directly integrate out heavy singlet fields extending the SM (here denoted N p ) using a seesaw mechanism of neutrino mass generation [14][15][16][17]. A seesaw mechanism for neutrino mass provides an explanation of the smallness of neutrino masses due to a hierarchy of scales. Such an extension of the SM is well described by an effective field theory approach for the same reason.
In this paper, we systematically develop the SMEFT implementation and matching of the seesaw model, integrating out the heavy N p states assuming a renormalizable ultraviolet (UV) extension to the SM. We examine the effect of higher dimensional operators in the SMEFT operator expansion, beyond the Weinberg operator, on the low energy neutrino mass matrix that results. We find by explicit calculation the tree level matching contributions to the SMEFT dimension seven operators.
The seesaw model has been studied many times in the past in an EFT context, see refs. [18][19][20][21][22][23][24][25][26][27][28]. Our results go beyond past work by reporting the complete matching for three generations of heavy singlet fields integrated out in sequence in the seesaw model for the first time up to dimension seven. Simultaneously we incorporate into this implementation of neutrino EFT the flavour space expansion of neutrino phenomenology previously developed in ref. [29]. We discuss how flavour space expansions can be used to relate the neutrino mass spectrum to the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix [30,31] phenomenology. We extend the results of ref. [29] using the fact that one can perturb neutrino phenomenology about the eigenvectors diagonalizing the lepton mass matrix in a general way, simply treating these vectors as an unknown basis of C 3 . We demonstrate the utility of the systematic expansion that can be constructed using this technique with a simple phenomenological example.
The method developed here can be used to study the growing data set on neutrino phenomenology. This can be done in a systematically improvable manner, using well defined expansions, in an effective field theory approach. This formalism is sufficiently general that it can accommodate flavour symmetries assumed in the UV sector, but is not limited to any such flavour symmetry requirement. This approach can also be extended to other UV models of neutrino mass generation in a straightforward manner. 1 The c superscript in eq. (1.1) corresponds to a charge conjugated Dirac four component spinor defined as ψ c = Cψ T with C = −iγ2 γ0 in the chiral basis we employ. The star superscript is reserved for the complex conjugation operation that is applied to scalar and vector quantities. Chiral projection and charge conjugation do not commute (see the discussion in ref. [10] for a review). c L denotes the doublet lepton field that is chirally projected and subsequently charge conjugated.

JHEP11(2017)088
2 Full theory for a minimal seesaw scenario We consider the full theory Lagrangian 2 as given by L = L SM + L Np . To fix our notation we define the SM Lagrangian (L SM ) as Here the fermion fields summed over are where T A are the SU(3) generators, t I = τ I /2 are the SU(2) generators, and y is the U(1) hypercharge generator. Flavour indicies are suppressed in eq. (2.1), restoring the flavour indicies one has for example: where the charge conjugate of N R is N c R . Following ref. [18], we define a field satisfying the Majorana condition N p = N c p in its mass eigenstate basis, with all Majorana phases θ p for each real mass eigenstate shifted into the effective couplings [18], The corresponding Lagrangian is defined as The ω p β = {x β , y β , z β } are each complex vectors in flavour space that have absorbed the Majorana phases. The invariants constructed from these vectors will allow a flavour space expansion as we discuss below. N p is a four component spinor satisfying the Majorana condition, not a two component Weyl spinor. We use greek letters such as β, κ for the label of a flavour vector in the heavy singlet field mass eigenbasis.

Equations of motion of the seesaw theory
We integrate out each N p in sequence, and utilize the Equation of Motion (EOM) to reduce to an operator basis. The EOM include the N p states still present in the spectrum. The relevant modifications of the SM EOM are and Note the last term in the EOM due to varying the fields in the L 5 operator in the SMEFT. 3 Finally the EOM for the N p are The usage of the EOM consistently drops a N p field when it is integrated out of the theory.
Matching onto the SMEFT is defined by requiring that the Wilson coefficients in the higher dimensional operators reproduce the low energy, or infrared (IR), limit of the full theory. For example, consider the IR limit where s 2 m 2 p for N p carrying four momenta s µ as illustrated in figure 1. The N p s-channel propagator is expanded in this limit as Note that we adopt a conventional normalization of the Wilson coefficient of the dimension five operator of the form

L (5) matching
Integrating out the heaviest N p state, denoted N 1 , the matching onto the leading L (5) operators is given by The matrix (x β ) T x κ /m 1 is complex with only one eigenvalue, as only N 1 was integrated out coupled to the complex flavour vector x β . The notation x T β x κ is an outer product of the complex vectors. Integrating out the remaining two lighter N p states in sequence gives where c β κ = (ω p β ) T ω p κ /m p and the flavour index p is summed over. Contracting the SU(2) L indicies of Q 5 and taking a matrix element where the Higgs field is taken as the background field value gives We define the mass eigenstate neutrino fields with prime superscripts. These field are related by the unitary rotation matricies (denoted U ) to the weak eigenstates used so far by Changing to the mass eigenstate basis we find The physical low energy neutrino masses m p ν at leading order in the SMEFT expansion in v/m p are then given as 4 The L (6) matching follows directly and we find (3.10) The operators Q H each with flavour indicies βκ are defined as in ref. [5] with the notation φ exchanged for H for the Higgs field. 5 Here we have reduced the operators to the Warsaw basis form using the EOM and combining terms into Hermitian derivatives defined as We have used the fact that Hermitian operators generate real eigenvalues, and hence the matching coefficient in eq. (3.10) is real. The derivative on the lepton doublet field has been reduced out using the EOM and using the fact thatH † H = 0. For previous results on dimension six matching comparable to the terms in eq. (3.10) see refs. [18,21,22,35]. Our results are distinct from past works in the SU(2) L field dependence. As eq. (2.6) contains the N p fields still in the spectrum when integrating out the heavy Majorana mass eigenstates in sequence, the following terms are also generated. Integrating out N 1 gives Here the notation a · b applied to complex flavour vectors {x, y, z} is a Hermitian inner product, see the appendix for details on the flavour space algebra.
Due to the presence of the Majorana mass scale in the EOM the following contributions to L (3.14) 5 Explicitly these operators are given by Q Dimension seven operators come about due to the expansion of a propagator, such as eq. (3.2), to third order, and from the contraction of the local contact operators present in L (6) N 2,3 once N 1,2 are integrated out in time ordered products. We follow the approach in refs. [1,5] of removing derivative operators in the basis. We define the short hand notation to aid in presenting the resultsC Using the Higgs EOM in eq. (2.6) on the results of the tree level exchange of N p expanded to third order, one finds the terms Q LH and the remaining operator notation for L (7) is defined in ref. [1]. 6 Eq. (3.16) vanishes when the Higgs takes on its background expectation value. This leads to a vanishing of the contributions to the low energy neutrino mass matrix from this sum of terms. Applying the Higgs EOM and reducing the direct matching contributions into field strengths of the SM fields leads to the L (7) operators 7 Here we have used Fierz relations and the EOM to reduce to this basis, utilizing refs. [33,37,38]. It is also important to include the effect of L 5 in determining the EOM for the lepton fields, as this contribution leads to a matching contribution to Q H of the form

JHEP11(2017)088
This contribution perturbs the neturino mass matrix, as we discuss below. The operators in L N 2,3 when N 1 is integrated out are given by In addition, when N 2 is integrated out in sequence the additional matching contributions to the operators involving N 3 are L (7) We have checked the L (5,6,7) matching results with multiple matrix elements to avoid any potential matching ambiguities. We also note the N p mass matrix gets perturbed after integrating out N 1 or N 2 . We have determined these corrections, but as they are dimension eight in the SMEFT they are neglected here.

Perturbation and non-perturbation of the neutrino mass matrix
At tree level if the N p states are integrated out simultaneously or not, the low energy neutrino mass matrix is perturbed due to L (7) matchings. The nature of the perturbations are however reflective of the orientations of the heavy singlet fields in flavour space, as well as their mass spectrum. It is interesting that a number of effects that would perturb the low energy neutrino mass matrix cancel out. For example, the terms in eq. (3.16) cancel in the limit that the Higgs takes on its vacuum expectation value, as previously mentioned. Integration by parts and EOM manipulations can be used to see this result in the complete basis, when considering the matching onto the operator Q H . This operator does lead to a contribution JHEP11(2017)088 to the neutrino mass matrix when the Higgs takes on its vacuum expectation value The dependence on this operator in the expansion of the propagator to third order can be seen to vanish integrated by parts, while also using eq. (2.8). One finds which cancels the corresponding Q H term in eq. (3.17). No additional terms that contribute to the neutrino mass matrix result from the manipulations in the previous equation; these manipulations also cancel the terms in eq. (3.13). Alternatively, one can integrate out N 2,3 using the interactions in eq. (3.13). Doing so, one finds a contribution to L (7) that directly cancels the Q H dependence in eq. (3.17). It is important to include L (6) N 2 ,N 3 and L (7) N 2 ,N 3 when defining the matching onto the theory to sub-leading order for this reason. Use of the EOM, and integration by parts on the N p states still in the spectrum when N 1 is integrated out leads to ambiguities if the full Lagrangian is not specified.
The fact that a subset of contributions to L (7) related to the expansion of the propagator does not lead to a perturbation of the neutrino mass matrix at tree level can also been understood intuitively. To obtain H † H times Q 5 requires two extra insertions of the coupling of the N p states to the SM fields. In eq. (2.5) this coupling is always accompanied by the light SM field so that no local operator is obtained in the heavy N p limit expanding the propagator in Feynman diagrams, as illustrated in figure 2. This viual argument is limited, as this fact is not preserved when reducing the operators obtained in the expansion of the propagator by the EOM. This is another example of the fact that EOM effects in a field theory do not have a trivial Feynman diagram interpretation.
The detailed nature of the neutrino mass matrix perturbations do change if the states are integrated out simultaneously or not, as we discuss below.  (6) and L (4) The limited argument in the previous section also does not forbid perturbations of the neutrino mass matrix due to integrating out the N p states in sequence. Directly expanding out the propagator at tree level to third order, a L (7) matching contribution comes about due to integrating out the heaviest N p mass eigenstate, and subsequently integrating out the lighter N p mass eigenstates. This always occurs as the N p cannot be indistinguishable and generate three distinct eigenvalues of the low energy neutrino mass matrix in the UV scenario we consider. (Different masses of the N p states alone still lead to only one low energy eigenvalue of the neutrino mass matrix, only with a different normalization.) These contribution also match onto L (7) and lead to additional effects perturbing the neutrino mass matrix. The action in the EFT generated when the heaviest N p state is integrated out has a time ordered product contribution of the form

Time ordered products of L
Reproducing the IR limit in the SMEFT give the following matching contributions Note that the terms in the second line are generated when consistently retaining EOM terms (including the N p that remain in the spectrum) to reduce the matching contributions to a minimal basis. Again, for the heavy states to be indistinguishible and integrated out simultaneously, they would have to be oriented in flavour space in the same manner, and have identical masses. Note that these perturbations to the neutrino mass matrix are proportional to these differences in flavour space and the multiple mass scales.

The neutrino mass matrix up to L (7)
The contributions to the neutrino mass matrix in the weak eigenbasis up to L (7) are given by 8 As the mass matrix is perturbed due to corrections at L (7) which are suppressed by O(v 4 /m 3 p ) and of order O(ω 4 ). As such it is established that these corrections can be neglected until perturbations of the Wilson coefficient in L (5) is pushed to relative order O(ω 2 v 2 /m 2 p ) compared to leading effects captured by Q 5 . Radiative corrections to Q 5 are generally larger than the non-pertubative corrections due to L 7 and must be incorporated for phenomenological studies as well if these corrections are to be considered.

Flavour space expansion for the eeesaw
The expansion that results when integrating out the heavy singlet states in sequence is not the only expansion present in lower energy Neutrino phenomenology. The usual matching that was developed in the previous sections leads to small perturbations on the neutrino mass matrix. A larger effect for phenomenology is expanding the Wilson coefficient of the Weinberg operator systematically due to the perturbations of integrating out the N 1,2,3 states. In the remainder of this work, we incorporate and improve on results of ref. [29] to develop perturbations of the U (ν, L) matricies, assuming a seesaw origin of neutrino mass. We use the SMEFT treatment of the seesaw model developed in the previous sections. The idea is to link perturbations of the PMNS matrix to perturbations of the neutrino mass generation mechanism. 9 A key point underlying this approach is Majorana mass terms, unlike Dirac mass terms, originate in bi-linears of the same field operators. As such, the complex mass matrix is diagonalized by a single rotation matrix of the field ν L introduced through ν p L = U (ν, L) pr ν L,r . For this reason any expansion of the neutrino mass matrix is more directly tied to an expansion of the unitary rotation matricies U (ν, L).
With the results of the previous section, the Flavour Space Expansion (FSE) of ref. [29] is now on a firmer theoretical footing. For example, the heaviest neutrino in the low energy theory is generically linked to integrating out the lightest singlet field (denoted in this work N 3 ). In ref. [29] the neutrino mass matrix is generated by first integrating out the lightest JHEP11(2017)088 singlet state, and then integrating out the heavier N p in sequence. Although this can be done, it is conceptually more clear to integrate out the three N p states by removing the heaviest state first, as done here, and subsequently perturb the low energy neutrino mass matrix after the lighter N p states are removed in sequence. Doing so the usual SMEFT expansion is present, clarifying the impact of the multiple expansions present on low energy neutrino phenomenology.

Developing the FSE
The FSE is distinct from the double expansion in (v/m p ) n , and (E/m p ) n that dictates the relative size of contributions in the SMEFT operator expansion. This eigenvector perturbation formalism can always be implemented in a type one seesaw model. However, there is no guarantee that the FSE will be quickly convergent, and therefore predictive, as it depends upon unknown UV physics parameters.
The basic expectation is that in seesaw models the FSE will be perturbative [29]. The reason is that the matrix M pr is expected to be approximately uniform in entries in the N p interaction eigenbasis to the SM states, as the N p do not carry (known SM) quantum numbers. Diagonalizing the corresponding mass matrix, any hierarchy in the couplings of the N R states to the SM states is washed out rotating the the mass eigenbasis, and the magnitude of the ω p are drawn together. For this reason it is expected that on general grounds. Here the notation refers to the Euclidean norm of the complex vectors in flavour space. The FSE is of the form The utility of the FSE depends upon with similar conditions for integrating out the state of mass m 1 . By construction the SMEFT matching has been formulated so that m 3 m 2 < 1. Using the Cauchy-Schwarz equality a · b = ∆ ab a b with ∆ ab < 1 the FSE depends upon Considering eq. (5.1), which directly follows from the quantum numbers of the N p states when tuning is avoided, it is expected that the FSE is present and convergent. In what follows we assume this is the case. Assuming the FSE exists, the results of ref. [29] follow directly, and can be expanded upon in the following way. To establish notation we define

JHEP11(2017)088
The rotation matrix is decomposed in eigenvectors such that U (ν, L) = ( ρ c , ρ b , ρ a ) with ρ i a column vector with ρ i = 1, and ρ i = 1. These eigenvectors are such that We choose the orthonormal eigenvector basis at leading order to be given by [29] 10 With this convention choice, the mass of the heaviest neutrino is given by m a = | z| 2 v 2 /2 m 3 at leading order and without loss of generality. The lighter neutrinos are introduced as perturbations [29]. This can also be done without loss of generality. If the FSE used to introduce these effects is a convergent expansion with small higher order terms depends upon the UV parameters in the seesaw model. The perturbations to the eigenvectors and eigenvalues are given in ref. [29]. We also define the eigenvectors U (ν, which include the perturbations of the eigenvectors to obtain the full complex mass matrix at dimension five in the SMEFT. Note that all eigenvectors in this discussion are normalized to ensure unitarity of the PMNS matrix order by order in the FSE. 11 The PNMS matrix is defined in direct analogy to the CKM matrix as The rotation matrix U (e, L) is introduced to diagonalize the lepton mass matrix Defining the orthonormal (column) eigenvectors of the lepton rotation matrix U † (e, L) as σ i with U † (e, L) = ( σ 1 , σ 2 , σ 3 ) T we have As we are assuming Majorana neutrino masses in a seesaw model, this matrix can be compared to the standard parameterization for unitary matricies. Define with the convention choice c ij = cos θ ij , s ij = sin θ ij and angles θ ij = (0, π/2). Here δ = (0, 2π), α 21 and α 31 are CP violating phases. Each entry in eq. (5.10) is a Hermitian inner product characterized by two parameters, naively leading to eighteen parameters. Comparing to eq. (5.13) which is a general low energy parameterization in terms of six parameters (three moduli angles and three phases) makes clear that there is a redundancy of description in this naive interpretation. However, the eigenvectors sets making up the rotation matricies have to be orthogonal to lead to three masses for the charged leptons and neutrinos. As such the third vector is not independent in its flavour space orientation, although it can carry a relative phase. This leads to nine parameters in each case. Using the relation This expression for ρ a then defines z/ z at leading order in the FSE. Further, without loss of generality z = 1 at leading order. The σ i are a set of orthonormal eigenvectors for the unitary matrix U (e, L). These vectors form a basis for the field C 3 , as they diagonalize M † e M e , a Hermitian positive matrix also defined over the field C 3 . We can expand the unknown complex flavour vector y into this orthonormal basis. Using the orthogonality and normalization properties of the basis vectors of this space, and a general parameterization of these vectors, then allows the use of the systematic EFT expansion, without the rotation matrix U (e, L) being chosen to have a fixed form. We can always define a flavour vector such that with A , B , C ⊂ C. The vectors σ i can be parameterized as discussed in the appendix. These vectors satisfy the complex algebra σ i × σ j = ijk σ k and σ i × σ j = ijk σ k without loss of generality, and we note that the σ i are projectable onto σ i . Solving the general system of equations is straightforward, if tedious. As an example of the utility of this formalism we examine and falsify a simple case. We show that a UV scenario where the second heavy state integrated out couples to the SM as does not satisfy eqs. (5.15)-(5.17) and eqs. (5.7) simultaneously in the limit s 13 → 0. 12 This simple example suffices for our purpose of demonstrating how to perturb in the unknown 12 As s 2 13 0.02 for δm > 0 or δm < 0 [13] the limit considered is experimentally motivated.

JHEP11(2017)088
This introduces ten parameters into the parameterization of this matrix, instead of the usual nine parameters for a unitary 3 × 3 matrix. We use the redundancy in one phase introduced to establish the algebra of the eigenvectors we wish to perturb in. From this general parameterization we have One can then directly determine the complex algebra σ i × σ j = ijk σ k is present when the phase convention choice α 1 + α 2 + α 3 + β 1 + β 2 + β 3 = 2 π n , n ⊂ Z is made. A phase choice of this form is allowed, and reduces the number of free parameters in the parameterization of the unitary matrix U (e, L) to nine. It follows directly that σ i × σ j = ijk σ k in general. It is also required to know the projection coefficients of the σ i onto the basis of vectors σ i to perform the eigenvector perturbations in a general way. They can be derived directly using the definition of the Hermitian inner product and recalling σ i · σ j = ( σ j · σ i ) .
To simplify the intermediate steps of the calculation involving σ i and σ i it can be convenient to re-phase the charged lepton field to make the eigenvalues of M † e M e positive and remove the α i from eq. (A.3) without physical effect. 13 The β i phases in eq. One can also choose a parameterization of U (e, L) where these β i intermediate unphysical phases vanish. This similarity transformation leaves the eigenvalues of M † e M e invariant but does not leave the eigenvectors invariant in general. Choosing this phase convention fixes a general class of σ i to perturb around as a basis for C 3 . As physically observable effects due to U PMNS only come about due to the relationship between the eigenvectors sets σ i and ρ i this can be done as a convention choice.

JHEP11(2017)088
B Operator basis of ref. [ Table 1. The operator basis of ref. [1] matched onto in this work. Here the spinors are in four component notation and C = −iγ 2 γ 0 in the chiral basis we employ.
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.