Abstract

We consider extensions of the Dirac equation with mass terms and . The corresponding Hamiltonians are Hermitian and pseudo-Hermitian ( Hermitian), respectively. The fundamental spinor solutions for all generalized Dirac equations are found in the helicity basis and brought into concise analytic form. We postulate that the time-ordered product of field operators should yield the Feynman propagator ( prescription), and we also postulate that the tardyonic as well as tachyonic Dirac equations should have a smooth massless limit. These postulates lead to sum rules that connect the form of the fundamental field anticommutators with the tensor sums of the fundamental plane-wave eigenspinors and the projectors over positive-energy and negative-energy states. In the massless case, the sum rules are fulfilled by two egregiously simple, distinguished functional forms. The first sum rule remains valid in the case of a tardyonic theory and leads to the canonical massive Dirac field. The second sum rule is valid for a tachyonic mass term and leads to a natural suppression of the right-handed helicity states for tachyonic particles and left-handed helicity states for tachyonic spin-1/2 antiparticles. When applied to neutrinos, the theory contains a free tachyonic mass parameter. Tachyons are known to be repulsed by gravity. We discuss a possible role of a tachyonic neutrino as a contribution to the accelerated expansion of the Universe “dark energy.”

1. Introduction

1.1. Generalized Dirac Equations: Mass Terms and Dispersion Relations

Dirac is often quoted saying in some of his talks that the equation that carries his name [1, 2] is “more intelligent than its inventor.” Of course, it needs to be added that it was Dirac himself who found most of the additional insight. Here, we are concerned with extensions of the Dirac equation which contain both tardyonic and tachyonic mass terms. Tardyonic (subluminal) mass terms lead to dispersion relations of the form , whereas tachyonic mass terms lead to superluminal dispersion relations of the form , where is the energy and is the momentum. The generalized, matrix-valued mass term enters the Dirac equation in the form . The are matrices that fulfill the relations , where we choose the space-time metric as . The denotes the partial derivative with respect to the space-time coordinate . It is quite surprising that a systematic presentation of the solutions of the generalized Dirac equations , in the helicity basis [3], has not been recorded in the literature to the best of our knowledge. While the following discussion is somewhat technical, we believe that it will be beneficial to give their explicit form, in order to fix ideas for the following discussion.

We set and use the Dirac matrices in the standard representation and define . For the ordinary Dirac theory, one has (one should say more precisely ) with a real mass , The dispersion relation is . The corresponding Dirac Hamiltonian reads Extensions of the Dirac equation with pseudoscalar mass terms that contain the fifth current have been introduced in the literature. In [4], it is shown that for a mass term of the form , the fermion propagator may obtain nontrivial gradient corrections already at the first order in derivative expansion, for a position-dependent mass. In that case, the fermion self-energy may contribute to a conceivable explanation for -violation during electroweak baryogenesis, as pointed out in [4]. We thus study the following generalized form of the tardyonic (subluminal) Dirac equation: The dispersion relation is . The Hermitian tardyonic Hamiltonian operator reads as We may indicate a further motivation for our study; namely, the unitarity of the matrix implies the existence of useful relations [5] for the even powers obtained upon expanding a one-loop amplitude, formulated with a mass term , in powers of . This implies that a better understanding of the tardyonic equation with two mass terms could be of much more general interest.

It has not escaped our attention that the chiral transformation connects the two Hamiltonians and for and , but it is computationally easier and more instructive to consider the real and imaginary parts of the mass term separately.

Within a systematic approach to generalized Dirac equations with pseudoscalar mass terms, we also consider tachyonic (superluminal) mass terms of the form which induce a superluminal dispersion relation . The corresponding generalized Dirac equation has been named the “tachyonic Dirac equation” and reads as follows [610]: The corresponding Hamiltonian reads The relation has been given in [9, 10]. However, it is much more instructive to observe that is Hermitian; that is, . The concept of Hermiticity is known in lattice theory [11, 12] and is otherwise called pseudo-Hermiticity [1323].

An obvious generalization of the tachyonic case contains an imaginary mass and a mass term, The dispersion relation is . The corresponding Hamiltonian reads and is Hermitian, . For , (9) has been discussed in [24, 25].

It is our goal here to present the fundamental eigenspinors corresponding to the plane-wave solutions of (2), (4), (7), and (9) in a unified and systematic manner. Furthermore, we discuss the second-quantized versions of the fermionic theories described by the generalized Dirac equations. Anticipating part of the results, we may point out that the massless Dirac equation “interpolates” between the tardyonic equations (2) and (4) and the tachyonic equations (7) and (9). For zero mass, helicity and chirality are equal. Helicity and chirality “depart” from each other in very specific directions, when the tardyonic and tachyonic mass terms are “switched on,” as we shall discuss in the following.

1.2. Tachyonic Dirac Equation and Neutrinos: Possible Connections

The tachyonic generalized Dirac equations (7) and (9) describe the motion of superluminal particles, which may either be important for astrophysical studies (neutrinos) or for artificially generated environments such as honeycomb photonic lattices in which pertinent dispersion relations become practically important [26]. The existence of superluminal particles would not falsify Einstein’s theory of special relativity [27], which according to common wisdom is based on the following postulates. (i) The principle of relativity states that the laws of physics are the same for all observers in uniform motion relative to one another. (ii) The speed of light in a vacuum is the same for all observers, regardless of their relative motion or of the motion of the source of the light. Predictions of relativity theory regarding the relativity of simultaneity, time dilation, and length contraction would not change if superluminal particles did exist. Furthermore, as shown by Sudarshan et al. [2831] and Feinberg [32, 33], the existence of tachyons, which are superluminal particles fulfilling a Lorentz-invariant dispersion relation , is fully compatible with special relativity and Lorentz invariance. According to special relativity, it is forbidden to accelerate a particle “through” the light barrier (because for ), but a genuinely superluminal particle remains superluminal upon Lorentz transformation. Significant problems are encountered when one attempts to quantize the tachyonic theories, but again, as shown in [9], these problems may not be as serious as previously thought. In particular, the so-called reinterpretation of solutions propagating into the past according to the Feynman prescription [31] is a cornerstone of modern field theory. Furthermore, it has been shown in [9] that tachyonic particles can be localized and equal-time anticommutators of the spin-1/2 tachyonic field involve an unfiltered Dirac- (see equation (37) of [9]).

Despite these arguments, we can say that, from the point of view of fundamental symmetries, accepting a superluminal neutrino would be equivalent to “an ugly duckling.” Adding to the difficulties, we notice that recent experimental claims regarding the conceivable observation of highly superluminal neutrinos have turned out to be false. One may point out that a relative deviation with at , as claimed by some recent experimental collaborations, would correspond to a negative neutrino mass square in the order of , if one assumes a Lorentz-invariant dispersion relation [34]. Still, there is at present no conclusive answer regarding the conceivable superluminality of at least one neutrino flavor [3539], and it is intriguing that all available direct measurements of the neutrino mass square have resulted in negative expectation values, still compatible with zero within experimental uncertainty, whereas published experimental best estimates for the neutrino speed [4043] have been superluminal, again still compatible with the speed of light within experimental error. The recent ICARUS result [43] is consistent with this trend [4450]; the best estimate for the neutrino velocity is superluminal, but the deviation from is statistically insignificant. The OPERA collaboration [51] has indicated a preliminary, revised result of . Neither subluminal nor superluminal propagation velocities are excluded based on the available experimental data. The “ugly duckling of a superluminal neutrino” is not beautiful; if we are to consider accepting it, then we should be able to hope that the emergence of at least one “added insight” should be the result of this operation.

Before we discuss the possible emergence of these benefits, let us include some historic remarks. According to reliable sources (Professor M. Fink from the University of Austin (Texas) was engaged in discussions with Professor J. A. Wheeler, who found the notion of a nonvanishing neutrino mass so unappealing that he discouraged experimentalists from undertaking any effort to measure the neutrino mass, based on arguments described in Section 1 (M. Fink, private communication, 2012.)), Professor J. A. Wheeler, in his later years at the University of Austin (Texas), used to argue that the neutrino has to be massless, necessarily, and that in his opinion, it could only be a massless Weyl particle with definite helicity (and chirality). Recently presented arguments [52] regarding the possibility of overtaking a subluminal, left-handed neutrino, looking back and seeing a right-handed neutrino, were supposedly already used by Wheeler in order to dispel the conceivable existence of a neutrino mass term. If we assume that right-handed neutrinos are much more massive than their left-handed counterparts, then this leads to a paradox unless one assumes that all light, right-handed neutrinos are sterile. (The problem with a right-handed sterile massive neutrino is that for massive neutrinos, chirality and helicity are different; hence a coupling of the form no longer vanishes for massive Dirac neutrinos if one uses the canonical eigenstates of the massive Dirac equation. One therefore has to invoke additional exotic mechanisms in order to ensure the “sterility” of the right-handed Dirac neutrinos.) Wheeler also disliked (Professor M. Fink from the University of Austin (Texas) was engaged in discussions with Professor J. A. Wheeler, who found the notion of a nonvanishing neutrino mass so unappealing that he discouraged experimentalists from undertaking any effort to measure the neutrino mass, based on arguments described in Section  1 (M. Fink, private communication, 2012.)) the notion of a Majorana neutrino, arguing that the charge conjugation invariance condition imposed on the Majorana particle precludes the existence of plane-wave solutions to the Majorana equation and maximally violates lepton number. Again, these arguments (We are grateful to Professor M. Fink for the clarification and confirmation.) are in full agreement with those recently given in [52].

The original standard model thus called for manifestly massless neutrinos. The commonly accepted observation of neutrino oscillations precludes the possibility that all three generations of neutrino mass eigenstates are massless. Lepton number conservation is based on the global gauge symmetry , applied simultaneously to all lepton fields. A Majorana neutrino would destroy lepton number as a global symmetry but solve the “autobahn paradox,” because a Majorana neutrino would be equal to its own antiparticle and thus, looking back, the right-handed neutrino state would consist of the same particle = antiparticle.

On the other hand, if we assume that the neutrino is described by the tachyonic Dirac equation, then the following statements are valid.(i)Statement 1: we can properly assign lepton number and use plane-wave eigenstates for incoming and outgoing particles, while allowing for nonvanishing mass terms and thus mass square differences among the neutrino mass (not flavor) eigenstates. (ii)Statement 2: there is a natural resolution for the “autobahn paradox” because a left-handed spacelike neutrino always remains spacelike upon Lorentz transformation and cannot be overtaken. (iii)Statement 3: the right-handed particle and left-handed antiparticle states are suppressed due to negative Fock-space norm.(iv)Statement 4: at least qualitatively, tachyonic neutrinos could yield an explanation for a repulsive force on intergalactic distance scales as they are repulsed, like all tachyons, by gravitational interactions (dark energy).

Pauli [53] postulated the existence of neutrinos, on the basis of the conservation of angular momentum and energy, and also introduced pseudo-Hermitian operators [13]. Here, we describe conceivable connections of neutrino physics and pseudo-Hermitian operators. Final clarification can only come from experiment. When in 1956, Reines and Cowan [54] discovered the electron neutrino, two and a half years before Pauli’s death, Pauli replied [55] by telegram “Thanks for message. Everything comes to him who knows how to wait. Pauli.” In defense of the tachyonic hypothesis, we would like to stress that a tachyonic Dirac neutrino would allow us to retain lepton number conservation as a symmetry of nature. We would thus like to write up these thoughts in the current paper, with attention to detail. We should point out that our approach fully conserves Lorentz invariance, in contrast to the extensions of the Standard Model based on Lorentz-violating terms which can otherwise lead to superluminal propagation (see Tables  11 and 13 of [56]).

Units with are used throughout the paper. The organization is as follows: in Section 2, we discuss massless and tardyonic theories. Tachyonic extensions of the Dirac equation are discussed in Section 3. Connections of the fundamental tensor sums over the eigenspinors with the derivation of the time-ordered propagator are analyzed in Section 4. A candidate for dark energy is presented in Section 5. Conclusions are reserved for Section 6.

2. Generalized Dirac Equations: Massless and Tardyonic Theories

2.1. Massless Dirac Theory

The massless Dirac equation and the massless Dirac Hamiltonian read as We note that is both Hermitian as well as Hermitian; that is, . The dispersion relation is . With , we seek positive-energy and negative-energy solutions of the form where denotes a quantum number which is equal to the helicity for positive-energy states and equal to the negative of the helicity for negative-energy states. With , we have . In the massless limit, the solutions to the Dirac equation are given as (see Chapter 2 of [57])The well-known helicity spinors are recalled as These fulfill the fundamental relations , as well as and , where and the sum over is over the values . The sums over the fundamental bispinors and fulfill the following sum rules: as well as Sum rule I can be obtained by a quick explicit calculation, and sum rule II holds because in the massless limit, helicity equals chirality (positive sign for positive energy, negative sign for negative energy). We denote the Dirac adjoint by . One can easily check by an explicit calculation that and . We can thus introduce a factor under the summation over spins in (15) and replace one of the factors by a multiplication of the Dirac adjoint spinor from the right by the fifth current. The Lorentz-invariant normalization of the massless solutions vanishes; that is, .

Eigenstates of the massless Hamiltonian have to be eigenstates of the chirality operator because the chirality commutes with the Hamiltonian, in the sense that . Furthermore, where is the vector of spin matrices, and the helicity operator is identified as . Let be the eigenvalue of chirality and let be the eigenvalue of the helicity operator. Then, the eigenvalue of the Hamiltonian is . Since and , we easily recover the known fact that helicity equals chirality for positive energy, whereas the relation is reversed for negative-energy states (see also Chapter 2.4 of [57]). We are aware of the fact that the considerations reported in the current section partly refer to the literature, but we give them in some detail because they are essential for the following considerations.

2.2. Massive Dirac Theory

We start from the ordinary Dirac equation given in (2), which reads . In the helicity basis, the fundamental spinor solutions read as The algebraic relations that have to be fulfilled by the bispinor amplitudes and read as follows: The dispersion relation is . In the helicity basis, the solutions of (19) with a tardyonic mass term are easily written down, using the identity . With an appropriate normalization factor and after some algebraic simplification, the positive-energy solutions read as follows:The negative-energy eigenstates of the tardyonic equations in the helicity basis are given asThese solutions are consistent with those given in [3] and in Chapter  23 of [58], and the normalizations are () One can change the normalization according to The Lorentz-invariant normalization is equal to one for the fundamental positive-energy bispinors and equal to minus one for the fundamental negative-energy bispinors, A little algebra is sufficient to reproduce the following known sums over bispinors: In accordance with general wisdom of the tardyonic case, these do not involve a helicity-dependent prefactor. The sum rule (25) is of type I (see (15)).

2.3. Two Tardyonic Mass Terms

Inspired by the discussion in Section 1, we consider an equation with two tardyonic mass terms, which has already been indicated in (2) and reads . For the corresponding bispinors in the fundamental plane-wave solutions, this implies thatThe dispersion relation is . The fundamental positive-energy bispinors read as follows:The negative-energy eigenstates of the equation with two tardyonic mass terms are given as for negative helicity (positive chirality in the massless limit) and for positive helicity (negative chirality in the massless limit). In the massless limit (first , then , and then ), we again reproduce the massless solutions , , and . Of course, in the limit , one also has to expand the normalization denominators in powers of . For the solutions (27a), (27b), (28a), and (28b) reduce to the solutions of the ordinary Dirac equation in (20a), (20b), (21a), and (21b) and can be expanded in to yield corrections to the ordinary Dirac equation for small ; that is, . The states are normalized with respect to the condition In the normalization the positive-energy solutions acquire a “positive Lorentz-invariant norm,” whereas the negative-energy solutions have “negative Lorentz-invariant norm,” After some algebra, one can derive the following sums over bispinors: These are easily identified as the positive- and negative-energy projectors. The sum rules do not involve helicity-dependent prefactor and are of type I (see (15)).

The solutions (27a), (27b), (28a), and (28b) approach the massless solutions if one replaces first and then let . They are thus useful for systems where the mass is greater than . For , one would like to calculate solutions that approach the massless case for the sequence , then . These read as follows: In comparison to the solutions (27a), (27b), (28a), and (28b), they acquire a nontrivial phase factor.

3. Generalized Dirac Equations: Tachyonic Mass Terms

3.1. Tachyonic Dirac Equation

The tachyonic Dirac equation is given in (7) and reads . The fundamental bispinors entering the equations fulfill the equations Using and some algebra, the prefactors in the fundamental bispinors (for positive energy) take a very simple form,For negative energy, the solutions read as follows: The normalization condition is () One can change to Lorentz-invariant normalization by a multiplication with , The “calligraphic” spinors fulfill the following helicity-dependent normalizations: where we observe that is a good quantum number because the helicity operator commutes with the Hamiltonian (8). The sum rule fulfilled by the fundamental plane-wave spinors is of type II (see (16)). For the positive-energy spinors, we have where for the negative-energy spinors, the sum rule reads The expressions on the right-hand sides are the positive- and negative-energy projectors.

3.2. Two Tachyonic Mass Terms

We study the equation , as given in (9). The fundamental spinors, which we denote and , fulfill the following equations:

The positive-energy solutions are obtained using the identity . With , they read as follows:

The negative-energy solutions for the tachyonic equation with two mass terms are given asThe normalization condition is . We use a definition of the “calligraphic” spinors analogous to (38), In analogy to (40a) and (40b), a sum rule of type II (see (16)) is fulfilled by the fundamental plane-wave spinors, We thus obtain that the desired projectors onto positive- and negative-energy solutions for the Dirac equation with two tachyonic mass terms (9). The generalized equation is fully compatible with the Clifford-algebra-based approach recently described in [59].

4. Theorems for Generalized Dirac Fields

4.1. Spinor Sums and Time-Ordered Propagator

Our central postulate regarding the quantized fermionic theory is that the time-ordered vacuum expectation value of the field operators should yield the time-ordered (Feynman) propagator, which, in the momentum representation, is equal to the inverse of the Hamiltonian (upon multiplication with ). This postulate implies that under rather general assumptions regarding the mathematical form of the elementary field anticommutators, sum rules have to be fulfilled by the tensor sums over the fundamental spinor solutions. It is perhaps not surprising that these sum rules are precisely of the form investigated in Sections 2 and 3 of this paper.

For definiteness, we consider the solution of the tachyonic Dirac equation (Section 3.1). The generalization to other generalized Dirac equations is straightforward. We start from the field operator [9] where annihilates particles and creates antiparticles. The following anticommutators vanish:We assume the following general form for the nonvanishing anticommutators,with arbitrary and functions of the quantum numbers and . One might argue that since and must be dimensionless, they can depend only on the dimensionless arguments and , but that is a detail of the discussion which we do not pursue any further. Our only assumption concerns the fact that the field anticommutators should be diagonal in the helicity and wave vector quantum numbers, leading to the corresponding Kronecker and Dirac-’s.

We assume that the spin-matrix either constitutes a Lorentz scalar or a pseudoscalar quantity, which is a scalar under the proper orthochronous Lorentz group. The time-ordered product of field operators reads as This equation contains the same coefficient functions and that enter into (48a) and (48b). In order to proceed with the derivation of the propagator, we must postulate that the following sum rules hold: The sum rule (50) is crucial for the further steps in the derivation of the time-ordered propagator. Introducing a suitable complex integral representation for the step function, one obtains from (49), using (50), after a few steps which we do not discuss in further detail (see also equations (3.169) and (3.170) of [57]) The convention is that in any integrals , the component is set equal to in the integrand when it occurs in scalar products of the form , and so forth, but if the integral is over the full , then the integration interval is the full . With the convention which is adopted in many quantum field theoretical textbooks, including [57, 60], we finally obtain the result The tachyonic propagator is identified, under the integral sign, as The sum rules (40a) and (40b) imply that the derivation is valid for the choice in which case the relations given in (50) are fulfilled. Note that this observation does not imply that the choice (54) necessarily is the only one for which we are able to fulfill the postulates given in (50), but it is the only structurally simple choice that we have found.

For the egregiously simple choice (54), let us study the transition to the massless limit (16) in some further detail. Indeed, in the limit , the denominator of the spin sums in (40a) and (40b) vanishes, and a finite limit is obtained after multiplication with ,In order to compare the normalizations of the fundamental spinors in the massless limit, we calculate the following quantities: where the fundamental spinors and of the massless equation have been given in Section 2.1. Observing that in the massless limit, the identifications and , as implied by (56) and (57), show that the identities (55a) and (55b) precisely reduce to the sum rule (16) in the massless limit.

4.2. Generalized Field Anticommutators for Tardyonic and Tachyonic Fields

For definiteness, we have considered the case of the tachyonic Dirac field in the above derivation. The decisive observation is that the tachyonic choice, is consistent with both massive tachyonic fields discussed in Sections 3.1 and 3.2, whereas the tardyonic choice, yields the time-ordered propagator for both massive tardyonic fields discussed in Sections 2.2 and 2.3. The nonvanishing anticommutators for tardyons take the simple form (cf. (48a) and (48b)), Again, compared with (48a) and (48b), the nonvanishing anticommutators for tachyons take the simple form, With these universal choices, the theory of the tardyonic and tachyonic spin-1/2 fields can be unified. The time-ordered propagator is given as We use the sum rules for the tensor sums over fundamental spinors given in (25) (for the tardyonic Dirac field), (32) (for the tardyonic Dirac field with two mass terms), (40a) and (40b) (for the tachyonic Dirac field), and (45) (for the Dirac field of imaginary mass). Going through the exact same derivation as outlined above in between (49) and (53), we obtain the following results for the time-ordered propagators of tardyonic and tachyonic fields. For the tardyonic Dirac field (Section 2.2), one has For the tardyonic field with two mass terms (Section 2.3), the Feynman propagator is easily found as where and are infinitesimal imaginary parts. Both tardyonic mass terms acquire an infinitesimal negative imaginary part, and the prefactor from (46) in the field operator needs to be replaced by , where the tardyonic energy is . For the tachyonic Dirac field (Section  3.1 in [9]), one has the result given in (53). Finally, for the Dirac field with two tachyonic mass terms, we have In the latter case, the prefactor from (46) in the field operator needs to be replaced by , where the tachyonic energy is . For both tachyonic fields discussed here, the mass acquires an infinitesimal positive imaginary part, as manifest in the results given in (53) and (65).

4.3. Tachyonic Gordon Identities

It is useful to illustrate the derivation outlined above by exploring its connection to tachyonic Gordon identities. For definiteness, we again concentrate on the tachyonic Dirac equation discussed in Section 3.1. The matrix element of the vector current finds the following Gordon decomposition for positive-energy spinors: For negative-energy solutions, the identity reads as For , one hasThe matrix element of the axial current reads whereas for negative-energy solutions For , this simplifies toThe results (69a), (69b), (70a), and (70b) for the tachyonic axial vector current have a similar structure as the Gordon decomposition for the tardyonic vector current obtained with the ordinary Dirac equation (see equation of [57]). The role of the Dirac adjoint for the tardyonic case is taken over by the “chiral adjoint” for the tachyonic particle. Here, the designation “chiral adjoint” is inspired by the fact that transforms as a pseudoscalar under Lorentz transformations.

The structure of (68a), (68b), (70a), and (70b) is somewhat peculiar with regard to parity. In (68a), an apparent vector current on the left-hand side appears to transform into an axial current on the right-hand side, whereas in (70a), an apparent axial vector on the left-hand side of the equation becomes what appears to be a vector on the right-hand side. The reason lies in the more complicated behavior of the tachyonic Dirac equation under parity as investigated in [10]. Namely, the tachyonic Dirac equation (7) contains a term which transforms as a scalar under parity, as well as a term which transforms as a pseudoscalar, The mass term in the tachyonic Dirac equation is pseudoscalar and changes sign under parity. Indeed, in [10], the tachyonic Dirac equation has been shown to be separately invariant, and invariant, but not invariant, due to the change in the mass term.

In order to put this observation into perspective, we recall that the entries of the electromagnetic field strength tensor are composed of axial vector components (magnetic field), as well as vector components (electric field). The transformation properties of the electromagnetic field strength tensor under the proper orthochronous Lorentz group are nevertheless well defined.

The transformation (71b) can be interpreted as a transformation under parity. Thus, if we interpret the mass as a pseudoscalar quantity, then the right-hand sides of (68a) and (70a) transform as a vector and an axial vector, respectively. It is the parity noninvariance of the mass term in the tachyonic Dirac equation which leads to the somewhat peculiar structure of (68a) and (70a).

4.4. Helicity-Dependence and Gupta-Bleuler Condition

The anticommutator relations for tardyons given in (60) imply that both left-handed and right-handed helicity states, for both particles and antiparticles, have positive norm. However, the anticommutator relations for tachyons given in (61) imply that right-handed particle as well as left-handed antiparticle states acquire negative norm. This is shown in equations (31) and (32) of [9]. Indeed, for one-particle states , where is the normalization volume in coordinate space. The Fock-space norm is negative for .

For clarification, the corresponding Gupta-Bleuler condition should be indicated explicitly. In full analogy to the instructive discussion of the Gupta-Bleuler mechanism for the photon field, as given in full clarity in Chapter  (9b) of Schweber’s textbook [61], we select the positive- and negative-frequency component of the ()-component of the neutrino field operator,and postulate that it annihilates any physical Fock state of the tachyonic field, These relations automatically imply that the Gupta-Bleuler condition also is realized in terms of the expectation value but the condition (74) is stronger. We recall that the Gupta-Bleuler condition on the photon field reads , where is a Fock state of the photon field. As stressed in Schweber’s book [61] on page 246, the condition is not sufficient for the suppression of the longitudinal and scalar photons, but one must postulate that , where is the positive-frequency component of the photon field operator. Because the tachyonic fermion, unlike the photon, is not equal to its own antiparticle, we need two conditions, given in (74).

A crucial question now concerns the possibility of reversing the helicity-dependence, that is, the question of whether or not a different choice for the helicity-dependent factors in (48a) and (48b) exists that would imply negative norm for left-handed particles and right-handed antiparticles. A related question is whether other tachyonic Dirac Hamiltonians exist for which the choice instead of would fulfill our general postulate, namely, the sum rule (50). For reasons outlined in the following, we can ascertain that this is not the case; tachyonic spin-1/2 particles should always be left handed.

The arguments supporting this conclusion are as follows. First, if we assume that the tachyonic field fulfills a sum rule of type II which for massless fields is given in (16), then it is impossible to replace by in the sum rule because of the necessity to preserve a smooth massless limit. The considerations in the text following (16) imply that if one were to replace by in the massless case, then one would violate the sum rule (16). The second argument is obtained by explicit calculation. We have checked that if one replaces in the tachyonic and imaginary-mass Dirac equations (7) and (9), then the sum rules fulfilled by the corresponding fundamental spinors still contain the characteristic factor . For the imaginary-mass Dirac equation, this result is obtained in [25]. Intuitively, we can understand this result as follows: the mass in the denominators of the right-hand sides of (40a), (40b), and (45) is obtained as the modulus and does not change if we replace in the superluminal Dirac equation. The mass in the numerator of the right-hand sides of (40a), (40b), and (45) changes sign, but this is consistent with the obvious change in the functional form of the positive-energy and negative-energy projectors as we change the sign of the mass term. Again, this consideration supports the conclusion that we cannot invert the helicity-dependence by choosing a different Hamiltonian; the factor persists.

The third argument comes from the tachyonic Gordon identities discussed in Section 4.3. We use the Gordon identity (70a) and (70b) and the normalization (39) to calculate the bispinor trace (with a multiplied from the right) of the left-hand side of (40a), The bispinor trace of the right-hand side of (40a) is which shows the consistency of the bispinor sum (40a) with the Gordon decomposition (70a) and (70b). If we were to replace by , the two sides of the relation (76a) would differ by a minus sign. The bispinor trace of the left-hand side of (40b) is We have used the tachyonic Gordon decomposition for negative-energy states as given in (70a) and (70b), which differs from the positive-energy Gordon decomposition by a minus sign, but an additional minus sign is obtained from the Lorentz-invariant normalization of the negative-energy fundamental bispinors. From the right-hand side of (40b), we have which again is fully consistent, but only because we have a factor in the sum rule (40a) and (40b), which combines with the factor from the Lorentz-invariant normalization of the , to give as a final result in (77a) and (77b).

5. Physical Interpretation: From Tachyonic Neutrinos to Cosmology

5.1. Arguments for and against Tachyonic Neutrinos

In the absence of conclusive experimental evidence, the hypothesis of tachyonic neutrinos has been controversially discussed in the literature. Three main arguments [62] have been brought forward against tachyonic neutrinos. (1) They would require us to give up the notion of a Lorentz-invariant vacuum state, and even the vacuum would become unstable in the presence of tachyonic fields. (2) Given the tachyonic dispersion relation , the role of states with needs to be clarified. (3) The physical (probability!?) interpretation of the conserved Noether current of the free tachyonic Dirac equation has been called into question [62], and it has been argued that no consistent interpretation can be given because certain zero components of the conserved current were conjectured to vanish for all tachyonic momentum eigenstates [62].

A possible answer for question (1) has been proposed in [9]. Summarizing the argument, it has been concluded in [9] that one can solve the problem in two ways. (i) One can Lorentz transform the vacuum state and Lorentz transform all fundamental creation and annihilation operators of the fermion field (some of these will change from annihilators to creators upon transformation, due to the spacelike nature of the tachyons). (ii) One keeps a Lorentz-invariant vacuum state, and only transform the space-time arguments and of the field operators, keeping all creation operators as creators and annihilation operators as annihilators. The amplitudes, cross sections, and so forth, obtained using approach (ii), then depend on scalar product of four vectors which are equal to the result obtained by first calculating the process in the original Lorentz frame and then performing the Lorentz transformation into the moving frame.

The conjecture regarding an expansion of a “false” vacuum in the presence of tachyons can be traced to the fact that most of the tachyonic theories discussed so far in the literature are scalar [2833]. Indeed, scalar tachyons have a problem with instability, because of the structure of the mass term in relation to the field Hamiltonians, which changes sign in a tachyonic theory, suggesting that the field energy can be lowered by creating tachyons. The problem does not occur in tachyonic spin-1/2 theories because a linear, not quadratic, mass term enters the field Lagrangian and Hamiltonian. Provided one reinterprets the spin-1/2 antiparticle solutions in the usual way (negative energy for propagation into the past becomes positive energy for propagation into the future), it then becomes immediately clear that the vacuum energy cannot be lowered upon spin-1/2 tachyon antitachyon pair production.

A solution to problem (2) has also been proposed in [9]. Namely, the energies with (for ) find a natural interpretation in terms of complex resonance and antiresonance energies, which describe unstable states which decay in time. Particle resonances are damped for propagation into the future, and antiparticle antiresonances are damped for propagation into the past, as they should be [9]. The occurrence of momentum eigenstates with real energies, and resonances with complex resonance energies, is a well-known phenomenon all across physics (e.g., in atomic physics, the ground-state energy of the helium atom is strictly real, whereas autoionizing resonances in the three-particle system have a manifestly complex resonance energy).

Another “myth” which should be refuted concerns a conceivable “runaway reaction” where a moving tachyon releases an arbitrarily large amount of energy, as it loses energy and accelerates, given its classical energy-velocity relation . According to this relation, a tachyon indeed accelerates as its energy is lowered and becomes commensurate with the invariant mass square. (At high energy, a tachyon approaches the light cone, though.) An infinitely fast tachyon takes the role of a tardyon at rest [28]; the explicit eigenstates have been indicated in [9]. However, energy conservation holds, and one cannot mix the notion of energy increase by acceleration, which only holds for tardyonic particles, in order to “convert” a tachyon losing energy by acceleration into a tardyon that gains energy in the same process. The “runaway reaction” is impossible; only a finite amount of energy is released as the tachyon accelerates and the energy goes from to zero. Energy conservation holds for tachyons, even if they accelerate, somewhat counterintuitively, when losing energy.

An answer to question (3) has not yet been provided so far in the literature to the best of our knowledge. Here, we aim to provide a possible physical interpretation for the conserved current and scalar product and also show where certain arguments presented originally by Hughes and Stephenson in their research paper [62] entitled “against tachyonic neutrinos” become inconsistent. The Lagrangian density of the tachyonic Dirac particle reads, in first quantization [6], This is equivalent, up to partial integration, to the symmetric form [10] where we suppress the space-time argument . The nonsymmetric form (78) clearly exhibits the presence of the “chiral adjoint” in the Lagrangian. The conserved current is The zero component assumes the following values in plane-wave eigenstates, according to (56) and (57):We conclude that is positive for left-handed particle states () and right-handed for antiparticle states (likewise, ; the helicity is equal to for antiparticles). For , the “axial norm” of the momentum eigenstates given in (82a) and (82b) vanishes, but this is a very special case. Namely, for tachyons, implies for the momentum and thus corresponds to an “infinitely fast” tachyon (), which remains “infinitely fast” upon Lorentz transformation. It would thus be wrong to conclude, as done in [62] in the text following of [62], that all tachyonic states have zero axial norm, because this would correspond to a forbidden generalization of a result which holds asymptotically, for , to all finite values of the tachyonic velocity . Our (82a) and (82b) gives the explicit result for any finite value of the energy .

The spatial integral of the zero component of the conserved current, is precisely equal (up to a sign) to the scalar product which is conserved under the time-evolution by the Hermitian (pseudo-Hermitian) Hamiltonian . This scalar product is not positive definite, as already noticed in the work of Pauli [13], and precisely corresponds to the scalar product introduced by Pauli in equation of [13], where in the notation of [13] we have . Similar observations have been made in equations and of [62]. According to [63, 64], one could otherwise define a so-called operator, which is not equal to the charge conjugation operator and “remedies” the problem of negative norm attained by some states under the norm, leading to a redefined, positive-definite scalar product. However, the negative norm finds a rather natural interpretation, and a redefinition of the scalar product therefore is not required.

For the ordinary Dirac equation, the conserved current is . Its timelike component is , where the subscript reminds us of the electron. The latter can be interpreted as a positive-definite probability density which is conserved under the time evolution generated by the ordinary Hermitian Dirac Hamiltonian .

The tachyonic Dirac current is obtained from the ordinary Dirac by the replacement . We have seen that the scalar product for the tachyonic Dirac Hamiltonian is equal to an integral over the timelike component of the conserved Noether current of the Dirac equation and is not positive definite. In order to put this observation into perspective, it is instructive to recall that for the Klein-Gordon equation , the zero component of the conserved current is not positive definite, either. Therefore, the zero component of the Klein-Gordon current cannot be interpreted as a probability density but must be interpreted as a charge density, which is positive for particles and negative for antiparticles.

This interpretation is not available for the zero component of the Noether current of the tachyonic Dirac equation, because the equation is not charge conjugation invariant and is primarily proposed to describe neutrinos [10]. However, we can come closer to a physical interpretation of the timelike component of the Noether current of the tachyonic equation if we compare the interaction Lagrangian of quantum electrodynamics to the weak interaction of a neutrino and a “heavy photon,” that is, a boson, where is the Weinberg angle. The axial vector part of is where the ordering of the matrices is important. For the interaction with the timelike component of the vector potential, we have the expressions For quantum electrodynamics (QED), we interpret as the probability density, and this suggests an interpretation of the expression as an “axial probability density” or “axial interaction density” of the neutrino field with the timelike component of the boson. According to (82a) and (82b), the “axial interaction density” is positive for the physically allowed states (left-handed particle and right-handed antiparticle states) and negative for the physically forbidden states (right-handed particle and left-handed antiparticle states). This consideration is independent of the suppression mechanism for the states of “wrong” helicity which, in second quantization, proceeds due to negative norm (see the discussion following (72)).

In general, the physical interpretation of tachyonic theories has been discussed by Feinberg [32, 33] and Saudarshan et al. [2831]. The reinterpretation principle is a cornerstone of the theory. The only physically sensible quantities in a quantum theory are transition amplitudes. If the time ordering of superluminal events changes upon Lorentz transformation, then one reinterprets the amplitude as connecting two space-time events whose coordinates are transformed according to the Lorentz transformation, so that the only physically sensible quantity (transition amplitude) simply connects two events, one of which happens before the other [31]. This may seem counterintuitive at first, but it is perhaps a little less counterintuitive if we take into account that the accepted formulation of quantum field theory is based on the Feynman propagator and on the reinterpretation principle for the “advanced” part of the Feynman Green function, which propagates antiparticle solutions into the past. In order to avoid problems with regard to causality, the canonical quantum field theory of subluminal particles has to be supplemented by a reinterpretation principle, just like the tachyonic theory. The physical interpretation of the Noether current is not affected by this consideration.

Finally, let us point out that even without invoking reinterpretation, superluminal propagation can be compatible with causality if we postulate that the tachyonic mass is so small that the superluminality is within the limits set forth by the uncertainty relation. With an energy with and , we have (in units with ). Thus, for an infinitesimal tachyonic mass parameter which does not exceed , causality is preserved within the limits set by the uncertainty principle. Quantum limitations, the role of unstable modes and quantum tunneling in superluminal propagation have been discussed in [6568].

5.2. Gupta-Bleuler Condition and Seesaw Mechanism

The commonly accepted mechanism for the suppression of right-handed neutrino and left-handed antineutrino states is the seesaw mechanism [69]. After integrating the heavy degrees of freedom (the sterile neutrino), it contains a nonrenormalizable dimension five operator and has a hierarchy problem: Namely, the neutrino masses are inversely proportional to the grand unification (GUT) scale and sensitively depend on fine-tuning of . Let us consider equation (18) of [69] and hypothetically consider a formulation that would result if the physically observable neutrino were right handed. Then, we could reformulate equation (18) of [69] as where With the expectation value of the Higgs field, instead of the left-handed Majorana neutrino masses, the right-handed ones would be small, The seesaw mechanism is not unique in suppressing a definite helicity of the neutrino. It is unique provided one formulates it in terms of the left-handed fermion fields, but it could be formulated with inverted helicities if the helicity of the observed neutrinos at low energy were different. In the latter case, the right-handed neutrino mass would be small, instead of the left-handed one. On the other hand, the seesaw mechanism has the distinct advantage that it is not necessary to assume a superluminal character of the neutrino.

The mechanism discussed here in the text following (72) is definite in making a prediction regarding the suppression of right-handed neutrino states, as explained by three independent arguments in Section 4.4. However, we have to assume a superluminal neutrino. While an interacting superluminal field theory is problematic, it is perhaps not as problematic as previously thought (see [70] and Section  4 of [9]).

Final clarification can only come from experiment. The seesaw mechanism is compatible with a Majorana neutrino. The tachyonic Dirac equation implies that the neutrino cannot be equal to its antiparticle, because it does not allow charge conjugation invariant solutions. It is only but not invariant. Experimental evidence for neutrinoless double beta decay is disputed [71], and direct measurements of the neutrino mass square currently exclude neither positive nor negative values [4450]. The generally accepted seesaw mechanism implies that neutrino masses are generated by a nonrenormalizable interaction with a concomitant hierarchy problem, and the mechanism in itself could be reformulated with opposite helicities. It is compatible with a Majorana neutrino. By contrast, a tachyonic neutrino is “automatically” left handed and not equal to its own antiparticle. It is described by a Hermitian Hamiltonian and is plagued with the conceptual difficulties associated with (ever so slightly) superluminal propagation. This means that it is experimentally possible to test the models.

5.3. Tachyonic Neutrinos as a Candidate for Dark Energy

The formulation of a gravitational interaction of a spin-1/2 particle is nontrivial in the quantized formalism. Brill and Wheeler [72] performed the pioneering steps in this direction. In order to formulate the gravitational coupling of a Dirac particle, one has to formulate generalized Dirac matrices , which fulfill anticommutation relations compatible with the local metric of curved space-time. Based on the Christoffel symbols , one formulates the Christoffel affine connection matrices in spinor space and calculates the covariant derivative as follows [7276]: The above formula for is valid in the case of a diagonal metric such as the Schwarzschild metric considered in [74]. For a general space-time metric, and with “West-Coast” conventions for the local vierbein (“East-Coast” conventions were used in [72]), the result reads as where is the spin matrix. The and coefficients transform the Dirac matrices to the local vierbein, With the covariant derivative , the gravitationally coupled Dirac equation reads as [7274] In the case of a tachyonic Dirac particle, it has to be reformulated as follows: The space-time coordinate-dependent matrix can be defined as where is the determinant of the metric and is the local tensor, while is the totally antisymmetric Levi-Cività tensor. The last identity in (96) is valid for a diagonal metric .

Let us first discuss (94) very briefly. A projection onto the upper and lower radial components and in a gravitational field can be found in equations (19) and (20) of [74]. For vanishing electrostatic potential , equation (20) of [74] is invariant under the replacement and . So, if is an eigenvalue of the gravitationally coupled Dirac equation, so is . Invoking reinterpretation and replacing for antiparticles, we find that the spectrum of the gravitationally coupled Dirac Hamiltonian is the same for particles and antiparticles. Therefore, the formalism makes the unique prediction that tardyonic antiparticles, like tardyonic particles, are attracted by a gravitational field. In passing, we note that the often cited motivation for the investigation of trapped antihydrogen and its interaction with the gravitational field therefore is faced with a unique theoretical prediction. Antiparticles are attracted by gravitation as much as particles are.

It has been confirmed within the last two decades [7780] that the Universe expands more rapidly on large distance scales than compatible with the attractive gravitional attraction of the matter density in the Universe. Coupling to a scalar field “quintessence” is usually invoked in order to explain the expansion of the expansion rate of the Universe [81]. As the scalar quintessence field “rolls down its potential,” it accelerates the expansion rate of the Universe. Because of a self-attracting property (the quintessence field energy can be lowered by increasing the local density of the quintessence field), quintessence has positive energy density but negative pressure. This property is necessary in order to constitute a candidate for dark energy [81]. In the standard model of cosmology, dark energy accounts for about 73% of the total mass-energy of the Universe. The quintessence field thereby acts like a time-dependent cosmological constant.

We can thus conclude that most of the energy in the Universe actually is not gravitationally attractive; that is, that gravity can repel. In the following, we shall present qualitative arguments which suggest that tachyonic neutrinos may play a role in the expansion of the Universe and may contribute to “dark energy.” Conceivable connections of tachyonic physics and dark energy have been explored in the literature, employing either scalar fields [82] or Lorentz-violating mechanisms [83]. In order to provide an alternative explanation for dark energy, it is necessary to invoke a mechanism that leads to a repulsive gravitational force on intergalactic distance scales in the Universe. It is not fully surprising that tachyonic neutrinos may provide for such an alternative mechanism. Namely, both on the classical level [84], and on the level of quantum theory (95), tachyonic particles are repulsed by gravitational fields.

On the classical level, this is seen as follows. We start from the familiar equation of motion of a particle with mass in curved space-time (see equation (86’) of [84]), which is a geodesic, With a suitably redefined proper time , the zero geodesic for a tachyon reads as (see equation (86’) of [84]) The force exerted on a tardyon reads, according to equation (79a) of [84], However, there is a sign change for a tachyon (see equation (79b) of [84]), corresponding to the change in the energy-velocity relation for tardyon versus tachyon, which leads to gravitational repulsion for tachyons. Another intuitive way to understand gravitational tachyonic repulsion is to observe that the quantity gets smaller when the position-dependent mass gets smaller, whereas the tachyonic energy gets larger when the position-dependent mass gets smaller.

As neutrinos are ubiquitous within the cosmos, one now needs to evaluate their conceivable contribution to the repulsive force on intergalactic distance scales. A detailed discussion is beyond the scope of the current paper, but we can formulate some initial considerations and order-of-magnitude estimates. One can identify, a priori, two sources of neutrinos which may become important: (i) thermalized neutrinos which decouple in the early Universe [8588] and (ii) nonthermalized, high-energy neutrinos which are continuously generated through various cosmic processes, such as nuclear fusion in stars like the sun, supernovae, and high-energy cosmic background.

It is generally assumed that today’s cosmic background (CMB) radiation is accompanied by a neutrino background, which at some point (low energy) decoupled from the other particles. The average energy of the background neutrinos is generally assumed to be of the same order-of-magnitude as compared to the cosmic microwave background (CMB), and the currently accepted value for the temperature of the neutrino background is times that of the electromagnetic background radiation, that is, about 1.9 K. The Universe, according to this hypothesis, would be filled with a sea of nonrelativistic background neutrinos. This picture changes when we assume that the neutrino is tachyonic. In the earliest stages of the Universe, tachyonic neutrinos are of course highly energetic, with . As they lose energy, they become faster, because (in the classical theory). The energy vanishes as , within quantum theory. For , according to [9], the spectrum of the tachyonic Hamiltonian contains unstable resonances and antiresonances, with a purely imaginary resonance energy. Because the real part of the resonance energy vanishes, the unstable states form a “background” of fluctuating quantum states in the Universe. Particle resonances () and antiparticle antiresonances () are damped for propagation into the future and past, respectively (see [9]). The latter case amounts to propagation into the future if one invokes reinterpretation for the antiparticle solutions [9].

A very course-grained estimate regarding the role of the fluctuating resonance and antiresonance states can be given as follows. We start from the grand canonical ensemble where is the Boltzmann constant, is the normalization volume, is the chemical potential, and is the energy of a tachyon as a function of . The multiplicity factor takes into account the particle and antiparticle resonances and is supplemented here in comparison to the derivation presented in [89]. Furthermore, we have assumed an isotropic energy dependence . The pressure is found as For , the contribution of resonances and antiresonances with to the pressure and energy density is easily evaluated according to equations (27) and (28) of [89]. For the pressure, we have whereas the energy density is The equation of state fulfilled by the zero-temperature limit of the tachyonic Dirac sea is . This is the required equation of state for a “vacuum” energy density that describes a nonvanishing cosmological constant [90]. If , then there is no net energy gain upon pulling on a “piston” which contains the Universe (an illustrative discussion can be found in [91]).

Both results in (103) and (104) are imaginary. This raises the pertinent question of how to incorporate the imaginary pressure and energy density of the neutrino field into the evolution equations of the Universe. In giving rise to unstable resonance states, the tachyonic Dirac field, if it exists, would be different from any other known fundamental quantum fields [57]. An exhaustive mathematical description would require an extension of scattering theory to complex energy and momentum exchanges, applied to low-energy collisions with neutrinos, and describe the continuous decay and repopulation of the unstable tachyonic states in the Universe. A detailed discussion is beyond the scope of the current paper. However, in a first approximation, we can argue that an established technique for the mathematical treatment of complex resonances entails complex scaling (for reviews, see [92, 93]), and we thus explore a complex scaling ansatz to the solution of the evolution equations in the following.

The imaginary results in (103) and (104) describe the quantum fluctuations due to the unstable states. We thus consider equations (8)–(11) of [90], which relate the accelerated expansion rate and the Hubble constant (the dot denotes differentiation with respect to time) to the energy densities of various cosmologically relevant quantities,Here, is the matter density in the Universe and is the energy density associated with the curvature of the Universe (), while is the energy density corresponding to the cosmological constant. The mass density is and . We aim to solve (105a) and (105b) using a complex variable ansatz and . Using (105a) and (105b), with a value  km s−1 Mpc−1 for the Hubble constant, one easily reproduces the accepted value for the critical mass density of  kg cm−3 (in SI units) which corresponds to a critical energy density of  J m−3 (we have set ). In natural units, the result converts to  GeV4.

We courageously set in (105a) and (105b) (flat Universe) and match the modulus of the (complex) accelerated expansion against the currently accepted values of (see [94]) and . According to (105b), the matching is performed by courageously assuming that the bulk of the cosmological constant is given by the complex energy density (104) and equating the complex modulus as follows: The solution is , and so our estimate for the neutrino mass reads as If this treatment is valid, then the mass necessarily has to be the mass of the “heaviest” neutrino mass eigenstate, that is, the one with the largest modulus of the tachyonic mass . Heavier eigenstates would lead to a larger cosmological constant, because is proportional to the fourth power of the neutrino mass. In terms of a conversion to flavor eigenstates, the heaviest mass eigenstates might well be the one closest to the electron neutrino and electron antineutrino lepton flavor eigenstates [95].

We have interpreted the quantities and given in (103) and (104) as contributions to a cosmological constant which describes the evolution of the Universe in complex space and time directions and . We should supplement the solutions for the complex rotation angles for space and time, which read as and , within our complex scaling approach. The evolution of the energy density with the aging of the Universe is given by (see equation of [90]) We therefore obtain a time-independent energy density , which is consistent with the fundamental character of the fluctuating unstable resonances of the neutrino field. We reemphasize that the mass is tachyonic and enters the gravitationally coupled tachyonic Dirac equation (95). The given mass value is not currently excluded by any terrestrial experiment, and it is consistent with the observed time spread of the arrival times of neutrinos from the supernova SN1987A (see [96]). The hypothesis that very light elementary spin-1/2 particles could be tachyonic has been pursued elsewhere [97].

6. Conclusions

In the current investigation, we present the fundamental solutions of generalized Dirac equations in the helicity basis in a systematic and unified manner. Of particular importance are the Dirac equation with two tardyonic mass terms and two tachyonic mass terms . Let us summarize the main results. We have discussed the ordinary (tardyonic) Dirac equation in Section 2.2, a tardyonic Dirac equation (with two mass terms) in Section 2.3, and two tachyonic Dirac equations in Sections 3.1 and 3.2. We give the fundamental eigenspinors that enter the plane-wave solutions of all of these equations (see (20a), (20b), (21a), (21b) (27a), (27b), (28a), (28b), (35a), (35b), (36b), (36a), (42a), (42b), (43a) and (43b)). The eigenspinors are obtained using projector techniques as outlined in Chap. 2 of [57]. For the “normal” Dirac equation (see (20a), (20b), (21a), and (21b)), our results are consistent with [3] and Chapter 23 of [58]. For the generalized Dirac equations, the solutions have not yet appeared in the literature in the compact form given in the current paper, to the best of our knowledge. For the ordinary (tardyonic) Dirac equation, for the tachyonic Dirac equation, and for the imaginary-mass Dirac equation, the prefactors are brought into compact analytic form (see (20a), (20b), (21a), (21b), (35a), (35b), (36a), (36b), (42a), (42b), (43a) and (43b)).

Finally, in Sections 4.1 and 4.2, we find that the tardyonic and tachyonic theories can be unified on the basis of the structurally simple anticommutator relations given in (60) and (61), which are independent of the magnitude of the mass terms. As outlined in Section 4.1, the coefficient functions and in the postulated form of the anticommutators (48a) and (48b) enter the tensor sums over the fundamental spinors in (50). For the egregiously simple choices indicated in (60) and (61), which are consistent with a smooth massless limit, the tensor sums over the fundamental eigenspinors yield the positive-energy and negative-energy projectors, for both tardyonic and tachyonic eigenspinors. Consistency with the massless limit requires the presence of the factor in the tensor sums over the eigenspinors for the tachyonic equations. This fact is verified, on the basis of tachyonic Gordon identities and related considerations, in Sections 4.3 and 4.4. The presence of the factor in the fundamental tachyonic field anticommutators in (61) implies the suppression of right-handed particle and left-handed antiparticle states, due to negative norm, as shown in (72).

Finally, in Section 5, we observe that since tachyons are repelled by gravity, it might be worthwhile to investigate their conceivable role in the mechanism(s) responsible for the accelerated expansion of the Universe “dark energy.” Furthermore, the tachyonic resonances and antiresonance energies might play a role in the sum over states that enter the thermodynamic potentials of a free tachyonic fermionic gas in the low-temperature limit. If we consider the states with imaginary energy to be unstable, fluctuating states, then it is intuitively obvious that they might contribute to a fluctuating energy density and pressure on large distance scales in the Universe. As described in an approach in Section 5.3, an order-of-magnitude calculation based on a complex scaling transformation of the time evolution of the Universe leads to a tachyonic neutrino mass which is not excluded at present by any terrestrial experiment. We believe that it might be worthwhile to explore the physical consequences of the “ugly duckling” (tachyonic neutrino) somewhat further. It allows us to retain, among other things, lepton number conservation as a symmetry of nature.

Acknowledgments

The authors acknowledge helpful conversations with C. Hirata, R. J. Hill, W. Rodejohann, P. J. Mohr, and I. Nándori. This work was supported by the NSF (Grant PHY-1068547) and by the National Institute of Standards and Technology (Precision Measurement Grant).