The Standard Model of particle physics with Diracian neutrino sector

The minimally extended standard model of particle physics contains three right handed or sterile neutrinos, coupled to the active ones by a Dirac mass matrix and mutually by a Majorana mass matrix. In the pseudo-Dirac case, the Majorana terms are small and maximal mixing of active and sterile states occurs, which is generally excluded for solar neutrinos. In a"Diracian"limit, the physical masses become pairwise degenerate and the neutrinos attain a Dirac signature. Members of a pair do not oscillate mutually so that their mixing can be undone, and the standard neutrino model follows as a limit. While two Majorana phases become physical Dirac phases and three extra mass parameters occur, a better description of data is offered. Oscillation problems are worked out in vacuum and in matter. With lepton number -1 assigned to the sterile neutrinos, the model still violates lepton number conservation and allows very feeble neutrinoless double beta decay. It supports a sterile neutrino interpretation of Earth-traversing ultra high energy events detected by ANITA.


Introduction
Thus far the Large Hadron Collider (LHC) has not produced evidence for physics beyond the standard model (BSM).But the neutrino sector must involve BSM because neutrinos have mass.Indeed, the 2015 Noble prize in physics was awarded to T. Kajita and A. B. McDonald "for the discovery of neutrino oscillations which show that neutrinos have mass" [1].
The standard neutrino model (SνM) with its three Majorana neutrinos has measured values for the mass-squared differences, the mixing angles θ 12 , θ 23 and θ 13 and the weak Dirac phase δ.But the absolute mass scale, the order of the hierarchy, normal or inverted, and the Majorana phases are unknown.There is stress in the fit to the standard solar model [2]; there is a reactor neutrino anomaly [3,4]; MiniBooNE finds 4.5σ evidence for a sterile neutrino [5], while MINOS/MINOS+ does not [6].At present, there is no definite conclusion about the existence of an eV sterile neutrino [7].
There is also input from cosmology.From the lensing of background galaxies by the large, reasonably relaxed galaxy clusters Abell 1689 [8][9][10] and Abell 1835 [11] there is indication for 3 active and 3 sterile neutrinos with common mass of 1.5-1.9eV, which act as the cluster dark matter.We shall not dwell here into the many questions raised by and counter-evidence to that possibility, but refer to the discussion and cited articles in these references.Be it as it may, the 3+3 case puts forward to consider the minimal extension of the standard model (SM) in the neutrino sector.By default, this accepts all SM physics without extension in the Higgs, gauge, quark and charged lepton sectors.Gauge invariance then forbids the presence of a 'left handed' Majorana mass matrix between the left handed active neutrinos, so that there must be a Dirac mass matrix to give them mass.As such a term mixes left and right handed fields, this presupposes the existence of 3 right handed neutrinos, also called sterile, i. e., not involved in elementary particle processes [12].For that reason, they are allowed to have a mutual 'right handed' Majorana mass matrix.In order to make up for half of the cluster dark matter, sterile neutrinos have to be generated in the early cosmos by oscillation of active ones.This is only possible when the Dirac mass matrix is accompanied by a non-trivial right handed Majorana mass matrix.
In the pseudo-Dirac limit, the right handed Majorana masses are much smaller than the eigenvalues of the Dirac mass matrix.The maximal mixing of the resulting pseudo-Dirac neutrinos implies that in principle half of the emitted solar neutrinos has become sterile here on Earth, and thus unobservable 1 ; this is ruled out by the standard solar model [2].Hence the pseudo-Dirac case is often considered to be ruled out.We intend to show, however, that there is a way out of this conundrum, so as to faithfully include neutrino mass in the SM without changing its high energy sector.
While excellent studies such as [12][13][14] discuss the theory for general number N s of sterile neutrinos, we shall work out the case N s = 3 in a nontrivial limit where the 6 Majorana neutrinos combine into 3 Dirac neutrinos so that the maximal mixing is harmless and can be circumvented.We call them Diracian neutrinos, i. e., Dirac neutrinos in a model with both Dirac and Majorana masses.In section 2 we treat the theory and in section 3 we consider various applications.We close with a summary.

The Lagrangian for active plus sterile neutrinos
In this section we concentrate on the neutrino sector of the SM.For completeness we present the full Lagrangian in Appendix B.

Active neutrinos only
1 See section 2.4 for details.
We start from the SM Lagrangian where the e, µ and τ fields are diagonal in the mass basis.Left handed neutrinos and right handed antineutrinos exist are called "active neutrinos" since they participate in the weak interactions 2 .Additional neutrinos are not involved in them and called sterile.If only active ones exist, they are Majorana particles.Their mass term involves the quantized left handed fermionic flavor fields ν eL , ν µL , ν τ L , where C is the charge conjugation matrix, T denotes transposition, † Hermitian conjugation, and h.c.Hermitian conjugated terms.M M L is called the left handed Majorana mass matrix.In the SM gauge invariance forces M M L to vanish [12]; if it is present, it must originate from high energy BSM, such as Weinberg's dimension-5 operator.Considering new physics only in the neutrino sector, we neglect M M L .

The Dirac and Majorana mass matrices
In absence of M M L , the only possibility to give mass to the active neutrinos is by a Dirac mass matrix.Since that involves products of left and right handed fields, this presupposes the existence of N s ≥ 3 sterile neutrinos, that must be right handed and represented by quantized fermionic fields ν iR , where the Dirac mass matrix M D is a complex 3 × N s matrix.The sterile fields do not enter the weak interactions; they are singlets under the U(1) Y ×SU(2) L × SU(3) C gauge groups of the SM and affect neither gauge invariance, anomalies nor renormalization.Hence they preserve its full functioning while accounting for neutrino masses.Moreover, the sterile fields may have a mutual mass term like Eq. (1), where the right handed Majorana mass matrix M M R is symmetric and complex valued, and where ν c iR is the charge conjugate of ν iR , While ν iR is a right handed field, ν c iR is left handed (see Appendix A for properties of γ and C matrices).The kinetic term has a common form for all species [12], Left and right handedness refers to the chirality; see Appendix A .
where the slash denotes contraction with γ matrices, and the partial derivatives acting as

The general mass matrix for 3 sterile neutrinos
Though the number of right handed neutrinos is not fixed in principle, the case N s = 3 has, if not a practical value [8][9][10][11], at least an esthetic one: for each left handed neutrino there is a right handed one, in the way it occurs for charged leptons and quarks.The three families of active left and sterile right handed neutrinos have the flavor 3 vectors3 With the combined left handed flavor vector the above mass Lagrangians combine into In general, the mass matrix consists of four 3 × 3 blocks, As stated, we take M M L = 0.In the ("standard", "pure" or "trivial") Dirac limit also M M R = 0.For pseudo-Dirac neutrinos M M R will be small with respect to M D , or, more precisely, small with respect to the variation in the eigenvalues of M D .Though we consider general M D , we are inspired by the case of galaxy cluster lensing where it has nearly equal eigenvalues with central value 1.5-1.9eV [8][9][10][11]; in section 3.1 we shall show that the entries of M M R then typically lie well below 1 meV.

Intermezzo: One neutrino family
In case of one family the flavor vector is The entries of Eq. ( 10) are scalars, so Its eigenvalues are The physical masses are their absolute values [12].For μ nonnegative, this leads to The corresponding eigenvectors are For small μ these are 45 • rotations, i. e., maximal mixing of the active and sterile basis vectors.Formally we may undo the rotations over 45 • , by considering where the approximations are to first order in μ, the pseudo Dirac regime.The first vector, e ã, has its main weight on the first component, so it is mainly active, which we indicate by the tilde on a.The second one, e s, is mainly sterile.But unless μ = 0, the masses m 1,2 are different, so that e ã and e s are not eigenvectors and have no physical meaning.In fact, the mass squares have the difference An initially active state, with momentum p will at time t have oscillated into where E 1,2 = p 2 + m 2 1,2 .The occurrence probability is where for p ≫ m In practice there will not be a pure initial state but some wave packet [12].For t ≫ /∆E the cosine in Eq. ( 19) will average out, so that the fraction of observable neutrinos is approximately 1 2 .In plain terms: for t large enough, half of the neutrinos are sterile and thus unobservable.For the solar neutrino problem the one-family approximation happens to work quite well [15] and the detection rates are well established.Hence for the pseudo Dirac model it would mean that twice as many neutrinos should be emitted as in the standard solar model.The corresponding doubling of heat generated by nuclear reactions is ruled out by the measurements of the solar luminosity, so the case is rarely discussed.
Only in the pure Dirac case, i. e., with Majorana mass μ = 0, the oscillations will not take place, since m 1,2 = m and ∆E = 0.When starting from an initial active state ν a (0), it now equals e ã, and this can be taken as eigenstate.The sterile state will merely be a spectator, "just sitting there and wasting its time".This can be generalized to three families.If one would follow the Franciscan William of Ockham (Occam's razor), it would be preferable for active neutrinos to be Majorana rather than Dirac with unobservable right handed partners.
The SνM differs from the neutrino sector in the SM by accounting for finite masses of its 3 Majorana neutrinos.Below we discuss a "Diracian" setup in which the sterile fields become physical, namely partly active, and the active fields partly sterile, even though the mass eigenstates have Dirac signature in vacuum.

Diagonalization of the Dirac mass matrix
We return to the 3 family case and its total mass matrix Eq. ( 10) with M M L = 0. We notice that any 3 × 3 unitary matrix U can be decomposed as a product of 5 standard ones, The diagonal matrix D M is called the Majorana phase matrix.Likewise we denote the diagonal phase matrix 4 D ′ by diag(e iη ′ 1 , e iη ′ 2 , e iη ′ 3 ).The matrix U D is the product of where c i = cos θ i , s i = sin θ i where the angles θ i are termed in standard notation θ 1 = θ 23 , θ 2 = θ 13 and θ 3 = θ 12 .The Dirac phase δ is also called weak CP violation phase.The complex valued Dirac mass matrix M D can be diagonalized by two unitary matrices of the form (21), viz.
The result reads with the real positive mi 5,6 .We identify U D L with the PMNS mixing matrix U D = U 1 U 2 U 3 and D M L with the Majorana matrix D M employed in literature.4 For U in Eq. ( 21) only 5 of the η i and η ′ i are needed; this can be seen by factoring out e iη1 from D M and e iη ′ 1 from D ′ and setting η 1 → η 1 − η ′ 1 .Both sides of Eq. ( 21) thus involve 9 free parameters.

5
We denote Dirac mass eigenvalues by mi to distinguish them from the physical masses m i , the eigenvalues in absolute value of the total mass matrix. 6 While the left hand side of Eq. ( 23) has 9 complex or 18 real parameters, the right hand side has 9 + 3 + 9; but since M d is diagonal, the diagonal matrices and D M * R only act as a product.Hence it is allowed to fix D M R before solving D M L , see below Eq. ( 28).The number of parameters available for the diagonalization is then still 18.
To connect the transformation (23) to M DM , we introduce the 6 × 6 unitary matrix and define, using that New active and sterile fields express (8) as With these steps the right handed Majorana mass matrix transforms into 21), one can, however, use the phases in D M R to make the off-diagonal elements of M N real and nonnegative 6,7 .We denote the diagonal elements of the Majorana matrix M N by μi , that may still be complex, and the real positive off-diagonal elements by µ i .The right handed Majorana mass matrix M N then takes the form so that the total mass matrix M reads Except in the pure Dirac limit where µ i = μi = 0, the n L are not rotations of mass eigenstates. 7 Actually, for n lepton families there are 1 2 n(n − 1) independent complex valued off-diagonal elements and n Majorana phases, so making all off-diagonal elements real and nonnegative is possible for n = 3 or 2.

Diracian limit
For reasons explained above, we wish to achieve pairwise degeneracies in the masses.The standard Dirac limit, just taking µ i and μi → 0, is a trivial way to achieve this; we shall, however, need finite values for them and design the more subtle "Diracian" limit.
To start, we notice that the eigenvalues of the mass matrix (29) follow from det(M − λI) = 0, where The criterion to get pairwise degeneracies in the eigenvalues (up to signs), is simply that the odd powers in λ vanish.Let us denote and express the μi in a common dimensionless parameter ū through The relations i m2 i ∆i = i ∆i = 0 make the coefficients of λ 5 and λ 1 of Eq. (31) vanish, respectively.To condense further notation, we express the µ i into dimensionless non-negative parameters u i , For normal ordering of the mi (notice that these are Dirac masses, not the physical masses), m1 < m2 < m3 implies ∆1 < 0, ∆2 > 0, ∆3 < 0, hence ∆1 ∆2 ∆3 > 0; this is also the case for the inverted ordering m3 < m1 < m2 whence ∆1 > 0, ∆2 < 0, ∆3 < 0. It thus holds that Equating the λ 3 coefficient of Eq. ( 31) to zero requires This cubic equation has the solutions for n = −1, 0, 1 and positive or negative We restrict ourselves to real solutions; there is always one.Then the matrix M is real valued.All solutions are real when D > 0, which occurs in particular when the u i are small, i.e., in the pseudo-Dirac case.Then there exist the large solutions n = ±1 with ū ≈ ±1, which in both cases leads to the eigenvalues λ ± i ≈ ±M i for i = 1, 2, 3.For small u i the n = 0 solution has a small ū and μi , viz.
The Diracian limit, defined by Eqs. ( 33), ( 34) and (37), reduces Eq. ( 31) to a cubic polynomial in λ 2 , Its analytical roots are intricate, but they are easily calculated numerically.Denoting them as m 2 i , the squares of the physical masses, the eigenvalues of 3 it holds that m 1 m 2 m 3 = m1 m2 m3 .The eigenvectors are set by M jk e and they are real and orthonormal.They can be expressed as with orthonormal e (iã) and e (is) for i = 1, 2, 3.For small µ i and μi the e (iã) and e (is) read to first order 8 e (1ã) = 1, 0, 0, μ1 4 m1 With the first three components of these vectors relating to active neutrinos and the last three to sterile ones, it is seen that for small µ i and μi the e (iã) are mainly active and the e (is) mainly sterile, which we indicate by the tildes.The e (2i−1) and e (2i) are 45 • rotations of the e (iã) and e (is) , which is maximal mixing.In the standard Dirac limit, it is customary to work with Dirac states and not with Majorana states.Likewise, in our Diracian limit the mass degeneracies allow the rotations to be circumvented by working with the e (iã) and e (is) themselves.Indeed, there holds the exact decomposition These steps allow to retrieve the standard Dirac expressions in the limit where the Majorana masses µ i , μi vanish, whence the e (iã) and e (is) become purely active and purely sterile states, respectively.
8 Notice that the i and i + 3 components of e (iã) and e (is) stem with the one-family case Eq. ( 15).
For neutrino oscillation probabilities in vacuum (see subsections 2.4 and 3.2) one needs the eigenvalues m 2 i and eigenvectors e (iã) and e (is) of M 2 , In terms of n L defined above Eq.( 27) and related there to the flavor states (8), the fields for the mass eigenstates are Here j = 1, 2, 3 label active fields and j = 4, 5, 6 sterile ones.Hence the fields ν i ãL annihilate chiral left handed, mainly active neutrinos and create similar right handed antineutrinos, while the ν i sR annihilate chiral right handed, mainly sterile neutrinos and create similar left handed antineutrinos.
The mass term of the 6 Majorana fields now takes the form of 3 Dirac terms, because fermion fields anticommute and left and right handed fields are orthogonal.The here introduced Dirac fields, combine left and right handed chiral fields, as usual.They are the mass eigenstates.In this basis the Dirac-Majorana neutrino Lagrangian is a sum of Dirac terms,

Charged and neutral current
Neutrinos also enter the currents coupled to the W and Z gauge bosons, which are part of the covariant derivatives in the Lagrangian, see Eq. (110) below.The W boson couples to the charged weak current.On the flavor basis it reads with g the weak coupling constant.The neutral weak current reads on the flavor basis with θ w the weak or Weinberg angle.
To express these in the mass eigenstates, we define A and S as matrices consisting of the active components of the 6-component eigenvectors e (iã) and e (is) , respectively, and, likewise, A s and S s for the sterile components From (42) we read off that for small µ i and μi From the orthonormality of the eigenvectors it follows that the real valued 3 × 3 matrices A and S satisfy the unitarity relation while From Eqs. ( 7), ( 8), ( 21), ( 22), ( 24) and (27), and denoting As shown below Eq. (50), the diagonal phase matrix D ′ L can be absorbed in the fields.Inverting Eq. (45) leaves Eq. (48) invariant and expresses the flavor eigenstates as superpositions of mass eigenstates ν mL = (ν ãL , ν c sR ).In vector notation, and using ν c sR = −ν T sR C † , one has Here U DM = U D D M , with U D is the standard PMNS matrix, see (21), while D M = diag(e iη 1 , e iη 2 , e iη 3 ) is the Majorana matrix of the 3-neutrino problem; its η i are Dirac phases now 9 .We also introduced The sterile field ν c sR = (ν 1 c sR , ν 2 c sR , ν 3 c sR ) similar to (56) reads The only current knowledge of the involved matrix elements lies in (54).

9
A word on nomenclature: The Majorana phases in the matrix D M stem from the 3+0 SνM, without sterile neutrinos.While they become physical Dirac phases in the 3+3 DMνSM, there appear no true 3+3 Majorana phases, so we propose to keep this name for them.Hence the DMνSM has 3 physical phases: 1 Dirac and 2 "Majorana" phases.They all appear in the CP-invariance breaking part of the neutrino oscillation probabilities, see Eq. ( 81).
The flavor eigenstates can also be written as single sums over mass eigenstates, with the 3 × 6 PMNS matrix U having elements the latter deriving from (55), while U † U = 1 6×6 , because U represents the 3 active rows of a unitary 6 × 6 matrix which also involves A s and S s .Hence the GIM theorem that J µ N C has the same form on flavor and mass basis, does not hold [12].
Inserted in the currents the relations (56) yield

Lepton number for sterile neutrinos
There is an ambiguity in defining the lepton number of the sterile neutrinos.The lepton number of neutrinos is investigated by making the transformation This leaves the kinetic terms invariant and for the standard choice L s = L a = 1 also the Dirac mass terms (2).Only the Majorana mass terms (1) and (3) will vary by factors e ±2 iφ : they violate lepton number conservation by ∆L = ±2.This approach connects the lepton number L a = 1 of active neutrinos also to sterile neutrinos, hence L νs = −1 for sterile antineutrinos (charge conjugated sterile ones).This assigns lepton number +1 to the components j = 1, 2, 3 of N f L of Eq. ( 8) and n L of Eq. ( 26), but −1 to the components j = 4, 5, 6.Then the mixing (45), or its reverse (56), (58), enforced by the nonvanishing right handed Majorana mass matrix, makes it impossible to consistently connect a lepton number to the particles connected to the mass eigenstates ν i ã and ν i s.
The opposite choice L a = 1, L s = −1 circumvents this problem for general models with active and sterile neutrinos.According to (8), ( 45) and (56) the lepton number L a = 1 of ν a particles is consistent with L ã = 1 of a ν ã particle and L s = −1 of ν s particles.This choice is henceforward consistent with (58).The benefit of this convention is that in pion and neutron decay both channels π − → µ + νãL , π − → µ + ν sR , and n → p + e + νẽ , n → p + e + ν s, respectively, satisfy lepton number conservation.
The minor price to pay is that now both the Dirac mass term and the Majorana terms violate lepton number conservation by two units, so that the unsolvable problem of lepton number violation remains unsolved.Indeed, as we shall discuss below, the Majorana mass terms still allow for neutrinoless double β decay, where a nucleus decays by emitting two electrons (or two positrons) but no (anti)neutrinos.

Estimates for the Dirac and Majorana masses
For later use we present the eigenvalues up to third order in µ 1,2,3 and μ1,2,3 , viz.
With ε ijk the Levi-Civita symbol, there hold the exact relations Let us investigate Eq. ( 64) for normal ordering.The effects of the Majorana masses are anticipated to occur at the level of ∆m 2 sol .We fix m2 and express m1,3 in d 1,3 as so that With m ≈ m2 we set also With d 1,3 and Θ 3 fixed we can determine µ 1,2 or, equivalently, u 1,2 .Imposing ∆m 2 sol ≪ ∆m 2 atm ≪ m2 , we deduce from (64c) and from the difference of (64a) and (64b) that These are expressions of order unity and exact to leading order in d 1 , d 3 − 1 and Θ 3 .Hence the typical scale is µ 1,2 ∼ ∆m 2 atm ∆m 2 sol / m = 4.3 10 −4 eV 2 / m and µ 3 ∼ 8 10 −5 eV 2 / m.In the eigenvectors (42) the coefficients are typically rather small.Hence the mixing matrix S will essentially involve elements ±Θ 1,2,3 with Θ 3 introduced in (68).If one of the Θ i dominates but is still small, it can be seen as a mixing angle.
In particular Θ 3 ∼ 0.1 is possible, which is relevant for the ANITA events to be discussed below.

Neutrino oscillations in vacuum
We consider neutrino oscillations in the plane wave approximation.See Ref. [17] for an excellent discussion of its merits.In the notation (59), (60) an initially pure active state vector reads It evolves after time t into with the Lorentz invariant phase at |r| ≈ ct given by This result can be motivated for a wave packet [17].The amplitude for transition into active state β is where we used that φ 2i−1 = φ 2i .The transition probability after time t may be expressed as two terms, where the first one represents the standard "Dirac" result and where the Majorana masses add the "Majorana" expression The fact that 3 β=1 P M να→ν β ≤ 0 reflects oscillation into sterile states.While the Majorana phases η i cancel in P D να→ν β as usual, they remain present in P M να→ν β .This occurs in the DMνSM because they are upgraded to physical Dirac phases, see footnote 9. P DM να→ν β has the schematic δ-dependence 1 + cos δ + sin δ + cos 2δ + sin 2δ, but the cos 2δ, sin 2δ terms are turned into cos δ, sin δ terms by taking η 3 → η ′ 3 = δ + η 3 .This is equivalent to replacing the U D of ( 21) by wherein δ enters only in the schematic form 1+e iδ .Absolute Majorana phases have no physical meaning; indeed, the expression (78) involves them only through their differences.The CP violation effect, can be read off from the above by switching α ↔ β.The terms with sin(2η i −η j −η k ), cos(2η i −η j −η k ) with j = k and j = k from (78) cancel, leaving the dependence on δ, the η i and t of the form Choosing η 1 = 0 this vanishes only for the trivial values δ, η 2 , η ,3 equal to 0 or π.It confirms that in the DMνSM two of the Majorana phases η i of the SνM are physical Dirac phases 9 .

Neutrino oscillations in matter
Relativistic neutrinos have energy Neutrino oscillations in vacuum are ruled by the Hamiltonian H 0 = E, which reads on the flavor basis The q term leads to qδ αβ and can be omitted, as it plays no role for the eigenfunctions.For propagation in matter one adds the matter potential.The charged and neutral currents induce the scalar potentials involving the electron number density n e = n e − − n e + and the neutron number density n n , and yielding The potential of the active neutrinos is diagonal on the flavor basis, while the sterile ones do not sense any.This results in the total matter potential on the flavor basis From Eq. ( 25) and its real, symmetric nature it follows that Hence the m 2 i , the eigenvalues of M 2 , arise from The total matter Hamiltonian therefore reads on the flavor basis The factor D ′ L in the decomposition (21) for U L , also drops out from V a m since V a m is diagonal, hence it can be totally omitted.Eq. ( 88) can be expressed as First solving the eigenmodes of H m and then going to the flavor basis allows to evaluate the effects of oscillations on the active neutrinos without having knowledge of the undetermined matrix U R .The eigenfunctions do not alter upon subtracting ( m2 /2q) 1 6×6 from H m , after which all elements are small.Inside matter the Diracian properties are lost, there are just 6 Majorana states with different masses.While the neutral current potential V N C can be omitted in the limit M N → 0, this is not allowed in general.In matter one has real potentials V α , α = e, µ, τ .Due to V CC the matrix V a m is complex hermitian.The hermitian H m has 6 different positive eigenvalues but complex valued eigenmodes.
Inside the Sun the neutrino transport is dominated by a mostly electron-neutrino mode [15]; in the 3+3 model this is represented by two nearly degenerate, nearly maximally mixed Majorana modes.The resonance condition in the standard solar model now splits up as a condition for each of them.
The so-called solar abundance problem stems from the inconsistency between the standard solar model parameterized by the best description of the photosphere and the one parameterized to optimize agreement with helioseismic data sensitive to interior composition [18].The biggest deviations in the solar composition are of relative order 1% and occur at ∼ 0.7R ⊙ .Our modified resonance conditions offer hope for an improved description of the data.

Pion decay
One of the simplest elementary particle reactions is It describes a negatively charged π − particle, π − = (dū), consisting of a down quark (charge −e/3) and an anti-up quark (charge −2e/3), decaying into a W − boson (charge −e), which in its turn decays into a muon (charge −e) and an muon-antineutrino (charge 0).Related reactions are π − → e + νe , π + → µ + + ν µ and π + → e + + ν e .In the DMνSM the current 91) by decays with any of the 6 mass eigenstates (ν i m ) R emitted.They can be grouped as The ν i ãL and the charge conjugated ν i c sR fields have identical chiral structure (up to a phase factor, see Ref. [12], eqs (2.139) vs. (2.356)),differing only by their creation and annihilation operators.Hence all decay channels involve the standard chiral factors, and a new factor, the sum over final neutrino states, for α = β = µ.(In this equality we employed Eq. ( 55)).So charged pions decay in the DMνSM at the same rate as in the SM.Neutral pion decay does not involve neutrinos, so it is also not modified.

Neutron decay
A neutron n = (ddu) consists of 2 down quarks and one up quark, and a proton p = (duu) of 1 down and 2 up quarks.Neutron decay n → p + e + νe involves a transition from a down quark to an up quark producing a virtual W − boson, which decays into an electron and an electron antineutrino.As coded in the charged current J µ † CC of (49), it occurs in the DMνSM in two channels, n → p + e + νã and n → p + e + ν s.Both decay channels involve the standard chiral factors, and a new factor, the sum over final neutrino states . This is the α = β = e element of Eq. ( 93), so it is equal to unity.Hence the neutron lifetime in the DMνSM stems with the one in the SνM.
The main decay channel is n → p + e + νã .With our convention L ν s = −1, also the channel n → p + e + ν s conserves the lepton number.The latter occurs at a slower rate due to the small term SS T in (93).We could not convince ourselves that it would be ineffective in beam experiments and hence be capable to explain the neutron decay anomaly between beam and bottle measurements [19].

Muon decay
With the neutron decay going into two channels, muon decay µ → e + νeR + ν µL goes into four, with rates of leading schematic order 1, SS T , SS T and (SS T ) 2 , respectively, adding up to the SM result.

Neutrinoless double β-decay
In a simultaneous double neutron decay (double β-decay) the emission of two electrons involves the schematic neutrino terms ν i 2 ãL + ν i c 2 sR + ν i ãL ν i c sR , corresponding to the emission of two mostly active antineutrinos, two mostly sterile neutrinos, or one of each.With L ν ã = 1 and L ν s = −1, all three channels conserve the lepton number.
In the standard neutrino model also neutrinoless double β-decay is possible.Then only the ν i 2 ãL term occurs, subject to the Majorana condition ν i c ãL = ν i ãL ; with A ij → δ ij and S ij → 0, it yields an amplitude proportional to m ee = ei e 2 iη i m i for small m i .The GERDA search puts a bound |m ee | ≤ 0.15 − 0.33 eV [20].Does this rule out the DMνSM for m ∼ 2 eV? Not, as we show now.
In our situation with Diracian neutrinos neither ν i 2 ãL nor ν i c 2 sR contributes, but neutrinoless double-β decay does arise from the ν ãL ν c sR terms.All spinor terms are again as in the SνM.The only change occurs in the effective mass, which now reads It involves cancellations, since S is nearly asymmetric while A and m/m are close to the identity matrix.But the cancellations are maximal.From the definitions (51) we can go back to the six eigenvectors e (i) of Eqs. ( 40) and (41).Recalling that λ 2i = m i and λ 2i−1 = −m i , (i = 1, 2, 3), it follows that From (30) it is seen that M jk = 0 for j, k = 1, 2, 3, because we neglected the left handed Majorana mass matrix M M L .In general m ee = M DM ee = (M M L ) ee [13].Hence m ee = 0 for this leading order diagram.Nevertheless, neutrinoless double β-decay, involving lepton number violation ∆ L = 2, is not forbidden in the DMνSM.It occurs in the m 3 i /q 2 correction to m i in (95) stemming from the internal propagator m i /(q 2 − m 2 i ) = m i /q 2 + m 3 i /q 4 + • • • .But the suppression factor m 2 i /q 2 makes its measurement impractical for realistic q ∼ MeV − GeV.Loop effects may fare better, but are also tiny.If a finite m ee is established, it points at new high-energy physics.
In conclusion, the non-detection of neutrinoless double β-decay is compatible with the DMνSM.

Small twin-oscillation
If the degeneracy of the solar twin modes is slightly lifted by non-cancellation of the last two terms in each line of Eq. (63), one gets From the standard solar model we know that oscillations should not occur underway to Earth [2], so that |∆ m2 1 | 10 −12 eV 2 .For the supernova SN1987A at distance of 51.4 kpc the absence of twin-oscillations even implies that |∆ m2 1 | 10 −22 eV 2 ; the alternative is that twin-oscillation did take place, and only half of the emitted neutrinos arriving here on Earth were active and could be detected.The implied doubling of power emitted in neutrinos then requires an adjustment of the SN1987A explosion model.

The similarly defined ∆ m2
2 and ∆ m2 3 may be larger.Either of them may describe the MiniBooNE anomaly [5] (disputed by MINOS [6] and still debated [7]) with ∆m 2 ≈ 0.04 eV 2 .But a more elegant approach hereto is to keep the 3 Diracian neutrinos and add a fourth sterile one.

Sterile neutrino creation in the early Universe
The creation of sterile neutrinos in cosmology is an important process based on loss of coherence in oscillation processes.It is well studied, see e.g.[14], and is important when sterile neutrinos are to make up half of the cluster dark matter [8][9][10][11].
The charged currents (61a,b) allow the creation of sterile neutrinos out of active ones via e + + ν e → W + → e + + νs , and creation of sterile antineutrinos out of active ones in the process e + νe → W − → e + ν s , that is, by four-Fermi processes with virtual W ± exchange.The first term in Eq. (61c) describes interaction of active neutrinos with Z, the second of sterile ones, and the last two the exchange of active versus sterile neutrinos, and vice versa.In particular the creation of sterile neutrinos out of active ones is possible in two channels via the four-Fermi process ν α + να → ν s + νs with the exchange of a virtual Z boson.All these processes conserve lepton number.As is seen from the sterile component in the flavor eigenstate (56) or from the charged and neutral currents (61), to achieve the sterile neutrino creation a finite matrix S is needed.Hence it does not occur in the standard Dirac limit where both Majorana mass matrices M M L and M M R vanish.

Muon g − 2 anomaly
The gyromagnetic factor of the muon is g µ = 2(1 + a µ ).Dirac theory yields g µ = 2 and a µ is the anomaly due to quantum effects.The leading term is Schwinger's famous result, a µ is known up to its 9th digit, but there is a ∼ 3.5σ discrepancy between measurement and prediction, a exp µ − a SM µ = 288(63)(43) 10 −11 [16], where the first error is statistical and the second systematic.Our interest lies in the contribution of neutrinos, which occurs in a simple triangle diagram with virtual W bosons. Ref. [21] presents the result for an arbitrary number of sterile neutrinos.For neutrino masses well below M W it reads The unitarity relation (60), viz.(UU † ) αβ = δ αβ , is valid even beyond our N s = 3 case [12].So the DMνSM reproduces the one-loop outcome of the SM, as well as the dominant two-loop electroweak contributions of Ref. [22].

ANITA detection of UHE cosmic neutrino events
Scattering of ultra high energy (UHE) cosmic rays on cosmic microwave background photons puts the GZK limit on their maximal energy [23,24] and acts as a source for EeV (10 18 eV) (anti)neutrinos via the creation and decay of charged pions [25], as considered in section 3.4.
The Antarctic Impulsive Transient Antenna (ANITA) is a balloon experiment at the South Pole that detects the radio pulses emitted when UHE cosmic neutrinos interact with the Antarctic ice sheet.In a set of 30 cosmic ray events, ANITA has discovered an upward going event with energy E ∼ 0.6 EeV [26] and one with ∼ 0.56 EeV [27].Both are consistent with the cascade caused by a τ lepton created beneath the ice surface 10 .But the SM connects a relatively large neutrino-nucleon cross section σ ∼ G 2 F m N E to an UHE neutrino [29], so that the probability for it to traverse a large path L through the Earth to reach Antartica is P T ∼ exp(−n σL) is small, where n is the effective density of nuclei.The numbers are P T ∼ 4 10 −6 and 2 10 −8 , respectively [26,27].In Ref. [30] it is pointed out that a sterile neutrino with a smaller cross section, viz.σ → σ sin 2 Θ, with Θ the mixing angle with respect to active neutrinos, may be involved.To explain the events and relate them to the detections at IceCube, AUGER and Super-Kamiokande, these authors fix Θ at 0.1 [30].
For typical models a sterile neutrino with such a large mixing angle should have been discovered already.In the DMνSM the situation is different, however.It contains the reaction ν s → τ + W + , where the W + is quickly lost locally but the τ escapes and decays while creating a shower.In the approximation where one mixing angle dominates, an emitted electron neutrino will have components of strength sin 2 Θ on mostly sterile states.Inside the Earth they are scattered less, and, given that they enter the other side of the Earth, can be measured in the τ -flavor mode with modified probability PT ∼ sin 2 Θ exp(−n σL sin 2 Θ) and modified flux Φ sterile /Φ active ∼ sin 4 Θ exp(−n σL sin 2 Θ).The estimates of section 3.1 show that the value Θ = 0.1 is reasonable for the component Θ 3 of the mixing matrix S, see in particular the expression for Θ 3 in Eq. (71).Hence the DMνSM supports the sterile-neutrino interpretation of ANITA events.
For determining the DMνSM parameters, it seems worth to include UHE data.

Summary and outlook
Since the maximal mixing of pseudo Dirac neutrinos runs into observational problems, neutrino mass is often supposed to stem from a high-energy sector beyond the standard model (BSM), for instance by the seesaw mechanism [12,14].We show that the mixing effects can be suppressed in the Dirac-Majorana neutrino standard model (DMνSM ), the minimal extension of the SM with 3 sterile neutrinos (3+3 model) with both a Dirac and a right handed Majorana mass matrix.Indeed, to have the 6 physical masses condense in 3 degenerate pairs poses 3 conditions, which leaves 3 Dirac and 3 Majorana masses free.In this Diracian limit the neutrino mass eigenstates act as Dirac particles like the other fermions in the SM.There is no change in the pion, neutron and muon decay, nor on the muon g − 2 problem.Compared to the general case, less mixing occurs since members of the same Dirac pair undergo no mutual oscillation.For small Majorana masses the left handed mass eigenstate is still mostly active, and the right handed one still mostly sterile.A flavor eigenstate has a component on mass eigenstates with a mostly sterile character.With mixing angles up to 0.2 -0.3, this allows to explain the ANITA ultra high energy events.Hence for determining the DMνSM parameters, it is natural to include UHE data. 10See [28] for a possible explanation due to sub-surface reflection in the ice for a downward event.
In the Diracian limit the model keeps some of its Majorana properties.Neutrino oscillations in matter involve the usual 6 nondegenerate Majorana states.Lepton number is not conserved.Neutrinoless double-β decay remains possible, be it at an impractically small rate.Sterile neutrino generation in the early cosmos is possible at temperatures in the few MeV range.
It is interesting to investigate whether processes involving the neutral current can further test the model.They are relevant e.g. in nonresonant sterile neutrino production in the early universe.
By connecting lepton number L = 1 to (mostly) active neutrinos but L = −1 to (mostly) sterile neutrinos, neutron decay and double β-decay conserve the lepton number, while lepton number violation is restricted to feeble neutrinoless double-β decay.Should that be observed, it would invalidate our assumption of negligible left handed Majorana mass matrix, and prove the presence of BSM physics in the high energy sector.
The SM has 19 parameters while 6 neutrino parameters are established and 2 anticipated 11 .The DMνSM adds 3 further Majorana masses.In the limit where they vanish, the sterile partners decouple and the standard neutrino model emerges.The extra Majorana masses and Dirac role of the "Majorana" phases may alleviate some of the tensions in solar, reactor and other neutrino problems.
From a philosophical point of view, we do not consider the values of the Dirac and Majorana masses and phases as problematic properties in urgent need of an explanation, but rather as further mysteries of the standard model.
Up to (110) included, this represents the SM itself: the first line contains the U(1) Y , SU(2) L and SU(3) C gauge fields, respectively, and the Higgs kinetic and potential energy.The second line contains the kinetic terms for three families of quarks, charged leptons and active neutrinos; the third line lists the quark couplings to the Higgs field with 3 × 3 Yukawa matrices Y U,D and the charged lepton couplings to the Higgs field with Yukawa matrix Y ℓ .Eq. (111) exhibits the kinetic term for 3 sterile neutrinos and the Yukawa couplings between the active leptons, the Higgs field and the sterile neutrinos with a 3 × 3 Yukawa matrix Y ν .Eq. ( 111) also contains the right handed Majorana mass terms of Eq. (3).
The Diracian limit of the main text refers a special form for M M R .The quark doublets in Eq. ( 110) contain the left handed up, down, charm, strange, top and bottom quarks, The lepton doublets contain the left handed electron, muon and tau (tau lepton, tauon), and their active neutrinos: the left handed electron, mu and tau neutrino, The right handed quark and lepton singlets are grouped as The covariant derivatives D µ and / D = γ µ D µ in Eq. (110) contain currents from gauge fields, except for / Dν R = / ∂ν R , since ν R is gauge invariant.Indeed, the weak currents J CC from Eq. (49) and J N C from Eq. (50) arise as parts of L L i / DL L .The strong currents from Q L i / DQ L do not involve the neutrino sector.In the unitary gauge the normal and conjugated Higgs doublets read, respectively, where v = −µ 2 /λ is the vacuum expectation value and H the dynamical Higgs field.
After spontaneous symmetry breaking the 3 × 3 Dirac mass matrices for the Down = {d, s, b} and Up = {u, c, t} quarks are M D q = vY D / √ 2 and M U q = vY U / √ 2; for the charged leptons M ℓ = vY ℓ / √ 2 and the Dirac mass matrix for the active neutrinos is M ν = vY ν / √ 2 (it is denoted as M D in the main text).From unitary transformations of the fields it follows that the matrices M D q and M ℓ can be taken diagonal with the respective particle masses as entries.Next, the diagonalization of M U q is performed with the CKM mixing matrix and of M D = M ν with the PMNS mixing matrix of Eq. ( 21).