On the Majorana equation - Relations between its complex two-component and real four-component eigenfunctions

We first derive without recourse to the Dirac equation the two-component Majorana equation with a mass term by a direct linearization of the relativistic dispersion relation of a massive particle. Thereby, we make only use of the complex conjugation operator and the Pauli spin matrices, corresponding to the irreducible representation of the Lorentz group. Then we derive the complex two-component eigenfunctions of the Majorana equation and the related quantum fields in a concise way, by exploiting the so called chirality conjugation operator that involves the spin-flip operator. Subsequently, the four-component spinor solutions of the real Majorana equation are derived, and their intrinsic relations with the spinors of the complex two-component version of the Majorana equation are revealed and discussed extensively.


Introduction
In this paper we first present a new derivation [1] of the two-component Majorana equation including a mass term by a direct linearization of the relativistic dispersion relation of a massive particle, in a way similar to that used originally by Dirac [2]. We subsequently derive from the two-component complex the four-component real version of the Majorana equation [3]. We thus obtain some new results and look at the problem from a fresh theoretical perspective, thus expanding the established theoretical framework as described in modern monographs [4,5] on this subject. One intention of this work is to obtain the Majorana equation "completely on its own rather than as an afterthought when treating the Dirac equation", as it was phrased by Case [6], who first reformulated half a century ago the Majorana theory of the neutrino. Another main objective is to address the key question of whether Majorana particles are their own antiparticles, which they are according to common textbook wisdom. We will show that a real four-component or complex two-component spinor wavefunction does not imply this conclusion necessarily, but following Mannheim [7] we shall consistently define antiparticles as being associated with the negative root of the dispersion relation and particles with the positive one.
The recent paper by Pal [8] titled "Dirac, Majorana, and Weyl fermions" discussed the definitions of and connections between these three fields, and their intrinsic symmetries. We consider the present paper partly as a complement to his largely tutorial work, yet here with emphasis on the two-component Majorana equation for a massive fermion, which is treated without recourse to the Dirac equation [1]. Thereby we only make use of the complex conjugation operator and the Pauli spin matrices, corresponding to the irreducible representation of the Lorentz group. The Lorentz-invariant complex conjugation [8] involves the here introduced spin-flip operator. Its connection to chiral symmetry is also discussed. Finally, we also show that Dirac's equation in the real Majorana representation is a direct consequence of the two-component complex Majorona equation as obtained below.
We derive in much detail the eigenfunctions of the Majorana equation and the related quantum fields, yet other than in the paper by Case [6] in a direct way, whereby we will exploit the spin-flip operator, or as we also may call it perhaps more adequately, the chirality conjugation operator. Pal [8] gave in the introduction to his paper a good motivation, which we fully agree with, for dealing again with this theoretical subject. Our main motivation is to contribute to the ongoing discussion of whether massive neutrinos are Dirac or Majorana fermions, and to better understand the latter theoretically in terms of the twocomponent theory and its connection to the four-component real version. In this respect we also complement the old pedagogical review on Majorana masses by Mannheim [7].

New derivation of the Majorana equation and its symmetries
In this section a new derivation of the two-component Majorana equation is presented, by linearizing the standard dispersion relation of a massive relativistic particle. Following Dirac [2] historically, one can derive his equation in a straightforward way from the relativistic energy and momentum relation, E = m 2 + p 2 , which in fact is constitutive for all field theories of massive particles. Here the four-momentum of the particle is denoted as p µ = (E, p), and the units are such that the speed of light c = 1, and the particle mass is m. The contravariant four-momentum operator in space-time is denoted by P µ = (P 0 , P), and its covariant form is given by the differential operator P µ = (P 0 , −P) = i∂ µ = i(∂/∂t, ∂/∂x) (conventionally we use units of = 1). The energy-momentum relation is usually written in manifestly covariant form as mass-shell condition, p µ p µ = m 2 . Dirac's equation results from a its linearization, which is achieved by introducing four new operators (usually represented by matrices), namely the matrix three-vector α and scalar β, yielding a Hamiltonian in the linear algebraic form As shown in any text book (see, e.g., [13,11]), in terms of matrices one requires at least four dimensions to represent the operators α and β. Furthermore, to satisfy the mass-shell requirement they must algebraically obey the (indicated by the symbol {, }) anti-commutator, and also have the property that β 2 = 1 and (α · p) 2 = p 2 . This is readily ensured, for example in the chiral representation, if we decompose the 4 × 4 matrix as follows, α =ᾱσ, withᾱ 2 = 1. Thus, α can concisely be expressed in terms of the spin operator σ andᾱ. The 2 × 2 Pauli matrices are known to have the required property (σ · p) 2 = p 2 . Similarly, we can write β =βσ 0 . Finally, we quote explicitly the two real matricesᾱ andβ (and a related one calledγ =βᾱ), which are going to be used later and are defined as The three Pauli matrices [9], to which we add here the 2 × 2 unit matrix denoted as σ 0 , have their standard form, but are quoted again because of their importance to what follows. They read and build the three-vector σ = (σ x , σ y , σ z ). We may also introduce the two four-vector forms of the Pauli matrices. They are defined according to Jehle [10] by the four-vectors σ µ ± = (σ 0 , ±σ) and obey σ µ ± = σ ∓µ . Making use of the above matrices, we can obtain Dirac's equation in standard chiral form [1,11] and the classical Hamiltonian representation.
Returning to the basic equation (1), the question may then be asked whether there are other than pure matrix representations of the operators α and β. In a recent paper [1] it was shown that the answer is yes. As this subject is not common knowledge, we shall repeat here some of the key algebra, making use of the operator of complex conjugation C, which transforms any complex number z into its complex conjugate z * . This appears naturally in the symmetry operations of time inversion and charge conjugation of Dirac's equation [11]. Therefore let us define an important operator, which is not of pure algebraic nature but involves this complex-conjugation operator C. We call (in partial nomenclature analogy to β) this operator τ and define it appropriately as τ = σ y C, which differs from the definition in Marsch's paper [1] by an unimportant phase factor i. The importance of τ was first recognized by Case [6]. It is antihermitian conjugate, τ † = −τ , and other than β obeys the relation τ 2 = −1.
Obviously, the operation of τ on the spin vector σ leads to its inversion, i.e., when using the commutation rules of the spin operators, we can show that the operation τ στ −1 = −σ yields a spin flip. Also remember from equation (4) that σ * = σ T , where the superscript T indicates the transposed matrix. As τ flips the spin, we have τ σ µ ± = σ µ ∓ τ , and τ i + iτ = 0, because of the action of C. Therefore, τ also anti-commutes with the momentum four-vector operator P µ , and thus we have τ P 0 + P 0 τ = 0 and τ P + Pτ = 0.
Making use of the above properties of the spin flip operator τ = σ y C, we can now still use a two-component matrix representation given by the Pauli matrices (4). Therefore, by help of τ we can go a decisive mathematical step beyond pure matrix algebra and define the linear energy operator in analogy to the equation (1). The operators involved obey the anti-commutation rule, and also have the property that (σ · p) 2 = p 2 , but τ 2 = −1. We stress again that for the space-time operators {P, τ } = 0 and {P 0 , τ } = 0. When putting P 0 = H in the equation (5), squaring it and multiplying out, the above key features of τ must be exploited. Thus we retain, when inserting the differential operators explicitly, the Klein-Gordon equation. As a consequence of (5), and without recourse to the Dirac equation, we obtain directly a linear wave equation from (5), which is named after Weyl [12] without the mass term and after Majorana with the mass term. It involves only the Pauli matrix operators acting on a two-component spinor φ, but introduces the complications that are caused by the operator τ , and reads as follows: In (7) we have not indicated the 2 × 2-unit matrix explicitly. This equation is nothing but what is called nowadays the two-component Majorana equation [4,5], although it was never written down this way by Majorana himself. As shown by Case [6], this equation can be also derived from Dirac's equation in chiral form, or more generally if one imposes on it the condition of Lorentzcovariant complex conjugation, a procedure which is clearly described in a recent tutorial paper of Pal [8]. However, the view we take here is that equally well we may consider equation (5) as basic, in the spirit of the citation [6] quoted in the introduction.
Obviously, there is a second version of the Majorana equation (7), which is obtained by the operation of τ on it. Namely, when we apply τ from the left side we find that Consequently, if φ solves the equation (7) with the plus sign in front of σ, then χ = τ φ solves it with the minus sign, but with a minus sign also at the mass term. Effectively, this amounts to replacing σ by −σ in (7), also implying τ is exchanged by −τ . Apparently, the two equations (7) and (8) are closely connected twins. It is interesting to note that the two together can provide a special but important four-component Dirac spinor in the form which solves the Dirac equation. Moreover, charge (or better chirality) conjugation [11,13] is, in the chiral representation, given by the operator C =γτ , which when operating on ψ C ± reproduces it. Since C 2 = 1, this operator has only the eigenvalue ±1. Therefore, the spinors ψ C ± constructed from the solution of the Majorana equation are the two eigenfunctions with eigenvalue ±1 of the chirality conjugation operator C. Moreover, the chiral projection operators are given by P ± = 1/2(1±ᾱ), and thus the left-and right-chiral spinor is ψ R,L = P ± ψ. As {γ,ᾱ} = 0, we obtain CP ± = P ∓ C, and consequently, if one applies C separately to the right-or left-chiral fields, one obtains Thus this operator exchanges the chirality of the chiral projections of the complexself-conjugate spinor ψ C , and it transforms left-into right-handed states and vice versa, a virtue which is essentially associated with the spin-flip operator τ as key ingredient of C.
In conclusion, in this special situation the operator C acts as chirality conjugation or reversal. While being a composite of the two related and by themselves irreducible Majorana equations, the reducible chiral Dirac equation has as an intrinsic symmetry this chirality conjugation symmetry.
Consequently, in what follows we can interpret the spinor φ as representing a right-chiral field, and vice versa χ = τ φ as representing a left-chiral field, which is obtained from φ by application of the spin-flip operator. Apparently, it translates from one into the other irreducible representation of the Lorentz group.
Let us consider now the symmetries of the two-component complex Majorana equation, in particular the charge exchange C, parity P, and time reversal T operations. Here we expand the lucid discussion of Case [6]. Generally speaking the Majorana equation is invariant under the symmetry operation O, if the twocomponent spinor also fulfils that equation. When applying the operation O from the left and its inverse O −1 from the right, whereby the unit operator is given by the decomposition OO −1 = 1, we obtain the result We define as usually the time and space coordinate inversion operations T and P on a spinor φ by and also recall the complex conjugation operation C, which gives CiC −1 = −i, and yields Here the asterisk denotes again the complex conjugate number. With these preparations in mind, it is easy to see which operators provide the various symmetry operations, which are composed in the Table 1. To complete the operator algebra, it is important to note that the coordinate reversal operators T and P commute with τ and σ, respectively. Let us first consider in (12) the time reversal, O = T = τ T. Apparently, it does not affect the mass term with τ , and also leaves the first term as well as the kinetic term iσ unchanged. Therefore, φ T = σ y φ * (x, −t) solves the Majorana equation (7) as well. Conversely, the parity operation O = P = iτ P anticommutes with the mass term and inverts the sign of the first term in (12), and it also does not leave the momentum term invariant, since as x changes its sign so does the spatial derivative. Consequently, φ P = iσ y φ * (−x, t) does solve the Majorana equation (7). Finally, we consider the chirality conjugation, O = C = τ , which only changes the sign of the first term. Therefore, φ C = σ y φ * (x, t) does not solve the Majorana equation (7) but solves its conjugate version (8).
In conclusion, of the symmetry operations given in Table 1 time reversal and parity are obeyed, but chirality inversion is maximally broken.
Finally, as the second major topic of our paper we shall reformulate the two-component Majorana equation in its real form by decomposing the spinor φ into its real and imaginary part: φ = φ R + iφ I . Making use of the three real 2 × 2 matrices defined in equation (3) we can write: and similarly we obtain both of which can be combined, by introducing φ ± = φ R ± φ I , in a single equation: Note thatγ 2 = −1,ᾱ 2 = 1, andβ 2 = 1. Furthermore, these matrices all anticommute, and thus by squaring (18) one immediately retains the scalar Klein-Gordon equation reading which in fact was our starting point, when deriving a relativistic wave equation. From the two real two-component Majorana equations (18), when both are combined into a single one, the standard Dirac equation in the real Majorana representation follows immediately. Namely, we may arrange the two spinors φ ± into a single four-component real Dirac spinor and write ψ † = (φ † + , φ † − ). Then the coupled system of (18) transforms into a 4 × 4 real matrix differential equation which reads γ 0 0γ (20) Consequently, we may now introduce the subsequent real 4 × 4 Majorana matrices in their natural (as deduced from the Pauli matrices) representation They mutually anticommute and obey:γ µγν +γ νγµ = −2g µν . Thus the real Dirac equation in Majorana representation reads which can, with the help of the purely imaginary Dirac matrices γ µ = iγ µ , easily be brought into the standard form of the Dirac equation: Therefore, the four-component Dirac equation (22) is, in its real Majorana representation, a direct consequence of the basic two-component complex Majorana equation (7), which was derived here without invoking the Dirac equation in the first place. Subsequently, we shall also derive the four-component real eigenspinors of equation (22).

Eigenfunctions of the two-component Majorana equation with mass term
Returning to the basic equation (7), we here derive its eigenfunctions. In order to solve it, we make the usual plane-wave ansatz, φ( . But note that we require the plane wave and its complex conjugate as well, because of the operator τ in (7). The twocomponent spinors u and v can not be assumed to be complex conjugated, i.e. the wave function cannot be expected to be real. The resulting linked eigenspinor equations are By insertion of the first into the second equation, or vice versa, the relativistic dispersion relation is obtained which yields the two eigenvalues which are obtained from the requirement that for nontrivial solutions of the spinors u and v to exist the determinant associated with the eigenvalue problem (26) must vanish. The negative root in (27) cannot be neglected, since as usually it is related to antiparticles. We can solve (25) for v and insert it back into the above ansatz for φ(x, t) to obtain finally the solutions of (7) and (8) in the concise forms To validate these solutions by direct differentiation, careful attention must be paid to the anticommutation rules between τ and i, respectively σ. We may also solve (25) for u instead, and then it insert it back into the above ansatz for φ(x, t).
We are free to choose for the eigenspinor u the two standard spin up and down eigenfunctions: u † 1 = (1, 0) and u † 2 = (0, 1), and similarly for v, but there is a better and more adequate choice if p is nonzero, as discussed below. Similar solutions like that of equations (28) and (29) are obtained for the negative energy root in (24), yielding the antiparticle wavefunctions. Superposition of all the Fourier modes and their summation over the momentum variable p leads finally to the general Majorana fields, whose quantization then follows from the canonical rules [2][3][4][5] and are given below. Above we obtained the formal solutions of the Majorana equation, which we reiterate here by introducing the Majorana operator M as follows The solution φ obeys Mφ = 0. Note, however, that the four-momentum operator P µ = i∂ µ = i(∂/∂t, ∂/∂x) does not commute with M, and neither does the spin-flip operator τ . Therefore, the φ and χ of the equations (28) to (29) are not eigenfunctions of any of these operators. However, the helicity operator, σ · P, does commute with M. Consequently, to chose the eigenfunctions of the helicity operator in φ for the free functions u or v is most convenient. The eigenvalue equation of the helicity operator in Fourier space reads (σ ·p)u ± (p) = ±u ± (p).
The two eigenvectors depend on the momentum unit vectorp = p/p and are given by in which the half angles of θ and φ appear. The eigenvectors for the samep are orthogonal to each other and normalized to a modulus of unity, and they obey the relation u ± (−p) = u ∓ (p). This is a consequence of the eigenvalue equation (31), which implies that u ± (p) is an eigenvector of the helicity operator, corresponding to a right-handed, respectively left-handed, screw with respect to the momentum direction. According to (32) we have u † ± (p)u ± (p) = 1 and u † ∓ (p)u ± (p) = 0. The dagger denotes as usually the transposed (denoted by the superscript T ) and complex conjugated vector, respectively matrix. The scalar product between two vectors (spinors) v and w is just the sum over the products of their two components, i.e. vw simply stands for v 1 w 1 + v 2 w 2 .
We emphasize that the spin-flip operator τ = σ y C, when operating on the above eigenspinors, leads to τ u ± (p) = ±iu ∓ (p), i.e., it connects the eigenfunctions of opposite helicity, and turns out to be a quite useful relation in the subsequent considerations. By its application, we can write the two possible associated eigenfunctions after (28) as where the advantage of the helicity eigenfunctions becomes obvious. Namely, by use of (31) we obtain the (now normalized to unity) eigenfunction, which is a mixed state involving both helicities in a symmetric fashion. Note that for vanishing mass, m = 0 corresponding to the Weyl equation, only the positive helicity remains, and thus the wavefunction becomes purely right-handed. Its left-handed version is then obtained by applying the spin-flip operator τ on (35) like in equation (29). Operation of M on this φ ± validates that it solves the Majorana equation, i.e. Mφ ± = 0. A simple form of the eigenspinor is obtained for a particle at rest, i.e. p = 0, which yields The helicity eigenvectors in this case can be chosen to be the standard ones obtained for the angles φ = 0 and θ = 0 according to (32).
Concerning the adequate choice of the eigenspinor u 1,2 in (28) for the two possible eigenvalues E 1,2 (p) = ± m 2 + p 2 , we recall that the helicity operator does commute with M, and generally we should consider a mixture or superposition of both helicities. But this is indeed already implied in the solution (35). Therefore, as we can only have two linearly independent eigenvectors to define the basis, we assume that u can be decomposed such that these two eigenvectors are defined by u 1,2 (p) = a 1,2 (p)u ± (p), with some complex amplitude a 1,2 (p) of module unity, |a 1,2 | 2 = 1, to ensure normalization. The natural but arbitrary association with the sign of the energy (27) is to take the positive sign for positive (right-handed) helicity and the negative sign for negative (left-handed) helicity. Thus we may conventionally refer to particles (index 1) and antiparticles (index 2). Upon insertion of this ansatz in (28) we obtain a yet more general expression for the two related eigenfunctions as follows where the obvious momentum arguments in u and a have been suppressed for the sake of lucidity. Apparently, the previous wavefunction (35) is retained for a 1,2 (p) = 1, and E 1,2 = E. It turns out to be convenient to introduce the two real quantities the squares of which add up to unity, ε 2 1 + ε 2 2 = 1. A useful property of the epsilons is the obvious relation (E ± p)ε ∓ = mε ± . From now on E(p) = m 2 + p 2 is used for the positive root obtained from (27), and the argument p in the epsilons will subsequently be omitted for simplicity. Then we can write the two eigenfunctions (with a new name) concisely as (42) To make the antiparticle wavefunction having the same standard plane wave phase like that of the particle, we invert its momentum and note again that the polarization vectors change according to u ± (−p) = u ∓ (p). This inversion is permitted as we will later sum over all momenta, and thus p is just a mute index. Consequently, we take φ P (p) =φ 1 (p) for the particle, but φ A (p) =φ 2 (−p) for the antiparticle, and we redefine the amplitudes as a(p) = a 1 (p) for the particle, but b(p) = ia * 2 (−p), respectively, b * (p) = −ia 2 (−p) for the antiparticle. The new wavefunctions finally read as follows Both eigenfunctions are normalized to unity, yielding the relation φ † P,A φ P,A = 1, which is obtained by using ε 2 + + ε 2 − = 1, exploiting the normalization and orthogonality of the helicity eigenvectors, and noting that a * a = 1 and b * b = 1. They form a complete eigenvector basis of the Majorana operator M and obey Mφ P,A = 0. However, we emphasize again that they are neither eigenfunctions of the four-momentum operator P µ = (P 0 , P), nor of the helicity operator or the spin-reversal operator τ .
By adding up the the two contributions we retain the most general solution of the Majorana equation (7) as the two-component spinor in the form, This superposition of the particle and antiparticle eigenspinors yields the full Majorana wavefunction in the form

The Majorana quantum field
At this point the transition from (46) to a quantum field is quite natural and obvious. We just have to replace the amplitudes in (46) by the canonical anticommuting fermion operators obeying Of course, all possible anticommutators between either a pair of creation, respectively destruction, operators are zero. Mutually, between the as and bs the anticommutators vanish, as they should since these two degrees of freedom are independent. The operator a † (p) creates, and vice versa a(p) annihilates the plane-wave state of a particle of positive helicity, with momentum p and energy E = m 2 + p 2 . The b-operators do the same, yet for the related antiparticle of the opposite negative helicity. We finally get, by keeping explicitly all momentum arguments, the following field operator for any given momentum p in the concise form which assumes the standard form (similar to what is known from the Dirac equation). The quantum field operators corresponding to φ are denoted by a capital Φ. We can now sum up over all momenta to obtain the full Majorana Above we introduced above as abbreviations two operators, which correspond to the creation, respectively annihilation, polarization-vector operators (indicated by an underscore) for a given mixed helicity state as follows: and similarly They obey the standard anticommutation rule, whereby the inner product between the original helicity eigen-vector is used, as well as the relation ε 2 + + ε 2 − = 1. Thus we obtain, by implying the standard scalar product between the complex two-component spinors involved, the result and correspondingly Again all possible mixed anticommutators between the c and d operators vanish. For the two particle number operators we obtain from the above equations (51) and (52) the linear combinations This combined coordinate transformation, involving the original eigenspinors of the helicity operator and their corresponding state operators (classically amplitudes), has a physically important interpretation. As stressed before, the helicity operator does commute with M, but one requires a mixed helicity state to represent its eigenfunctions. This mixing depends on the ratio of the momentum, p, to the total energy, E(p), of a particle. Only if m = 0, we have pure left-or right-handed states, otherwise the massive states are mixed. That MΦ p = 0, follows readily from the notion that where the auxiliary relation (E ± p)ε ∓ = mε ± has been used. We should mention here, that according to its definition the operation of τ on a creation or annihilation operator is defined such that τ a = a † τ , or similarly τ a † = aτ .
As an interlude, we may now consider the unmixed massless case m = 0, which gives us the Weyl field. Then ε + = 1 and ε − = 0. The resulting quantum field operator is given by Apparently, this field involves only the polarization vector for positive helicity. Therefore, in this case the right-chiral Weyl field operator annihilates righthanded particles and creates left-handed antiparticles. Operating with the spin flip operator τ , according to τ aτ −1 = a † and equation (33), on this field yields the left-chiral Weyl field, which reads Consequently, in the case of the left-chiral Weyl field the operator creates righthanded particles and annihilates left-handed antiparticles.
Similarly, we obtain for the massive right-chiral Majorana quantum field and by application of τ on it the left-chiral Majorana quantum field c(p)).
Ultimately, we are interested in the expectation value of a given hermitian operator O for the quantum field Φ, for example the Majorana fields (61) and (62) or the Weyl fields of (59) and (60), which are all not hermitian fields. For the massless Weyl fields, the related plane waves are also eigenfunctions of the four-momentum operator P µ . In contrast, for the massive Majorana field this is not true, and therefore only the expectation value of P µ can be evaluated. The same comment applies to the kinetic helicity operator, which yet does commute with M. This interesting feature was discussed at length and clarified by Mannheim [7], who concluded that for the massive Majorana field theory "second quantization is necessary a priori in order to produce a sensible particle interpretation". Therefore, we can only calculate average expectation values, which are defined by volume integrals of the real part of the binary form This formula shall be now applied to calculate the expectation value of the energy-momentum four-vector, the helicity and number operator of the rightchiral Majorana quantum field Φ M R (x, t) of equation (61). Upon its insertion into (63), and after some lengthy but standard calculations (similar to what is usually done with the Dirac equation, e.g. see the textbooks [11] or [13]) we obtain the concise result where conventionally an infinite constant, − p E(p), has been discarded as irrelevant zero-point energy. Let us consider the helicity density at a given p of the Majorana operator (49). We find and thus by using these relations and by summing up all contributions of equation (65), we obtain the total helicity operator H M R , which is (within an unimportant constant) given by It involves the difference of the number operators of the particles and antiparticles, which therefore are expected to contribute oppositely to the net helicity. Interestingly, we can go through the same procedure with the left-chiral field (62) and end up with exactly the same expectation values as calculated above. One gets a similar result for the Weyl quantum field, by putting either m = 0 in the Majorana field or by using directly equations (59) and (60). In other words it suffices to consider mathematically the original and constitutive equation (5). Chirality conjugation, i.e. the replacement of the spin σ by its negative inverted vector −σ, does not yield new physical information concerning the above expectation values (64) and (66).

Four-component eigenspinors of the Majorana equation
Following the derivation of the real Majorana equation (22), its eigenfunctions can, by construction according to equation (20), easily be derived. Starting point is the complex two-component original wave function (46), in which the wave amplitudes a and b correspond after the previous section to particles, respectively antiparticles, and thus transform into the related annihilation operators for the quantum fields. The four-component eigenfunction is therefore given by the spinor where the real and imaginary parts are taken from (46). This spinor can be written as a sum of the particle and antiparticle components and reads These two contributions can be expressed, in terms of the complex four-component polarization spinors to be defined below, separately as follows: Apparently, ψ P,A = ψ * P,A , and thus the wave functions are real. The polarization spinorsα andβ can be constructed from the two-component eigenfunctions, given in equations (32), of the helicity operator (31). Thus they can be concisely written asα Use has again been made of the symbols ε ± as defined in equation (40). Using their properties and the orthogonality of u ± , one can readily show thatα †α = 1 andβ †β = 1, respectively,α †β = 0 =β †α .
We may now insert the spinors (69) and (70) into the real Majorana equation (22) or complex Dirac equation (23), to validate that they are solutions of those equations. With the covariant four-momentum p µ = (E, −p), one finds that the real and imaginary parts ofα =α R + iα I (respectivelyβ) must obey the coupled equations:γ µ p µαI + mα R = 0, and similarlyγ µ p µαR − mα I = 0, which can be combined to yield the complex solution of the Dirac equation as The same equation must hold true for the polarization spinorβ, and this is indeed the case. The latter equation can explicitly be written in matrix form as Mα = 0, with the non-hermitian matrix which is the 4 × 4-matrix analogue of the Majorana operator M in equation (30). The validation of (75) requires some lengthy algebraic calculations that shall not be done here. Nontrivial solutions ofα to exist requires that the determinant of M vanishes. The corresponding polynomial in the real variable E is of fourth order, yet only yields the two real roots already quoted in equation (27), corresponding to particles and antiparticles. For the sake of completeness we finally give the full four-component polarization spinor α, which in terms of the angles of the momentum unit vector p = (sin θ sin φ, sin θ cos φ, cos θ) and p reads as follows: Similarly, the polarization spinorβ can be derived explicitly from (72). This derivation essentially concludes the section on the solution of the real Majorana equation with mass term. We have presented its eigenspinor for particles ψ P in equation (69) and for antiparticles ψ A in (70). The related quantum fields are readily obtained by interpreting a as annihilation operator, and by replacing a * by a † and interpreting it as the particle creation operator, both together obeying the usual fermion anticommutation rule (47). The resulting field operator Ψ P is not hermitian, though, asα † is not equal toα * . The same comments apply to the antiparticle quantum field described by the creation operator b † and polarization spinorβ * .

Summary and conclusions
The principal goal of this paper was to rederive and discuss the two-component and four-component Majorana equations completely on their own rather than as a spinoff of the chiral Dirac equation. We obtained these equations including a mass term in a new way by a direct linearization of the relativistic dispersion relation of a massive particle. Thereby we only made use of the complex conjugation operator and the Pauli spin matrices, corresponding to the irreducible representation of the Lorentz group. We then calculated the eigenfunctions of the complex two-component Majorana equation and its related quantum field in a transparent way, exploiting the spin-flip or related chirality conjugation operator. Subsequently, we showed the four-component Dirac equation in its real Majorana representation to be the natural outcome of the genuine, twocomponent and complex Majorana equation (7). As the analysis of this version of the Majorana equation clearly reveals, the two associated Majorana particles (with creation operator a † and similarly b † ) represent linearly independent states and are not their own antiparticles, when being defined like in reference [7] in the sense of having positive or negative energies (27). We recall here that the original eigenfunctions (46) and (77) result from a superposition of different states, and thus do neither describe a pure helicity state nor do they have a well defined four-momentum. These distinct features only emerge by evaluating the expectation values (64) and (66) of the quantum fields. They do in fact describe particles and antiparticles with opposite helicities, consistently with the standard perception of neutrinos and antineutrinos.
The symmetries of the two-component Majorana equation were also analysed. It was shown to obey time reversal and parity, yet apparently violates chirality conjugation by construction. As the Majorana quantum field is uncharged, the common term "charge conjugation" does not appear adequate here, but better Lorentz-covariant complex conjugation [8]. Even more appropriately, while stressing clearly the physical meaning, one may speak of "chirality conjugation", which refers to the two conjugate versions of the Majorana equation having opposite chirality. So, the Majorana equation (7), while breaking chiral symmetry by definition, is as basic as Dirac's equation, and when considered as a quantum field can describe massive and uncharged right-handed or lefthanded particles or antiparticles, corresponding to massive neutrinos. In our opinion, the approach starting from scratch with equation (5) gives us new insights into the nature of the Majorana equation, and shows that it deserves to be considered in its own right.