Preprint
Review

This version is not peer-reviewed.

A Fermion Primer

Submitted:

12 August 2026

Posted:

13 August 2026

You are already at the latest version

Abstract
This is a review paper intended to illustrate a few critical issues concerning Dirac, Weyl and Majorana fermions and their differences. The first part consists in basic introductions to fermions, and more in detail to fermions in 4d, focusing in particular in what differentiate the three type of fermions: chirality, helicity, mass, field equations, properties under discrete symmetry transformations, Actions, Observables. On this basis we tackle a series of challenging and sometime controversial problems: the non-existence of Weyl fermion propagators, and the ways to circumvent it; the regularizations for Weyl fermion amplitudes, in particular the appropriateness of using the PV regularization; the definition of a functional integral for Weyl fermion, the difference between Weyl and massless Majorana fermions and the difference between the Dirac and Majorana mass terms. Most notably we come to the conclusion that the trick of replacing the non-existing Weyl fermion propagator with a massless Dirac one, although it may yield in some cases correct results, is flawed at the very origin by a logical loophole. We then show which is the correct way to proceed in this case. We consider also the topic of applying the Wick rotation at the classical Action level for Weyl fermions and conclude that this procedure leads to a nonequivalent theory. Finally, although anomalies are not the central focus here, we have deemed it useful to summarily review the relation, and its cohomological basis, that exists between the lack of a Weyl fermion propagator (non-invertibility of the Weyl-Dirac operator) and the appearance of dangerous anomalies in the theory.
Keywords: 
;  ;  

1. Introduction

Fermions present surprising and unfamiliar aspects, with respect to standard macroscopic objects. When we take a turn of 360o we see the same world as before turning. For a fermion this is not the case: in our anthropic language we say that it has to turn twice, i.e. 720o, in order to return to the initial situation. A fermion sees a world twice as large as our world. We say that fermions have half-integer spin, at variance with the other set of elementary particles, the bosons, which have integer spin. The consequence of this peculiarity is that the way we represent fermions is drastically different from the way we represent bosons. In order to represent fermions we use mathematical objects called spinors. The latter are called that way because they belong to representations of a spin group (see below), which has precisely to do with the revolving property alluded to above. Spin one-half spinors in 4d can be represented by two or four complex numbers, usually ordered in a column vector. As a consequence the action on them is represented by matrices, the Pauli or gamma matrices. The spinor wave functions or spinor fields are thus represented by columns of two or four functions. There are three distinct types of spinor wave-functions or fields: Weyl, Majorana and Dirac. The equations of motion are matrices of linear differential operator acting on them. The non-standard mathematics of spinors brings into the game some problems which do not seem to have been uniformly elaborated in the literature, or, at least, are resolved in ambiguous ways. This is a review where we discuss such issues and the way we have resolved them. We utilize in part material taken or ri-elaborated from our past research [1,2,3].
A first problem arises concerning the mass term for spinors, which can be vanishing or non-vanishing. Massive spinors are either Dirac or Majorana. They are bispinors with four components. In the Majorana case they can be cast in real form. Massless spinors can be of three types: Weyl spinors with two complex components or chiral bispinors; self-conjugate Majorana bispinors; Dirac bispinors. The Majorana mass term leads inevitably to a self-conjugate Majorana bispinor, for a Weyl spinor has perforce zero mass the two chiral components being decoupled (the mass terms does in fact couple them, as is well known). If a spinor has mass it cannot be Weyl, and viceversa. This is what we try to clarify in the simplest possible way, also to remove the debris generated by the illegal identification of massless Majorana and Weyl fermions.
The second problem arises when one quantizes a theory containing Weyl fermions. Quantization requires the existence of the corresponding propagator. A propagator is the inverse of the kinetic operator, and for Weyl fermions formally it does not exist because the kinetic operator contains a chiral projector. In perturbative field theory one resolves the puzzle by assigning to Weyl fermions the same propagator as for massless Dirac fermions (which does exist). As will be commented upon below this is an arbitrary way of proceeding and one does not have any warranty that the result so obtained makes any sense. In fact experience shows that this trick may work as long as one does not meet critical issues such as anomalies (for Weyl fermions anomalies are the symptoms that a proper propagator does not exist). We will argue below that in such cases this procedure is based on a logical paradox. Nevertheless this may work in perturbative field theory because the fermionic chiral nature is preserved by vertices which contain the chiral projector. An explanation will be provided for the cases in which this trick leads to correct results. In general it has to be assumed that the perturbative approach (Feynman diagrams) applied to Weyl fermion, in which the propagator employed is the Dirac one, is not reliable. For nonperturbative approaches (such as the heat kernel one) the same trick may lead astray because it changes the nature of the problem as will be explained below.
Additional problems and ambiguities arise when one comes to the effective Action of a Weyl fermion theory. One can compute it perturbatively or non-perturbatively. In any case one faces the issue of regularization. In the perturbative case one needs propagators and vertices, which, as we have just seen, becomes a formidable obstacle for Weyl fermions. The tentative remedy is to use the Dirac propagator instead of the missing Weyl one. Once this choice is made one faces the issue of regularization. We have several choices: dimensional [4,5,6] Pauli-Villars [7,8] or cutoff regularizaton - below we consider these three in a simple enough example. Here we make clear a fundamental point: a PV regularized propagator is not the inverse of any kinetic term of whatsoever classical local Action, independently of its being Dirac or Majorana (whose propagators coincide). A somewhat similar problem arises with the technique of Wick rotation: one may be tented to do the Wick rotation in the original classical Minkowski theory and formulate an equivalent problem in a Euclidean spacetime. However we will show that the so obtained Euclidean theory, provided it exists, has very little to do with the original Minkowski one. In simple words, one cannot regularize the original local form of the Action. The regularization process is something separate, which has to do with the mathematical (distribution theoretical) properties of Feynman diagrams. In the case of Weyl fermions this process is particularly delicate, one must make sure that throughout it chirality is preserved. Careful not to be victim of a semantic trap: it is standard lore that in the regularization process a symmetry may be broken, but chirality is not a symmetry, it is the defining status of the fields in the theory. If `chirality is broken’ it means that one, without realizing it, has landed in a foreign theory. For this reason it is of upmost importance how we define the relevant Feynman diagrams or the relevant heat kernel operator. From this point of view PV regularization must be particularly under scrutiny. It needs masses and, as anticipated above, massive fermions can only be Dirac or Majorana, which carry both chiralities. PV can be therefore applied only as a mathematical algorithm to regularize chiral fermion integrals (see an example below), where masses are simple regularizing parameters, without any reference to actual fields and particles. In general a good amount of skepticism is in order when using a perturbative method based on Feynman diagrams for a theory involving Weyl fermions. Better methods are available.
The paper is organized as follows. Section 2 is devoted to a general introduction of fermions as representations of a Clifford algebra and a spin group. This Section can be skipped by those who are not specifically interested in it without prejudice for the comprehension of the subsequent Sections. Section 3 contains a very basic introduction to Weyl, Dirac and Majorana spinors in 4d1. Weyl fermions are the basis for the construction also of Dirac and Majorana ones. Classical solutions, symmetries and canonical quantization, with observables and causal Green’s functions, are explicitly derived and exhibited. In Section 4 a summary is given of Dirac, Majorana and Weyl fermions in 4D with their properties in the language more appropriate for quantum field theory. Section 5 is devoted to explicit examples of regularizations (dimensional, Pauli-Villars and cutoff) for Weyl fermions. The subject of Section 6 are Dirac and Weyl propagators, together with the severe problems carried along by the latter, of which both false and true solutions are discussed. Two more critical issues are also discussed in this Section: the difference between Weyl and massless Majorana fermions, and the Dirac and Majorana mass terms. Section 7 is devoted to Euclidean fermion field theories and, in particular, to the impossibility to write down a real Action for Weyl fermions. Finally in Section 8 the Atiyah-Singer family’s index theorem is presented and its relevance for the invertibility of the Dirac-Weyl operator (=existence of propagator) is highlighted.
Notation. We use a metric g μ ν with mostly - signature. The gamma matrices satisfy { γ μ , γ ν } = 2 g μ ν and
γ μ = γ 0 γ μ γ 0
The generators of the Lorentz group are Σ μ ν = 1 4 [ γ μ , γ ν ] . The charge conjugation operator C is defined to satisfy
γ μ T = C 1 γ μ C , C C * = 1 , C C = 1
The chiral matrix γ 5 = i γ 0 γ 1 γ 2 γ 3 has the properties
γ 5 = γ 5 , ( γ 5 ) 2 = 1 , C 1 γ 5 C = γ 5 T
In the sequel, when more convenient, we shall use a more explicit representation of the gamma matrices, see beginning of section 3.

2. Clifford Algebras and Spinors

At the basis of all the elaborations on spinor fields in quantum field theories is the mathematical concept of Clifford algebra and its representations. Spinors are in fact representations both of a Clifford algebra and of a spin group, which is in turn a subgroup of the Clifford algebra. This Section is a short introduction to this subject and to its connection with spin and spinor fields (for more complete accounts see, for instance,[11,12]). The uninterested reader can skip it without prejudice. The only essential thing to be kept in mind for the sequel is the table at the end.
Given a vector space V with a quadratic form q ( v ) consider the tensor algebra T ( V ) = r = 0 r V . The quotient by the ideal I q ( V ) formed by all the elements of the type v v q ( v ) 1 defines the Clifford algebra
C ( V , q ) = T ( V ) / I q ( V )
It can be defined alternatively as the algebra generated by V (and 1), subject to the relation
v · w + w · v = 2 g ( v , w )
for u , w V , where 2 g ( v , w ) = q ( v + w ) q ( v ) q ( w ) . The dot in (2) is the product
v · w = v w + g ( v , w )
If the quadratic form is the one defined on the basis elements { e i } , i = 1 , n of V by g ( e i , e j ) = η i j , where η i j is the diagonal matrix with r entries equal to 1 and s entries equal to -1, with r + s = n , the relation (2) becomes
e i · e j + e j · e i = 2 η i j
One easily recognizes the form of the γ -matrix algebra (although the definition (2) is more general): in fact the field theory γ -matrices are a representation of the Clifford algebra. When we refer specifically to the canonical form (4), the Clifford algebra will be denoted C r , s . These algebras have been classified in terms of fields and real or complex 2 × 2 matrix rings, M ( 2 , R ) and M ( 2 , C ) , respectively. For instance
C 3 , 0 = M ( 2 , C ) , C 3 , 1 = M ( 2 , R ) M ( 2 , R ) , etc .
So far we have considered real Clifford algebras. Complex Clifford algebra (which are the complexification of real ones) are denoted C C ( n ) . Their classification is rather simple
C C ( 2 k ) = M ( 2 k , C ) , C C ( 2 k + 1 ) = M ( 2 k , C ) M ( 2 k , C )
The automorphism of C ( V , q ) induced by v v ˜ = v for v V , defines a decomposition into even and odd subspaces
C ( V , q ) = C 0 ( V , q ) C 1 ( V , q )
and gives C ( V , q ) a Z 2 -graded algebra structure: C i ( V , q ) · C j ( V , q ) C i + j ( V , q ) , i , j = 0 , 1 .
Inside a Clifford algebra one can define various groups. The largest one is the group of invertible elements
C r , s × = { a C r , s | a 1 }
An important subgroup is the Clifford-Lipschitz group
C r , s = { a C r , s × | a v a 1 V , v V = R r , s }
This definition is based on the adjoint action in V. We can extend it to the full Clifford algebra
σ a ( x ) = a x a 1 , a C r , s , x C r , s
σ a is an automorphism of C r , s . A more instrumental operation is the twisted adjoint
σ ˜ a ( x ) = a ˜ x a 1
where a ˜ has been defined before.
Let us define the group
O ( V , q ) = { t G L ( V ) : t * q = q }
where t * q is the pullback of q by t : V V . I.e. O ( V , q ) is the group of linear invertible transformations of V that leave the quadratic form q unchanged. And let us define the subgroup
S O ( V , q ) = { t O ( V , q ) : det ( t ) = 1 }
As it turns out, we have the correspondences
σ ˜ ( C r , s ) = O ( r , s ) , σ ˜ ( C r , s 0 ) = S O ( r , s )
where C r , s 0 = C r , s C r , s 0 is the even Clifford subalgebra.
Now let us come to the spin groups. We can define a norm N ( a ) for any element a C ( V , q ) , such that N ( v ) = q ( v ) for any v V . Then the spin group of V endowed with a quadratic form q is
S p i n ( V , q ) = { a C r , s 0 | N ( a ) = ± 1 }
The restricted spin group S p i n + ( V , q ) is identified by the condition N ( a ) = 1 .
One can prove that the following sequence is exact
0 Z 2 S p i n ( V , q ) σ ˜ S O ( V , q ) 0
where the twisted adjoint map σ ˜ , for any v , w V , can also be represented as
σ ˜ v ( w ) = w 2 q ( v , w ) q ( v ) v
By the above exact sequence S p i n ( V , q ) is a double covering of S O ( V , q ) .
A K ( K = C or R ) representation of C ( V , q ) is a K -algebra homomorphism
C ( V , q ) H o m K ( W , W )
of C ( V , q ) into the linear transformations of a finite dimensional vector space W. W is called a K -module of C ( V , q ) . Clearly this defines also representations of S p i n ( V , q ) and S O ( V , q ) as groups contained in C ( V , q ) . In particular C ( V , q ) can be itself a representation space for the adjoint or twisted adjoint action defined above. Of course this will be a representation only of C × ( V , q ) or of its subgroups S p i n ( V , q ) and S O ( V , q ) .
We are now ready to define spinors. Among various definitions that can be met in the literature we choose the following one:
  • Given a vector space V = R r , s with canonical quadratic form q and the associated Clifford algebra C ( V , q ) , a spinor is any element of an irreducible representation space of the associated spin group S p i n ( V , q ) .
A classification of spin representations can be found in the literature. Since, in the problems considered in this paper, we are interested in complex representations, we limit ourselves to them. They are as follows
  • n = 2 k , C 2 k 1 C 2 k 1
  • n = 2 k + 1 , C 2 k
( C is the field of complex numbers).
The meaning of the decomposition in two irreducible representations when n = 2 k is related to the volume element. In the complex case the latter is
ω C = i k e 1 e 2 k ,
where e i is an orthonormal basis of V = R r , s . It has the property ω C 2 = 1 , k . The two irreducible representations correspond to the eigenvalues ± 1 of ω C and are identified by the projectors π ± = 1 2 ( 1 ± ω C ) .
The element ω C is evidently an alias of the chirality operator in quantum field theory.

2.1. Spinor Representations in Even Dimension

The properties of spinors in even dimensional spaces are summarized in the following Table 1:
where L , R denotes left-handed or right-handed spinors; the numbers in the columns with C or R tag denote the number of complex or real components according to whether they are in column C or R , respectively. `Self’ means self-conjugate and `Mutual’ means that the L and R representations are conjugate to each other.
After this general and abstract introduction we focus in the sequel on 4d spinors and fields.

3. Spinors in 4 D

In this Section we shall analyze in depth the two basic forms of spinors: namely, the Weyl 2-component spinors, their chiral bispinor representations and the corresponding dynamical properties, as well as the charge self-conjugated Majorana bispinors. As we shall see, the longstanding well-known Dirac bispinors will be obtained as complex combinations of Majorana bispinors and in terms of a local Bogoliubov transformation. For the sake of clarity and completeness, we would like to specify that in the sequel we shall employ the following (more detailed) conventions and notations.
The Pauli spin matrices are the three Hermitean 2 × 2 matrices
σ 1 = 0 1 1 0 σ 2 = 0 i i 0 σ 3 = 1 0 0 1
which satisfy
σ j σ k = δ j k + i ε j k σ ( j , k , = 1 , 2 , 3 )
where ε 123 = 1 = ε 123 so that
[ σ j , σ k ] = 2 i ε j k σ { σ i , σ j } = 2 δ i j
The Dirac matrices in the Weyl, or spinorial, or even chiral representation are given by
γ 0 = 0 1 1 0 γ k = 0 σ k σ k 0 γ 5 = 1 0 0 1
with the Hermitean conjugation property
γ μ = γ 0 γ μ γ 0 { γ μ , γ ν } = 2 g μ ν I
so that
β = γ 0 = γ 0 = 0 1 1 0 α k = γ 0 γ k = σ k 0 0 σ k ( k = 1 , 2 , 3 )
We denote by ψ L the 2-component spinor belonging to the irreducible bidimensional left representation D ( 1 2 , 0 ) of the Lorentz group, while ψ R belongs to its nonequivalent right counterpart D ( 0 , 1 2 ) . The same holds true for the spin-states u L ( p ) and u R ( p ) et cetera. The four-component spinors, or bispinors, ψ ( x ) do belong to the reducible tetradimensional representation D ( 1 2 , 0 ) D ( 0 , 1 2 )
ψ ( x ) = ψ L ( x ) ψ R ( x )
We denote the rank-two chiral projectors with P ± 1 2 ( 1 ± γ 5 ) , in such a manner that for any four-component bispinor ψ ( x )
ψ ( x ) = P ψ ( x ) = ψ L ( x ) 0
ψ + ( x ) = P + ψ ( x ) = 0 ψ R ( x )
and the chiral spin-states u ± ( p )  et cetera are defined in the very same way.

3.1. Weyl Spinor Fields

Let us first analyze the massless and chiral Weyl spinor fields and their specific properties at the classical and quantum level.
Consider a massless Dirac bispinor, i.e., any complex and Grassmann valued bispinor in the Minkowski space, which is a solution of the differential equation
/ ψ ( x ) = 0 ψ ( x ) = ψ L ( x ) ψ R ( x )
In order to obtain the most general solution of the above equation, one has to sensibly move from the conventional approach and formalism that has been developed for the celebrated Dirac equation
i / ψ ( x ) = M ψ ( x )
In the chiral representation for the gamma matrices we can build up a very convenient set of spin-states as follows. Consider the matrix α ν γ 0 γ ν = ( I , γ ) with matrix elements
α ν = σ ν 0 0 σ ˜ ν
where we have set σ ν = ( I , σ k ) with k = 1 , 2 , 3 , while σ ˜ ν = ( I , σ k ) . Then the massless Dirac equation can also be cast in the form
α · ψ ( x ) = ( 0 + α k k ) ψ ( x ) = 0 σ · ψ L ( x ) = 0 σ ¯ · ψ R ( x ) = 0
in such a manner that it can be rewritten à la Schrödinger
i ψ t = H 0 ψ
where the massless Dirac Hamiltonian operator reads, in physical units,
H 0 = c α k ( i k ) = c α k p ^ k
Solving the massless Dirac equation is nothing but finding the general solutions of the pair of decoupled Weyl equations (14).

3.1.1. Spin-States

We can search the plane waves solutions of eq. (14) in the form
ψ ( x ) = Γ ( p ) exp { i p μ x μ } = Γ ( p ) e i p · x p μ = p ˜ μ = ( p 0 , p )
where Γ ( p ) are the so called spin-states associated to the plane wave solutions of the massless Dirac-Weyl equation (14). Now, if we set H 0 ( p ) α k p k = α · p , from the commutation relation [ γ 5 , H 0 ] = 0 we see that we can always select a basis in the 4d complex linear space of the bispinor spin-states according to their chirality χ = 1 and energy (or frequency) p 0 = ± | p | pair of eigenvalues.
Thus, in order to realize a basis, we have to search among states of opposite chirality and frequency. Then, in so doing, we will be able to have at hand a complete and orthogonal quartet of spin-states for the 4d massless bispinor space. The construction goes as follows. Explicit evaluation of the relevant spin-matrices yields for p 0 = | p | =
σ · p = σ μ p μ = I + σ k p k = + σ · p = + p z p x i p y p x + i p y p z
σ ˜ · p = σ μ p ˜ μ = I σ k p k = σ · p = p z p x + i p y p x i p y + p z
γ ν p ν p / = 0 σ ˜ · p σ · p 0 = 0 0 p z p x + i p y 0 0 p x i p y + p z + p z p x i p y 0 0 p x + i p y p z 0 0 γ ν p ˜ ν p ˜ / = 0 σ · p σ ˜ · p 0 = 0 0 + p z p x i p y 0 0 p x + i p y p z p z p x + i p y 0 0 p x i p y + p z 0 0
so that we find
σ · p σ ˜ · p = σ ˜ · p σ · p = p 2 = 2 p 2 = p / 2 = p ˜ / 2 = 0
( σ · p ) 2 = 2 σ · p ( σ ˜ · p ) 2 = 2 σ ˜ · p
Notice that p ˜ ν = ( , p ) = p ν is the dual light-cone tetra-vector which satisfies
p ˜ 2 = p 2 = 0 p · p ˜ = 2 p 2 = 2 2
Moreover it is useful to define the on-shell ( p 0 = ) Weyl-Dirac 2d projectors
Π ( p ) p ˜ / p / 4 2 = 1 2 σ · p 0 0 σ ˜ · p = 1 2 1 α · p = α · p 2
Π ˜ ( p ) p / p ˜ / 4 2 = 1 2 σ ˜ · p 0 0 σ · p = 1 2 1 + α · p = α · p ˜ 2
which fulfill by construction
Π = Π Π 2 = Π Tr Π = 2
Π ˜ = Π ˜ Π ˜ 2 = Π ˜ Tr Π ˜ = 2
Π Π ˜ = 0 = Π ˜ Π [ γ 5 , Π ] = [ Π ˜ , γ 5 ] = 0
γ 0 Π γ 0 = Π ˜ γ 0 Π ˜ γ 0 = Π
It means that Π and Π ˜ are projectors onto 2d orthogonal chirally invariant though parity exchanging spaces. The momentum space Weyl-Dirac equation (14) for the positive energy spin-states evidently reads Π ( p ) Γ ( , p ) = 0 , the general solution of which being Γ ( , p ) = Π ˜ ( p ) ξ , where ξ is any arbitrary constant bispinor. Hence, without loss of generality, we can set, for example, with p ν = ( , p ) = ( , p )
u ( p ) = 1 p z α ν p ˜ ν ξ = u L ( p ) 0 0 = 1 p z p z p x i p y 0 0 ξ = 1 0 0 0 u L ( p ) = σ · p ˜ p z 1 0
The above spin-state is a positive frequency solution of equation (14) and exhibits negative chirality
H 0 ( p ) u ( p ) = u ( p ) γ 5 u ( p ) = u ( p )
Π u ( p ) = 0 Π ˜ u ( p ) = u ( p )
Moreover, the 2-component left-handed Weyl spinor u L ( p ) does fulfill
u L ( p ) u L ( p ) = 2 u L ( p ) u L ( p ) = 2 P ˜ L ( p 0 = )
where the rank-two square matrix
P ˜ L σ ˜ · p 2 = σ · p ˜ 2
does satisfy for p 0 =
σ · p P ˜ L = 0 P ˜ L = P ˜ L P ˜ L 2 = P ˜ L Tr P ˜ L = 1
so that it corresponds to the one dimensional projector on positive energy and left-handed spin-states.
In order to set up a right-handed positive energy spin-state, we can involve the charge conjugated of the above constructed spin-state: namely,
u + ( p ) i γ 2 u * ( p ) = u c ( p )
As a matter of fact, it is easy to check that the chiral bispinor u + ( p ) satisfies the distinguished features
γ 5 u + ( p ) = u + ( p ) H 0 ( p ) u + ( p ) = u + ( p )
Proof. We immediately find
γ 5 u + ( p ) = i γ 5 γ 2 u * ( p ) = i γ 2 [ γ 5 u ( p ) ] * = i γ 2 u * ( p ) = u + ( p )
Furthermore we can write
H 0 ( p ) u + ( p ) = H 0 ( p ) i γ 2 u * ( p ) = H 0 * ( p ) i γ 2 u ( p ) * H 0 * ( p ) γ 2 = ( α 1 p x + α 3 p z α 2 p y ) γ 2 = γ 2 H 0 ( p ) H 0 ( p ) u + ( p ) = i γ 2 u ( p ) * = i γ 2 u * ( p ) = u + ( p )
that proves the above statements (28).
On the other side, it is immediate to realize that a right-handed spin-state with negative energy is provided by
v + ( p ) γ 0 u ( p )
which obviously fulfills
γ 5 v + ( p ) = v + ( p ) H 0 ( p ) v + ( p ) = v + ( p )
its charge conjugated being left-handed, i.e.,
v ( p ) i γ 2 v + * ( p ) = v + c ( p )
and satisfying
γ 5 v ( p ) = v ( p ) H 0 ( p ) v ( p ) = v ( p )
The ultimate meaning of the above construction is that, once we have fixed the initial left-handed spin-state of positive energy, the rest of the spin-state orthogonal basis can be fully determined by the charge conjugation and parity discrete symmetry operations.
It turns out that the above defined quartet of massless and chiral bispinor spin-states does fulfill orthogonality and closure relations. As a matter of fact, spin-states of opposite chirality and/or opposite frequency are orthogonal
u ± ( p ) u ( q ) = 0 = v ± ( p ) v ( q ) u ± ( p ) v ( q ) = 0 u ± ( p ) v ± ( p ) = 0 = v ( p ) u ( p )
Moreover we have
u ± ( p ) u ± ( p ) = v ± ( p ) v ± ( p ) = 2
that exhausts the orthogonality relations and furthermore we get
( 2 ) 1 [ u ± ( p ) u ± ( p ) + v ± ( p ) v ± ( p ) ] = P ±
with of course P ± = 1 2 ( I ± γ 5 ) which represent the closure or completeness relations for the massless and chiral bispinor spin-states basis quartet.
The corresponding plane-wave functions, which are normal-mode solutions of the massless Dirac equation, read
u ± , p ( x ) = [ ( 2 π ) 3 2 ] 1 2 u ± ( p ) e i t + i p · x v ± , p ( x ) = [ ( 2 π ) 3 2 ] 1 2 v ± ( p ) e i t i p · x
and fulfill in turn orthonormality and closure relation with respect to the usual Poincaré invariant inner product: for instance
d x u ı , p ( x ) v j , p ( x ) = 0 d x u ı , p ( x ) u j , p ( x ) = δ ı j δ ( p p ) = d x v ı , p ( x ) v j , p ( x ) ( ı , j = + , p , p R 3 )
et cetera.
We recall that, according to the general principles of quantum field theory, the plane wave functions v ± , p ( x ) of negative energy or frequency are associated to the anti-particles of the corresponding spin 1 2 charged quantum field.
In the Standard Model of the fundamental interactions the neutrino and anti-neutrino particles are supposed to be massless2 and charge conjugated one to each other, the charge being identified with the so called lepton number. Moreover, the statu quo of the present phenomenological frame [10] does require that neutrino particles are always left-handed and with a negative helicity, while anti-neutrinos always appear to right-handed and with a positive helicity. Sic stantibus rebus it turns out to be natural to associate u ( p ) , u , p ( x ) to the left-handed spin-state and normal mode, respectively, of the Standard Model neutrino Weyl field, while v + ( p ) , v + , p ( x ) will refer to the corresponding negative frequency spin-state and normal modes, which will be related to the anti-neutrino particles of opposite chirality and helicity, as far as their masses can eventually be neglected. It is worthwhile to stress once more that this mandatory strong correlation between chirality and helicity does actually occur only for the Weyl spinor field and never for any other kind of relativistic spin 1 2 wave fields.
The remaining pair of chiral spin-states and normal-modes u + , v would correspond to hypothetical right-handed neutrinos and related left-handed anti-neutrino spin-states and normal modes, that have never been observed so far, although their existence does represent, if any, one of the most intriguing and exciting actual possibilities to go beyond the Standard Model.

3.1.2. Symmetries

The classical Lagrange density for the massless Dirac field can be written in different forms: namely,
L 0 = ψ ( x ) i α · ψ ( x ) 1 2 ψ ( x ) i α · ψ ( x )
1 2 ψ L ( x ) i σ · ψ L ( x ) + 1 2 ψ R ( x ) i σ ¯ · ψ R ( x )
up to irrelevant tetra-divergences, where the spinor fields are understood as Grassmann valued functions on the Minkowski space, whereas
ψ ( x ) = ψ L ( x ) ψ R ( x ) α ν = σ ν 0 0 σ ¯ ν
The Euler-Lagrange field equations are
α · ψ ( x ) = 0 σ · ψ L ( x ) = σ ¯ · ψ R ( x ) = 0
the general plane-wave solutions of which have been obtained and discussed in the previous Section. The canonical energy-momentum tensor reads
T μ ν = ψ ( x ) α μ i ν ψ ( x ) = ψ L ( x ) σ μ i ν ψ L ( x ) + ψ R ( x ) σ ¯ μ i ν ψ R ( x )
which is not symmetric and drives to the conserved energy-momentum vector
P 0 = d x ψ ( x ) i 0 ψ ( x ) = d x ψ ( x ) H 0 ψ ( x )
= d x ψ L ( x ) H L ψ L ( x ) + ψ R ( x ) H R ψ R ( x )
P = d x ψ ( x ) ( i ) ψ ( x )
Notice that the energy-momentum tensor is trace-less g μ ν T μ ν = 0 , owing to the validity of the field equations. The canonical total angular momentum density tensor can always be written in the form
M λ μ ν = x μ T λ ν x ν T λ μ + S λ μ ν
S λ μ ν = 1 2 ψ ( x ) γ 0 { γ λ , σ μ ν } ψ ( x )
so that from Noether theorem related to the Lorentz symmetry we get
T ν μ T μ ν = λ S λ μ ν
which notoriously entails the helicity conservation in any direction along which the energy-momentum tensor appears to be symmetric. In order to determine the helicity assignements we can proceed as follows. Since the generators of the rotation group for both the irreducible Weyl spinor representations are S j k = 1 2 ε j k σ it follows that the corresponding Pauli-Lubanski operator reads
W 0 = P · S = p · 1 2 σ
W j = 1 2 σ j 1 2 i ε j k p k σ
which entails
g μ ν W μ W ν = 1 4 p j p k σ j σ k 1 4 σ j i ε j k p k σ σ j i ε j r s p r σ s = 1 2 2 + 1 4 δ k r δ s δ k s δ r p k p r σ σ s = 0
where the 2 × 2 identity matrix has been understood. Then, thanks to the light-like nature of the Pauli-Lubanski operator, we can identify the helicity operator to be
h W 0 = | W | = 1 2 n ^ · σ n ^ = p /
together with the Weyl Hamiltonian operators
H 0 = H L 0 0 H R H R = c p · σ = H L
in such a manner that we find H L = c h , H R = c h and consequently
H L u L ( p ) = c u L ( p ) h u L ( p ) = 1 2 u L ( p ) H L v L ( p ) = c v L ( p ) h v L ( p ) = 1 2 v L ( p ) H R u R ( p ) = c u R ( p ) h u R ( p ) = 1 2 u R ( p ) H R v R ( p ) = c v R ( p ) h v R ( p ) = 1 2 v R ( p )
which means that :
  • neutrino spin-states carry negative helicity, while the anti-neutrino spin-states carry positive helicity, i.e. γ 5 u ( p ) = γ 5 v + ( p ) = 1 .
  • The situation is reversed for the remaining pair of the spin-states. Thus, positive helicity is exhibited by positive frequency particles - the would-be right neutrinos - while negative helicity by negative frequency anti-particles, i.e. γ 5 u + ( p ) = γ 5 v ( p ) = + 1 .
The Weyl Action integrals
S L = d 4 x 1 2 ψ L ( x ) i σ · ψ L ( x ) S R = d 4 x 1 2 ψ R ( x ) i σ ¯ · ψ R ( x )
and the massless Dirac Action integral, which is the sum of the left- and right-handed Weyl Action integrals, are invariant under scale or dilatation symmetry. The latter acts upon any inertial coordinate system and spinor field variables on the Minkowski space according to
x λ x λ = e ϱ x λ ( ϱ R )
Ψ ( x ) Ψ ( x ) = e 3 2 ϱ Ψ ( e ϱ x ) ( Ψ = ψ , ψ L , ψ R )
In order to obtain the infinitesimal operators of the scale transformations we consider an infinitesimal dilatation/contraction, viz.,
x μ = ( 1 δ ϱ + ) x μ ( | δ ϱ | 1 )
Ψ ( x ) = 1 + 3 2 δ ϱ + Ψ ( x + x δ ϱ + ) = 1 + 3 2 δ ϱ Ψ ( x ) + x · Ψ ( x ) δ ϱ + = Ψ ( x ) + 3 2 δ ϱ Ψ ( x ) + x · Ψ ( x ) δ ϱ +
and thereby
δ Ψ ( x ) = δ ϱ 3 2 + x · Ψ ( x ) = 3 2 δ ϱ Ψ ( x ) d Ψ ( x )
because δ x = x δ ϱ . Then we can write the infinitesimal total variation of the spinor fields in the form
Δ Ψ ( x ) = ( δ + d ) Ψ ( x ) = 3 2 δ ϱ Ψ ( x )
in accordance with the well-known general rule that connects local ( δ ) and total ( Δ ) variations of any classical relativistic wave field. Hence, we are allowed to identify
D Ψ 3 2 + x ·
with the generator of the local scale transformations on the classical spinor fields. Now, since Ψ ( x ) = e 3 2 ϱ Ψ ( x ) , it follows that the Lagrangian (37) does satisfy
L 0 ( x ) = e 4 ϱ L 0 ( x )
so that the Action integral for the massless Dirac spinor field appears to be invariant, viz.,
S = d 4 x L 0 ( x ) = d 4 x L 0 ( x )
The dilatation tetra-current can be readily obtained from Noether theorem and can be written in different though equivalent forms, viz.,
D μ ( x ) = δ L 0 δ μ ψ ( x ) D ψ ψ ( x ) L 0 ( x ) x μ = δ L 0 δ μ ψ ( x ) 3 2 ψ ( x ) + x ν T μ ν ( x ) = 3 2 ψ ( x ) i α μ ψ ( x ) + x ν T μ ν ( x )
so that
· D = g μ ν T μ ν = 0
Finally, it turns out that the massless Dirac Action integral turns out to be invariant under the U(1) groups of the ordinary and chiral phase changes,viz.,
ψ ( x ) ψ ( x ) = e i θ ψ ( x ) ψ ( x ) = e i γ 5 θ ψ ( x ) ( 0 θ < 2 π )
or equivalently
ψ L ( x ) ψ L ( x ) = e i θ ψ L ( x ) ψ L ( x ) = e i θ ψ L ( x ) ( 0 θ < 2 π )
ψ R ( x ) ψ R ( x ) = e i θ ψ R ( x ) ( 0 θ < 2 π )
in such a manner that Noether theorem drives to the internal currents
J μ ( x ) = ψ ( x ) α μ ψ ( x ) J 5 μ ( x ) = ψ ( x ) α μ γ 5 ψ ( x )
which satisfy the continuity equation so that the internal charges
Q = d x ψ L ( x ) ψ L ( x ) + ψ R ( x ) ψ R ( x )
Q 5 = d x ψ L ( x ) ψ L ( x ) + ψ R ( x ) ψ R ( x )
are conserved in time Q ˙ = Q ˙ 5 = 0 . To sum up, for a massless Dirac field we have seen until now the occurrence of a 13-parameter Lie group of space-time and internal symmetries leading to thirteen conserved charges. Actually it can be shown, as we shall see later on, that scale invariance implies in fact the invariance under a larger 15-parameter symmetry group called the conformal group, which includes the ten dimensional Poincaré group. Scale invariance implies conformal invariance if and only if there exists in the considered theory a symmetric traceless conserved energy-momentum tensor Θ μ ν , called the improved energy-momentum tensor, such that · D = g μ ν Θ μ ν = 0 . It turns out that this is true for all power-counting renormalizable massless models involving fields of spin 1 . Hence, the massless Dirac Action integral is actually invariant under a 17-parameter Lie group involving space-time as well as internal symmetries.
Concerning discrete symmetries, it turns out that the Action integral S = d 4 x L 0 ( x ) does respect C , P and T transformations separately. The charge conjugation C is an internal discrete symmetry which acts on the classical Dirac spinors according to
ψ ( x ) ψ c ( x ) = γ 2 ψ * ( x )
while parity P transforms yields ψ ( t , x ) = γ 0 ψ ( x ) , whereas time reversal T is provided by ψ ( t , x ) = Θ ψ ( x ) , Θ = γ 3 γ 1 . It is straightforward to verify by direct inspection that the Action integral is invariant. Notice that the above mentioned discrete transforms are not even well defined for chiral spinors, for they mix the chiral components and necessarily drive out of the specific chiral spinor spaces of the two irreducible nonequivalent Weyl representations of the Lorentz group D ( 1 2 , 0 ) and D ( 0 , 1 2 ) .

3.1.3. Quantum Theory

The setting up of a canonical quantum theory for the massless Dirac-Weyl spinor field is far from being straightforward and trivial. First of all it is not at all unique: there are several options available in the market, the best buy to be selected and picked up through the attainment to the general principles of quantum field theory as well as by the correspondence with the actual phenomenology. In order to be appropriate and applicable to the concrete high energy and particle Physics context, the operation of parity and charge conjugation must be actually well defined and truly implementable: this basic requirement clearly rules out a priori the possibility to formulate the quantum theory in terms of the 2-component Weyl spinor field ψ L and ψ R , as the afore mentioned operations are not supported by the 2d Weyl spinors. Hence, the only possible consistent and physically meaningful formulation is in terms of the chiral bispinors ψ .
The massless Dirac quantum field is the operator valued tempered distribution
ψ ( x ) = ψ ( x ) + ψ + ( x )
ψ ( ) ( x ) = d p [ ( 2 π ) 3 2 ] 1 2 c , p u ( p ) + c + , p u + ( p ) e i p · x i t
ψ ( + ) ( x ) = d p [ ( 2 π ) 3 2 ] 1 2 d , p v ( p ) + d + , p v + ( p ) e i t i p · x
ψ ( x ) = P ψ ( x ) = d p c , p u , p ( x ) + d , p v , p ( x )
ψ ( y ) = ψ ( y ) P = d p c , p u , p * ( y ) + d , p v , p * ( y )
which corresponds to the most general solution of the Weyl-Dirac equations (14), in which the creation and destruction operators satisfy the canonical anti-commutation relations
{ c , p , c , p } = { d , p , d , p } = δ ( p p )
all the other anti-commutators being null. It follows that, for example,
{ ψ ( x ) , ψ ( y ) } x 0 = y 0 = d p d p × { c , p , c , p } u , p ( x ) u , p ( y ) + { d , p , d , p } v , p ( x ) v , p ( y ) x 0 = y 0 = d p u , p ( x ) u , p ( y ) + v , p ( x ) v , p ( y ) x 0 = y 0 = d p ( 2 π ) 3 2 u ( p ) u ( p ) e i p · ( x y ) + v ( p ) v ( p ) e i p · ( x y ) = d p ( 2 π ) 3 2 u ( p ) u ( p ) + v ( p ) v ( p ) e i p · ( x y ) = d p ( 2 π ) 3 2 2 P e i p · ( x y ) = P δ ( x y )
and quite analogously
{ ψ + ( x ) , ψ + ( y ) } x 0 = y 0 = P + δ ( x y )
while all the remaining anti-commutators do vanish at arbitrary times, e.g. { ψ ( x ) , ψ + ( y ) } = 0  et cetera. It follows that the canonical equal time anticommutation relations for the massless Dirac quantum field (58) will be provided by
{ ψ ( x ) , ψ ( y ) } x 0 = y 0 = δ ( x y ) { ψ ( x ) , ψ ¯ ( y ) } x 0 = y 0 = γ 0 δ ( x y )
while the general canonical anticommutation relations take the standard usual form
{ ψ ( x ) , ψ ( y ) } = 0 = { ψ ¯ ( x ) , ψ ¯ ( y ) }
{ ψ ( x ) , ψ ¯ ( y ) } S 0 ( x y ) = γ ρ D 0 ( x y ) / x ρ
where D 0 ( x y ) is the well known massless Pauli-Jordan distribution, which fulfill the D’Alembert wave equation and endorses the microcausality property
γ 0 S 0 ( x y ) = 0 D 0 ( x y ) ( x y ) 2 < 0
Now the main point. A consistent and reliable canonical quantum theory for a free non interacting field must give rise to the causal Green’s function, or Feynman propagator, in the absence of which the perturbative expansion and the renormalization procedure cannot be implemented, so that any meaningful physical interpretation and predictions are prevented. The causal Green’s functions is the time-ordered product of field operators and provides the specific inversion of the kinetic operator which guarantees causality: namely, first a particle of the quantum field is created out of the vacuum and only later it is destroyed. Moreover, the causal Green’s function always admits the smooth transition to the Euclidean formulation, which is essential in order to evaluate the Feynman integrals in perturbation theory, in particular the divergent and regularized ones. As we shall see below, it turns out that only the canonical quantum theory of the massless Dirac spinor field can provide this featuring bonus: the Weyl chiral bispinor quantum fields () or any other combination of those ones but (58) can do the job. Let’s analyze this crucial issue in more detail.

3.1.4. Causal Green’s Functions

Consider the rank-two square matrix of the vacuum expectation value of the time ordered product of a pair of e.g. right-handed quantum Weyl fields: namely,
S c ( x y ) 0 | T ψ R ( x ) ψ R ( y ) | 0
After insertion of the normal-mode expansion we get
S c ( x y ) = θ ( x 0 y 0 ) ψ R ( x ) ψ R ( y ) 0 θ ( y 0 x 0 ) ψ R ( y ) ψ R ( x ) 0 = θ ( x 0 y 0 ) d p [ ( 2 π ) 3 2 ] 1 u R ( p ) u R ( p ) × exp i ( x 0 y 0 ) + i p · ( x y ) θ ( y 0 x 0 ) d p [ ( 2 π ) 3 2 ] 1 v R ( p ) v R ( p ) × exp i ( x 0 y 0 ) i p · ( x y ) = θ ( x 0 y 0 ) i · σ d p [ ( 2 π ) 3 2 ] 1 exp i ( x 0 y 0 ) + i p · ( x y ) + θ ( y 0 x 0 ) i · σ d p [ ( 2 π ) 3 2 ] 1 exp i ( x 0 y 0 ) i p · ( x y )
Now, if we recall the definitions of the massless Wightman distributions [15]
i D 0 ( ± ) ( x y ) ± D p exp ± i p · x p 0 = = ± d p ( 2 π ) 3 2 exp ± i ( x 0 y 0 ) i p · ( x y )
then we obtain
S c ( x y ) = θ ( x 0 y 0 ) i · σ ( i ) D 0 ( ) ( x y ) + θ ( y 0 x 0 ) i · σ i D 0 ( + ) ( x y ) = i · σ [ θ ( y 0 x 0 ) i D 0 ( + ) ( x y ) θ ( x 0 y 0 ) i D 0 ( ) ( x y ) ] δ ( x 0 y 0 ) [ D 0 ( ) ( x y ) + D 0 ( + ) ( x y ) ] = i x · σ D 0 c ( x y )
where D 0 ( ξ ) = D 0 ( ) ( ξ ) + D 0 ( + ) ( ξ ) is the massless Pauli-Jordan distribution [15], which fulfills D 0 ( 0 , ξ ) = 0 , while
D 0 c ( x y ) θ ( y 0 x 0 ) i D 0 ( + ) ( x y ) θ ( x 0 y 0 ) i D 0 ( ) ( x y )
is the Feynman propagator, or causal Green’s function, of a massless scalar field, its Fourier representation being
D 0 c ( x y ) = i ( 2 π ) 4 d 4 p e i p · ( x y ) p 2 + i ε
so that i D 0 c ( ξ ) acts as the unique fully causal inverse of the d’Alembert differential operator.
As a consequence we eventually find [ in this paragraph we use the notation σ ˜ μ = σ ¯ μ  ]
σ ¯ · x S c ( x y ) = δ ( 4 ) ( x y )
which means that the rank-two matrix-valued causal Green’s function S c is the fully causal inverse of the right-handed differential operator σ ¯ · . By following step-by-step the above derivation mutatis mutandis, it is not difficult to prove that after setting
S ¯ c ( x y ) 0 | T ψ L ( x ) ψ L ( y ) | 0
we get
S ¯ c ( x y ) = i x · σ ¯ D 0 c ( x y ) σ · x S ¯ c ( x y ) = δ ( 4 ) ( x y )
Consider now a change of the inertial reference frame, i.e. x μ = Λ μ ν μ ( x + a ) ν . We have
D 0 c ( x y ) = D 0 c ( x y )
Λ L 1 σ ¯ μ Λ R μ = Λ R σ ¯ μ Λ R μ = Λ μ κ μ σ ¯ κ Λ μ μ ρ ρ = σ ¯ μ μ
and thereby
σ ¯ μ μ = Λ L σ ¯ μ μ Λ R 1
It follows that the transformation law of the Weyl causal Green’s functions
S ¯ c ( x y ) = Λ L S ¯ c ( x y ) Λ R 1
S c ( x y ) = Λ R S c ( x y ) Λ L 1
does necessarily involve the S L ( 2 , C ) matrices of both IRREPs D ( 0 , 1 2 ) and D ( 1 2 , 0 ) . To better grasp this subtle point, let us rewrite e.g. the causal Green’s function in its Fourier integral representation in its most suggestive form, viz.,
S ¯ c ( x y ) = i ( 2 π ) 4 d 4 p p 0 σ · p p 0 2 p 2 + i ε e i p · ( x y ) = i ( 2 π ) 4 d 4 p e i p · ( x y ) p 0 σ · p ( p 0 + i ε ) ( p 0 + i ε ) = d 4 p ( 2 π ) 4 e i p · ( x y ) i p 0 + i ε i p 0 + i ε p 0 σ · p 2
It turns out that in the above off-shell expression the rank-two square matrix
W = p 0 σ · p 2
is no longer a projector on left-handed two-component spin states because
W 2 = p 0 2 2 p 0 σ · p + p 2 4 p 2 W Tr W = p 0 R
as it was true for the corresponding on-shell matrix with p 0 = . Hence, in the process of inversion for the Weyl differential operators, leading to the Weyl causal Green’s functions, the very notion of chirality is definitely and unavoidably lost, just owing to the nature of the time ordered products of the Weyl quantum fields.
In order to get a bona fide fully causal inversion for the kinetic differential operator of a massless spinor we have to turn back to Dirac bispinors. Consider in fact the quantity
S 0 c ( x y ) = 0 | T ψ ( x ) ψ ¯ ( y ) | 0 = 0 S ¯ c ( x y ) S c ( x y ) 0 = 0 i · σ ¯ i · σ 0 D 0 c ( x y )
so that
i / S 0 c ( x y ) = ( 1 ) 0 0 D 0 c ( x y ) = i δ ( 4 ) ( x y )
which shows, as expected, that the massless Schwinger propagator S 0 ( x y ) is the unique inverse of the massless Dirac differential operator / , that fulfills causality forward and backward in time, its Fourier representation being
S 0 c ( x y ) = i ( 2 π ) 4 d 4 p p / p 2 + i ε e i p · ( x y )

3.1.5. Observables

As we have seen above, the consistent quantum theory of the massless Weyl-Dirac free field necessarily involves four types of quanta, although the phenomenological constraints seem to indicate that only two of them are physical. The 1-particle state associated to a real left-handed neutrino of momentum p and a negative helicity is given by c , p | 0 , while a real right-handed anti-neutrino of momentum q and positive helicity is provided by d + , q | 0 . The plane wave function that describes an incoming neutrino in a scattering experiment or decay process is therefore
0 | ψ ( x ) c , p | 0 = u , p ( x )
Quite analogously, an incoming right-handed anti-neutrino of momentum q and positive helicity will be described by the plane wave function
0 | ψ ¯ ( x ) d + , q | 0 = v ¯ + , q ( x )
The outgoing states for the very same physical particles will be clearly described by the plane wave functions
0 | c , p ψ ¯ ( x ) | 0 = u ¯ , p ( x ) 0 | d + , q ψ ( x ) | 0 = v + , q ( x )
All the remaining 1-particle state are not physical, as far as we know from experiments, in such a manner that the Standard Model gauge interactions will forbid the transitions from physical massless left-handed neutrinos and right-handed anti-neutrinos to unphysical Weyl particles, as it does.
Let us come now to the evaluation of the Observables - or better to say invariant and conserved charges - in the massless spinor quantum theory. For example, since we have
ψ ± ( y ) i μ ψ ± ( x ) = : ψ ± ( y ) i μ ψ ± ( x ) : d p p μ v ± , p * ( y ) v ± , p ( x )
with p 0 = , one can easily obtain
P μ = ι = , + d p p μ c ι , p c ι , p + d ι , p d ι , p p 0 = = P L μ + P R μ
In order to set up the spin operator, consider for instance a left Weyl field ψ ( t , z ) travelling along the O z axis so that
T 12 ( t , z ) T 21 ( t , z ) = 0 = ν S ν 21 ( t , z )
Σ = d z : S 021 ( t , z ) : Σ ˙ = 0
Σ = 1 2 d z : ψ ( 0 , z ) Σ 3 ψ ( 0 , z ) : = 1 2 d z : ψ L ( 0 , z ) σ 3 ψ L ( 0 , z ) :
Now, since for a 1d motion we have for = p z sgn ( p z )
ψ L ( t , z ) = d p z [ 4 π ] 1 2 × c , p z u L ( p z ) e i t + i z p z + d , p z v L ( p z ) e i t i z p z u L ( p z ) = σ 3 p z p z 1 0 v L ( p z ) = + σ 3 p z p z 0 1
and thereby
σ 3 u L ( p z ) = u L ( p z ) , σ 3 v L ( p z ) = v L ( p z ) ,
σ 3 ψ L ( t , z ) = d p z [ 4 π ] 1 2 × c , p z u L ( p z ) e i t + i z p z d , p z v L ( p z ) e i t i z p z
This entails that the general form of the spin operator for a one dimensional motion along e.g. the O z axis becomes
Σ z = 1 2 d p z c , p z c , p z d + , p z d + , p z c + , p z c + , p z + d , p z d , p z
which endorses the fact that neutrino particles are of negative helicity, while the antineutrinos exhibit positive helicity, owing to the opposite verse of the wave vector, the situation being reversed for the hypothetical right-handed neutrinos of positive chirality.
Turning to the internal symmetries and charges we find
Q = d x : ψ ( x ) ψ ( x ) + ψ + ( x ) ψ + ( x ) : Q + Q + Q = d x d p : c , p u , p * ( x ) + d , p v , p * ( x ) × d p c , p u , p ( x ) + d , p v , p ( x ) :
= d p c , p c , p d , p d , p
Q 5 = Q + Q
Finally, the dilatation charge actually reads
D = d x : 3 2 i ψ ( x ) ψ ( x ) + x ν T 0 ν ( x ) : = x 0 P 0 + i 3 2 Q d x x k : ψ ( x ) k ψ ( x ) :
We find
0 D = P 0 + d x x k : i 0 ψ ( x ) k ψ ( x ) : d x x k : ψ ( x ) k i 0 ψ ( x ) : = P 0 + d x x k : ψ ( x ) H 0 k ψ ( x ) : d x x k : ψ ( x ) k H 0 ψ ( x ) : = P 0
Since we have the well-known relation
0 D = i [ D , P 0 ] = P 0
it follows that we can write, as expected,
D ˙ = d D d t = 0 D + i [ P 0 , D ] = 0
so that the dilatation charge is globally time independent.
From the above analysis and correspondences, it follows that the 1-particle states c ( p ) | 0 and d + ( p ) | 0 of the Fock space of the massless Dirac spinor field do respectively represent the physical particles with the following quantum number assignements: namely,
i) opposite chirality χ = 1 ;
ii) the same energy-momentum p ν = ( , p ) ;
iii) opposite internal charge, called lepton number, = ± 1 ;
iv) opposite helicity h = 1 .

3.1.6. Discrete Symmetries

Turning to the discrete symmetries beside the Nöther theorem, we start with the charge conjugation unitary transformation of the quantum theory for a massless Dirac spinor field, which amounts to the simultaneous exchange between chirality and particle-antiparticle operators, up to a phase factor, viz.,
C c , p C 1 = e i θ d ± , p C d , p C 1 = e i θ c ± , p
C 2 = I C = C = C 1 0 θ < 2 π
that yields
C ψ ( x ) C = e i θ d p d ± , p u , p ( x ) + c ± , p v , p ( x )
Now, from the charge conjugation relations for the plane wave functions, viz.,
u , p ( x ) = i γ 2 v ± , p * ( x ) v , p ( x ) = i γ 2 u ± , p * ( x )
we can write
C ψ ( x ) C ψ c ( x ) = e i θ d p d ± , p u , p ( x ) + c ± , p v , p ( x ) = i γ 2 e i θ d p c ± , p u ± , p * ( x ) + d ± , p v ± , p * ( x ) = γ 2 ψ ± ( x ) for θ = π / 2
Now, since we have ψ ( x ) = ψ ( x ) + ψ + ( x ) , we come to the result
ψ c ( x ) = C ψ ( x ) C = γ 2 ψ ( x )
which is nothing but the quantum generalization of the charge conjugation rule for classical spinor fields, i.e. ψ c ( x ) = γ 2 ψ * ( x ) , as expected.
Thus the charge conjugated of a particle of e.g. negative chirality, negative helicity and positive lepton number is an anti-particle of opposite chirality, positive helicity and negative lepton number. It turns out that the charge conjugation transformation necessarily involves both chiralities, in order to be implemented. Hence the C transform is well defined only for 4-component bispinors and not for 2-component Weyl spinors belonging to the two irreducible 2d representations.
The parity transformation is a discrete space-time symmetry of the Action integral S = d 4 x L 0 ( x ) that can be implemented at the quantum field theory level as
P c , p P = e i η c ± , p P d , p P = e i θ d ± , p 0 η < 2 π 0 θ < 2 π
After insertion of the normal mode expansion of the massless Dirac quantum field, then we find with x μ = x ˜ μ = ( x 0 , x ) = x μ
ψ ( x ) = P ψ ( x ˜ ) P = d p P c , p P u , p ( x ˜ ) + P d , p P v , p ( x ˜ ) = d p e i η c ± , p u , p ( x ˜ ) + e i θ d ± , p v , p ( x ˜ ) = d p e i η c ± , p u , p ( x ˜ ) + e i θ d ± , p v , p ( x ˜ )
Now we have
u , p ( x ˜ ) = [ ( 2 π ) 3 2 ] 1 2 u ( p ) e i t + i p · x = γ 0 u ± , p ( x ) v , p ( x ˜ ) = [ ( 2 π ) 3 2 ] 1 2 v ( p ) e i t i p · x = γ 0 v ± , p ( x )
so that, after choosing
η = k π θ = 2 k π ( k Z )
we definitely obtain the usual transformation law
ψ ( x ) = P ψ ( x ˜ ) P = γ 0 ψ ± ( x )
ψ ( x ) = P ψ ( x ˜ ) P = γ 0 ψ ( x )
P 2 = I P = P = P 1
Once again, it is worthwhile to remark that even the parity transform is not a well defined operation for a two-component Weyl fermion, because it necessarily involves both chiralities, i.e. both irreducible 2d representations of the proper Lorentz group.
In quantum field theory the time reversal operator T corresponds to an antilinear and antiunitary transformation [13,14] which acts as follows: namely,
[ P 0 , T ] = 0 { P , T } = 0 { M ρ σ , T } = 0
where, as usual, ( P μ , M ρ σ ) are the ten generators of the Poincaré group. It follows therefrom that the antilinear and antiunitary time reversal operation does reverse the signs of all particle momenta and spin angular momenta, the same being true for antiparticles. In this regard we need to explain what is meant by spin-flip [14]. Consider for example the spin state
u ( p ) = u L ( p ) 0 0 = 1 p z p z p x i p y 0 0 u L ( p ) = σ · p ˜ p z 1 0
This positive energy, left-handed spin-state of negative chirality does involve the spin-up constant spinor 1 0 and to our purpose it will be better labelled as u ( p , ) . Its corresponding spin reversed state u ( p , ) has the very same structure, but for the constant spinor that will be turned to the spin-down 0 1 . Thus we can denote the pair of normalized spin flip states according to
u ( p , ) = 1 p z p z p x i p y 0 0
u * ( p , ) = 1 p z p x + i p y p z 0 0
and quite analogous relationships for the remaining spin states. Then we have
T c , p T 1 = e i η c , p T d , p T 1 = e i θ d , p ψ ( x ) = T ψ ( x 0 , x ) T 1
= d p e i η c , p u , p * ( x ˜ ) + e i θ d , p v , p * ( x ˜ )
= d p e i η c , p u , p * ( x ˜ ) + e i θ d , p v , p * ( x ˜ )
the complex conjugation being due to the antilinear and antiunitary nature of the T operation, while the spin-flip operation in the spinor plane wave functions is understood. Consider now, for instance,
T ψ ( ) ( t , x ) T 1 = d p [ ( 2 π ) 3 2 ] 1 2 × e i η c , p u * ( p ) + e i η + c + , p u + * ( p ) e i t + i p · x
It can be readily seen that the right-hand side of the above equation becomes a local expression for the quantum field at the instant t, notably Θ ψ ( ) ( t , x ) , if a 4 × 4 matrix can be found such that
e i η u * ( p ) = Θ u ( p )
From the explicit form (94-) it turns out that the above relations hold true for η = 0 . The othogonality and closure relations for the spin-states require the unitarity property Θ Θ = I = Θ Θ . Moreover, from the spin-state equations we obtain
α · p u ( p ) = u ( p ) α * · p u * ( p ) = u * ( p ) α * · p Θ u ( p ) = Θ u ( p ) Θ α * Θ = α
The unitary solution of the very last equation is unique except for an irrelevant phase factor. In the Weyl representation of the Clifford algebra, the real matrix
Θ = γ 1 γ 3 = i Σ 2
is a solution with the important property Θ * Θ + I = 0 , which does not depend on the specific representation of the Dirac matrices. Thus we definitely find
ψ ± ( x ) = T ψ ± ( t , x ) T 1 = Θ ψ ± ( t , x ) ( x = x ˜ )
which shows that the time reversal is a chirality preserving transformation. Moreover, time reversal is well defined even for two component Weyl spinors, since the spin-flip operation does not mix left and right components of a massless Dirac bispinor. This means that also the CP and CPT transformations are well defined for two-component Weyl spinors. In summary, the quantum theory which arises from the classical Action integral
S 0 d 4 x L 0 ( x ) = d 4 x ψ ( x ) i α · ψ ( x )
is invariant under the discrete C , P , T transformation separately, while the Weyl Action integrals
S L = d 4 x 1 2 ψ L ( x ) i σ · ψ L ( x ) S R = d 4 x 1 2 ψ R ( x ) i σ ¯ · ψ R ( x )
and related quantum theories are invariant under T and CP transforms, which are well defined and good symmetries for left and right-handed Weyl spinors.

3.2. Majorana Spinor Field

One can write a relativistic invariant field equation for a massive 2-component spinor field: to this aim, let us start from a left Weyl spinor ψ L D ( 1 2 , 0 ) that transforms according to the S L ( 2 , C ) matrix Λ L . Call such a 2-component spinor field χ a ( x ) ( a = 1 , 2 ) .

3.2.1. Definitions and Lagrangians

Let us consider the Weyl spinor wave field as a classical anti-commuting field, i.e. a Graßmann valued left Weyl spinor field function over the Minkowski space which satisfies
{ χ a ( x ) , χ b ( y ) } = 0 x , y M | a , b = 1 , 2
together with the complex conjugation rule
( χ 1 χ 2 ) * = χ 2 * χ 1 * = χ 1 * χ 2 *
so as to imitate the Hermitean conjugation of quantum fields. A Majorana classical spinor field is a self-conjugated bispinor, that can be constructed, for example, out of the left-handed spinor χ a ( x ) ( a = 1 , 2 ) as follows: namely,
χ M ( x ) = χ ( x ) σ 2 χ * ( x ) = χ M c ( x )
the charge conjugation rule for any classical bispinors ψ being defined by the general relationship
ψ c ( x ) = e i θ γ 2 ψ * ( x ) ( 0 θ < 2 π )
which is a discrete internal - i.e. space-time point independent - symmetry transformation. Here below we shall suitably choose θ = 0 . The Majorana bispinor has a right-handed lower Weyl spinor component σ 2 χ * D ( 0 , 1 2 ) , albeit functional dependent, due to the charge self-conjugation constraint, in such a manner that χ M possesses both chiralities and polarizations, at variance with its left-handed Weyl building spinor χ ( x ) . There is another kind of self-conjugated Majorana bispinor, which can be set-up out of a right-handed Weyl building spinor φ D ( 0 , 1 2 )
φ M ( x ) = σ 2 φ * ( x ) φ ( x ) = φ M c ( x )
From the Majorana self-conjugated bispinors, one can readily construct the most general Poincaré invariant and power counting renormalizable Lagrangian. For instance, by starting from the bispinor χ M ( x ) we have
L M = 1 4 χ M ( x ) α μ i μ χ M c ( x ) 1 2 m χ ¯ M ( x ) χ M c ( x )
where α ν = γ 0 γ ν , while the employed notation reminds us that the upper and lower components of a Majorana bispinor can never be treated as functionally independent, even formally, due to the presence of the self-conjugation constraint. It follows that the Majorana mass term can be written in the two equivalent forms
L M m = 1 2 m χ ( x ) σ 2 χ ( x ) + χ ( x ) σ 2 χ * ( x ) = i m χ 1 ( x ) χ 2 ( x ) + χ 1 * ( x ) χ 2 * ( x ) = i m χ 2 * ( x ) χ 1 * ( x ) + χ 2 ( x ) χ 1 ( x ) = L M m *
whence it is clear that the corresponding integral d 4 x L M m ( x ) is not invariant under the internal U ( 1 ) phase transformation
χ ( x ) χ ( x ) = χ ( x ) e i θ ( 0 θ < 2 π )
Concerning the kinetic term, from the relations
χ M = χ χ σ 2 χ = ( χ * )
together with
α ν = γ 0 γ ν = σ ν 0 0 σ ¯ ν σ μ = σ ¯ μ = ( 1 , σ )
α ν γ 2 = 0 σ ν σ 2 σ ¯ ν σ 2 0
one can obtain the Majorana kinetic term in the 2-component formalism
1 4 χ M ( x ) α μ γ 2 i μ χ M * ( x ) = 1 4 χ ( x ) σ μ i μ χ ( x ) + 1 4 χ ( x ) σ 2 σ ¯ μ σ 2 i μ χ * ( x ) = 1 4 χ ( x ) σ μ i μ χ ( x ) + 1 4 χ ( x ) ( σ ¯ μ ) * i μ χ * ( x ) = 1 2 χ ( x ) σ μ i μ χ ( x ) + 1 2 χ ( x ) ( σ ¯ μ ) * i μ χ * ( x ) = χ ( x ) σ μ i μ χ ( x ) 1 2 i μ χ ( x ) σ μ χ ( x )
which coincides with the Lagrangian for a left-handed Weyl spinor, as expected, where the above complex conjugation rule (102) for Grassmann valued field has been used.
It turns out that even the Majorana Action integral
S M = 1 4 d 4 x χ M ( x ) α μ γ 2 i μ χ M * ( x ) 2 m χ M ( x ) α 2 χ M * ( x )
= d 4 x χ ( x ) σ μ i μ χ ( x ) + 1 2 m χ ( x ) σ 2 χ * ( x ) + χ ( x ) σ 2 χ ( x )
is not invariant under the overall phase transformation χ M ( x ) = e i θ χ M ( x ) of the Majorana bispinor, just like the previously discussed mass term. Hence, as it will be further endorsed after the transition to the Majorana representation of the Dirac matrices, there is no invariant scalar charge or lepton number for a Majorana massive spinor. According to the general Noether theorem for a Lie group of global internal symmetry transformations
μ J a μ = δ L / δ ω a ( a = 1 , 2 , , s )
where ω a are the parameters of the internal symmetry Lie group of dimension s, we find
J M μ = ( δ L M / δ μ χ ) i χ = χ σ μ χ δ L M / δ θ = i m χ σ 2 χ χ σ 2 χ *
so that
μ J M μ = μ ( χ σ μ χ ) = i m ( χ σ 2 χ * χ σ 2 χ )
We recall that the canonical dimensions of any kind of spinor ψ in physical units are [ ψ ] = eV cm 1 . It follows that the Noether charge is the lepton number, a pure number , which can be identified as follows: namely,
( t ) = 1 c d x χ ( t , x ) χ ( t , x )
which satisfies the evolution equation
d d t = i m 2 d x χ ( x ) σ 2 χ ( x ) χ ( x ) σ 2 χ * ( x )
that leads to the lepton number conservation only in the massless case.
Actually, as it will be discussed and clarified in the sequel, there is a relic continuous U(1) symmetry only for Majorana massless spinors, which drives to the existence of a conserved pseudo-scalar charge, the meaning of which will be better focused further on.
Another key observation is now in order. Starting with a left-handed Weyl 2-component spinor χ ( x ) D ( 1 2 , 0 ) one can also build up a chiral left-handed bispinor
ν L ( x ) = P L χ M ( x ) = χ ( x ) 0 P R ν L = 1 2 ( 1 + γ 5 ) ν L = 0
together with its right-handed counterpart
ν L c ( x ) = γ 2 ν L * ( x ) = 0 σ 2 χ * ( x ) ν R ( x ) P L ν R = 1 2 ( 1 γ 5 ) ν R = 0
which still involves only the very same left-handed Weyl 2-component spinor χ ( x ) D ( 1 2 , 0 ) . With the help of those two chiral left- and right-handed bispinors, one can suitably describe a massive neutrino, once we employ the related Lagrangian density
L L = 1 2 ν ¯ L ( x ) γ μ i μ ν L ( x ) 1 2 m ν ¯ L ( x ) ν R ( x ) + c . c .
= 1 2 χ ( x ) σ μ i μ χ ( x ) + 1 2 m χ ( x ) σ 2 χ * ( x ) + c . c .
However, it turns out that the above Lagrangian does actually coincides with the Majorana Lagrangian in (104) in such a manner that the field equations for a chiral left-handed massive bispinor ν L ( x ) = P L χ M ( x ) are nothing but the very same as the massive Majorana spinor χ ( x ) , the general solutions of which will be fully discussed in the sequel, after the transition to the Majorana representation.
The full Majorana Lagrangian in the 2-component formalism reads
L ^ M = 1 2 χ ( x ) σ μ i μ χ ( x ) + 1 2 m χ ( x ) σ 2 χ ( x ) + c . c .
so that the Euler-Lagrange field equation becomes
i σ μ μ χ ( x ) + m σ 2 χ * ( x ) = 0
which is known as the Majorana field equation. Multiplying to the left by σ 2 and taking complex conjugation yields
i σ 2 0 χ * ( x ) i σ 2 σ k * k χ * ( x ) + m χ ( x ) = 0
Remembering that we have σ 2 σ σ 2 = σ and that σ ¯ μ = ( 1 , σ ) we come to the equivalent form of the Majorana wave equation, viz.,
i σ ¯ μ σ 2 μ χ * ( x ) + m χ ( x ) = 0
Now, if we act from the left with the operator i σ ¯ ν ν to equation (113) and use equation (114) we find
σ ¯ ν σ μ ν μ χ ( x ) + m 2 χ ( x ) = ( + m 2 ) χ ( x ) = 0
which means that the left-handed spinor χ D ( 1 2 , 0 ) , i.e. the building block of the self-conjugated Majorana bispinor, is actually solution of the Klein-Gordon wave equation. Moreover, by understanding the Majorana bispinor to be defined by the self-conjugation constraint χ M ( x ) = χ M c ( x ) , it is easy to check that the pair of coupled wave equations (113) and (114) is equivalent to the single bispinor wave equation
( α ν i ν β m ) χ M ( x ) = 0
where use has been made of the Dirac notation β = γ 0 , while the Majorana Lagrangian can be recast in a further 4-component form
L M = 1 4 χ M ( x ) α μ i μ χ M c ( x ) 1 2 m χ M ( x ) β χ M c ( x )
It is immediate to verify that the bispinor form (115) of the field equations does coincide with the pair of functionally dependent and equivalent forms (113) and (114) of the Majorana spinor wave equation. Notice that the Majorana self-conjugated bispinor transforms under the Poincaré group according to
χ M ( x ) = Λ 1 2 χ M c ( x ) Λ 1 2 = Λ L 0 0 Λ R
with x = Λ ( x + a ) . Hence it follows that the generators of the Lorentz group acting upon any kind of bispinor are provided by the usual expressions
S 0 k = i 2 γ 0 γ k = 2 i σ k 0 0 i σ k
S j k = i 4 [ γ j , γ k ] = 1 2 ε j k σ 0 0 σ = 1 2 ε j k Σ
in such a manner that we obtain the infinitesimal variation
Δ χ M ( x ) = 1 2 i σ μ ν χ M ( x ) δ ω μ ν σ μ ν = i 4 [ γ μ , γ ν ]

3.2.2. The Majorana Real Representation

Let us now search the general complex solutions of the Majorana linear wave equation (115) for a self-conjugated bispinor χ M ( x ) . Once obtained, we can readily get the general complex solution for the massive chiral left-handed neutrino ν L ( x ) after a projection on the upper component of χ M ( x ) . To implement this task, it appears utmost convenient to employ the so called Majorana representation [16]. It turns out that, by definition, the Majorana bispinor χ M ( x ) = χ M c ( x ) must fulfill the self-conjugation constraint, which linearly relates the lower spinor component to the complex conjugate of the upper spinor component. Then, a representation must exist which makes the Majorana bispinor real, with the previously introduced two independent complex functions χ a C ( a = 1 , 2 ) replaced by the four real functions ψ M , α R ( α = 1 , 2 , 3 , 4 ) . To obtain this real representation, we notice that for any complex bispinor ψ we find
ψ c = γ 2 ψ * ψ * = γ 2 ψ c
in such a manner that for the self-conjugated complex bispinor χ M we can write
χ M * ( x ) = γ 2 χ M ( x )
Then a linear transformation to a real bispinor field ψ M = ψ M * can be found by noting that
χ M = S ψ M χ M * = S * ψ M = S * S 1 χ M γ 2 + S * S 1 = 0
The explicit form of the unitary transformation matrix can be readily found to be
S = 2 2 1 σ 2 σ 2 1 = 2 2 ( 1 + γ 2 ) S = S 1 = 2 2 ( 1 γ 2 )
or even more explicitly
S = 2 2 1 0 0 i 0 1 i 0 0 i 1 0 i 0 0 1
From the above relation ψ M = S 1 χ M = ψ M * one can immediately obtain the correspondence rule between the complex and real forms of the self-conjugated Majorana bispinor: namely,
ψ M 1 = 2 e χ 1 ψ M 2 = 2 e χ 2 ψ M 3 = 2 m χ 2 ψ M 4 = 2 m χ 1
χ 1 = ( ψ M 1 i ψ M 4 ) / 2 χ 2 = ( ψ M 2 + i ψ M 3 ) / 2
Thus we can suitably make use of the so called Majorana representation for the gamma matrices, which is purely imaginary, and provided by a similarity transformation acting on the γ matrices in the Weyl representation, viz.,
γ M μ S γ μ S
that yields
γ M 0 = σ 2 0 0 σ 2 γ M 5 = 0 σ 2 σ 2 0
γ M 1 = i σ 3 0 0 i σ 3 γ M 2 = 0 σ 2 σ 2 0 γ M 3 = i σ 1 0 0 i σ 1
which satisfy by direct inspection the Clifford algebra
{ γ M μ , γ M ν } = 2 g μ ν { γ M ν , γ M 5 } = 0
γ M 0 = β M γ M k = γ M k γ M 5 = γ M 5
together with
γ M ν = γ M ν * γ M 5 = γ M 5 *
The result is that, at the place of a complex self-conjugated bispinor, which has been constructed out of a left-handed Weyl spinor χ M ( x ) , one can safely and more suitably employ a real Majorana bispinor: namely,
χ M ( x ) = χ M c ( x ) ψ M ( x ) = S χ M ( x ) = ψ M * ( x )
a quite analogous construction being obviously there, had we started from a right-handed Weyl spinor φ D ( 0 , 1 2 ) . Then, the Majorana Lagrangian and the ensuing Majorana wave equation take the manifestly real forms
L M = 1 4 ψ M ( x ) α M ν ν ψ M ( x ) + 1 2 i m ψ M ( x ) β M ψ M ( x )
( i / M m ) ψ M ( x ) = 0 ψ M ( x ) = ψ M * ( x )
α M ν = γ M 0 γ M ν α M 0 = I β M γ M 0
with the real symmetric matrices
α M 1 = σ 1 0 0 σ 1 α M 2 = 0 1 1 0 α M 3 = σ 3 0 0 σ 3
which manifestly fulfill α M = α M * = α M . It turns out that, from the manifestly real form of the Majorana Lagrangian, the only relic internal symmetry of the Majorana’s Action integral - but in the massless case - is the discrete Z 2 symmetry, i.e. ψ M ( x ) ψ M ( x ) .
Singular massless limit. Actually, it turns out that only in the massless casethere exists a further accidental invariance of the Action integral under the U(1) internal symmetry group
ψ M ( x ) ψ M ( x ) = exp ± i θ γ M 5 ψ M ( x ) ( 0 θ < 2 π )
the imaginary unit being convenient to keep the reality of the trasformed Majorana bispinor. From the Noether theorem we get the corresponding real tetra-current, which satisfies the continuity equation
j 5 μ ( x ) = 1 2 ψ M ( x ) α M μ i γ M 5 ψ M ( x ) · j 5 ( x ) = 0
as well as the ensuing conserved pseudo-scalar charge
± Q 5 = ± 1 2 d x ψ M ( t , x ) i γ M 5 ψ M ( t , x ) Q ˙ 5 = 0
the overall ± sign being conventional and irrelevant. In order to better understand the meaning of this quantity we notice that, if we turn to the complex form χ M ( x ) of the self-conjugated bispinor, we can immediately realize that the above U(1) internal symmetry of the massless case is nothing but the ordinary phase transformation for the independent building spinor χ ( x ) . As a matter of fact, the phase transformation χ ( x ) = e i θ χ ( x ) just induces the chiral transformation χ M ( x ) = e i θ γ 5 χ M ( x ) on the complex form of the Majorana bispinor. Notice that the quantity
ψ M i γ M 5 ψ M = ψ M ψ M 0 i σ 2 i σ 2 0 ψ M ψ M
does vanish for ordinary real functions ψ M , α R ( α = 1 , 2 , 3 , 4 ) while in the case of anti-commuting Grassmann-valued functions we get
± Q 5 = ± d x ψ M , 1 ( t , x ) ψ M , 4 ( t , x ) ψ M , 2 ( t , x ) ψ M , 3 ( t , x )
the integrated quantity being nothing but than i χ ( x ) χ ( x ) as expected.
A very important remark is now in order. The Majorana Lagrangian (121) and wave equation (122) have been derived starting from an original two-component left-handed spinor χ ( x ) or, equivalently, a self-conjugated bispinor χ M ( x ) = χ M c ( x ) . However, a closely related and still real Lagrange density, just leading to the very same field equation, can be immediately written even for any arbitraty complex bispinor Ψ M ( x ) : namely,
L M = 1 2 Ψ M ( x ) α M ν i ν Ψ M ( x ) m Ψ M ( x ) β M Ψ M ( x )
( i / M m ) Ψ M ( x ) = 0
At variance with (121) the Lagrange density (123) does enjoy the invariance under the global U(1) symmetry transformations Ψ M ( x ) = e i θ Ψ M ( x ) ( 0 θ < 2 π ) driving eventually to the charge or number conservation. It follows therefrom that the real Majorana wave equation will admit in general both real as well as complex solutions, with quite different physical properties and behaviors.

3.2.3. Solution of the Majorana Wave Equation

One can rewrite the Majorana wave equation (122) in the Schrödinger form and in physical units
i ψ M t = m c 2 β M ψ M c α M k i k ψ M = H M ψ M ( x )
whence we can read the Majorana one-particle Hamiltonian operator
H M = α M k p ^ k c + m c 2 β M ( p ^ k = i k )
To solve the Majorana linear and real wave equation (122) we turn to the Fourier representation
ψ M ( x ) = d 4 p ( 2 π ) 3 / 2 ψ ˜ M ( p ) exp { i p · x } = d 4 p ( 2 π ) 3 / 2 ψ ˜ M ( p ) exp { i p · x }
so that from (122) we get
( p / M m ) ψ ˜ M ( p ) = 0 ( p / M + m ) ψ ˜ M ( p ) = 0 p / M p ν γ M ν
which implies
ψ ˜ M , α ( p ) = ( p / M + m ) α β ϕ ˜ β ( p ) p 2 m 2 ϕ ˜ α ( p ) = 0 ( α = 1 , 2 , 3 , 4 ) ϕ ˜ α ( p ) = δ p 2 m 2 f α ( p )
and in turn
ψ ˜ M , α ( p ) = ( m p / M ) α β ϕ ˜ β ( p ) = ( m p / M ) α β f β ( p ) δ p 2 m 2
The momentum space Majorana equations (126) are a pair of algebraic equation
( p 0 β M p k γ M k m c ) ψ ˜ M ( p ) = 0 ( p 0 β M p k γ M k + m c ) ψ ˜ M ( p ) = 0 c p 0 = ± ω p = ± p 2 c 2 + m 2 c 4
which can be turned into a pair of energy eigenvalue problems
H ˜ M ( p ) ψ ˜ M ( p ) = ± ω p ψ ˜ M ( p ) c p k α M k + m c 2 β M H ˜ M ( p )
where H ˜ M ( p ) is the 1-particle Majorana Hamiltonian in the momentum space representation, which satisfies the property H ˜ M 2 ( p ) = ω p 2 .
In order to find explicit solutions in momentum space, let us first consider the following matrix-valued functions: namely,
E M ± ( p ) = m c ± p / M / 2 m c with c p 0 = ± ω p = ± p 2 c 2 + m 2 c 4
which do satisfy the following relationships
[ E M ± ( p ) ] * = E M ( p ) [ E M ± ( p ) ] = E M ± ( p )
[ E M ± ( p ) ] 2 = E M ± ( p ) E M ± ( p ) E M ( p ) = 0
tr E M ± ( p ) = 2 E M + ( p ) + E M ( p ) = I
The above matrix-valued functions can be fully understood as projectors, if we simultaneously perform the matrix Hermitean conjugation and the full inversion, or parity transform, on the function variable p p . The physical meaning of those projectors can be better featured and focused on the Majorana particle rest frame p = 0 , in which p / M = m c β M . To this concern, we recall the spin operators S k = 1 2 Σ k ( k = 1 , 2 , 3 ) in the Majorana representation where
Σ M , 1 = i γ M 2 γ M 3 = 0 i σ 3 i σ 3 0 Σ M , 3 = i γ M 1 γ M 2 = 0 i σ 1 i σ 1 0
Σ M , 2 = i γ M 3 γ M 1 = σ 2 0 0 σ 2
From the relationships
β M 2 = 1 = Σ M , 2 2 [ β M , Σ M , 2 ] = 0
it appears that one can actually set up a complete orthogonal tetrad of bispinors, which are common eigenvectors of the Hermitean commuting matrices β M and Σ M , 2 . It turns out that the spin states ξ r ( r = ± ) are the pair of degenerate eigenstates of the Majorana 1-particle Hamiltonian H ˜ M ( p ) in the particle rest frame p = 0 , with positive eigenvalue p 0 = m and with opposite spin projections on the O Y axis, viz.,
β M ξ ± = ξ ± Σ M , 2 1 ξ ± = 0 ξ ξ ± = 0
In the very same way we obtain
β M η ± = η ± Σ M , 2 1 η ± = 0 η η ± = 0
which means that the spin states η r ( r = ± ) are in fact the two degenerate eigenstates of the Majorana 1-particle Hamiltonian in the particle rest frame p = 0 , with negative eigenvalue p 0 = m and with opposite spin projections on the O Y axis. Altogether we can write
ξ r ξ s = 2 δ r s = η r η s ξ r η s = 0 r , s = ±
Notice that, from the relations { γ 5 , γ 0 } = 0 = [ γ 5 , Σ 2 ] , it follows that the chiral matrix γ 5 maps the positive energy space into the negative one, and viceversa, in the Majorana particle rest frame p = 0 . Hence, without loss of generality, we can always assume
γ M 5 ξ ± = η ± γ M 5 η ± = ξ ±
In general, the freedom in the choice of an orthogonal pair, in the positive and negative energy planes of the rest frame p = 0 , is provided by the S O ( 2 , R ) Abelian symmetry little group
F ± ( θ ) = exp { ± i Σ M , 2 θ } = I cos θ ± γ M 1 γ M 3 sin θ [ F ± ( θ ) ] 1 = [ F ± ( θ ) ] = F ± ( θ )
because β M + β M = 0 , γ M k = γ M k , in such a manner that
ξ r = s = + , [ F + ( θ ) ] r s ξ s η r = s = + , [ F ( ϑ ) ] r s η s
To sum up, one can set up the positive and negative energy spaces projectors ε ± , in the rest frame of the massive Majorana particle, which satisfy
ε + = 1 2 ( ξ + ξ + + ξ ξ ) ε = 1 2 ( η + η + + η η )
ε ± 2 = ε ± ε ± ε = 0 ε + + ε = 1 tr ε ± = 2
Now, in order to move from the rest frame to an arbitrary inertial frame we can proceed as follows: consider the so called spin-states of positive and negative energy, viz.,
u r ( p ) 2 m ( 2 ω p + 2 m ) 1 2 E M + ( p ) ξ r v r ( p ) 2 m ( 2 ω p + 2 m ) 1 2 E M ( p ) η r ( r = + , )
which again satisfy
γ M 5 u r ( p ) = v r ( p ) γ M 5 v r ( p ) = u r ( p ) ( r = + , )
The above spin-states do fulfill by construction
2 m E M ( p ) u r ( p ) = ( m p / M ) u r ( p ) = 0 2 m E M + ( p ) v r ( p ) = ( m + p / M ) v r ( p ) = 0 ( p 0 = ω p , r = + , )
and they are the two pair of degenerate eigenstates of the positive and negative energy projectors
E M + ( p ) u r ( p ) = u r ( p ) E M ( p ) v r ( p ) = v r ( p ) ( r = + , )
as well as of the Majorana 1-particle Hamiltonian in momentum space, viz.,
H ˜ M ( p ) u r ( p ) = ω p u r ( p ) H ˜ M ( p ) v r ( p ) = ω p v r ( p ) ( r = + , )
By direct inspection one can readily verify the orthogonality relations
u r ( p ) u s ( p ) = 2 ω p δ r s u ¯ r ( p ) u s ( p ) = 2 m δ r s v r ( p ) v s ( p ) = 2 ω p δ r s v ¯ r ( p ) v s ( p ) = 2 m δ r s ( r , s = + , )
together with
u r ( p ) v s ( p ) = u ¯ r ( p ) v s ( p ) = 0 ( r , s = + , )
For example we have
v r ( p ) v s ( p ) = ( 2 ω p + 2 m ) 1 η r ( m p / ˜ ) ( m p / ) η s = ω p ( ω p + m ) 1 η r ( ω p + p · γ β M m β M ) η s = 2 ω p η r η s = 2 ω p δ r s
where use has been made of the reations β M η s = η s and η r γ M η s 0 . Furthermore we get
u ¯ r ( p ) u s ( p ) = ( 2 ω p + 2 m ) 1 ξ r ( m + p / ˜ ) β M ( m + p / ) ξ s = ( 2 ω p + 2 m ) 1 ξ r ( m + p / ) 2 ξ s = 2 m ( 2 ω p + 2 m ) 1 ξ r ( m + p / ) ξ s = 2 m δ r s u ¯ r ( p ) v s ( p ) = ( 2 ω p + 2 m ) 1 ξ r ( m + p / ˜ ) β M ( m p / ) η s = ( 2 ω p + 2 m ) 1 ξ r ( m + p / ) ( m p / ) η s 0 u r ( p ) v s ( p ) = ( 2 ω p + 2 m ) 1 ξ r ( m + p / ˜ ) ( m p / ˜ ) η s 0
Finally, from the relation β M η ± = η ± we get for N = ( 2 ω p + 2 m ) 1 / 2
v ± ( p ) = N ( m p / M ) β M η ± = N ( ω p m β M + p k α M k ) η ± = N [ ω p H ˜ M ( p ) ] η ±
in such a manner that from the relationship H ˜ M 2 ( p ) = ω p 2 we can immediately obtain
H ˜ M ( p ) v r ( p ) = N ω p [ H ˜ M ( p ) ω p ] η r = ω p v r ( p ) ( r = + , )
as expected.
Moreover, it is possible to prove the closure relations
r = + , u r ( p ) u ¯ r ( p ) = p / + m r = + , v r ( p ) v ¯ r ( p ) = p / m ( p 0 = ω p )
Proof. Consider for instance the closure relation for the negative energy spin-states. It turns out that
r = + , v r ( p ) v ¯ r ( p ) v s ( p ) = 2 m v s ( p ) ( p / m ) v s ( p ) = 2 m E M ( p ) v s ( p ) = 2 m v s ( p )
r = + , v r ( p ) v ¯ r ( p ) u s ( p ) = 0 = ( p / m ) u s ( p )
which means that the two matrix-valued functions on both sides of the above equality do share the same rank two, the same degenerate eigenvalues ( 2 m , 0 ) and the corresponding very same eigenvectors. A quite analogous conclusion holds obviously true even concerning closure relation for the positive energy spin-states. It follows that the matrix-valued functions in the equality (137) are proportional. However, if we multiply by β M from the right both closure relations and take the trace, then we get the same result 4 ω p , which means that the proportionality constant is one, as required.
As a corollary of the closure relations, after direct substitution, we come to the nontrivial general identity
( m ± p / ) ε ± ( m ± p / ) = ( p / ± m ) ( ω p + m )
Concerning the Majorana spinor wave functions, it turns out that a complete orthonormal set of complex plane wave bispinor solutions of the Majorana wave equation (122) is provided by
u p , r ( x ) [ ( 2 π ) 3 2 ω p ] 1 2 u r ( p ) exp { i ω p t + i p · x } v p , r ( x ) [ ( 2 π ) 3 2 ω p ] 1 2 v r ( p ) exp { i ω p t i p · x } ( r = + , )
which fulfill the energy eigenvalue problem
H M u p , r ( x ) = ω p u p , r ( x ) H M v p , r ( x ) = ω p v p , r ( x )
as well as the orthogonality and closure relations: namely,
d x u p , r ( x ) u q , s ( x ) = δ r s δ ( p q ) = d x v p , r ( x ) v q , s ( x )
i r = + , d p u p , r ( x ) u ¯ p , r ( y ) = ( i / M + m ) D ( ) ( x y ) i r = + , d p v p , r ( x ) v ¯ p , r ( y ) = ( i / M + m ) D ( + ) ( x y )
where D ( ± ) ( x y ) are the well known tempered distribution
D ( ± ) ( x y ) = ± 1 i ( 2 π ) 3 d 4 p exp { ± i p · ( x y ) } δ ( p 2 m 2 ) θ ( p 0 )
Proof. By direct replacement we find
i r = + , d p u p , r ( x ) u ¯ p , r ( y ) = i r = + , d p [ ( 2 π ) 3 2 ω p ] 1 u r ( p ) u ¯ r ( p ) exp { i ω p ( x 0 y 0 ) + i p · ( x y ) } = i ( i γ M μ μ , x + m ) d p [ ( 2 π ) 3 2 ω p ] 1 exp { i ω p ( x 0 y 0 ) + i p · ( x y ) } = i ( 2 π ) 3 ( m + i / M ) d 4 p exp { i p · ( x y ) } δ ( p 2 m 2 ) θ ( p 0 ) = ( i / M + m ) D ( ) ( x y ) i r = + , d p v p , r ( x ) v ¯ p , r ( y ) = i r = + , d p [ ( 2 π ) 3 2 ω p ] 1 v r ( p ) v ¯ r ( p ) exp { i ω p ( x 0 y 0 ) i p · ( x y ) } = i ( i γ M μ μ , x + m ) d p [ ( 2 π ) 3 2 ω p ] 1 exp { i ω p ( x 0 y 0 ) i p · ( x y ) } = i ( 2 π ) 3 ( m + i / M ) d 4 p exp { i p · ( x y ) } δ ( p 2 m 2 ) θ ( p 0 ) = ( i / M + m ) D ( + ) ( x y )
Now a key observation is in order, for the genius of Ettore Majorana to be fully recognized and appreciated. It turns out that only in the Majorana representation of the purely imaginary Dirac matrices, two further complete orthonormal sets of positive and negative frequency solutions of the field equation (122) can be immediately obtained by complex conjugation of the previously obtained ones. As a matter of fact we find
v r * ( p ) 2 m ( 2 ω p + 2 m ) 1 2 E M + ( p ) η r * ( r = ± p 0 = ω p )
which manifestly satisfies
H ˜ M ( p ) v ± * ( p ) = ω p v ± * ( p )
owing to β M η ± * = η ± * and β M ξ ± * = ξ ± * . In a quite analogous way, the spin states
u r * ( p ) 2 m ( 2 ω p + 2 m ) 1 2 E M ( p ) ξ r * ( r = ± p 0 = ω p )
do satisfy
H ˜ M ( p ) u ± * ( p ) = ω p u ± * ( p )
Notice en passant that, e.g., γ M 5 u ± * ( p ) = v ± * ( p ) , whilst for the adjoint spin-states we find u ¯ ± * ( p ) = [ u ± ( p ) β M ] * = u ± ( p ) β M . It turns out that any spin-state is truly orthogonal to its complex conjugate because, e.g.,
u ¯ r ( p ) u s * ( p ) = ξ r ( m + p / ˜ ) β M ( m p / ) ξ s * ( 2 m + 2 ω p ) 1 0 ( r , s = + , )
while the closure relation can be directly obtained from the complex conjugation of eq.s (137)
r = + , u r * ( p ) u ¯ r * ( p ) = r = + , u r * ( p ) u r ( p ) β M = p / + m
It follows therefrom that the most general classical solutions of the real Majorana wave field equation (122) might appear in both real and complex forms: namely,
ψ M ( x ) = r = + , d p a p , r u p , r ( x ) + c . c . = ψ M * ( x )
ψ ^ M ( x ) = r = + , d p b p , r v p , r * ( x ) + c . c . = ψ ^ M * ( x )
Ψ D ( x ) = r = + , d p [ c p , r u p , r ( x ) + d p , r * v p , r ( x ) ]
where all the coefficients a p , ± , b p , ± , c p , ± , d p , ± are Grassmann valued complex numbers. As already noticed, the dynamics and the symmetry properties related to the above real and complex solutions of the Majorana real wave field equations (122) and (124) will be quite different and will be analyzed in the sequel, after the transition to the quantum theory. The complex solutions of the Majorana real wave equation are usually called Dirac bispinors, or spinors, whilst the real solutions are the self-conjugated Majorana bispinors, or spinors.

3.2.4. Canonical Quantum Theory

The transition to the quantum theory is achieved by passing from the complex Grassmann valued c-number coefficients a p , ± , b p , ± , c p , ± , d p , ± to the corresponding creation-annihilation operators, which satisfy the canonical anti-commutation relations, viz.,
{ a p , r , a q , s } = δ ( p q ) δ r s { b p , r , b q , s } = δ ( p q ) δ r s
{ c p , r , c q , s } = δ ( p q ) δ r s { d p , r , d q , s } = δ ( p q ) δ r s
all the other anti-commutators being null. In so doing, it appears that the Majorana quantum fields, which arise out of the real classical solutions ψ M ( x ) and ψ ^ M ( x ) , will become Hermitean and neutral quantum fields, whilst that one arising from Ψ M ( x ) will represent a charged quantum field. Consider in fact the neutral, Hermitean, self-conjugated quantum spinor field
ψ M ( x ) = p , r a p , r u p , r ( x ) + a p , r u p , r * ( x ) = [ ψ M ( x ) ] ψ M c ( x )
ψ ¯ M ( x ) = ψ M ( x ) β M = p , r a p , r u ¯ p , r ( x ) + a p , r u p , r ( x ) β M = ψ M ( x ) β M
in which we have introduced the shorthand notation
p , r d p r = + ,
Here the dagger symbol means simultaneously full Hermitean conjugation in the Fock space, as well as in the complex bispinor part. Let us first evaluate its charge or lepton number
d x ρ ( x ) = d x : ψ M ( x ) ψ M ( x ) : = d x p , r q , s : a p , r u p , r ( x ) + a p , r u p , r ( x ) a q , s u q , s ( x ) + a q , s u q , s * ( x ) : = p , r q , s δ ( p q ) δ r s : a p , r a q , s + a p , r a q , s : = p , r a p , r a p , r a p , r a p , r 0
which endorses the neutral nature of the Hermitean and self-conjugate Majorana spinor quantum field. Needless to say, the above equality is a formal one and deserves a regularization procedure - such as the point-splitting e.g. - but it is eventually and essentially true, thanks to the normal ordering prescription and the canonical anticommutation relations.
Let us determine the canonical anticommutators for the Hermitean Majorana quantum field, which is an operator valued tempered distribution. In order to be as clear as possible, it is better to recall that the Majorana spinor quantum field is a Hermitean operator valued tempered distribution, which also lies in the four dimensional bispinor space ψ M , α ( x ) ( α = 1 , 2 , 3 , 4 ) . More precisely, it turns out that the Majorana bispinor quantum field do actually live in a two dimensional subspace of the four dimensional linear bispinor space, as it satisfies the constraints
E M ψ M ( x ) = 0 E M + ψ M ( x ) = ψ M c ( x )
E M ± ( m ± i / M ) / 2 m
thanks to the basic relations (132) and (133). It is worthwhile to remark that the rank-two differential operators E M ± are real in the Majorana representation, owing to the purely imaginary nature of the gamma matrices, so that
E M ± * = ( m i γ M ν * ν ) / 2 m = E M ±
Then we can write
{ ψ M ( x ) , ψ ¯ M ( y ) } = { ψ M ( x ) , ψ M ( y ) } β M = p , r q , s a p , r u p , r ( x ) + a p , r u p , r * ( x ) , a q , s u q , s ( y ) + a q , s u q , s ( y ) β M = p , r q , s { a p , r , a q , s } u p , r ( x ) u ¯ q , s ( y ) + { a p , r , a q , s } u p , r * ( x ) u q , s ( y ) β M = p , r u p , r ( x ) u p , r ( y ) + u p , r * ( x ) u p , r ( y ) β M = d p [ ( 2 π ) 3 2 ω p ] 1 ( p / + m ) e i p · ( x y ) + ( p / m ) e i p · ( x y ) p 0 = ω p = ( m + i / M ) d p [ ( 2 π ) 3 2 ω p ] 1 e i p · ( x y ) e i p · ( x y ) p 0 = ω p = ( i ) ( m + i / M ) D ( x y ) S M ( x y )
where we have set / M γ M μ μ , while use has been made of the closure relation (143) whereas D ( x y ) is the well known Pauli-Jordan tempered distribution
D ( x y ) = i ( 2 π ) 3 d 4 k exp { i k · ( x y ) } δ ( k 2 m 2 ) sgn ( k 0 )
As a matter of fact, the canonical spinor anticommutator { ψ ( x ) , ψ ¯ ( y ) } = S ( x y ) is generally and uniquely determined by 1. Lorentz covariance 2. the field equation ( m i / ) S ( x y ) = 0 and 3. the general spinor microcausality property
γ 0 S ( x y ) = 0 D ( x y ) ( x y ) 2 < 0 S ( 0 , x y ) = γ 0 δ ( x y )
It is important to emphasize that the above mentioned noteworthy featured properties do not depend at all upon the specific Dirac matrices representations, so that they keep holding true for both Dirac and Majorana spinor quantum fields, in spite of their deep difference in the internal symmetries, the same conclusion being also valid, mutatis mutandis, for the off-shell causal Green’s functions or Feynman propagators. However, it turns out that only for the Majorana spinor we get S M ( x y ) + S M * ( x y ) = 0 , thanks to the Majorana spinor self-conjugation and the purely imaginary nature of the Dirac matrices in the Majorana representation.
A further neutral, self-conjugate, Hermitean quantum Majorana spinor can be obtained out of the classical solution (145), viz.,
ψ ^ M ( x ) = p , r b p , r v p , r * ( x ) + b p , r v p , r ( x ) = [ ψ ^ M ( x ) ] ψ ^ M c ( x )
ψ ^ M ( x ) β M = q , s b q , s v q , s ( x ) β M + b q , s v ¯ q , s ( x )
which again satisfies the canonical anticommutation relations
{ ψ ^ M ( x ) , ψ ^ M ( y ) } β M = S M ( x y ) with { ψ ^ M ( x ) , ψ ¯ M ( y ) } = 0
Consider now the non-Hermitean combination
Ψ D ( x ) = [ ψ M ( x ) + i ψ ^ M ( x ) ] / 2
which still satisfies the Majorana wave equation (122). Its normal modes expansion reads
Ψ D ( x ) = 2 1 / 2 p , r a p , r u p , r ( x ) + i b p , r v p , r * ( x ) + a p , r u p , r * ( x ) + i b p , r v p , r ( x ) = p , r c p , r w p , r ( x ) + d p , r z p , r ( x )
whence the identity follows
a p , r u p , r ( x ) + i b p , r v p , r * ( x ) = 2 c p , r w p , r ( x )
a p , r u p , r * ( x ) + i b p , r v p , r ( x ) = 2 d p , r z p , r ( x )
which is nothing but a local Bogoliubov transformation, as we shall check here below. As a matter of fact, after taking inner products with some suitable spin-state wave functions, it is possible to prove that the creation/annihilation operators c p , r , c p , r , d q , s , d q , s do still satisfy the canonical anticommutation relations, which endorses the local nature of the above Bogoliubov transformation.
Proof. Consider the inner product
d x u q , s ( x ) a p , r u p , r ( x ) + i b p , r v p , r * ( x ) = 2 c p , r d x u q , s ( x ) w p , r ( x )
We find
d x u q , s ( x ) v p , r * ( x ) = ( 2 ω p ) 1 u s ( p ) v r * ( p ) δ ( p q ) = 1 2 ξ s η r * δ ( p q ) ( 2 ω p ) 1 u s ( p ) v r * ( p ) = ( 2 ω p ) 1 ( 2 m + 2 ω p ) 1 ξ s ( m + ω p β M p k γ M k ) ( m + ω p β M + p γ M ) η r * = ( 2 ω p ) 1 ( 2 m + 2 ω p ) 1 ξ s [ 2 ω p 2 + 2 m ω p β M + ω p β M p γ M p k γ M k ω p β M ] η r * = ( 2 ω p ) 1 ( 2 m + 2 ω p ) 1 ξ s [ 2 ω p 2 + 2 m ω p + ω p p γ M p k γ M k ω p ] η r * = 1 2 ξ s η r *
because ξ s β M = ξ s , β M η r * = η r * and ξ s γ M k η r * = ξ s { β M , γ M k } η r * 0 . The two orthogonal pair ξ ± , η ± * of the positive energy bispinors in the Majorana particle rest frame p = 0 are evidently connected by the S O ( 2 , R ) transform (128) so that 1 2 ξ s η r * = 1 2 ξ s [ F + ( θ ) ] r t ξ t f r s ( θ ) in which f r s 1 ( θ ) = f s r ( θ ) = f r s ( θ ) ( r , s = + , ) and the summation over repeated indexes is understood. Then we obtain
a p , r + i b p , t f t r ( θ ) = 2 c p , t f t r ( θ ¯ ) a q , s i b q , t f t s ( θ ) = 2 c q , t f t s ( θ ¯ )
2 c p , r = a p , t f t r ( θ ¯ ) + i b p , t f t r ( θ θ ¯ ) 2 c q , s = a q , u f u s ( θ ¯ ) i b q , u f u s ( θ θ ¯ )
{ c p , r , c q , s } = 1 2 { a p , t f t r ( θ ¯ ) + i b p , t f t r ( θ θ ¯ ) , a q , u f u s ( θ ¯ ) i b q , u f u s ( θ θ ¯ ) } = δ ( p q ) δ r s
Analogous procedures drive to show the full validity of the canonical commutation relation { d p , r , d q , s } = δ ( p q ) δ r s , as well as the vanishing of all the remaining ones.
It follows that a Dirac spinor, which represents a charged and non-Hermitean quantum spinor field, can always be expressed as a complex combination of two independent neutral, self-conjugated and Hermitean Majorana quantum spinor fields, i.e.,
ψ D ( x ) = [ ψ M ( x ) + i ψ ^ M ( x ) ] / 2
From the similarity transformation rule γ μ = S γ M μ S 1 we can go back to the Weyl representation of the gamma matrices and thereby
ψ D ( x ) = p , r c p , r S w p , r ( x ) S 1 + d p , r S z p , r ( x ) S 1
The normal modes - in the Weyl representation of the gamma matrices - are provided by
w p , r ( x ) [ ( 2 π ) 3 2 ω p ] 1 2 w r ( p ) exp { i ω p t + i p · x } z p , r ( x ) [ ( 2 π ) 3 2 ω p ] 1 2 z r ( p ) exp { i ω p t i p · x } ( r = + , )
where
w r ( p ) = ( 2 m + 2 ω p ) 1 / 2 ( m + p / ) ς r z r ( p ) = ( 2 m + 2 ω p ) 1 / 2 ( m p / ) ζ r
where
γ 0 ς ± = ς ± γ 0 ζ ± = ζ ± ( Σ 2 1 ) ς ± = 0 ( Σ 2 ± 1 ) ζ ± = 0
with Σ 2 = i γ 3 γ 1 . The above spin states do satisfy z ± ( p ) = i γ 2 w ± * ( p ) and w ± ( p ) = i γ 2 z ± * ( p ) .
Proof. After setting N = ( 2 m + 2 ω p ) 1 / 2 we find
i γ 2 w ± * ( p ) = i N γ 2 ( m + γ 0 ω p + γ 1 p x + γ 3 p z γ 2 p y ) ς ± * = N ( m p / ) i γ 2 ς ± * γ 0 i γ 2 ς ± * = i γ 2 ς ± * Σ 2 i γ 2 ς ± * = i γ 2 ς ± * i γ 2 ς ± * ζ ± z ± ( p ) = i γ 2 w ± * ( p )
Charge conjugation is conventionally defined as the operation in which particles and antiparticles are interchanged, up to an arbitrary overall phase factor. It follows thereby that if we set
C c p , r C 1 = e i η d p , r C d p , s C 1 = e i η c p , s
r , s = + , p R 3
then the repeated application of the charge conjugation operation yields
C C c p , r C 1 C 1 = e i η C d p , r C 1 = c p , r
so that we can always assume the following properties: namely,
C 2 = I C = C = C 1
It turns out that the above defined spin-states do fulfill the remarkable relationship
z ± ( p ) = i γ 2 w ± * ( p ) w ± ( p ) = i γ 2 z ± * ( p )
Hence, the charge conjugation transformation law can be rewritten as
ψ D c ( x ) = p , r e i η d p , r w p , r ( x ) + c p , r z p , r ( x ) = e i ( η π / 2 ) p , r d p , r γ 2 z p , r * ( x ) + c p , r γ 2 w p , r * ( x ) = e i ( η π / 2 ) γ 2 ψ D ( x )
Thus, if we choose η = π / 2 then we obtain the transformation rule
ψ D c ( x ) = γ 2 ψ D ( x )
which precisely corresponds to the classical transformation ψ c ( x ) = γ 2 ψ * ( x ) .

3.2.5. Observables

Let us now analyze and discuss the observables of a Majorana quantum field and the related properties of the Majorana particles, if eventually detected to any concern in some forthcoming experiment.
First of all, it is worthwhile to notice that from the orthogonality relations
u p , r , u q , s = d x u ¯ p , r ( x ) β M u q , s ( x ) = δ r s δ ( p q ) u p , r , u q , s * = d x u ¯ p , r ( x ) β M u q , s * ( x ) = 0 u p , r * , u q , s * = d x u p , r ( x ) β M u q , s * ( x ) = δ r s δ ( p q )
one can easily obtain all the observable quantities involving the Majorana massive quantum spinor field. For example the energy-momentum tetra-vector takes the form
P μ = i 4 d x : ψ ¯ M ( x ) β M μ ψ M ( x ) : = i 4 d x : q , s a q , s u ¯ q , s ( x ) + a q , s u ¯ q , s * ( x ) × β M μ p , r a p , r u p , r ( x ) + a p , r u p , r * ( x ) : = 1 2 p , r p μ : a p , r a p , r a p , r a p , r : = p , r p μ a p , r a p , r ( p 0 = ω p )
where the normal ordering prescription has been customarily adopted.
Concerning the total angular momentum operator, owing to the rotation invariance of the theory and without loss of generality, one can always arrange the Observer reference frame in such a manner that the direction of the motion of a Majorana particle lies along the O y axis. In such a circumstances, for example, if x ψ M = z ψ M = 0 that implies ψ M ( t , y ) = ψ M ( t , 0 , y , 0 ) , then we obtain for the energy-momentum tensor
2 i T 13 ( t , y ) = ψ ¯ M ( t , y ) γ M 1 z ψ M ( t , y ) = 0
2 i T 31 ( t , z ) = ψ ¯ M ( t , y ) γ M 3 x ψ M ( t , y ) = 0
and thereby
μ M 13 μ = μ S 13 μ = 0
whence it follows that the helicity - the projection of the spin angular momentum on the direction of the motion - is conserved in time. After insertion of the normal modes expansion one gets
h = d y : 1 2 ψ M ( t , y ) Σ M , 2 ψ M ( t , y ) : = d y : p , r a p , r u p , r ( t , y ) + a p , r u p , r * ( t , y ) × 1 2 Σ M , 2 q , s a q , s u q , s * ( t , y ) + a q , s u q , s ( t , y ) :
in which we have set
p = ( 0 , p , 0 ) q = ( 0 , q , 0 ) ω p = p 2 + m 2
u p , r ( t , y ) = [ 4 π ω p ] 1 / 2 u r ( p ) exp { i p y i t ω p } ( r = + , )
the normalization being now consistent with the occurrence that the spinor plane waves are independent of the transverse spatial coordinates x = ( x 1 , x 3 ) . From the commutation relation
[ ω p γ M 0 p γ M 2 , Σ M , 2 ] = 0
together with the definition
u ± ( p ) ( 2 ω p + 2 m ) 1 / 2 m + ω p γ M 0 p γ M 2 ξ ±
it can be readily derived that
( Σ M , 2 1 ) ξ ± = 0 ( Σ M , 2 1 ) u ± ( p ) = 0
which yields in turn
h = 1 2 d p a p , + a p , + a p , a p ,
It follows therefrom that the 1-particle states a p , ± | 0 do represent charge self-conjugated, or neutral, Majorana massive particles with energy-momentum p μ = ( ω p , p ) and two positive/negative helicity states, that has to be eventually detected with equal probability.3
Atiah-Singer index. We recall that, only for massless self-conjugated spinors, there exists an accidental U(1) internal symmetry for the Majorana quantum spinor field, the corresponding conserved pseudo-scalar charge being expressed by
Q 5 = ± 1 2 d x : ψ ¯ M ( t , x ) β M γ M 5 ψ ( t , x ) : Q ˙ 5 = 0
the overall sign being conventional and irrelevant. After insertion of the normal mode expansion we get
Q 5 = ± d x : q , s a q , s υ ¯ q , s ( x ) + a q , s υ ¯ q , s * ( x ) × 1 2 β M γ M 5 p , r a p , r υ p , r ( x ) + a p , r υ p , r * ( x )
where we have denoted by
υ r ( p ) = p / M υ r / 2 ( | p | , r = + , ) υ + = 1 i 1 i υ = i 1 i 1 υ r , p ( x ) = [ ( 2 π ) 3 2 ] 1 2 υ r ( p ) exp { i t + i p · x }
the spin-states and plane wave-functions respectively of the massless self-conjugated spinor field. The spin-states satisfy the conventional orthogonality and closure relations as they correspond to the regular limit lim m 0 w r ( p ) = υ r ( p ) . Now we have
γ M 5 = 0 0 0 i 0 0 i 0 0 i 0 0 i 0 0 0
so that
γ M 5 υ ± = ± υ ± γ M 5 υ ± ( p ) = υ ± ( p )
which means that the spin-states in the massless case are eigenstates of the chiral matrix in the Majorana representation.
It turns out that in the 1d case, i.e. the dynamics along e.g. the O y axis, which is the relevant situation for the construction of the helicity operator, we actually have
υ ± ( p y ) = ( β M γ M 2 p y ) υ ± / 2 [ = p y sgn ( p y ) ]
in such a manner that the above spin-states are common eigenstates of the chiral and spin matrices because
[ β M γ M 2 p y , Σ M 2 ] = 0 [ Σ M 2 , γ M 5 ] = 0
Now we obtain
d x υ r , p ( , x ) γ M 5 υ s , q * ( t , x ) = 1 2 υ r ( p ) γ M 5 υ s * ( p ) exp { 2 i t } υ r ( p ) γ M 5 υ s * ( p ) = p ˜ 2 2 υ r γ M 5 υ s * 0 ( r , s = + , p ˜ ν = p ν )
Moreover we find
d x υ r , p ( , x ) γ M 5 υ s , q ( t , x ) = 1 2 υ r γ M 5 υ s = 1 2 δ r s ( δ r δ s + ) ( r , s = + , )
whence it follows, after a suitable choice of the overall sign, that the above introduced pseudo-scalar charge Q 5 is nothing but the quantum counterpart of the Atiah-Singer index for a Majorana spinor field
Q 5 = ν M d p a p , + a p , + a p , a p ,
i.e. the difference between the number of quanta with positive and negative helicity.
We can therefore conclude that the 1-particle states a p , ± | 0 do represent neutral, i.e. charge self-conjugated, spin 1 2 Majorana particles, massive or massless, with energy-momentum p μ = ( ω p , p ) and two equally probable positive/negative helicity states.
From the Majorana Lagrangian and related wave equation, one can immediately realize that the special distributions for the Majorana neutral spinor field are the very same as for a Dirac charged quantum field: namely,
ψ M ( x ) , ψ ¯ M ( y ) = S ( x y ; m ) 0 | T ψ M ( x ) ψ ¯ M ( y ) | 0 = S F ( x y ; m )
the massless limit being smooth.

4. Dirac, Majorana and Weyl Fermions in 4D

In this Section we briefly summarize and highlight the main properties of the Dirac, Weyl and Majorana fermion fields that we have thoroughly discussed at length in the previous Sections.
Let us start from a few basic definitions and properties of spinors on a 4d Minkowski space. A 4-component Dirac fermion ψ under a Lorentz transformation transforms as
ψ ( x ) ψ ( x ) = exp 1 2 λ μ ν Σ μ ν ψ ( x ) ,
for x μ = Λ μ ν x ν . Here λ μ ν + λ ν μ = 0 are six real canonical coordinates for the Lorentz group, Σ μ ν are the generators in the 4d reducible representation of Dirac bispinors, while Λ μ ν are the Lorentz matrices in the irreducible vector representation D ( 1 2 , 1 2 ) . The invariant kinetic Lagrangian for a free Dirac field is
i ψ ¯ γ μ μ ψ = i ψ ¯ γ · ψ
where ψ ¯ = ψ γ 0 . A Dirac fermion admits a Lorentz invariant mass term m ψ ¯ ψ .
A Dirac bispinor can be seen as the direct sum of two Weyl spinors4
ψ L = P L ψ , ψ R = P R ψ , where P L = 1 γ 5 2 , P R = 1 + γ 5 2
with opposite chiralities
γ 5 ψ L = ψ L , γ 5 ψ R = ψ R .
A left-handed Weyl fermion admits a Lagrangian kinetic term
i ( ψ L , γ · ψ L ) ψ ¯ L γ μ μ ψ L
but not a mass term, because ( ψ L , ψ L ) = 0 , since γ 5 γ 0 + γ 0 γ 5 = 0 . Hence, a Weyl fermion is massless and this property is protected by chirality conservation.
For classical Majorana spinors we need the notion of Lorentz covariant charge conjugate spinor ψ c
ψ c = C ψ C 1 = η C γ 2 ψ *
where η C is a phase which, for simplicity, in the sequel we set equal to 1. It is not hard to show that if ψ transforms like (169), then
ψ c ( x ) exp 1 2 λ μ ν Σ μ ν ψ c ( x ) .
Therefore one can impose on ψ the condition
ψ = ψ c
because both sides transform in the same way. By definition, a spinor satisfying (174) is a Majorana spinor. As we have seen before it admits both kinetic and mass term.
The Lorentz group spinor representation properties are as follows: the chiral matrix γ 5 commutes with Lorentz transformations matrix exp 1 2 λ μ ν Σ μ ν and so do P L and P R . This means that the Dirac representation is reducible. Multiplying the spinors by P L and P R selects the Weyl irreducible representations. To state it more precisely, Weyl representations are irreducible representations of the group S L ( 2 , C ) , which is the covering group of the proper ortochronous Lorentz group. They are usually denoted ( 1 2 , 0 ) and ( 0 , 1 2 ) in the S U ( 2 ) × S U ( 2 ) notation of the S L ( 2 , C ) irreducible representations. As we have seen in (173), Lorentz transformations commute also with the charge conjugation operation (172). This also implies that Dirac spinors are reducible and suggests another possible reduction: by imposing (174) we single out another irreducible representation, the Majorana one. The Majorana representation is the minimal irreducible representation of a (one out of eight) covering of the complete Lorentz group [17,18]. It is evident, and well-known, that the Majorana and Weyl representations in 4d are incompatible.
Let us consider next the charge conjugation and parity discrete symmetries and recall the relevant properties of a Weyl spinor. We have
C ψ L C 1 = P L C ψ C 1 = P L ψ c = ( ψ c ) L
The charge conjugate of a Majorana field is, by definition, itself. While the Action of a Majorana field is invariant under charge conjugation, for a Weyl spinor we have [  γ μ μ γ ·  ]
C ψ L ¯ i γ · ψ L C 1 = ( ψ c ) L ¯ i γ · ψ L c = ψ R ¯ i γ · ψ R
i.e. a Weyl fermion is, so to say, maximally non-invariant.
The parity operation is defined by
P ψ L ( t , x ) P 1 = η P γ 0 ψ R ( t , x )
where η P is a phase, which we set to 1. In terms of the Action we have
P ψ L ¯ i γ · ψ L P 1 = ψ R ¯ i γ · ψ R
For a Majorana spinor the Action is parity invariant.
This also suggests a useful representation for a Majorana fermion. Let ψ R = P R ψ be a generic Weyl spinor: as seen above, we have P R ψ R = ψ R and it is easy to prove that P L ( ψ R ) c = ( ψ R ) c , i.e. ( ψ R ) c is left-handed. Therefore the sum ψ M = ψ R + ( ψ R ) c is a Majorana spinor because it satisfies (174) and any Majorana fermion can be represented in this way.
Considering next the CP discrete symmetry. From the above it follows that the Action of a Majorana spinor is obviously invariant under it. On the other hand for a Weyl spinor we have
CP ψ L ( t , x ) ( CP ) 1 = γ 0 P R ψ c ( t , x ) = γ 0 ( ψ c ) R ( t , x )
Applying now the CP transformation to the Action for a Weyl spinor, one gets
CP ψ L ¯ i γ · ψ L ( CP ) 1 = ( ψ c ) R ¯ ( t , x ) i γ · ( ψ c ) R ( t , x ) = ( ψ c ) R ¯ ( t , x ) i γ · ( ψ c ) R ( t , x )
But one can prove as well that
( ψ c ) R ¯ ( t , x ) i γ · ( ψ c ) R ( t , x ) = ψ L ¯ ( x ) i γ · ψ L ( x )
It follows therefore that the Action for a Weyl spinor is CP invariant. Moreover, as we have seen in the previous Section, it is also, separately, T invariant and so CPT invariant.
The transformation properties of the Weyl and Majorana spinor fields are summarized in the following table:
Majorana Weyl P : P ψ ( t , x ) P 1 = γ 0 ψ ( t , x ) P ψ L ( t , x ) P 1 = γ 0 ψ R ( t , x ) C : C ψ C 1 = γ 2 ψ * = ψ C ψ L C 1 = P L ψ c = ( ψ c ) L CP : CP ψ ( t , x ) ( CP ) 1 = ψ ( t , x ) CP ψ L ( t , x ) ( CP ) 1 = γ 0 ( ψ c ) R ( t , x )
In the conventional Standard Model of the fundamental interactions the neutrino and antineutrino fields are treated as massless chiral Weyl bispinors: according to our conventions and notations we can set
ψ ( x ) = ν L ( x ) + ν R ( x ) = d p c , p u , p ( x ) + d + , p v + , p ( x ) + ν R ( x )
where u , v + are fixed and independent left- and right-handed spinors that, without loss of generality, we can identify with those ones of eq. (36). It is worthwhile to stress that, in spite of being a spin 1 2 field, the left-handed neutrino field necessarily exhibits positive frequency plane waves with only negative helicity and negative frequency plane waves with solely positive helicity. Contrary to some popular belief, this characteristic striking constraint for the Weyl spinors ψ is not merely due to their masslessness, because Majorana and Dirac massless spinors do behave quite differently.
We can interpret the Fourier expansion (183) of the Standard Model neutrino field as follows: c , p annihilates a left-handed particle with negative helicity, while d + , p creates a right-handed anti-particle with positive helicity. Needless to say, it appears to be evident that in this game right-handed neutrinos ν R and left-handed anti-neutrinos play only the role of sterile and non-interacting spectators, this role being confined to provide a consistent causal Green’s function: it turns out that the latter have never been detected neither in laboratory experiments, nor in astrophysical or cosmological observations.

5. Regularizations for Weyl Spinors

The classical Lagrange density for a Weyl (left) spinor in the four component formalism
ψ ( x ) = χ L ( x ) = χ ( x ) 0
reads
K ( x ) = ψ ¯ ( x ) i / ψ ( x ) = χ L ( x ) α ν i ν χ L ( x )
It follows that the corresponding matrix valued Weyl differential operator
w L α ν i ν P L [ P L = 1 2 ( I γ 5 ) ]
is singular and does not possess any rank-four inverse. After minimal coupling with a real massless vector field A μ ( x ) we come to the classical Lagrangian
L = χ L α ν i ν χ L + g A ν χ L α ν χ L 1 4 F μ ν F μ ν
where 0 g < 1 is a non-negative small parameter, while F μ ν = μ A ν ν A μ . It turns out that the classical Action
S = d 4 x L
is invariant under the Poincaré group, as well as under the internal U(1) phase transformations χ L ( x ) e i g θ χ L ( x ) . The Action is invariant under the so called scale or dilatation transformations, viz.,
x μ = e ϱ x μ χ L ( x ) = e 3 2 ϱ χ L ( e ϱ x ) A μ ( x ) = e ϱ A μ ( e ϱ x )
with ϱ R , as well as with respect to the local phase or gauge transformations
χ L ( x ) = e i g θ ( x ) χ L ( x ) A ν ( x ) = A ν ( x ) + ν θ ( x )
which amounts to the ordinary U(1) phase transform in the limit of constant phase. It follows therefrom that there are twelve conserved charges in this model at the classical level and, in particular, owing to scale and gauge invariance, no mass term is allowed for both spinor and vector fields. The question naturally arises if those symmetries hold true after the transition to the quantum theory and, in particular, if they are protected against loop radiative corrections within the perturbative approach. Now, in order to develop perturbation theory, one has to face the problem of the lack of an inverse for both the Weyl and gauge fields, owing to chirality and gauge invariance. In order to solve this problem, it is expedient to add to the Lagrangian non-interacting terms, which are fully decoupled from any physical quantity. These terms break chirality and gauge invariance, albeit in a harmless way, just to allow the setting up of a Feynman propagator, or causal Green’s functions, for both the Weyl and gauge quantum fields. The simplest choice, which keeps Poincaré and internal U(1) phase change symmetries, is provided by
L = φ R α ν i ν φ R 1 2 ( · A ) 2
where
ψ ( x ) = φ R ( x ) = 0 φ ( x )
is a left-chirality breaking right-handed Weyl spinor field. Notice en passant that the modified Lagrangian L + L does exhibit the further U(1) internal symmetry under the so called chiral phase transformations
ψ ( x ) = ( cos θ + i g sin θ γ 5 ) ψ ( x ) ψ ( x ) = χ ( x ) φ ( x )
so that the modified theory involves another conserved charge at the classical level. From the modified Lagrange density we get the Feynman propagators for the massless Dirac field ψ ( x ) , as well as for the massless vector field in the so called Feynman gauge: namely,
S ( p ) = i p / p 2 + i ε D μ ν ( k ) = i g μ ν k 2 + i ε
and the vertex i g γ ν P L , with k + p q = 0 , which involves a vector particle of momentum k and a Weyl pair of particle and anti-particle of momenta p and q respectively and of opposite helicity.5
The lowest order 1-loop correction to the Weyl kinetic term p / P L is formally provided by the Feynman rules in the Minkowski space, viz.,
Σ 2 ( p / ) = i g 2 d 4 ( 2 π ) 4 γ μ D μ ν ( p ) S ( ) γ ν P L
By naïve power counting the above 1-loop integral turns out to be UV divergent. Hence a regularization procedure is mandatory to give a meaning and evaluate the radiative correction Σ 2 ( p / ) to the Weyl kinetic operator. Here in the sequel we shall examine in detail the dimensional, Pauli-Villars and UV cut-off regularizations.

5.0.1. Dimensional Regularization

In a 2 ω dimensional space-time the radiative correction to the Weyl kinetic term takes the form
reg Σ 2 ( p / ) = i g 2 μ 2 ϵ d 2 ω ( 2 π ) 2 ω D μ ν ( ) γ μ S ( + p ) γ ν P L
where ϵ = 2 ω > 0 is the shift with respect to the Minkowski space. Since the above expression is traceless and has the canonical engineering dimension of a mass in natural units, it is quite apparent that the latter cannot generate any mass term, which would be proportional to the unit matrix. Hence, mass is forbidden and it is convenient to set and evaluate
reg Σ 2 ( p / ) f ( p 2 ) p / P L tr [ p / reg Σ 2 ( p / ) ] = 1 2 2 ω p 2 f ( p 2 ) tr [ p / reg Σ 2 ( p / ) ] = g 2 μ 2 ϵ ( 2 π ) 2 ω d 2 ω ( i ) tr p / γ λ / γ λ P L ( p ) 2 + i ε 2 + i ε
Concerning dimensional regularization, we collect here below the definitions and key properties for the 2 ω × 2 ω γ matrices in a 2 ω dimensional space-time with a Minkowski signature
γ μ = γ ¯ μ μ = 0 , 1 , 2 , 3 γ ^ μ μ = 4 , , 2 ω 4 γ ¯ μ , γ ¯ ν = 2 g ¯ μ ν I γ ^ μ , γ ^ ν = 2 g ^ μ ν I γ ¯ μ , γ ^ ν = 0 g ¯ = diag ( + , , , ) g ^ = I ^ γ 5 i γ ¯ 0 γ ¯ 1 γ ¯ 2 γ ¯ 3 γ 5 2 = I γ ¯ μ , γ 5 = 0 = γ ^ μ , γ 5
Taking all the above listed equations into account, it is not difficult to check the following trace formulas: namely,
tr γ μ γ ν = g μ ν tr I = 2 ω g μ ν 2 ω tr γ κ γ λ γ μ γ ν = g κ λ g μ ν g κ μ g λ ν + g κ ν g λ μ tr γ ¯ κ γ ¯ λ γ ^ μ γ ^ ν = 2 ω g ¯ κ λ g ^ μ ν tr γ 5 γ ¯ μ γ ¯ λ γ ¯ ρ γ ¯ ν = i 2 ω ϵ μ λ ρ ν
Traces involving an odd number of Dirac matrices do vanish. Then we get
tr p / γ λ / γ λ P L = 2 ω ( ϵ 1 ) p ·
and thereby
p 2 f ( p 2 ) = i g 2 μ 2 ϵ ϵ 1 ( 2 π ) 2 ω 2 p · d 2 ω ( p ) 2 + i ε 2 + i ε
Turning to the Feynman parametric representation we obtain
p 2 f ( p 2 ) = i g 2 μ 2 ϵ ϵ 1 ( 2 π ) 2 ω 0 1 d x 2 p · d 2 ω 2 2 x p · + x p 2 + i ε 2
Completing the square in the denominator and after shifting the momentum x p , dropping the linear term in in the numerator owing to symmetric integration, we have
f ( p 2 ) = 2 i g 2 μ 2 ϵ ϵ 1 ( 2 π ) 2 ω 0 1 d x x d 2 ω 2 + x ( 1 x ) p 2 + i ε 2
When p 2 < 0 one can perform the Wick rotation and readily get the result
f ( p 2 ) = 2 g 2 μ 2 ϵ ϵ 1 ( 4 π ) ω 0 1 d x x 0 d τ τ ϵ 1 e τ x ( 1 x ) p E 2 = 2 g 4 π 2 [ Γ ( ϵ ) Γ ( 1 + ϵ ) ] 4 π μ 2 p 2 ϵ 0 1 d x x 1 ϵ ( 1 x ) ϵ = 2 g 4 π 2 [ Γ ( ϵ ) Γ ( 1 + ϵ ) ] 4 π μ 2 p 2 ϵ B ( 2 ϵ , 1 ϵ )
Expansion around ϵ = 0 yields
f ( p 2 ) = g 4 π 2 1 ϵ C 1 + 1 + ϵ ln 4 π μ 2 p 2 + × 1 ϵ [ ψ ( 2 ) + ψ ( 1 ) 2 ψ ( 3 ) ] + = g 4 π 2 1 ϵ + 1 + 3 C + ln 4 π μ 2 p 2 + evanescent

5.0.2. Pauli-Villars Regularization

Let us now repeat the very same calculation for the left Weyl spinor self-energy in the Pauli-Villars regularization. The latter is simply implemented by the following replacement of the massless Dirac propagator
reg Σ 2 ( p / ) = i g 2 d 4 ( 2 π ) 4 γ μ D μ ν ( p ) s = 0 S C s S ( , M s ) γ ν P L
where M 0 = 0 , C 0 = 1 while { M s λ s M | λ s 1 ( s = 1 , 2 , , S ) } is a collection of very large auxiliary masses. The set of constants C s will be suitably selected, as we shall see in the sequel, in such a manner to obtain a specific and mathematically meaningful form for the ultraviolet divergences that will manifest themselves in the limit λ s . Since we have
S ( , M s ) = i ( / + M s ) 2 M s 2 + i ε ( s = 1 , 2 , , S )
it is now convenient to set
reg Σ 2 ( p / ) f ( p 2 ) p / P L M υ ( p 2 )
in such a manner that we can write
tr reg Σ 2 ( p / ) = 4 M υ ( p 2 ) = i g 2 d 4 ( 2 π ) 4 s = 1 S C s 8 M s [ ( p ) 2 + i ε ] 2 M s 2 + i ε tr p / reg Σ 2 ( p / ) = 2 p 2 f ( p 2 ) = i g 2 d 4 ( 2 π ) 4 s = 0 S C s 4 p · [ ( p ) 2 + i ε ] 2 M s 2 + i ε
Once again, if we take profit of the Feynman parametric formula we obtain
υ ( p 2 ) = d 4 ( 2 π ) 4 s = 1 S C s 0 1 d x 2 i g 2 λ s 2 2 x p · + x p 2 ( 1 x ) M s 2 + i ε 2 = 0 1 d x d 4 ( 2 π ) 4 s = 1 S 2 i g 2 C s λ s 2 + x ( 1 x ) p 2 ( 1 x ) M s 2 + i ε 2 = 2 g 2 0 1 d x d 4 E ( 2 π ) 4 s = 1 S C s λ s E 2 + x ( 1 x ) p E 2 + ( 1 x ) M s 2 2
In order to proceed further on it is convenient to consider the generating integral
I n ( z E ) ( 1 ) n d 4 E ( 2 π ) 4 s = 1 S C s λ s exp { i E · z E } ( E 2 + Δ s 2 ) n
where Δ s 2 = ( 1 x ) [ x p E 2 + M 2 λ s 2 ] > 0 , s = 1 , 2 , , S . Explicit calculation yields
I n ( z E ) = ( 1 ) n 8 π 2 Γ ( n ) s = 1 S C s λ s 2 Δ s | z E | 2 n K 2 n ( Δ s | z E | )
where z E = ( z , z 4 ) while | z E | = z 2 + z 4 2 , whence we immediately obtain
I 2 ( z E ) = 1 8 π 2 s = 1 S C s λ s K 0 ( z s ) ( z s = Δ s | z E | )
where K 0 is the modified Bessel function of the third kind and of order zero, the series representation of which is provided by
K 0 ( z ) = k = 0 1 ( k ! ) 2 z 2 2 k ψ ( k + 1 ) ln z 2 = ln 2 z C + O ( z 2 ln z )
Then we can write
1 4 M tr reg Σ 2 ( p / ) = υ ( p 2 ) = 2 g 2 0 1 d x lim z E 0 I 2 ( z E )
whence it is clear that the limit exists, so that the regularization works, iff the Pauli-Villars condition is fulfilled, viz.,
s = 1 S C s λ s = 0
which provides in turn υ = 0 , i.e., no mass term is allowed in the 1-loop Weyl kinetic term, both in dimensional and Pauli-Villars regularizations, as expected.
Let us turn now to the other form factor
tr p / reg Σ 2 ( p / ) = 2 p 2 f ( p 2 ) = i g 2 d 4 ( 2 π ) 4 s = 0 S C s 4 p · [ ( p ) 2 + i ε ] 2 M s 2 + i ε = d 4 ( 2 π ) 4 s = 0 S C s 0 1 d x 4 i g 2 p · 2 2 x p · + x p 2 ( 1 x ) M s 2 + i ε 2 = 4 i g 2 0 1 d x d 4 ( 2 π ) 4 s = 0 S C s x p 2 2 + x ( 1 x ) p 2 ( 1 x ) M s 2 + i ε 2 = 2 g 2 p 2 0 1 d x x d 4 E ( 2 π ) 4 s = 0 S C s E 2 + x ( 1 x ) p E 2 + ( 1 x ) M s 2 2
whence we obtain
f ( p 2 ) = g 2 0 1 d x x d 4 E ( 2 π ) 4 s = 0 S C s E 2 + Δ s 2 2 = g 2 0 1 d x x d 4 E ( 2 π ) 4 s = 0 S C s 0 d τ τ exp τ E 2 τ Δ s 2 = g 2 16 π 2 0 1 d x x 0 d τ τ s = 0 S C s exp τ ( 1 x ) [ x p E 2 + M 2 λ s 2 ] = g 2 16 π 2 0 1 d x x 0 d τ τ exp τ x ( 1 x ) p 2 × s = 0 S C s 1 n = 1 τ n n ! ( 1 x ) n M 2 n λ s 2 n
Integrability requires the further Pauli-Villars condition
s = 1 S C s = 1
and consequently
f ( p 2 ) = g 2 16 π 2 0 1 d x x n = 1 1 n ! ( 1 x ) n M 2 n s = 1 S C s λ s 2 n × 0 d τ τ τ n exp τ x ( 1 x ) p 2 = g 2 16 π 2 s = 1 S C s 0 1 d x x n = 1 1 n λ s 2 M 2 x p 2 n = g 2 16 π 2 s = 1 S C s 0 1 d x x ln 1 + λ s 2 M 2 x p 2 = g 2 16 π 2 0 1 d x x s = 1 S C s ln λ s 2 x p 2 M 2 ln x p 2 M 2 = g 2 16 π 2 1 2 ln p 2 M 2 0 1 d x x ln x + 1 2 s = 1 S C s ln λ s 2 + 0 1 d x x s = 1 S C s ln 1 x p 2 M s 2 = g 4 π 2 s = 1 S C s ln λ s + 1 4 + 1 2 ln M 2 p 2 + evanescent g 4 π 2 1 ϵ + 1 + 3 C + ln 4 π μ 2 p 2 + evanescent
where, in the last line, a comparison with the result in dimensional regularization has been exhibited. It follows that by identifying
1 ϵ = s = 1 S C s ln λ s
in which the two Pauli-Villars conditions6 hold true, viz.,
s = 1 S C s = 1 s = 1 S C s λ s = 0
the coefficient of the divergent part turns out to be the same for both regulators, as expected.

5.0.3. Cut-Off Regularization

Let us repeat once more the one loop calculation of the Weyl spinor self-energy with a physical very large cut-off regulator, e.g. K = l P 1 with l P = G N / c 3 the Planck length. Of course, once again, no mass term is allowed for the trace of an odd number of gamma matrices is null. Moreover, for the previously introduced coefficient of the Weyl kinetic term we find with μ = ( 0 , )
p 2 f ( p 2 ) = i g 2 ( 2 π ) 4 2 p · θ K 2 2 d 4 ( p ) 2 + i ε 2 + i ε = i g 2 ( 2 π ) 4 0 1 d x 2 p · θ K 2 2 d 4 2 2 x p · + x p 2 + i ε 2 = i g 2 ( 2 π ) 4 0 1 d x d θ K 2 2 × lim η 0 d d η d 0 2 p 0 0 2 p · 2 2 x p · + x p 2 η + i ε
After translation of the integration variable x p , for a suitably large cut-off K and space-like momentum p 2 < 0 we get
p 2 f ( p 2 ) = i g 2 ( 2 π ) 4 0 1 d x d θ K 2 ( + x p ) 2 × lim η 0 d d η d 0 2 p 0 0 2 p · + 2 x p 2 2 + x ( 1 x ) p 2 η + i ε = i g 2 ( 2 π ) 4 0 1 d x d θ 4 K 2 2 × lim η 0 d d η d 0 2 x p 2 0 2 2 + x ( 1 x ) p 2 η + i ε = g 2 p 2 ( 2 π ) 3 0 1 d x x d θ 4 K 2 2 × lim η 0 d d η 2 x ( 1 x ) p 2 η 1 / 2 f ( p 2 ) = g 2 4 π 2 0 1 d x x 0 2 K d 2 2 x ( 1 x ) p 2 3 / 2 = g 2 4 π 2 0 1 d x x 0 1 d y y 2 y 2 x ( 1 x ) p 2 4 K 2 3 / 2
Consider u = a + y 2 , a = x ( 1 x ) ( p 2 / 4 K 2 ) > 0 , then from GR2.272 4. p. 86 we get the elementary result
0 1 d y y 2 u 3 = y u + ln ( y + u ) 0 1 = 1 a + 1 + ln 1 + a + 1 ln a
For a very large cut-off K we obtain
f ( p 2 ) = g 2 4 π 2 0 1 d x x 1 + ln 2 + 1 2 ln 4 K 2 p 2 1 2 ln x ( 1 x ) + = g 4 π 2 ln 4 K 2 p 2 + ln 4 + evanescent
To sum up, we have verified that the 1-loop correction to the (left) Weyl spinor self-energy has the general form, which is universal, i.e. regularization independent: namely,
reg Σ 2 ( p / ) f ( p 2 ) p / P L f ( p 2 ) : = g 4 π 2 1 ϵ + 1 + 3 C + ln 4 π μ 2 p 2 ( DR ) : = g 4 π 2 s = 1 S C s ln λ s + 1 4 + 1 2 ln M 2 p 2 ( PV ) : = g 4 π 2 ln 4 K 2 p 2 + ln 4 ( CUT OFF )
Some final remarks and comments are in order.
1.
In the present model of a left Weyl spinor minimally coupled with a gauge vector potential, no mass term can be generated by the radiative corrections in any regularization scheme. The left-handed part of the classical kinetic term does renormalize, while its right-handed part does not undergo any radiative correction and keeps on being free. The latter has to be necessarily introduced in order to define a Feynman propagator for the massless spinor field, like the gauge fixing term to invert the kinetic term of the gauge potential. The (one loop) renormalized Lagrangian for a Weyl fermion minimally coupled with a gauge vector potential has the universal - i.e. regularization independent - form
L 1 loop = L = χ L α ν i ν χ L + g A ν χ L α ν χ L 1 4 F μ ν F μ ν + φ R α ν i ν φ R 1 2 ( · A ) 2 ( Z 3 1 ) 1 4 F μ ν F μ ν + ( Z 2 1 ) χ L α ν i ν χ L + ( Z 1 1 ) g A ν χ L α ν χ L ( Z 2 1 ) = g 4 π 2 1 ϵ + F 2 ( ϵ , p 2 / μ 2 ) + = g 4 π 2 s = 1 S C s ln λ s + F ˜ 2 ( λ s , p 2 / M 2 ) + = g 4 π 2 ln 4 K 2 p 2 + ln 4 + F ^ 2 ( K 2 / p 2 ) +
where conventional notations have been employed. Notice that, as usual, the arbitrary finite parts F 2 , F ˜ 2 , F ^ 2 of the countertems are analytic for ϵ 0 and λ s , K , respectively, and have to be uniquely fixed by the renormalization prescription.
2.
The interaction definitely preserves left chirality and scale invariance of the counterterms in the transition from the classical to the (perturbative) quantum theory: no mass coupling between the left-handed (interacting) Weyl spinor χ L and right-handed (free) Weyl spinor φ R can be generated by radiative loop corrections.
3.
While the cut-off and dimensional regularized theory does admit a local formulation in D = 4 or D = 2 ω space-time dimensions, there is no such local formulation for the Pauli-Villars regularization. The reason is that the PV spinor propagator
s = 0 S C s S ( , M s )
where M 0 = 0 , C 0 = 1 while { M s λ s M | λ s 1 ( s = 1 , 2 , , S ) } , cannot be the inverse of any local differential operator of the Calderon-Zygmund type. Hence, there is no local Action involving a bilinear spinor term that can produce, after a suitable causal inversion, the Pauli-Villars regularized spinor propagator.
One such pitfall (that is confusing the inverse of a sum of differential operators with the sum of their inverses) is not rare in the literature and can be found in several approaches nominally utilizing a PV regularization. They have nothing to do with the Pauli-Villars regularization in the present calculation, notwithstanding their misleading denomination.
Now the main point concerning the chiral anomaly and the unitarity of the collision operator. As we have seen in details, in the present model of a left-handed Weyl spinor interacting with a gauge vector potential, all the symmetries of the classical Action do survive the divergent and finite parts of the counter-terms, for any choice of the regulators and renormalization prescription. As it is well-known since a long time, this is not so for the famous chiral left-handed triangle diagram, which turns out to be UV finite and yields
Γ ˜ μ 1 μ 2 μ 3 ( L ) ( k 1 , k 2 , k 3 ) 1 loop = g 3 12 π 2 ε σ μ 1 μ 2 μ 3 k 1 σ k 2 σ k 3 2 I 12 ( k 1 , k 2 , k 3 ) + cyclic permutations + g 3 4 π 2 ε σ τ μ 1 μ 2 k 1 σ k 2 τ k 3 μ 3 I 12 ( k 1 , k 2 , k 3 ) + cyclic permutations ( k 1 + k 2 + k 3 = 0 ) I r s ( k 1 , k 2 , k 3 ) = 0 1 d x 1 0 1 d x 2 0 1 d x 3 2 x r x s δ ( 1 x 1 x 2 x 3 ) k 1 2 x 2 x 3 + k 2 2 x 3 x 1 + k 3 2 x 1 x 2 ( r , s = 1 , 2 , 3 )
The one loop chiral left-handed triangular amplitude does satisfy the anomalous Ward identity
k 1 σ Γ ˜ σ μ 2 μ 3 ( L ) ( k 1 , k 2 , k 3 ) = g 3 12 π 2 ε σ τ μ 2 μ 3 k 2 σ k 3 τ
and related cyclic permutations.
The above expression for the 1-loop triangular amplitude is finite, it does not depend at all upon regulators, or counterterms or even renormalization prescriptions and it does not undergo any higher loop corrections (Adler-Bardeen theorem).However, it turns out that, owing to the anomalous Ward identity, the left-handed Weyl current does no longer satisfy the continuity equation, at variance with the classical case, viz.,
ν χ L α ν χ L = 1 3 g 4 π 2 F * μ ν F μ ν F * μ ν = 1 2 ε μ ν ρ σ F ρ σ
This finite and universal quantum breaking of the U(1) symmetry under phase transformations on the left Weyl spinor field does actually put in jeopardy the unitarity issue. As a matter of fact the Lorentz invariant quantum theory of a gauge vector field unavoidably involves a Fock space of states with indefinite norm. Now, in order to select a physical Hilbert subspace of the Fock space a subsidiary condition is necessary. In the Abelian case, when the fermion current satisfies the continuity equation then the equations of motion lead to ( · A ) = 0 , so that a subspace of states of non-negative norm can be selected through the auxiliary condition
· A ( ) ( x ) | phys = 0
A ( ) ( x ) being the destruction, positive frequency part of a d’Alembert quantum field. Conversely, in the present chiral model we find
g 1 ( · A ) = ν χ L α ν χ L = 1 3 g 4 π 2 F * μ ν F μ ν 0
in such a manner that nobody knows how to select a physical subspace of states with non negative norm, if any, where a unitary restriction of the collision operator S could be defined.

6. Weyl and Dirac Propagators

Equipped with the previous definitions, clarifications and results, in the rest of the paper we cope with some critical issues one comes across in quantum field theories involving Weyl fermions. We start with the fermion functional integral.
We denote by D the Dirac operator: namely, the massless matrix-valued differential operator applied in general to Dirac spinors on the 4d curved space with Minkowski signature ( + , , , )
D = i + V
where V μ is any real vector potential, valued in a some Lie algebra. We understand a spin connection when in presence of a non-trivial background metric. The Dirac fermions are represented in the four component formalism. The functional integral, i.e. the effective Action for a quantum Dirac spinor in the presence of a classical background vector potential
Z [ V ] = D ψ D ψ ¯ e i d 4 x g ψ ¯ D ψ
is formally understood as the determinant of D :   det ( D ) = i λ i . Concretely, the latter can be operationally defined in two alternative ways: either in perturbation theory, i.e. as the sum of an infinite number of 1-loop Feynman diagrams, some of which containing UV divergences by naive power counting, or in a non-perturbative approach, i.e. as the suitably regularized infinite product of the eigenvalues of D by means of the analytic continuation tool. In the perturbative approach, in order to assign a meaning to a finite number of UV divergent 1-loop diagrams by naive power counting, one needs a regulator (dimensional, cutoff, PV or others). In the non-perturbative framework the complex power construction and the analytic continuation tool, if available, may provide by themselves the necessary to set up for the infinite product of eigenvalues of a normal operator, without need of any further regulator, in such a manner to produce the Heat-Kernel expansion and related topics [22,23].
Any variation of (186) with respect to a parameter contained in the Dirac operator requires the existence of an inverse of the kinetic operator. Formally
δ det ( D ) = det ( D ) tr ( D 1 δ D )
In other words the existence of the inverse of D is essential for the definition itself of the functional integral. It turns out that, in the case of Dirac fermions, an inverse of D does exist and, if full causality is required in forwards and backwards time evolution on e.g. Minkowski space, it is the Feynman propagator or Schwinger distribution S , which is unique and characterized by the well-known Feynman prescription, which satisfies
D x S ( x y ) = δ ( x y ) , D S = 1 ,
The second equation is a synthetic notation7.
An example is the way we extract the trace anomaly from the functional integral. It the response of the latter under a Weyl (or even a scale) transform δ ω g μ ν = 2 ω g μ ν :
δ ω log Z = d 4 x ω ( x ) g μ ν ( x ) T μ ν ( x )
which is nothing but a particular case of (187). Here g μ ν ( x ) T μ ν ( x ) is the quantum trace of the energy-momentum tensor. Again, the latter can be calculated in various ways with perturbative or non-perturbative methods, that is with the Feynman diagram technique or with the so-called analytic functional method, respectively. The latter include the Schwinger’s proper-time method [21], the heat kernel method [22,23], the Seeley-DeWitt [19,24,25] and the zeta-function regularization [26].
The results concerning the trace anomaly for Dirac (and Majorana) fermions are well established. It is worth recalling that they rely on the possibility to turn to the Euclidean formulation, which is allowed thanks to the above mentioned Feynman prescription to get the Schwinger causal Green’s function or Feynman propagator for the quantum Dirac spinor field. In the absence of such a transition to an Euclidean formulation, the whole construction becomes mathematically meaningless and unreliable.
For Weyl fermions things drastically change. Let us denote a left-handed Weyl fermion by ψ L = P L ψ , where P L = 1 2 ( 1 γ 5 ) . Then, for example, the classical Action on the 4d Minkowski space reads
S L = d 4 x ψ ¯ L D ψ L
The Dirac operator, acting on left-handed spinors maps them to right-handed ones. Hence, the eigenvalue problem itself is not even well posed, so that the Weyl determinant cannot be defined at all. This is reflected in the fact that the inverse of
D L = D P L = P R D
does not exist, since it is the product of an invertible operator times a projector. The full propagator of a Weyl fermion does not exist in this naive form (this problem can be circumvented in a more sophisticated approach, see below). It is incorrect to pretend that the propagator is S L = S P R = P L S . First because such an inverse does not exist, second because, even formally,
D L S L = P R , and S L D L = P L
The inverse of the Weyl kinetic operator is not the inverse of the Dirac operator multiplied by a chiral projector. Therefore the propagator for a Weyl fermion is not the Feynman propagator for a Dirac fermion multiplied by the same projector. The impossibility to define eigenvalues of the kinetic operator or the lack of an inverse for the chiral Weyl kinetic term is a very severe problem. It is clear that the corresponding usual formulas for the chiral Weyl quantum theory do not exist at all.
How can we cope with these problem? In perturbation theory, in order to construct Feynman diagrams, one usually employs the ordinary free Feynman propagator for Dirac fermions. The use of a free Dirac propagator is formally justified, because one can add a free right-handed fermion to allow the inversion of the kinetic operator, as was done above. This is an easy way out, but it immediately raises several questions: in this way one explicitly breaks gauge and Lorentz invariance, and, if a background metric is minimally coupled to the Weyl fermion also diffeomorphisms and local Lorentz invariances are explicitly broken; moreover how can one preserve in this way the notion of chirality, which is a distinctive feature of Weyl spinors? There is no intrinsic way to guarantee that it is preserved all the way through the nontrivial steps of the calculations. In any case one must check a posteriori that diffeomorphism invariance is satisfied8.
The reason why one does such an expedient is because one hopes that the information about chirality is preserved by the fermion-boson-fermion vertex in (190), which contains the P L projector. As it has been shown in the example of Section 4, this method is believed to work correctly as long as it does not come across anomalies. Therefore absence of anomalies is a warranty that we are not introducing unwanted (opposite-handed) degrees of freedom. But this precisely cast doubts on the validity of results when anomalies are present. How can we trust an anomaly calculation when we know that the calculation is reliable only when anomalies are absent? We face a logical ouroboros9. We cannot trust our calculations precisely when we wish to compute anomalies. It sounds very much like the liar’s paradox: we are asking true information to somebody we know may be a liar. This logical short-circuit is not surprising because the procedure we have outlined has a bag at the very beginning: when we add a free right-handed fermion to the left-handed one, the partition function of which we aim to compute, the resulting Dirac operator has well defined eigenvalues, as peculiar as they may be. Hence, the use of the path integral definition as a product of eigenvalues is justified, but here lies the contradiction: this is not what we are looking for, because what we want to compute is the determinant of an operator whose eigenvalues are not even defined.
When we pass to nonperturbative approaches it is not possible to simply replace the missing Weyl propagator with the Dirac one. The reason is that the latter does not contain the correct information about chirality. To follow closely the perturbative approach we should consider the inverse of
i + V P L
The operator (193), as it stands, just treats the two chiralities asymmetrically. We will return to this point later on, but what matters here is to remark that the differential operator (193) is not the inverse of (185) multiplied by P L . If one pretends to replace the full Weyl propagator with the full Dirac propagator, one loses any information concerning the chirality, and trying to recover it by multiplying the Dirac propagator by P L , rather than a method, is a desperate deception.
Even though (193) can be used to compute consistent gauge anomalies yielding satisfactory results (for reasons that will be explained further on), it shares the same drawbacks as the just above outlined perturbative method. One should always and firmly keep in mind that all these approaches are based on the above manipulations, where we assume that the inverse of the relevant kinetic operator exists and if, by any chance, the kinetic operator is not invertible we pretend that it is and go on undeterred with the formalism. This attitude, based on the time-honored faith on practical calculations, may however hide unpleasant implications such as the logical loophole pointed out above. The inevitable conclusion is that all the results (even ours own) obtained with such procedure - that understands the invertibility of the kinetic operator when this is not the case - are intrinsically contradictory and should be considered with a good amount of skepticism.
There is however a way to avoid the logical trap pointed out before. A little thought shows what is critical in the previous procedures: there is a discontinuous jump between say (193) and the operator (191) which we should be computing. It is somewhat like studying a function only at the point where it is singular. There is no continuous parameter that connects one to the other. An elegant way out to bypass the difficulties with Weyl fermions is precisely following this suggestion; it consists in enlarging the space of potentials, by adding a spectator axial vector field, as was first done by W.A.Bardeen [27]. After all the most practical way to overcome the liar’s paradox is to ask many people.
Limiting ourselves for simplicity to the Abelian gauge case the Action integral for a Dirac fermion coupled to a vector V μ and an axial potential A μ is
S [ V , A ] = d 4 x i ψ ¯ i V i A γ 5 ψ
It is invariant under the following vector and axial-vector transformations
δ α V μ = μ α , δ α A μ = 0 δ β V μ = 0 , δ β A μ = μ β
One can now do perturbative quantization of this theory without any difficulty since the fermion is Dirac and its propagator is the standard one. Simultaneously the definition of the functional integral as the product of eigenvalues makes sense because the eigenvalue problem is well defined for it, as complicated as these eigenvalues may be. After carrying out all the necessary calculations one can extract the results relevant for the Weyl fermion case by taking the limit V V / 2 , A V / 2 . For in this limit in (194) one recovers the kinetic operator (193). This yields the bonafide correct results for a theory of Weyl fermions.
This method can be easily extended to non-Abelian potentials. The calculation of the effective Action generated out of (194) can be carried out along the traditional lines. For the perturbative approach everything needed for a correct computation is well defined: the propagator is the one for Dirac fermions, the vertices are of two types: for fermion-fermion-vector potential and for fermion-fermion-axial potential. The approach is the appropriate and traditional one for an operator that has a well defined spectrum of eigenvalues, like the Dirac operator discussed at the beginning of this Section. No poetic license is needed and any logical loophole is avoided. Now it is also evident why the very common (and basically contradictory) approach based on the simple replacement of the (missing) Weyl propagator with the Dirac one, may yield correct results. It is because the calculation steps are the same as in the (194) case as long as one imitates the substitution A μ = 0 for the covariant case and V V / 2 , A V / 2 for the consistent case at the right points: in most cases a risky but often successful approach.
The extension to calculations with a non-trivial background metric requires a further step. This is easily understandable by looking at the kinetic operator (193) which is supposed to be applied to Dirac fermions. In this case we must add a spin connection ω μ , i.e. V is replaced by V + ω , and γ μ = γ a e a μ where e a μ is the inverse tetrad. It is immediately evident that the corresponding theory is not invariant under diffeomorphisms. If, on the other hand, we replace V P L with ω + V P L diffeomorphisms are preserved in the classical Action, but we loose the effect the chirality can have on the metric dependence. The surplace can be overcome by introducing a second spectator metric. The corresponding construction is conceptually simple but requires a laborious extension of the formalism. We give up the idea of presenting it here, and refer to the literature. It was first formulated in [28]. The interested reader can found an extensive review with applications in [29]. One thing is however worth being reminded: in this formalism a central role is played by diffeomorphism invariance. There are unfortunately in the literature examples of calculations where such invariance has simply been explicitly broken, thus compromising the validity of the results.

6.1. Weyl and Massless Majorana Fermions

The second criticality is related to the (wrong) identification between Weyl and Majorana massless spinors. If this identification were possible the problem discussed in the previous Section would not exist, because we could simply replace the Weyl kinetic operator with the massless Majorana one, and the missing Weyl propagator with the existing Majorana propagator (identical to the Dirac one). We have already explained in detail above that Weyl and Majorana spinors belong to different representations of the Lorentz group. A Weyl spinor is a complex spinor with definite chirality (an eigenstate of γ 5 ) which corresponds to its helicity, constituting a minimal representation of the Lorentz group. A Majorana spinor is a self-conjugate bispinor, which can be cast in real form with a suitable choice of the gamma matrices. A Majorana spinor can be represented as a superposition of two Weyl spinors, conjugate to each other and thus with two opposite chirality. Therefore a Majorana spinor has undefined chirality, so that the relation with its helicity is also undefined. Moreover a parity operation maps the Majorana Action into itself, while it maps the Weyl Action (171) into the same Action for the opposite chirality. The same holds for the charge conjugation operator.
There is no possibility of confusing a Weyl spinor with a Majorana massless spinor. Therefore one may wonder how such misunderstanding may have arisen and spread. It certainly comes from a superficial knowledge of the mathematical properties of spinors in 4d, in which a false analogy with 2d, where Weyl and Majorana spinors coexist, may have played a role. But, above all, we suspect, the common properties of the two types of spinors may be at the origin of the misunderstanding. Beside being both massless they have other properties in common. For instance we can decompose any (massless) Dirac fermion into the sum of two Weyl spinors of opposite chirality or into the sum of two independent Majorana spinors. In this sense in model building it is equivalent to use Weyl or Majorana, provided the end result consists of Dirac spinors. Another reason why they are sometimes considered the same object may be due to the fact that we can establish a one-to-one correspondence between the components of a Weyl spinor and those of a Majorana spinor in such a way that the Lagrangian, in two-component notation, looks the same. If, in the chiral representation, we write ψ L as ω 0 , where ω is a two component spinor, then (171) above becomes
i ω σ ¯ μ μ ω
which has the same form (up to an overall factor) as a massless Majorana Action (see (103)). But the Action is not everything in a theory, it must be accompanied by a set of specifications. Even though numerically the Actions coincide, the way they respond to a variation of the Weyl and Majorana fields is different. One leads to the Weyl, the other to the Majorana equation of motion. The delicate issue is precisely this: when we take the variation of an Action with respect to a field in order to extract the equations of motion, we must make sure that the variation respects the symmetries and the properties that are expected in the equations of motion10. If we wish the eom to preserve chirality we must use variations that preserve chirality, i.e. must be eigenfunctions of γ 5 . If instead we wish the eom to transform in the Majorana representation we have to use variations that transform suitably, i.e. must be eigenfunctions of the charge conjugation operation. If we do so we obtain two different results, which are irreducible to each other, no matter what Action we use.
At this point it is perhaps not useless to clarify an issue concerning the already mentioned U(1) continuous symmetry of Weyl fermions. The latter is often confused with an axial R symmetry of Majorana fermions and used to justify the identification of Weyl and massless Majorana fermions. To start with, let us consider a free massless Dirac fermion ψ . Its free Action is clearly invariant under the transformation δ ψ = i ( α + γ 5 β ) ψ , where α and β are real numbers. This symmetry can be gauged by minimally coupling ψ to a vector potential V μ and an axial potential A μ , in the combination V μ + γ 5 A μ , so that α and β become arbitrary real functions. For convenience let us choose the Majorana representation for gamma matrices, so that all of them, including γ 5 , are imaginary. If we now impose ψ to be a Majorana fermion, its four components will be real and only the symmetry parametrized by β makes sense in the Action (let us call it β symmetry). If instead we impose ψ to be Weyl, say ψ = ψ L , then, since γ 5 ψ L = ψ L , the symmetry transformation will be δ ψ L = i ( α β ) ψ L .
This may be the origin of the confusion, because it looks like we can merge the two parameters α and β into a unique one and identify it with the β of the Majorana axial β symmetry. However this is not possible because for a right handed Weyl fermion the symmetry transformation is δ ψ R = i ( α + β ) ψ R . Forgetting β , the Majorana fermion does not transform. Forgetting α , both Weyl and Majorana fermions transform, but the Weyl fermions transform with opposite signs for opposite chiralities. In terms of anomalies, it is well-known that the axial R Majorana symmetry is anomalous: this is the well-known covariant or ABJ anomaly d 4 x β F A F A ( F A is the curvature of A). For a left (right) Weyl fermion the symmetry with parameter α β ( α + β ) is anomalous. This is a consistent anomaly, which, in the Abelian case we are considering, coincides in form with the covariant anomaly (but not in the non-Abelian case), although with a different coefficient and with opposite signs for opposite chiralities. Since a Dirac fermion can be regarded as a sum of two Weyl fermions with opposite chiralities, we see that the anomalies triggered by the α transformation cancel out, while the anomalies triggered by β add up. This is consistent with the well-known fact that for a Dirac fermion the U ( 1 ) α symmetry is not anomalous, while the axial β symmetry is anomalous and corresponds to twice the anomaly of a Majorana fermion. As we see the symmetries and anomalies of Majorana fermions are different from the symmetries and anomalies of Weyl fermions.
No, there is no room for confusing massless Majorana spinors with chiral Weyl spinors. Actually, a Majorana spinor describes neutral spin 1/2 objects - not yet detected in Nature - and consequently there is no phase transformation (U(1) continuous symmetry) involving self-conjugated Majorana spinors, independently of the presence or not of a mass term. Hence, e.g., its particle states do not admit antiparticles of opposite charge, simply because charge does not exist at all for charge self-conjugated spinors (actually, this was the surprising discovery of Ettore Majorana, after the appearance of the Dirac equation and the positron detection). The general solution of the wave field equations for a free Majorana spinor always entails the presence of two polarization states with opposite helicity. On the contrary, it is well known that a chiral Weyl spinor, describing massless quarks and leptons (before the electroweak symmetry breaking11) in the Standard Model, admits only one polarization or helicity state. It always involves antiparticles of opposite helicity and it always carries a conserved internal quantum number such as the lepton number, which is opposite for particles and antiparticles.
Notwithstanding their common properties, the difference between Weyl and massless Majorana spinors becomes crucial in certain instances. One of these is the functional measure in the path integral. The functional measure is well defined for Dirac or Majorana fermions, while it is undefined for Weyl fermions. Here the problem of course arises in nonperturbative calculations. K. Fujikawa invented a method [30] to compute ABJ anomalies (i.e. anomalies that appear in theories with Dirac or Majorana fermions). It is a heuristic method based upon the calculation of the symmetry variation of the functional measure. It works very well for Dirac and Majorana fermions. But the functional measure is computed relying on the decomposition of the fermion fields into a complete set of eigenvectors of the relevant Dirac operator. We know that the eigenvalue problem is not even defined in the case of Weyl fermions, thus it is mandatory to resort to other methods.
There are other more flexible methods, beside Fujikawa’s one, which in addition are mathematically well founded. These are the functional method alluded to above, which we denote with the term functional or `heat kernel’ methods. They translate into a mathematically reliable language the idea of determinant of the kinetic operator. Since a formulation of these methods with a linear Dirac operator is not known (yet), one resorts to a quadratic elliptic operator (for instance the square of the Dirac operator) and bases all calculations on the square root of the resulting path integral. For the method to be valid, however, such a quadratic operator must preserve all the properties of the linear one. Now, it is not difficult to concoct an elliptic operator out of the Weyl-Dirac one, but it is impossible to do it while preserving all the critical information, i.e. while preserving the chirality of the model and its classical symmetries12.
There are in fact insurmountable difficulties. As a matter of fact, first of all, in order to turn the square of the Dirac operator into an elliptic normal differential and matrix-valued operators, the transition to the Euclidean formulation is mandatory. The latter is absolutely legitimate and viable for Dirac fermions, whilst it does not exist at all for Weyl fermions. The very reason is deep and sharp: the Euclidean 4d symmetry group is O ( 4 , R ) which is locally isomorphic to the direct product O ( 3 , R ) × O ( 3 , R ) , in such a manner that any item in this theory must be invariant with respect to the exchange of any spin representation of the two identical and equivalent orthogonal trivial factors, no room being left to the very concept of chirality, which requires two non-equivalent irreducible 2d representations of the Lorentz group.
The solution to these series of problems exists without straining any fundamental principle or lore. But before coming to that we have to clarify another concept, that of Dirac mass and Majorana mass.

6.2. Dirac and Majorana Mass Terms

Another particularly confusing concept seems to be connected with mass terms for fermions. Let us start from a basic element, the mass term ψ ¯ ψ for a Dirac spinor, which is uncontroversial. It can be also rewritten by projecting the latter into its chiral components
ψ ¯ ψ = ψ L ¯ ψ R + ψ R ¯ ψ L
If ψ is a Majorana spinor this can be written
ψ ¯ ψ c = ( ψ c ) L ¯ ψ R + ψ R ¯ ( ψ c ) L ,
which is therefore well defined and Lorentz invariant by construction. Now, using the Lorentz covariant conjugate we can rewrite (198) as
( ψ L ) T C 1 ψ L + ψ L C ( ψ L ) * ,
which is expressed solely in terms of ψ L . The above expression (199) for the mass term may create the illusion that there exists a mass term also for Weyl fermions. But this is not the case. If we add this term to the kinetic term (171) we obtain an Action which gives rise to a field equation that is not Lorentz covariant: the kinetic and mass terms in the field equation belong to two different representations. To be some more explicit, a massive Dirac equation of motion for a Weyl fermion would be
i γ μ μ ψ L m ψ L = 0 ,
but this equation breaks Lorentz covariance because the first piece transforms according to a right-handed representation while the second piece according to a left-handed one, in such a manner that it cannot be derived from a Lorentz invariant Lagrange density.
Instead of the second term in the LHS of (192) one could use m C ψ L ¯ T , which has the right Lorentz properties, but the corresponding Lagrangian term would not be self-adjoint and one would be forced to introduce the adjoint term and end up again with (199). The reason of all this is, of course, that (199) is not expressible in the same canonical form as (171), while it is an allowed mass term for a Majorana fermion. This implies, in particular, that there does not exist such a thing as a “massive Weyl propagator”, that is a massive propagator involving only one chirality. This is something that, in particular, renders the use of the Pauli-Villars regularization problematic. In fact it would be natural to regularize the theory by means of massive Weyl ghost fields. But this is impossible. Sometimes in place of the missing massive Weyl propagator a Dirac or Majorana propagator is used. But this is a leap in the dark: since in these propagators both left and right chiralities appear, one has no warranty that chirality is preserved throughout the calculations and thus that the result one obtains is the correct one. A minimal precaution, in such a case, would be to crosscheck the results obtained with this regularization by comparing them with those obtained with others. In Sectiom 5, for instance, we have shown that the masslessness result of the Weyl kinetic term is preserved at one loop with three different regularization, including a PV one. These results, which are clearly visible in the four component formalism used so far, are much less recognizable in the two-component formalism.
It turns out that, on the one hand, the renowned and famed Goldhaber et al. series of experiments [10] unequivocally proves that neutrinos are always left-handed and of negative helicity, while anti-neutrinos are distinguished and charge conjugated, since they always appear to be right-handed and of positive helicity: it means, beyond any reasonable doubt, that neutrinos and anti-neutrinos are truly Weyl fermions. On the other hand, soon after the discovery of the neutrino flavor oscillations it has been promptly suggested to understand the latter in terms of a Majorana mass term which provides, as we have seen (109), a breakdown of the lepton number conservation and an open door towards flavor changes. Unfortunately, as we have thoroughly discussed and elucidated in the present review, the hypothesis of a Majorana massive neutrino badly and strongly contradicts the Goldhaber et al. long ago well established and solid results: it is impossible to simultaneously keep a full barrel and a drunk wife. The deep reasons to explain and understand neutrino oscillations are elsewhere but in the Majorana mass, at least if we believe in the principles of quantum field theory and the Standard Model construction and if we trust in [10].
Finally, going to nonperturbative approaches, we have already noted that the fermion functional integral measures are different in the two cases. This is a rather heuristic way of formulating the issue. One can express the same idea and calculations in terms of more rigorous functional or heat-kernel like methods. In fact in the literature one meets various different formulations of this problem, which oscillate between rigorous and heuristic approaches. Anyhow, the main central object is the quadratic operator that should represent the square of the kinetic operator, the Dirac operator in our case. Since, as we have already remarked, a version of the heat-kernel involving the linear Dirac operator is not available, one considers a quadratic version, i.e. a functional integral in which the linear Dirac operator is replaced by its square, and then takes the square root of it. This is the logical extension to continuous operators of the determinant of a matrix represented by the product of its eigenvalues. The square of the matrix is represented by the product of the squared eigenvalues.
Once this procedure has been accepted, the most important move is the choice of the quadratic operator. Let us explain this point in some detail. In the case of a Dirac fermion, by applying twice the Dirac operator makes sense because, for instance, for a gauge transformation δ ψ = i α ψ with α = α a T a , where T a are the Hermitean generators of a Lie algebra, we have
δ D ψ = i α D ψ
Thus it makes sense to apply twice D to ψ , because * * D * * 2 ψ under a gauge transformation has the same properties as D ψ . This was the meaning of our statement in the previous Section whereby the quadratic operator must respect the symmetries and features of the original kinetic operator. In many cases one meets in the literature this is not the case.
Let us consider a simple example. The functional integral for a free Dirac fermion (170) D = i + V (where V denotes any potential), is the (suitably regularized) product of its eigenvalues. This is true for a massless Majorana fermion, while for a Weyl fermion it is not so straightforward. Since the Dirac operator anticommutes with γ 5 , it maps a left-handed spinor to a right-handed one. Therefore the eigenvalue problem is not even defined for D L = D P L . We may replace the looked for det D L with det D L D L 1 2 . If D is self-adjoint (after a Wick rotaton) this reduces to det D 2 P L 1 2 . But this is deceivingly simple. It is true that the symmetries are preserved but the most important structure, i.e. the interaction vertex, is lost. Even if the quadratic operator contains a chiral projector the latter is totally ineffective. It is analogous to stating that the inverse of the Dirac-Weyl operator is the inverse of the Dirac operator multiplied by chiral projector. To protect the interaction vertex from disappearing one should use (193) as basic kinetic operator. But, as already noticed, (193) breaks gauge and diffeomorphism covariance. The puzzle can be solved with the already mentioned Bardeen’s approach. The real important point to be retained from this discussion is that the squared Dirac operator one is obliged to use in the functional methods must preserve the same symmetries as the linear Dirac operator which appear in the defining Action.
In most of the nonperturbative functional methods, after turning to a Euclidean setting, one uses the dimensional regularization or a regularization based on analytic continuation such as the zeta function one. At variance with these standard treatments, among the heat-kernel- like approaches there have been also attempts to use masses as regulators, much as in the PV regularization method. We have already remarked that the presence of massive terms introduce in general two chiralities into the game, so that the PV regularization does not seem the most suitable method for calculations in which the chirality must be preserved. In the perturbative case, however, an educated use of PV can lead to the correct result. But, due to its original sin on which it is based, the result must always be confirmed with other regularizations.
The nonperturbative use of masses as regulators consists in choosing instead of a massless quadratic operator, the square of a massive operator, which can only be a Dirac-Majorana massive operator. This method is sometime sold as PV. Now, first of all, as already noticed several times, the quadratic operator inevitably breaks the conservation of chirality, since it involves both chiralities in a balanced form and excludes any possibility to capture, for instance, consistent chiral anomalies. Moreover, the perturbative 1-loop corrections to the classical Weyl Action, as obtained from the original true Pauli-Villars regularization method, by no means can generate any dynamical mass term, neither Dirac nor Majorana, since a PV regularization preserves gauge invariance and U(1) internal chiral phase transformation invariance, so that it protects classical scale invariance and forbids any mass term at the quantum level. Beside other less formal reserves, it is roughly inappropriate to call this type of approach a PV regularization.
To sum up, the use of a Majorana-Dirac operator inevitably drives to a chirally symmetric operator, which cannot produce any left-right asymmetry, i.e. in particular no chiral split anomaly. i.e. with opposite sign for opposite chiralities. Such methods, in particular, are unable to intercept consistent chiral anomalies.

7. Euclidean Fermion Field Theories

As explained in the introduction, the technique of the Wick rotation to regularize the divergent integrals appearing in Feynman diagrams cannot be applied to the original classical Action. Below we explain why, but to do this job we need a short introduction to the Euclidean field theory for fermions, with the warning that the problem of the transition from a Minkowski fermion theory to a Euclidean counterpart [31,32,33,34,35,36] does not admit a unique solution. Below we outline a simple approach.
Let us start from the theory of a massive Dirac fermion and its spinor causal Green’s function, or Schwinger function13 . We can safely perform the natural replacements
i p 0 = p 4 i x 0 = x 4
so as to obtain a positive definite denominator in the Fourier transform, that is
S α β F ( i x 4 , x ) = i ( 2 π ) 4 d p d p 4 exp { i p 4 x 4 + i p j x j } p 4 2 + p 2 + M 2 ( γ 0 p 4 i γ k p k + i M ) α β
Next, we define the Euclidean Dirac Matrices
γ ¯ 4 γ 0 , γ ¯ k i γ k = i γ k ( k = 1 , 2 , 3 )
γ ¯ μ = γ ¯ k , γ ¯ 4 = γ ¯ μ , γ ¯ μ , γ ¯ ν = 2 δ μ ν
so that we can write
S α β F ( i x 4 , x ) = i ( 2 π ) 4 d 4 p E exp { i p E · x E } γ ¯ μ p E μ + i M p E 2 + M 2 α β = i ( 2 π ) 4 d 4 p E exp { i p E · x E } 1 p E i M α β S α β E ( x E )
where we understand
p = ( p 1 , p 2 , p 3 ) = ( p x , p y , p z ) = p E = ( p ¯ 1 , p ¯ 2 , p ¯ 3 )
together with
p 4 γ ¯ 4 + γ ¯ k p ¯ k = γ ¯ μ p E μ p / E / E γ ¯ μ E μ = γ ¯ k k + γ ¯ 4 4
The above move, however, is not cost-free. The Schwinger function (204) conflicts with hermiticity. The way out suggested by Osterwalder and Schrader, which ensures reflection positivity, is to assume that the proper classical variables for a Euclidean formulation of the Dirac spinor field theory, should be two distinct Euclidean bispinors, which we may denote again ψ E and ψ ¯ E , but obey
{ ψ E ( x ) , ψ E ( y ) } = { ψ ¯ E ( x ) , ψ ¯ E ( y ) } = { ψ E ( x ) , ψ ¯ E ( y ) } = 0
for all points x and y of the four-dimensional Euclidean space R 4 , while there is a non-trivial equal time canonical bracket between ψ E and the Hermitean conjugate of ψ ¯ E . The last of the (205) relations is crucial, for it implies that ψ ¯ E does not coincide with the Hermitean conjugate of ψ E times some matrix γ 4 . Thus, if we want to set up a meaningful Euclidean formulation for the Dirac spinor field theory, then we can treat ψ E and ψ ¯ E as totally independent classical Grassmann valued fields. This independence is the main novelty of the Euclidean fermion field theory, the rest of the construction being straightforward.
For instance, we use the definition of the Euclidean matrices γ ¯ μ to derive the O ( 4 ) transformation law for ψ E in the usual way. As a matter of fact we have the six Hermitean generators
Σ μ ν = i 4 γ ¯ μ , γ ¯ ν Σ ı ε ı j k Σ j k ( ı , j , k = 1 , 2 , 3 )
Σ 23 Σ 1 = 1 2 σ 1 0 0 σ 1 Σ 31 Σ 2 = 1 2 σ 2 0 0 σ 2
Σ 12 Σ 3 = 1 2 σ 3 0 0 σ 3 Σ k 4 Σ ¯ k = 1 2 σ k 0 0 σ k
s o ( 4 ) Lie algebra. The six O ( 4 , R ) generators satisfy the following Lie algebra commutation relations, viz.,
[ Σ j , Σ k ] = i ε j k Σ = [ Σ ¯ j , Σ ¯ k ] [ Σ j , Σ ¯ k ] = i ε j k Σ ¯
which can be suitably recast into the four-vector notations to yield
[ Σ μ ν , Σ κ λ ] = i δ μ κ Σ ν λ i δ λ μ Σ ν κ + i δ ν λ Σ μ κ i δ κ ν Σ μ λ
It follows that we have the following orthogonal transformation rules for the two-component spinors and four-component bispinors of the Euclidean formulation,
ψ ( x E ) = exp { ( i / 2 ) σ · ( α + η ) } ψ ( x E )
ψ ( x E ) = exp { ( i / 2 ) σ · ( α η ) } ψ ( x E )
with x E = R x E , R O ( 4 ) 14 while 0 | α | < 2 π and 0 | η | < 2 π , that yields
ψ E ( x E ) = ψ ( x E ) ψ ( x E )
ψ ( x E ) = exp ( i / 2 ) Σ μ ν θ μ ν ψ E ( x E )
θ k 4 = η k θ ı j = ε ı j k α k
It is crucial to gather that the upper and lower spinor components of a bispinor transform according to equivalent fundamental representations D 1 2 of SU(2), at variance with left and right Weyl spinors on the Minkowski space, so that the Euclidean bispinor ψ E ( x E ) belongs to the 4d reducible representation D 1 2 D 1 2 . Thus, the very concept of non-equivalent Weyl spinors of different chirality is lost in this Euclidean formulation of the spinor field theory. Moreover, ψ ¯ E transforms just like the Hermitean conjugated of ψ E , by definition. As a matter of fact, from the unitary property
U ( θ ) = exp { ( i / 2 ) Σ ¯ μ ν θ μ ν } U ( θ ) = U 1 ( θ ) = U ( θ )
we have
ψ ¯ E , ψ E = d 4 x E ψ ¯ E ( x E ) ψ E ( x E ) = d 4 x E ψ ¯ E ( x E ) U ( θ ) U ( θ ) ψ E ( x E ) = ( ψ ¯ E , ψ E )
with x E = R x E , R = R 1 . In this regard, let us focus on the conversion of the Minkowski space chiral matrix
γ 5 = i γ 0 γ 1 γ 2 γ 3 = γ 5 { γ 5 , γ μ } = 0
The standard definition in the Euclidean D dimensional Clifford algebra is
γ ¯ 5 γ ¯ 1 γ ¯ 2 γ ¯ 3 γ ¯ 4 = γ 5 { γ ¯ 5 , γ ¯ μ } = 0
in such a manner that the projector on the upper and lower components of a Euclidean bispinor become
P = 1 2 ( 1 + γ ¯ 5 ) P = 1 2 ( 1 γ ¯ 5 )
the upper and lower spinor components belonging to equivalent fundamental representations of SU(2).
The Euclidean Action integral for the bispinor field is given by
S E [ ψ E , ψ ¯ E ] = d 4 x E ψ ¯ E ( x E ) / E + M ψ E ( x E )
where ψ E is a column bispinor, whereas ψ ¯ E is an independent transposed or row bispinor. Here the overall sign as well as any overall factor are actually irrelevant and arbitrary, so that we could always absorb them into e.g. ψ ¯ E - remember that we are allowed to change ψ ¯ E and not ψ E or viceversa. Conversely, the lack of the factor i in front of the derivative term is not at all conventional: it is there just to ensure that the Euclidean fermion propagator, a.k.a. the 2-point Schwinger’s function, is proportional to ( i p / E M ) / ( p E 2 + M 2 ) ; if it were not for this i , then we would have tachyon poles after transition back to the Minkowski space.
It is worth to notice that the above spinor Euclidean Action integral can be obtained from the corresponding one in the Minkowski space, after the customary standard replacements (201) and (202)
x 4 = i x 0 , γ ¯ k i γ k ( k = 1 , 2 , 3 ) , γ ¯ 4 γ 0
ψ ( x ) ψ E x E ψ ¯ ( x ) e i θ ψ ¯ E ( x E ) ( 0 θ < 2 π )
in such a manner that we can always set
i S [ ψ , ψ ¯ ] S E [ ψ E , ψ ¯ E ] = d 4 x E ψ ¯ E ( x E ) / E + M ψ E ( x E )
Furthermore, the Euclidean Dirac operator / E + M is precisely the one that gives, according to the definition (204), the 2-point Schwinger function inversion formula
/ E + M α β S β η E ( x E ) = δ x E δ α η
Finally we summarize the basic formulas
S E [ ψ ¯ , ψ E ] = d 4 x E ψ ¯ ( x E ) / E + M ψ E ( x E ) S α β F ( i x 4 , x ) S α β E ( x E )
S α β E ( x E ) = d 4 p E ( 2 π ) 4 exp { i p E · x E } i p / E + i M α β / E + M α β S β η E ( x E ) = δ ( x E ) δ α η
It is important to highlight that the Euclidean Action integral for the bispinor field is by no means a real quantity.
What we have just presented is a possible solution for a Euclidean field theory obtained via a Wick rotation from a Minkowski field theory. The conclusion is that, as long as we consider the prescription here introduced, we end up with the doubling of the fermionic fields: the ψ ¯ in the Action has to be replaced with a bispinor independent of ψ , moreover we have to abandon the expectation of a real Action integral. Moreover, in this presentation there is no room for a Euclidean field theory of Weyl fermions. There are however more refined prescriptions, [37,38,39,40,41,42,43], which, for instance, supplement the Wick rotation with a further rotation of the spinor fields. In these new approaches one can impose hermiticity on the Action for a Euclidean Dirac spinor. It is also possible to define a Euclidean Action for Weyl spinors, but it remains impossible to comply with hermiticity: in the latter case the Action is not Hermitean and, as a consequence, the integrand in the path integral is complex.
The gist of this Section can be summarized with the following warning: when a field theory contains Weyl fermions one may be tempted to turn it into a classical Euclidean field theory and quantize this theory, instead of the original Minkowski one, with the aim of obtaining automatically Wick-rotated Feynman diagrams, already prepared for dimensional regularization. Nonetheless, this procedure does not make sense, because we have seen above that the Euclidean field theory one obtains in this way does represent a completely different physics, if any. This fact shares some close resemblance with the problem, illustrated above, of PV regularization: one cannot formulate a local field theory that automatically produces PV regularized Feynman diagrams, for the trivial reason that the inverse of the sum of kinetic operators is not the sum of their inverses.

8. The Family’S Index Theorem and Weyl Fermion Propagators

The existence or non-existence of a fermion propagator in a theory of Weyl fermions is related to the appearance of a special kind of anomalies which we call type O (O stands for obstructive). The latter are also referred to, somewhat neutrally, as consistent chiral anomalies. They are considered as a hazard for a theory, but they are only a symptom of a disease: the true disease is the lack of the Weyl fermion propagator. A fermion primer which does not cover this important correspondence would be largely incomplete; on the other hand an acceptable treatment of this issue would require a long digression in advanced mathematics, which is outside the scope of this paper. Therefore we end it with a summary we hope will be enough to appreciate the importance of the point and the rule that comes out of it.
For simplicity we limit here to the stage of gauge field theories. In a gauge field theory we have to do with spinors that transform according to definite representations of the gauge group G . So, we consider the tensor product of a spinor bundle S C ± with a vector bundle E corresponding to a representation ρ of the structure group G of P ( X , G ) : S C ± E S C ( E ) . The relevant connection will be the spin connection plus a gauge connection V = V μ d x μ valued in the representation ρ of the Lie algebra of G with antihermiten generators. The corresponding Dirac operator
D = D + i V
acts on the space of sections S C ( E ) , i.e. the space of spinor fields, and maps it to itself. The latter splits as
S C ( E ) = S C + ( E ) S C ( E )
where + and denote chiralities, and eq.(220) splits accordingly as
D = 0 D D + 0
The Atiyah-Singer index theorem reads:
I n d ( D + ) d i m ( k e r D + ) d i m ( k e r D ) = X c h ( E ) A ^ ( X )
k e r D ± denotes the space of zero modes of D ± in a specific spacetime manifold X. The symbol c h ( E ) indicates the Chern character of the E bundle, i.e. the rational characteristic class given in terms of the curvature F of V by
c h ( E ) = r + i 2 π tr F + i 2 2 ( 2 π ) 2 tr F 2 + +
where r is the dimension of the representation ρ . The symbol A ^ ( X ) denotes the A ^ -genus, which is the (rational) Pontryagin characteristic class of X. It can be expressed in terms of the Riemann curvature R as follows
A ^ ( X ) = 1 + 1 ( 4 π ) 2 1 12 tr R 2 + 1 ( 4 π ) 4 1 288 tr ( R 2 ) 2 + 1 360 tr R 4 +
This is however the index theorem for a fixed potential V. In gauge field theories we have to do with a family of potentials A , forming orbits of gauge transformations, in which the relevant space is the quotient Q = A G , the moduli space of the theory. In this case the k e r D ± , i.e. the spaces of zero modes of D ± is not fixed, but varies in dimension from point to point of Q , i.e. they are bundles over Q . Therefore the index in this case is the difference of two bundles. Differences of bundles are described by K theory. Using the latter Atiyah and Singer could write down the formula for the index of a family of operators.
In the gauge theory case just outlined in which the Dirac operator acts between two families of sections Γ ( S C + V ) and Γ ( S C V ) the index theorem can be expressed in the form
c h i n d ( D + ) = M c h ( V ) · A ^ ( T Q )
where c h is the Chern character, which for a vector bundle V is defined by
c h ( V ) = r a n k ( V ) + c 1 ( V ) + 1 2 ( c 1 ( V ) 2 2 c 2 ( V ) ) +
In (226), M is a generic spacetime manifold, Q A G is the orbit space of connections, T Q is its tangent space and V is the gauge bundle. The Chern character of the index measures the extent to which k e r D + differs from c o k e r D + = k e r D . It is intuitive that, as long as, i n d ( D + ) differs from 0, the inverse of D + cannot exists
The cohomology classes relevant to a gauge theory in 4d are of order 4 and 6. The order 6 ones, via the transgression formula, give rise to the consistent gauge anomalies. Therefore we see that consistent gauge anomalies are in one to one correspondence with cohomology classes that obstruct the invertibility of the Dirac-Weyl operator and for this reason are called obstructive (type O). But there are also the classes of order 4: they are proportional to tr ( R 2 ) and tr ( F 2 ) , where R is the Riemann curvature and F, as usual, the gauge curvature. They are known as Pontryagin and Chern classes, respectively. In field theory they may appear as odd parity trace anomalies. It is clear that, independently of the field theory interpretation, the corresponding obstructions must be canceled in order to permit the inversion of the Dirac-Weyl operator, i.e. the existence of the corresponding propagator.
Some specifications are in order. When the potential A is a fixed one and M a precise space-time manifold, like in the case of the ordinary index theorem (223), one can, disregarding the zero modes, still define a product of nonzero eigenvalues and interpret it as a partition function. But in the case considered here the zero modes vary from point to point of Q and M represents all possible manifolds, see [3,29], such an operation is generically impossible. The integration over M should not mislead the reader. Contrary to the case of the ordinary index theorem where the topology of spacetime plays a major role, the family’s index theorem is formulated in a universal form, in terms of classifying space and classifying bundles, meaning that it holds independently of the topology of the spacetime. This matches exactly the sense of perturbative approach to field theory where locality requires to make the calculations in a local patch, so that the topology of spacetime is irrelevant.
Another due specification is about the proof of the family’s index theorem: to the best of our knowledge it is available only in the Euclidean version. Following a well established and sensible tradition we considers its statement valid also in Minkowski spacetimes. This does not contradict what has been said in the previous Section: in the present Section we consider the problem of inverting the Dirac-Weyl kinetic operator, in the previous section the object was the construction of a canonical Action involving the same operator as kinetic operator. Both issues are studied in a Euclidean setting, but are different.
Once these clarifications are set forth, the family’s index theorem establishes an uncontroversial rule for perturbative calculations in local field theory. Any method or procedure that concerns the computations of effective Actions, conservation laws and anomalies must be in agreement with this theorem. The family’s index theorem provides the backbone for the calculations of this type.

9. Conclusion

The purpose of this paper was trying to clarify a few controversial or unclear issues concerning fermions, their definition and their use. For this reason we have started in Section 2 and 3 with a review of the their basic properties. In particular, Section 3 is an elementary and detailed construction of Weyl fermion solutions of the 4d Dirac equations, by means of which both Dirac and Majorana massless spinors have been constructed. The corresponding quantum fields have been introduced, their specific properties under discrete symmetry transformations have been established and their causal Green’s functions properly constructed. The general properties of Dirac, Weyl and Majorana fields have been summarized in Section 4. All this preparatory material has been brought to bear in the subsequent Sections. Section 5 is about the perturbative calculation of the effective Action of a Weyl fermion coupled to an Abelian gauge potential. This example exhibits one of the difficult issues considered in this review: the non-existence of the Weyl fermion propagator. This problem is tackled in the most popular way met in the literature: one adds to the original theory a free Weyl fermion with opposite chirality.
The example carried out in Section 5 with three different regularization (dimensional, PV and cutoff) shows that an accurate treatment lead to the same result for all three, as it should, but especially that no mass term can arise from quantization and that the added Weyl fermion of opposite chirality remain free after quantization, so that chirality of the original theory is preserved. The case considered in Section 5 has been then analyzed in general in Section 6. The main result of this analysis has been that the practice of replacing the Weyl propagator, which does not exist, with a Dirac propagator, is however flawed by an underlying logical loophole, so that this trick may accidentally lead correct results but does not have any probatory value. It has been shown that there are ways to circumvent this obstacle by enlarging the space of potentials, such as for instance adding an axial potential to the original vector one.
Section 6 has been devoted to dispelling a not infrequent identification of massless Majorana with Weyl fermions. It is argued that while these two types of fermions have a few characteristics in common, they are crucially different in other fundamental aspects: Weyl fermions are eigenstates of the chirality matrix, while Majorana fermions do not have any definite chirality and no definite helicity. The chirality issue is fundamental in the framework of anomalies: all anomalies with opposite sign for opposite chiralities (the dangerous ones) are absent for Dirac and Majorana fermions. For the very same reason Dirac and Majorana fermions can have a mass term, while for the Weyl ones a mass is forbidden. The difference between the Dirac and Majorana mass terms has been subsequently clarified. The main outcome of our analysis is that, in spite of being build-up by the very same two-component Weyl spinor, the chiral Weyl bispinor and the massless Majorana spinor do exhibit quite different and distinguished mathematical and physical properties and behavior.
There are two cases in which locality of the classical Action is put in jeopardy. One is when the PV regularization is inappropriately transferred from the realm of Feynman diagrams to that of local Actions by confusing the sum of inverses of differential operators (the PV ghost propagators) with the inverse of the sum of the corresponding kinetic differential operators to construct a local Action. This point has been insisted on throughout the paper. There is a similar issue related to Wick rotation. One may be tempted to turn to the Euclidean metric already in the classical Action, but it was shown in Section 7 that the theory one ends up with in this way is a physical system different from the starting one. The Wick rotation is an algorithm concerning the regularization of Feynman diagrams and related amplitudes, which cannot be transferred back mechanically to te original Action integral.
The final Section is devoted to the family’s index theorem. Although anomalies are not a central issue in this review, they nevertheless loom all through it. The anomalies called obstructive (the dangerous ones) pop up in chiral theories as a signal of a disease: the lack of chiral fermion propagators. The connection between anomalies and (non-existence of) chiral propagators is provided by the Atiyah-Singer index theorem. For this reason we have deemed it useful to give a summary review of it, certainly insufficient but hopefully enough to stir the reader’s interest.

References

  1. Bonora, L.; Soldati, R.; Zalel, S. (Eds.) L. Bonora, R. Soldati and S. Zalel, <italic>Dirac, Majorana, Weyl in 4d</italic>, Universe 6 (2020) 8, 111. [arXiv:2006.04546[hep-th]].
  2. Bonora, L.; Soldati, R. L. Bonora and R. Soldati, <italic>On the trace anomaly for Weyl fermions</italic> , arXiv:1909.11991 [hep-th]. [PubMed]
  3. Bonora, L. Fermions and Anomalies in Quantum Field Theories; Springer, 2023. [Google Scholar]
  4. P. Breitenlohner, D. Maison: Dimensional renormalization and the action principle. Comm. Math. Phys. 1977, 52, 11–38. [Google Scholar] [CrossRef]
  5. Breitenlohner, P.; Maison, D. Dimensionally renormalized Green’s functions for theories with massless particles . I Comm. Math. Phys. 1977, 52, 39–54. [Google Scholar] [CrossRef]
  6. Breitenlohner, P.; Maison, D. Dimensionally renormalized Green’s functions for theories with massless particles. II. Comm. Math. Phys. 1977, 52, 55–75. [Google Scholar] [CrossRef]
  7. Pauli, Wolfgang; Villars, Felix. On the Invariant Regularization in Relativistic Quantum Theory. Rev. Mod. Phys. 1949) [ Issue 3 July 1949, 21, 434–444. [Google Scholar] [CrossRef]
  8. Itzykson, Claude; Zuber, JeanBernard. Quantum Field Theory; McGraw-Hill: New York, 1980. [Google Scholar]
  9. Lounesto, P. Clifford Algebras and Spinors; Cambridge University Press, 2001. [Google Scholar]
  10. Goldhaber, Maurice; Grodzins, Lee; Sunyar, Andrew W. Helicity of Neutrinos . Phys. Rev. 1958, vol. 109(n. 3), 1015–1017. [Google Scholar] [CrossRef]
  11. Lawson, H. B., jr.; Michelsohn, M. -L. Spin geometry; Princeton Univ. Press: Princeton, 1989. [Google Scholar]
  12. Vaz, J., jr.; da Rocha, Ronaldo, jr. An introduction to Clifford Algebras and Spinors; Oxford University Press: Oxford, 2016. [Google Scholar]
  13. Merzbacher, E. Quantum Mechanics; John Wiley & Sons, 1961. [Google Scholar]
  14. Peskin, M. E.; Schroeder, D. V. An introduction to Quantum Field Theory; Westview Press, 1995. [Google Scholar]
  15. Bogoliubov, N. N.; Shirkov, D. V. Introduction to the theory of quantized fields; Interscience Publisher: New York, 1959. [Google Scholar]
  16. Ettore Majorana, Teoria simmetrica dell’elettrone e del positrone. In Il Nuovo Cimento; 1937; volume 14, pp. 171–184.
  17. Racah, G. On the symmetry of particles and antiparticles . Nuovo Cim. 1936, 14, 322. [Google Scholar]
  18. Cornwell, J. F. Group theory in physics; Academic Press: London, 1984; vol.II, p. Ch.17. [Google Scholar]
  19. DeWitt, B. S. The dynamical theory of groups and fields; Gordon & Breach, New York, 1965. [Google Scholar]
  20. DeWitt, B. S. Global approach to quantum field theory; Oxford University Press, 2003; vol. I and II. [Google Scholar]
  21. Schwinger, J. On gauge invariance and vacuum polarization . Phys. Rev. 1951, 82, 664–679. [Google Scholar] [CrossRef]
  22. Kirsten, K. Heat kernel asymptotics: more special case calculations . Nucl. Phys. B (Proc. Suppl.) 2002, 104, 119–126. [Google Scholar] [CrossRef]
  23. Vassilevich, D. V. Heat kernel expansion: user’s manual . Phys. Rep. 2003, 388, 279–360. [Google Scholar] [CrossRef]
  24. Seeley, R. T. Complex powers of an elliptic operator . In Proc. Sympos. Pure Math. 10; Amer. Math. Soc., 1967; pp. 288–307. [Google Scholar]
  25. Seeley, R. T. The resolvent of an elliptic boundary value problem . Am. J. Math. 1969, 91, 889–920. [Google Scholar] [CrossRef] [PubMed]
  26. Hawking, S. W. Zeta function regularization of path integrals in curved space-time . Commun. Math. Phys. 1977, 55, 133–162. [Google Scholar] [CrossRef]
  27. Bardeen, W. A. Anomalous Ward Identities in Spinor Field Theories . Phys. Rev. 1969, 184, 1848. [Google Scholar] [CrossRef]
  28. Bonora, L.; Cvitan, M.; Dominis Prester, P.; Giaccari, S.; Paulisic, M.; Stemberga, T. Axial gravity: a non-perturbative approach to split anomalies . Eur. Phys. J. C 2018, arXiv:1807.0124978, 652. [Google Scholar] [CrossRef]
  29. Bonora, L.; Giaccari, S. L. Bonora and S. Giaccari, <italic>Something Anomalies can tell about SM and Gravity</italic>, Symmetry 17 (2025) 273, arXiv:2412.07470 [hep-th]. [PubMed]
  30. Fujikawa, K. On the evaluation of chiral anomaly in gauge theories with γ5 couplings . Phys. Rev. 1984, D29 25. [Google Scholar]
  31. On the Euclidean structure of relativistic field theory . Proc. Nat. Ac. Sc. 1958, 44, 956. [CrossRef] [PubMed]
  32. Schwinger, J. Euclidean Quantum Electrodynamics . Phys. Rev. 1959, 115 117, 721 1407. [Google Scholar]
  33. Osterwalder, K.; Schrader, R. Feynman-kac formula for euclidean fermi and bose fields . Phys. Rev. Lett. 1972, 29, 1423. [Google Scholar] [CrossRef]
  34. Osterwalder, K.; Schrader, R. Euclidean Fermi fields and a Feynman-Kac formula for boson-fermion models . Helv.Phys. Acta 1973, 46, 277. [Google Scholar]
  35. Osterwalder, K.; Schrader, R. Axioms for Euclidean Green’s functions, Comm.Math.Phys. 31 (1973) 83 and. Axioms Euclidean Green’s functions.2, Comm.Math.Phys. 42 (1975) 281.
  36. Fro¨hlich, J.; Osterwalder, K. Is There a Euclidean Field Theory for Fermions . Helv. Phys. Acta 1974, 47, 781. [Google Scholar]
  37. Nicolai, H. A Possible constructive approach to (SUPER ϕ**3) in four-dimensions. 1. Euclidean formulation of the model . Nucl.Phys. B 1978, 140, 294–300. [Google Scholar]
  38. Mehta, M. R. Euclidean continuation of the Dirac fermion . Phys.Rev.Lett. 1990) 1983, 65. [Google Scholar] [CrossRef] [PubMed]
  39. Kupsch, J.; Thacker, W. D. Euclidean Majorana and Weyl Spinors . Fortsch.Phys. 1990, 38 35. [Google Scholar]
  40. Haba, Z.; Kupsch, J. Supersymmetry in euclidean quantum field theory . Fortsch.Phys. 1995, 43, 41–66. [Google Scholar] [CrossRef]
  41. van Nieuwenhuizen, P.; Waldron, A. On Euclidean spinors and Wick rotations . Phys. Lett. B389 1996, 29, e-print: 9608174 [hep-th. [Google Scholar]
  42. Wetterich, C. Spinors in euclidean field theory, complex structures and discrete symmetries, Nucl.Phys. B 852 (2011) 174-234, e-Print:1002.3556 [hep-th].
  43. Streater, R.S.; Wightman, A.S. PCT, Spin, Statistics, and all that; Princeton Univ. Press: Princeton, 2000. [Google Scholar]
1
In this paper we limit ourselves to the traditional classification of fermions. We do not consider the enlargements of the spinor concept due to Lounesto [9].
2
We recall that the neutrino and/or anti-neutrino mass has not yet been directly measured: nonetheless, the present experimental limit is severely stringent < 0.45 eV, i.e. a million times lighter than the lightest known particles.
3
Quite often in the literature it has been suggested that, were the neutrinos described by massive Majorana particles, then the lepton number conservation would be no longer true, so that the neutrino’s flavor oscillations could be accounted for and actually explained. However, it turns out that the renowned Goldhaber et al. experiments [10] clearly proved that neutrinos and antineutrinos are distinguished chiral particles and charge conjugated one to each other, and that their chirality and helicity are strongly correlated, in such a manner that neutrinos are merely left-handed with solely negative helicity, while antineutrinos are strictly right-handed and solely with positive helicity. These experimental facts strongly disagree with the hypothesis that neutrinos and antineutrinos could be described by massive Majorana charge self-conjugated, i.e. neutral, particles with two equivalent helicity states.
4
Here we use the ( L , R ) indexes to understand chiral Weyl bispinors with a nonvanishing (left,right)-handed 2-component Weyl spinors.
5
Customarily, the on-shell 1-particle states of a left Weyl spinor field are a left-handed particle with negative helicity 1 2 and a right-handed antiparticle of positive helicity 1 2 .
6
This means the the minimal choice is S = 2 .
7
For simplicity we understand factors of g , which are necessary, see [19,20], but inessential in this discussion.
8
In 4d we know that there is no room for consistent diffeomorphisms or local Lorentz anomalies: however, spurious (trivial) anomalies of this kind may arise, depending on the regularization procedure. Such spurious anomalies should be eliminated by subtracting counterterms and thus by modifying for consistency the effective Action.
9
A snake that eats its own tail, with the meaning of a self-sustaining argument that produces truth from nothing.
10
For instance, in gravity theories, the metric variation δ g μ ν is generic while not ceasing to be a symmetric tensor.
11
Thanks to chiral anomaly cancelation quarks and leptons pair up, except possibly for neutrinos, to form Dirac spinors which acquire masses after symmetry breaking.
12
Even the symmetry whose anomaly is being computed. The anomaly must result from the clash with the regulator, not because of the explicit breaking of the symmetry at the level of the quadratic operator.
13
The n point Green’s functions in a Euclidean spacetime are usually named n point Schwinger functions.
14
The orthogonal group O ( 4 ) of the rotations in the Euclidean space R 4 is a semi-simple Lie group O ( 4 ) isomorphic to O ( 3 ) × O ( 3 ) .
Table 1. Spinors in even dimension
Table 1. Spinors in even dimension
metric sign. Weyl, C Conjugacy Dirac, C Majorana-Weyl, R Majorana , R
      (2,0)    1L      1R Mutual       2    -            -    2
      (1,1)    1L      1R Self       2    1L         1R    2
      (4,0)    2L      2R Self       4    -            -    -
      (3,1)    2L      2R Mutual       4    -            -    4
      (6,0)    4L      4R Mutual       8    -            -    8
      (5,1)    4L      4R Self       8    -            -    -
      (8,0)    8L      8R Self       16    8L         8R    16
      (7,1)    8L      8R Mutual       16    -            -    -
      (10,0)    16L    16R Mutual       32    -            -    32
      (9,1)    16L    16R Self       32    16L         16R    32
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.