Calculation of the Decay Rate of Tachyonic Neutrinos against Charged-Lepton-Pair and Neutrino-Pair Cerenkov Radiation

We consider in detail the calculation of the decay rate of high-energy superluminal neutrinos against (charged) lepton pair Cerenkov radiation (LPCR), and neutrino pair Cerenkov radiation (NPCR), i.e., against the decay channels nu ->nu e+ e- and nu ->nu nubar nu. Under the hypothesis of a tachyonic nature of neutrinos, these decay channels put constraints on the lifetime of high-energy neutrinos for terrestrial experiments as well as on cosmic scales. For the oncoming neutrino, we use the Lorentz-covariant tachyonic relation E_nu = (p^2 - m_nu^2)^(1/2), where m_nu is the tachyonic mass parameter. We derive both threshold conditions as well as decay and energy loss rates, using the plane-wave fundamental bispinor solutions of the tachyonic Dirac equation. Various intricacies of rest frame versus lab frame calculations are highlighted. The results are compared to the observations of high-energy IceCube neutrinos of cosmological origin.


I. INTRODUCTION
We describe a calculation of the decay rate and energy loss rate of tachyonic (superluminal, "faster-than-light") neutrinos due to (charged) lepton pair Cerenkov radiation (LPCR) and neutrino pair Cerenkov radiation (NPCR). These two decay channels proceed via virtual Z 0 bosons. The processes are kinematically allowed for tachyonic (spacelike) neutrinos, and in the case of LPCR, above a certain energy threshold dependent on the neutrino mass. We base our treatment on a Lorentz-covariant theory of tachyonic (faster-than-light) spin-1/2 particles, i.e., on the tachyonic Dirac (not Majorana) equation [1][2][3][4][5]. Solutions of this equation [6][7][8][9][10] fulfill the Lorentz covariant dispersion relation E = ( k 2 − m 2 ν ) 1/2 , where E is the energy and k is the spatial momentum vector, while m ν is the tachyonic parameter, corresponding to a negative Lorentz-invariant mass square −m 2 ν . The quantity p µ p µ = E 2 − k 2 = −m 2 ν is Lorentz invariant. (Again, we shall assume here that neutrinos are Dirac particles and use the tachyonic Dirac equation [1][2][3][4][5] as a candidate for their physical description.) Tachyonic kinematics are somewhat counter-intuitive. For example, tachyons accelerate as they lose energy. For a subluminal ("tardyonic") particle, one can perform a Lorentz transformation into the rest frame where the spatial momentum k ′ of the particle vanishes. For a tachyonic particle, one can show that, starting from a state with real (as opposed to complex) energy k 2 − m 2 ν , that the Lorentz-transformed momentum always remains greater or equal than than m ν , i.e., k ′ ≥ m ν , and the Lorentz-transformed energy E ′ remains real [10]. One thus cannot possibly enter the rest frame where otherwise we would have k ′ = 0, and the energy would become complex. All that we can do for a tachyon is to transform into a frame where the Lorentz-transformed energy of the neutrino vanishes, i.e., we can enforce E ′ = 0, but not k ′ = 0. The latter frame constitutes a (distant) analogue of the "rest frame" of a tachyonic particle, where according to the classical dispersion relation, the fact that E = m ν / v 2 ν − 1 = 0 otherwise implies an infinite velocity v ν = ∞. All of these intricacies have to be taken into account in the calculation of threshold conditions and decay rates.
Here, we analyze the decay of energetic tachyonic neutrino via LPCR and NPCR. In the calculation of the decay and energy loss rates, we make extensive use of a recently developed formalism which expresses the solutions of the tachyonic Dirac equation in terms of helicity spinors [6][7][8][9][10]. Indeed, helicity remains a good quantum number for tachyonic solutions while the chirality operator does not commute with the tachyonic Dirac Hamiltonian, a fact which, among other things, leads to a natural explanation for the V − A structure of the weak leptonic current [11]. On a different issue, in particle physics, one usually carries out sums over the bispinor solutions using Casimir's trick [12], which is based on sum formulas that allow one to express the sum over the spin orientations of the spin-1/2 in a very concise, analytic form. For the tachyonic Dirac equation, the analogous sum formulas have recently been found [6,9], in the helicity basis.
A further complication arises because the time ordering along a space-like trajectory of a tachyonic neutrino is not unique. For a straight space-like trajectory with velocity v ν > c, it is possible to boost into a system with velocity u = c 2 /v ν , where the tachyonic particle assumes an infinite velocity, according to the velocity addition theorem v ′ = (v ν − u)/(1 − u v ν /c 2 ). Because v ν > c, we still have u = c 2 /v ν < c, which makes the boost permissible. A boost into any frame with velocity u ′ (with u < u ′ < c) will reverse the time ordering along a tachyonic trajectory. The time ordering problem for a tachyonic trajectory is connected with the problem that some fundamental tachyonic particle operators necessarily transform into tachyonic antiparticle operators upon Lorentz transformation [13][14][15][16][17][18][19][20][21]. For the decay of a tachyonic neutrino into an electron-positron pair, this consideration implies that one is at risk of picking up a contribution from neutrino-antineutrino annihilation when considering the decay of an incoming tachyonic neutrino. One can avoid this pitfall by introducing helicity projectors; these eliminate the spurious contribution from the annihilation channel. A clear exposition of the underlying formalism is one of the purposes of the current investigation.
The observation of highly energetic cosmic neutrinos by the IceCube collaboration [22][23][24] puts constraints on the superluminality of neutrinos because they need to "survive" the decay processes ν → ν e + e − and ν → ν ν ν. So, if the decay rate is otherwise sufficiently large in order to account for a substantial energy loss on interstellar time and distance scales, then one may relate the tachyonic threshold to a conceivable high-energy cutoff of the cosmic neutrino spectrum at a threshold energy E = E th ≈ 2 PeV [22][23][24]. Namely, in principle (see Ref. [25]), the tachyonic theory allows us to express the threshold energy as a function of the electron and neutrino masses, E th = f (m e , m ν ); a specific value of the threshold thus implies a definite value of m ν and also determines a numerical value for δ ν = v 2 ν /c 2 − 1, because of the dispersion relation E = m ν /δ 1/2 ν . However, all these conjectures crucially depend on the overall magnitude of the decay and energy loss rates: If these should turn out to be negligible on cosmic distance and time scales, then it will be impossible to relate the tachyonic mass parameter to the cutoff; hence, it is very important to have explicit results for the decay rates at hand.
include Lorentz-violating models [31][32][33][34][35][36][37][38] which lead to superluminality; such models have been applied to the analysis of astrophysical data [39,40]. Energy loss mechanisms due to LPCR have been subjected to alternative statistical analyses [41], and compared to other energy loss mechanisms, e.g., due to neutrino splitting [42]. Neutrino speed modifications have been linked to conceivable (local) variations in fundamental constants [43], and a connection of neutrino speed and neutrino flavor oscillations has been highlighted in Ref. [44]. Gravitational interactions have also been linked to neutrino speed modifications [45,46]. In terms of conceptual questions underlying both spinless as well as spin-1/2 tachyonic theories, including the stability of the vacuum, we refer to the discussion in Refs. [9,47]. A lengthy further discussion on the conceptual issues underlying the tachyonic model would otherwise be beyond the current paper, which already is quite verbose. We organize our investigations as follows. In Sec. II, we derive the energy threshold for LPCR as a function of the tachyonic mass parameter m ν . The derivation of the decay and energy loss rates due to lepton pair Cerenkov radiation is described in Sec. III. For neutrino pair Cerenkov radiation, formulas can be found in Sec. IV. Phenomenological consequences (IceCube data) are discussed in Sec. V. Units with = c = ǫ 0 = 1 are used throughout this paper.

A. Tachyonic Lepton Pair Threshold Based on a Space-Like Dispersion Relation
We consider the process shown in Fig. 1(a), which is lepton-pair Cerenkov radiation (LPCR). The threshold condition reads as where q is the four-momentum of the virtual Z 0 boson, while the incoming and outgoing neutrino momenta are p µ 1 = (E 1 , k 1 ) and p µ 3 = (E 3 , k 3 ). Threshold is reached when, depending on the geometry, the energy transfer from initial to final state is maximum, while the spatial momentum transfer is minimum. This implies that a larger spatial momentum transfer actually is disfavored from a point of view of pair production, because it leads to lesser values of q 2 . In other words, the greater the spatial momentum transfer, the smaller is the four-momentum transfer. Geometrically, we want the outgoing spatial momentum to be as close to the incoming spatial momentum as possible. At threshold, we can thus safely assume that the final neutrino state actually propagates into the same direction as the initial state.
Threshold is reached for a collinear geometry of maximum symmetry. The incoming and outgoing tachyonic particles are on the mass shell, i.e., E 1 = k 2 1 − m 2 ν and E 3 = k 2 3 − m 2 ν . The four-vector notation can thus be reduced to just two components, q = (E 3 , k 3 ) − (E 1 , k 1 ), and the momentum transfer q 2 carried by the Z 0 boson therefore reads as follows, (II.1) Electron-positron pair production threshold is reached at For minimum energy and momentum of the final neutrino state, we have E 3 = 0 and k 3 = m ν . Then, the condition (II.2) transforms into The energy of the tachyonic neutrino at threshold is given as where the latter approximation is valid for δ ≪ 1. One can rewrite this result, based on the tachyonic dispersion relation We note that the tachyonic threshold is a definite function of the mass parameters m ν and m e of the tachyonic neutrino and of the electron (positron), respectively.

B. Tachyonic Neutrino Pair Threshold Based on a Space-Like Dispersion Relation
In the previous section, we found that for electron-positron (charged lepton) pair production, threshold is reached for a collinear geometry. In order to investigate the presence or absence of a threshold for tachyonic pair production (i.e, with two outgoing tachyons), it is instructive to have a look at various geometries. Let us consider a pair of tachyons, both of them on the mass shell, E 2 − k 2 = −m 2 ν . The outgoing particles of the pair are labeled with the indices 2 and 4, as in Fig. 1. If we assume that the tachyons are emitted collinearly and with the same energy, then p µ = (E, k) = p µ 2 = (E 2 , k 2 ) = p µ 4 = (E 4 , k 4 ), and which is negative. For two neutrinos of different energy, emitted collinearly. i.e., with k 2 = k 2êz and k 4 = k 4êz ), one has For small tachyonic mass parameter m ν , a Taylor expansion yields In the limits k 2 → 0, k 4 → ∞ (or vice versa), q 2 assumes very large negative numerical values, demonstrating the absence of a lower threshold.
One might ask, however, if there is perhaps a higher cutoff for the allowed q 2 in relativistic tachyonic pair production kinematics. For the production of an anti-collinear pair, one has For large k 1 and k 2 , this expression assumes arbitrarily large positive numerical values. The only condition relevant to the allowed range of q 2 for tachyonic pair production thus is This result has important consequences for the calculation of neutrino-pair Cerenkov radiation (see Fig. 1(b)).

C. Tachyonic Maximum Momentum Transfer and Fermi Theory
One crucial question one might ask concerns the applicability of Fermi theory for the decay processes shown in Figs. 1(a) and (b), in the high-energy regime. The question is whether the condition q 2 ≪ M 2 Z , which ensures the applicability of Fermi theory, remains valid for a highly energetic, oncoming neutrino. Concerning this question, we first recall that, as already shown, threshold for pair production is reached for collinear geometry, i.e., when the final neutrino momentum k 3 is along the same direction as the initial state momentum k 1 . This implies that the maximum momentum transfer, for given energy E 1 of the incoming particle, also is reached for collinear geometry. Reducing space to one dimension, we find for the square q 2 of the momentum transfer, (II.11) For given k 1 , maximum four-momentum transfer is reached when the momentum of the outgoing particle is equal to k 3 = m ν , and thus E 3 = 0. This implies that The condition for using the effective Fermi theory for the virtual Z 0 boson exchange in Figs. 1(a) and (b) is q 2 ≪ M 2 Z , which in the high-energy limit can be reformulated as where we have conservatively estimated the neutrino mass parameter to be on the order of 1 eV. (In general, one estimates the neutrino masses to be of order (0.01 ÷ 0.05) eV, see Sect. 1 of Ref. [48] and the discussion around Eq. (14.21) of Ref. [49].) The condition E 1 ≪ 10 22 eV = 10 7 PeV is easily fulfilled by the most energetic neutrinos seen by the IceCube collaboration [22,23], which do not exceed ∼ 2 PeV in energy. Hence, we can safely assume Fermi theory to be valid in the entire range of incoming neutrino energies relevant for the current investigation.

D. "Rest" Frame of the Tachyon
Let us briefly analyze the role of the "rest frame" of the faster-than-light, incoming neutrino in the context of the tachyonic dispersion relation E = k 2 − m 2 ν . As is evident from a Minkowski diagram (see Fig. 2), the rest frame of the tachyonic "space-like" neutrino cannot be reached via a Lorentz transformation. By contrast, for a tachyon, it is possible to transform into a frame where the time interval (not the space interval!) swept on the tachyonic trajectory is zero, i.e., the tachyonic particle assumes an infinite velocity. This frame of infinite velocity, in some sense, constitutes the equivalent of the rest frame; namely, the incoming particle has zero time evolution (as opposed to zero space evolution), and thus infinite velocity. According to tachyonic kinematics, it then has zero energy. For illustration, we consider the boost into a frame with energy 0 < u < 1 (see Fig. 3), For a boost velocity u = E/k = k 2 − m 2 ν /k < 1, we have E ′ = 0, k ′ = m. One might be tempted to suggest that the decay rate calculation could be simplified in the tachyonic "rest" frame (with respect to the time, not space, i.e., E ′ = 0). However, in this frame, one cannot calculate the decay rate. This is most easily seen from an energy conservation condition. The oncoming neutrino energy vanishes for infinite velocity (E ′ = 0). Hence, the oncoming particle cannot provide the energy necessary to produce an electron-positron pair.
By contrast and for comparison, for a tardyonic (subluminal) particle, the Dirac "gap" between positive-energy and negative-energy states ensures that the energy of an oncoming, say, muon, is always bound by its rest mass m µ from Under a Lorentz transformation, the modulus of the momentum vector k ′ indeed never gets smaller than k ′ > mν (see the solid curve). The energy E ′ (indicated by long dashes), however, can go to zero, and in fact changes sign at the point where the modulus of the momentum vector just becomes equal to k ′ = mν which is the point of infinite velocity. (The constant curve k ′ = mν is indicated via short dashes.) When the energy E ′ changes sign, the propagation direction of the neutrino changes sign, and it moves in the negative x direction as opposed to the positive x direction. From the plot (solid curve), one might think that the momentum component along the boost does not change sign, but this is not physical. From the Minkowski diagram (see Fig. 2), one sees that the neutrino is still moving along the positive x axis, but with the time ordering of the start and end point interchanged.
below, even under a Lorentz transformation. Hence, the muon decay from rest, with E = m µ ≫ 2m e , is kinematically possible. Because the oncoming muon is timelike, the emitted virtual W boson can still carry enough momentum transfer q 2 > 0 in order to produce an electron, and an electron antineutrino. This is not the case for an oncoming tachyonic neutrino, whose energy is not bound from below, and can in fact vanish. When the oncoming neutrino energy vanishes, so does the decay rate.
Alternatively, one can observe that in its own rest frame (the "real rest frame" where the tachyon has a vanishing spatial momentum), the neutrino has the following properties, The energy becomes imaginary in its own rest frame. According to Fig. 3, the rest frame of a space-like, tachyonic neutrino cannot be reached via a Lorentz transformation, consistent with the purely real (rather than complex) quantities which enter Eq. (II.14). A further kinematic consideration is illustrative. Namely, according to Fig. 3, the energy of the tachyonic particle decreases as one "chases" it, then approaches zero and eventually flips sign at boost velocity u. For boost velocities beyond this point, the energy becomes negative, or alternatively, the time ordering of the start and end point of the trajectory of the tachyon reverses. A left-handed tachyonic neutrino, for boost velocities beyond u, would be seen as a right-handed antitachyon moving in the opposite spatial direction, for the moving observer. The spatial momentum k ′ , as seen from Fig. 3, always remains in the region k ′ ≥ m ν . The region with imaginary energies E ν = ±i m 2 ν − k 2 ν with k ν < m ν , never can be reached for an initial plane-wave tachyonic state with k ν > m ν , via a Lorentz transformation.
These considerations, together, afford an immediate explanation for the observation that the final result for the decay rate must necessarily vanish with the energy of the oncoming neutrino, and in fact, shall later be seen to contain the neutrino energy as a linear term. The calculation of the decay rate of the tachyonic neutrino needs to be done in the laboratory frame (lab frame).

E. Particle-Antiparticle Transformations and Tachyonic Decays
A few final considerations regarding the time ordering of tachyonic world lines and the calculation of the decay rate are in order. As already emphasized, decay rates are normally calculated most easily in the rest frame of the decaying particle. For tachyons, we cannot go into the true rest frame of the decaying particle, because the frame with k ′ = 0 cannot be reached for a tachyon, as already described. In the case of a tardyonic particle, there is an energy gap of twice the rest mass between the spectrum of positive-energy (particle) versus negative-energy (anti-particle) states. This energy gap vanishes for tachyons. A tardyonic oncoming particle state cannot transform into an incoming anti-particle state, irrespective of the Lorentz frame in which the process is observed. This implies that, e.g., for the decay of an oncoming muon into a muon neutrino, electron and electron anti-neutrino, there is no Lorentz frame in which the same process would be observed as a time-reversed process, i.e., the annihilation of an incoming muon antineutrino and an incoming muon, into a W boson, and the eventual production of an electron and an electron antineutrino.
Furthermore, it is known that the energy of a tachyonic particle may change sign upon a Lorentz transformation (see Fig. 3), so that particle trajectories may become anti-particle trajectories (with the time ordering of start and end points reversed). Indeed, the fact that some particle creation and annihilation operators transform into anti-particle operator upon a Lorentz transformation, has been mentioned as an important problematic aspect of tachyonic field theories [13][14][15][16][17][18][19][20][21], while possible re-interpretations have recently been discussed in Ref. [6].
For a tachyonic decay of an oncoming initial tachyonic neutrino into an electron-positron pair, and an energetically lower neutrino, this means the following. The interpretation of the process may depend on the Lorentz frame in which it is observed; tachyonic trajectories have no definite time ordering. (The only ordering in the tachyonic case concerns the helicity: A left-handed particle state will transform into a right-handed anti-particle state, and vice versa.) The decay of a highly energetic oncoming neutrino ("Big Bird", see Refs. [22,23]) into a energetically lesser one ("Tweety") via electron-positron pair production is interpreted equivalently as the annihilation of an incoming tachyonic neutrino, and an incoming tachyonic anti-neutrino, in specific Lorentz frames [see Fig. 2 In other Lorentz frames, it is even reinterpreted as the decay of an incoming highly energetic anti-neutrino, into a less energetic anti-neutrino and an electron-positron pair [see Fig. 2 We now consider the kinematics of the tachyonic decays displayed in Fig. 2 in detail. In Fig. 2(a), the world-line trajectories of the oncoming neutrino (joining space-time points labeled 1 and 2) and of the final zero-energy neutrino (joining space-time points labeled 2 and 3) are displayed. When "chasing" the decaying neutrino with a Lorentz boost, transforming the x and t axes into x ′ and t ′ , respectively, then from visual inspection, it is evident that the time ordering on the final decay product trajectory has reversed (t 2 = t 3 but t ′ 3 < t ′ 2 ). The decay product has turned into an incoming antineutrino, and the Lorentz-transformed process describes neutrino-antineutrino annihilation (into an electron-positron pair, but the world lines of the decay products are not displayed in Fig. 2). Physical reality has to be ascribed to both interpretations [13][14][15][16]. The observation of the moving ("primed") observer is equally valid. For the lab frame, this means that unless we have a counter-propagating beam of antineutrinos, the neutrino-antineutrino annihilation process does not contribute to the discussion of the "decay", which only converts the oncoming highly energetic neutrino into one with lesser energy.
Let us now consider Fig. 2(b). The incoming neutrino is chased by a "faster" Lorentz boost. The transformed axes become x ′′ and t ′′ , and the first, the decaying neutrino, now constitutes a zero-energy decay product, for the decay of an incoming antineutrino (time-ordered trajectory 3 → 2). As the boost velocity crosses the x ′′ and t ′′ axes, the decay "product" (from the point of view of the lab frame) has turned into an incoming, highly energetic, antineutrino, which in the triple-primed frame in Fig. 2(c), decays into an energetically lower antineutrino (from the point of view of the boosted frame). At the point where the x ′′ and t ′′ axes are crossed, the initial, incoming neutrino has transformed into an outgoing zero-energy neutrino or anti-neutrino state (the interpretation changes exactly at the point where the energy changes sign).
What do these considerations imply for the description of the tachyonic decay of a neutrino? We are working in the lab frame, and we need to calculate the process in the lab frame. Processes with incoming anti-neutrinos must be excluded from the integration, because they cannot contribute to the decay of an incoming neutrino. The interpretation of a process involving tachyons may depend on the Lorentz frame; for the calculation of the decay rate, only processes with incoming and outgoing positive-energy neutrino states may be considered, even if these states may transform into anti-particle states upon a Lorentz transformation. The final results are still Lorentz-invariant, as discussed below in Sec. V C.

A. Interaction Terms in Glashow-Weinberg-Salam Theory
In order to proceed to the calculation of the decay rate of the tachyonic, incoming neutrino, we briefly compile known Lagrangians from standard electroweak theory (see also Appendix A). We denote the weak coupling constant as g w . According to Chap. 12 of Ref. [50], quantum electrodynamics (QED) is described by the coupling of the electron to the photon, where θ W is the Weinberg angle and e and e describe the electron-positron field operators, while A µ is the electromagnetic field operator. Furthermore, the charged vector boson W ± interacts with a neutrino-electron current, where the addition of the Hermitian adjoint is necessary in order to include the W + boson. For the calculation of the muon decay, one needs the full Lagrangian given in Eq. (III.2), even twice, namely, once for the muon-muon-neutrino current, and a second time for the decay of the W into the electron and electron antineutrino, i.e., the same current is used in the electron and in the neutrino sector.
For the decay process of the tachyonic neutrino, one needs the coupling of the neutrino to the Z 0 boson, as well as the coupling of the left-and right-handed electron to the Z 0 boson, The latter form allows us to identify the vector-coupling and axial-vector coupling coefficient c V and c A . According to Eq. (12.237) of Ref. [50], the vacuum-expectation value v of the Higgs, the weak coupling constants g w and g ′ w , and the masses of the vector gauge bosons W ± , Z 0 and A, are related by These values match the experimental observations of M W = 80.385(15) GeV/c 2 and M Z = 91.1876 (21) GeV/c 2 . The matching with Fermi's effective coupling constant is given as Let us anticipate a certain consideration regarding the prefactors encountered in the calculation of invariant matrix elements, in the weak decay of the muon, and in the weak decay of a tachyonic neutrino. For the weak decay of the muon, one uses the Lagrangian (III.2), whose prefactors give a numerical factor 1/(2 √ 2) 2 = 1/8. For the weak decay by LPCR, we need to use the Lagrangians (III.3) and (III.4), whose combination results in a prefactor However, the weak decay of the tachyonic neutrino is mediated by a Z boson as opposed to a W boson, which results in a factor which is the same prefactor that we encounter in the invariant matrix element for the weak decay of the muon.

B. Degrees of Freedom in Three-Body Decay
Let us analyze the degrees of freedom in the phase-space of the final state, in three-body decay of a tachyonic neutrino into a less energetic neutrino, and a light fermion-antifermion pair. The momentum transfer is q 2 > (2m e ) 2 from the first fermion line. The decay rate is then obtained as an integral over the differential decay rate, This decay rate is 9 times differential, with 4 conservation conditions. We thus have 5 effective free variables. These can be assigned as follows: For the decay of a tachyonic neutrino via pair production, we may fix the three momentum components of the outgoing neutrino. Because both the incoming as well as the outgoing neutrino have to be on the mass shell, this fixes the four-vector q µ = p µ 1 − p µ 3 = (q 0 , q) completely. We can then go into the rest frame of the virtual Z 0 boson and argue that the decay must be completely symmetric there; i.e., the electron and positron should come out in directions exactly opposite of each other. This gives us two more degrees of freedom, namely, the polar and azimuthal angles of one of the outgoing fermions.
The three momenta of the outgoing neutrino and the two light fermion angles add up to the five effective degrees of freedom. So, once we have q µ = (q 0 , q), we have only two degrees of freedom left for the electron-positron pair.

C. Rationale of the Investigation
The rationale behind the calculations reported below can be summarized as follows. We shall approach the eventual calculation of the decay rate of a tachyonic neutrino due to electron-positron pair production in two steps.
• Step 1 (Complexities in the lab frame): As already emphasized, the tachyonic calculation, in which we are eventually interested, requires us to consider amplitudes in the lab frame, as opposed to the rest frame of the decaying particle. We thus need experience with calculations in the lab frame. The calculation of the muon decay rate is in principle very well known for the rest frame of the decaying particle. Here, we generalize the calculation to a muon decay rate calculation in the lab frame, where as we shall see, the allowed k 3 momenta (in the conventions used for Fig. 1) are inside an ellipsoid. Lorentz invariance of the integral over the allowed outgoing momenta is explicitly shown. • Step 2 (Decay of tachyonic, space-like particles): In the calculation of the decay rate of the tachyonic neutrino, we assume (in the spirit of Fig. 1) that both the incoming as well as the outgoing neutrinos are on the tachyonic mass shell, Under these circumstances, tachyonic decay is made possible exclusively due to the mass terms in the dispersion relations; hence we cannot ignore these terms. Furthermore, as we have already discussed, we need to remember that the region with | k 3 | < m ν actually is excluded from the region of allowed tachyonic momenta. We find that the physically allowed outgoing momenta are located inside the rotationally symmetric, shallow hull of a cupola-like structure, centered about the axis of the oncoming (decaying) neutrino (which we choose to be the positive z axis). As already anticipated in Sec. II E, we shall need to explicitly exclude from the calculation any processes related to neutrino-antineutrino annihilation. This necessity, in turn, makes the use of the explicit spinor solutions of the tachyonic Dirac equation [6,9] necessary.
In the calculation, we also need to overcome the pitfall connected with the time ordering of the tachyonic trajectories, anticipated in Sec. II E. We shall consider the weak decay of the muon, in the conventions of Fig. 4. The interaction terms from the electroweak standard model is used according to Eq. (III.2). For momentum transfers q 2 ≪ M 2 W , the effective four-fermion Lagrangian thus is where we use the matching (III.6). The Lorentz-invariant matrix element thus is (within the conventions of Fig. 4) Summing over the final spin states and averaging over the spin projections of the initial state leads to The procedure of integration in the rest frame of the decaying particle is discussed in Eq. (10.16) ff. of Ref. [12]. Furthermore, a mixed approach, where certain intermediate integrals are carried out covariantly, and only the final stages of the calculation are carried out in the rest frame of the decaying particle, is outlined in Chap. 7.2.2 of Ref. [51].
In the actual evaluation, in the conventions of Fig. 4, we keep the outgoing neutrino momentum as our final integration variable and write the decay rate in the lab frame as follows [for the expression in the lab frame, see Ref. [52]], We have used the following result, derived in Eq. (B.1), which is obtained for two outgoing particles with labels 2 and 4 which are on the electronic mass shell (III.14) By symmetry, we have J λρ (q) = J λρ (q). We have carried out the p 2 and p 4 integrals covariantly. Then, for the remaining integral over p 3 , we need the appropriate integration limits. We thus need to integrate assuming an incoming muon with energy E 1 in the positive z direction, with the incoming p 1 on the muon mass shell, p 2 1 = m 2 µ . The final p 3 describes the muon neutrino, so that within our approximations (p 3 ) 2 ≈ 0. The domain of the d 3 p 3 integration in Eq. (III.15) contains all four-momenta p 3 for which q 2 = (p 1 − p 2 ) 2 > 0.
In the rest frame of the decaying muon, the integration domain would consist of a sphere composed of vectors p 3 = (| k 3 |, k 3 ), with p 1 = (E 1 , 0) and | k 3 | ≤ m ν /2. By contrast, in the lab fame, we consider a muon moving up the z axis, with energy E 1 and wave vector k 1 , and an outgoing muon neutrino with energy E 3 and wave vector k 3 , The momentum transfer reads as follows, The allowed vectors k 3 are located inside a rotationally symmetric ellipsoid (see Fig. 5), which is centered at the point (0, 0, k z0 ) on the z axis. Let us denote by a the half axis of the ellipsoid in the radial direction ("away" from the z axis) and by b the half axis of the ellipsoid in the z direction (see Fig. 5). These half axes are given as follows, For the integration of the final phase-space in the expression (III.15), we need to calculate One uses the following parameterization (k x ≡ k 3x , k y ≡ k 3y , and k z ≡ k 3z ), k x = a ξ sin θ cos ϕ , k y = a ξ sin θ sin ϕ , k z = k z0 + b ξ cos θ . (III.20) The Jacobian is (III.21) Then, with u = cos θ and k 1 = χm µ , one has after the trivial integration over ϕ, After a somewhat tedious u integration, the result can be written in terms of the variable The decay rate is naturally written as Γ = Γ 1 + Γ 2 , where the integration domains are such that ξ 1 + χ 2 − χ assumes either of the values given in Eq. (III.23). Here, (III.23b) The two contributions evaluate to the expressions, so that This is the expected result for the muon decay width, with the 1/E 1 prefactor already included, which is here obtained directly by an explicit integration in the lab frame.
We have also verified [25] the results of Cohen and Glashow for superluminal neutrino decay, based on the noncovariant dispersion relation used in Ref. [31], via an independent calculation in the lab frame, as envisaged in Ref. [33]. The treatment described in Ref. [31] is based on a Lorentz-violating dispersion relation E = | k| v ν with v ν > 1, which constitutes a fundamentally different theoretical model as compared to the Lorentz-invariant tachyonic treatment presented here. In addition to the decay rate, we shall also consider the energy loss rate of an incoming neutrino beam in the lab frame, due to LPCR and NPCR. A remark is in order: The calculation of the energy loss per time of an incoming muon beam, to complement a corresponding calculation for the tachyonic neutrino beam, is not applicable, because the end product of the decay is not a less energetic muon, but the muon disappears from the beam altogether (see Fig. 4). 2) 2 is incoming along the positive z direction. The electron mass is set equal to unity, me = 1. The threshold condition therefore reads as q 2 = (p1 − p3) 2 ≥ 4. The boundaries of the volume of allowed k3 vectors, with Cartesian components kx = k3x, ky = k3y, and kz = k3z, are mainly concentrated in a narrow, rotationally symmetric cone about the z axis. Final states with | k3| < mν correspond to evanescent outgoing waves, lead to a complex-valued momentum transfer, and have to be excluded.

E. Step 2: Tachyonic Neutrino Decay (Covariant Dispersion Relation)
We calculate the decay width of the incoming tachyonic neutrino, in the lab frame, employing a relativistically covariant (tachyonic) dispersion relation, with both incoming as well as outgoing neutrinos on the tachyonic mass shell [E 1 = ( k 2 1 − m 2 ν ) 1/2 , and E 3 = ( k 2 1 − m 2 ν ) 1/2 , in the conventions of Fig. 1]. In the lab frame, in full accordance with Ref. [52], the decay rate simply is Here, spins refers to the specific way in which the average over the oncoming helicity states, and the outgoing helicities, needs to be carried out for tachyons. As explained in Sec. II D, we cannot reach the rest frame of the decaying particle by a Lorentz transformation, in the case of a tachyonic neutrino. Furthermore, as outlined in Sec. II C, a calculation with just the Fermi effective coupling constant actually is sufficient for the tachyonic case. For the relevant interaction terms, we use the same expression as in Sec. III E. We recall the coupling of the decaying neutrino to the Z 0 boson according to Eq. (III.3), For the coupling of the electron to the Z boson, we have according to Eq. (III.4), (III. 28) In view of the compensation mechanism given by Eq. (III.8), the effective four-fermion Lagrangian thus is given by The matrix element of the fundamental spinor solutions reads as follows, Here, the u T (p 1 ), and u T (p 3 ), constitute Dirac spinor solutions of the tachyonic Dirac equation. In the helicity basis [9,53], denoted by a subscript σ = ±, the tachyonic particle and antiparticle spinors are where the a σ ( k) are the fundamental helicity spinors (see p. 87 of Ref. [50]).
The properties of the tachyonic bispinor solutions differ somewhat from those of the "normal" tardyonic bispinors. The well-known sum formula for the tardyonic states is (σ denotes the spin orientation) σ u σ (p) ⊗ u σ (p) = ✁ p + m e . (III.32) For the tachyonic spin sums, one has the following sum rule for the positive-energy spinors, where we use p = (E, p) as the convention for the four-momentum andp = p/| p| is the unit vector in the p direction. Upon promotion to a four-vector, we havep µ = (0,p). The sum rule can thus be reformulated as where τ = (1, 0, 0, 0) is a time-like unit vector. In Refs. [6,9], it has been established that a consistent formulation of the tachyonic propagator is achieved when we postulate that the right-handed neutrino states, and the left-handed antineutrino states, acquire a negative Fock-space norm after quantization of the tachyonic spin-1/2 field. Hence, in order to consider the decay process of an oncoming, left-handed, positive-energy neutrino, we should consider the projection onto negative-helicity states, The squared and spin-summed matrix element is with c V ≈ 0, c A ≈ −1/2 (the last step implicitly defines the expression S). Here, the meaning of the notation spins becomes clear: We have summed over the spins of the outgoing electron-positron pair, while only one specific helicity is taken into account for the oncoming (decaying) neutrino. The Dirac γ traces in Eq. (III.36) give rise to a rather lengthy expression, which can be simplified somewhat because incoming and outgoing particles are on their respective mass shells, p 2 2 = p 2 4 = m 2 e , while on the tachyonic mass shell, we have p 2 1 = p 2 3 = −m 2 ν . Some other scalar products vanish, e.g., the scalar product of the time-like unit vector τ and the space-like unit vector, which is τ ·p = 0.
The result of the Dirac γ traces from Eq. (III.36) can then be inserted into Eq. (III.26), and the d 3 p 2 and d 3 p 4 integrals can be carried out using the following formulas, which we recall from Appendix B [see Eq. (B.1)], (III.38) After the d 3 p 2 and d 3 p 4 integrations, we are left with an expression of the form Both the expressions for S(p 1 , p 2 , p 3 , p 4 ) as well as F (p 1 , p 3 ) are too lengthy to be displayed in the context of the current paper. However, approximate formulas can be given, e.g., when the incoming energy E 1 is near threshold.
In order to obtain a better intuitive picture for the domain of allowed p 3 four-vectors, we have to analyze the tachyonic kinematics in some more detail. We calculate in the lab frame and assume that the oncoming neutrino has the energy-momentum four-vector The condition k 3 > m ν is naturally imposed for tachyonic kinematics (see Fig. 3). The squared four-momentum transfer then reads as where it is convenient to define u = cos θ. One may solve for the threshold angle cos θ th , For given E 1 and E 3 , all angles θ with cos θ > cos θ th , i.e., for θ < θ th , are permissible. Conversely, setting cos θ th = 1 and E 3 = 0, one may solve for E 1 and rederive the threshold condition (II.3). All of this implies that the domain of permissible k 3 vectors, near threshold, is centered about the z axis and forms a "cupola" of inner radius m ν (see Fig. 6). Within the kinematically allowed region, the tachyonic momentum transfer q 2 is plotted in Fig. 7.
For given E 1 , the widest opening angle θ = θ th is reached when E 3 becomes zero. One finds where the latter form is valid in the high-energy limit. Here, k th = 2m 2 e mν + m ν is the threshold momentum. It means that the produced pairs will be emitted in a very narrow cone, centered in the forward direction with respect to the decaying neutrino.
Maximum squared momentum transfer is reached at E 3 = 0 and θ = 0, confirming Eq. (II.12). Maximum outgoing energy E 3 is reached for minimum momentum transfer q 2 = q 2 min = 4m 2 e , with the final-state neutrino propagating into the positive z direction. Its energy is where the latter form is valid in the limit of large E 1 . This corresponds to the maximum allowed k 3 , (III.49) where, again, the latter form is valid in high-energy limit k 1 → ∞. A plot of the physically relevant range for q 2 is given in Fig. 7. Because of azimuthal symmetry, it is easily possible to find a convenient parameterization of the regime of allowed k 3 in spherical coordinates [using the parameterization given in Eq. (III. 43)]. We may finally express the integrand F from Eq. (III.41) in terms of the initial and final energies E 1 and E 3 of the decay process, and of the scattering angle θ (with u = cos θ). The decay rate given by Eq. (III.40) can thus be written as follows, Here, (k 3 ) min = m ν , while (k 3 ) max is given by Eq. (III.49). Furthermore, we have (E 3 ) min = 0, while (E 3 ) max is given by Eq. (III.48), and we have used the identity It is instructive to consider the double-differential energy loss d 2 E 1 , for a particle traveling at velocity v ν ≈ c (we restore factors of c for the moment), as it undergoes a decay with energy loss E 1 − E 3 , due to the energy-resolved decay rate (dΓ/dE) dE, in time dt = dx/c. It reads as follows, Now we revert to natural units with c = 1, divide both sides of the equation by dx and integrate over the energy loss. One obtains Hence, the energy loss rate is obtained as After a long, and somewhat tedious integration one finds the following expressions, which have been briefly indicated in Ref. [25], In the high-energy limit, one may rewrite the expressions as follows, The ratio of the energy loss rate to the decay rate is given as (III.57) Here, according to Eq. (II.4), the threshold energy is For all results given in Eqs. (III.55), (III.56) and (III.57), we have assumed that m e ≫ m ν . Interpolating formulas are given as . (III.59) These formulas interpolate between the regimes E 1 E th and E 1 ≫ E th given in Eq. (III.55).

A. Preliminary Steps
Having laid out the formalism in Sec. III, we can be brief in the current section. In the lab frame, again, the decay rate evaluates to Here, spins refers to the specific way in which the average over the oncoming helicity states, and the outgoing helicities, needs to be carried out for tachyons [25]. We use the Lagrangian (III.3). The effective four-fermion interaction thus is The matrix element M evaluates to in the notation for the tachyonic bispinors adopted previously. We use, again, the helicity-projected sum rule where τ = (1, 0, 0, 0) is a time-like unit vector. This leads to The squared and spin-summed matrix element for the tachyonic decay process thus is (IV.6) Again, we have chosen the convention to denote the by p 2 the momentum of the outgoing antiparticle.

B. Integration and Results
We now turn to the integration over the four-momenta of the outgoing particles. In the calculation, one may use the fact that the helicity projector is well approximated equal to the chirality projector for tachyonic particle in the high-energy limit (with the energy being significantly larger than the tachyonic mass). On the tachyonic mass shell, one has p 2 1 = p 2 2 = p 2 3 = p 2 4 = −m 2 ν . The trace over the Dirac γ matrices can be evaluated with standard computer algebra [54,55] The expressions for S(p 1 , p 2 , p 3 , p 4 ) as well as F (p 1 , p 3 ) are too lengthy to be displayed in the context of the current paper. We assume the same kinematics as in Eqs. (III.42) and (III.43). The integrations are done with under the conditions that all 0 < E 3 < E 1 , and all q 2 = (p 2 + p 4 ) 2 for the pair are allowed [see Eq. II.10], leading to in full analogy with Eq. (II.10). Finally, the energy loss rate is obtained in full analogy with Eq. (III.54), After rather tedious integration one finds the following expressions, Strictly speaking, these formulas are valid only for E 1 ≫ m ν , but this condition is easily fulfilled for all phenomenologically relevant neutrino energy, assuming that the modulus of neutrino masses |m ν | does not exceed 1 eV. We reemphasize that, unlike in Eq. (III.55), there is no further threshold condition. Parametrically, the results in Eqs. (III.55) and (IV.11) are of the same order-of-magnitude. Hence, neutrino pair emission is the dominant decay channel in the medium-energy domain, for an oncoming tachyonic neutrino flavor eigenstate.

A. Decay Processes on Cosmic Distance and Time Scales
If we assume that the tachyonic neutrino hypothesis is real, then a natural question to ask concerns the phenomenological consequences of the calculations outlined above. Parametrically, the decays by LPCR and NPCR described by Eqs. (III.55) and (IV.11) might set important limits on the observability of tachyonic neutrinos, provided the absolute magnitude of the decay energy loss rates are sufficiently large in order to induce a significant decay probability for neutrinos traveling across the Universe. This is because neutrinos registered by IceCube have to "survive" the possibility of energy loss by decay, and if they are tachyonic, then lepton and neutrino pair Cerenkov radiation processes become kinematically allowed.
Indeed, it is known that even very small Lorentz-violating parameters in a Lorentz-violating extension of the standard model may induce very significant energy loss processes at high energies [39,40]. This is because at high energies, small violations of the Lorentz symmetry correspond to very high virtualities of the particles (kinematic deviations from the mass shell), and therefore, the magnitude of the Lorentz-violating parameters is in fact severely constrained by the 37 neutrinos with E > 60 TeV which are believed to be of cosmological origin and which have been registered by the IceCube collaboration [22,23]. Meanwhile, preliminary evidence for a through-going muon depositing an energy of ≥ (2.6±0.3) PeV has been presented by some members of the IceCube collaboration [56]. The event could be interpreted in terms of a decay product of a neutrino of even higher energy [56]. (The difference of the energy deposited inside the detector and the neutrino energy, according to Fig. 4 of Ref. [23], is small.) If confirmed, this event would lead to even more restrictive bounds on the Lorentz-violating parameters.
The results given in Eqs. (IV.11) for the decay rate and energy loss rate due to NPCR are not subject to a threshold energy. Parametrically, they are of the same order-of-magnitude as those given for lepton pair Cerenkov radiation in Eq. (III.55), but the threshold energy is zero. Let us take as a typical cosmological distance 15 billion light years, L = 15 × 10 9 ly = 1.42 × 10 26 m , (V.1) and assume a (relative large) neutrino mass parameter of m 0 = 10 −2 eV. One obtains for the relative energy loss according to Eq. (IV.11a), Even at the large "Big Bird" energy of E ν ≈ 2 PeV, the relative energy loss over 15 billion light years does not exceed 5 × 10 −20 . Again, assuming that m ν = 10 −2 eV), one obtains for the decay rate the result Even for E ν ≈ 2 PeV, this means that the decay rate is only of order 10 −28 rad s , corresponding to a lifetime of ∼ 10 20 years, far exceeding the age of the Universe. Within the tachyonic model, quite surprisingly, both LPCR as well as NPCR are phenomenologically irrelevant, even for the highest-energy neutrinos registered by IceCube.

B. Neutrino Mass and Flavor Eigenstates
The tachyonic Dirac equation [1,6] reads as where the latter form holds for the plane-wave ansatz ψ(x) = u T (p) exp(−ip · x). The bispinor solutions u T (p) have been discussed at length in Refs. [6,9] and are used here in Eqs. (III.30) and (IV.3). They apply, first and foremost, to a mass eigenstate, with a definite tachyonic mass parameter. The Fermi couplings are universal among all neutrino flavors, and hence, the interaction Lagrangians used in our paper share this property. Our results (III.55) and (IV.11) for the decay and energy loss rates thus apply, at face value, to an incoming neutrino mass eigenstate. The results are thus relevant to the non-sterile neutrino flavor if at least one of the three observed non-sterile neutrino mass eigenstates is tachyonic.
We recall that the flavor eigenstates ν f are connected to the mass eigenstates m i by the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix, where U f i denote the elements of the flavor-mass mixing matrix. The decay and energy loss rates of the flavor eigenstates are given as For a slower-than light mass eigenstate i, one sets Γ(m i ) and dE dx (m i ) to zero. Here, just to be pedantic, we should point out that the calculation of NPCR in this case has to be modified to include all tachyonic mass eigenstates in the exit channel, conceivably modifying the overall results as much as by adding a multiplicative factor three (if all mass eigenstates are available in the exit channel of the tachyonic NPCR decay, see Fig. 1(b)).

C. Lorentz Invariance
The tachyonic dispersion relation E 2 − k 2 = −m 2 ν conserves Lorentz invariance. Hence, one might ask about the Lorentz invariance of our results, and in particular, about the Lorentz invariance of the threshold condition (II.5); finally, one might "chase" the high-energy neutrino, lowering its energy in the Lorentz-transformed, moving frame to a value below threshold. This question finds an answer in the subtleties of the tachyonic theory; we follow the discussion in Ref. [17]. Namely, upon a Lorentz transformation of the vacuum state, because there is no "energy mass gap" between the positive-and negative-energy states, some of the annihilation operators of quantized fields will turn into creation operators, and vice versa. This point is explained in detail around Eqs. (4.7)-(4.9) of Ref. [17]. (Incidentally, it is observed at the same place that the fundamental creation and annihilation operators of tachyonic fields have to be quantized according to fermionic statistics, which is another argument in favor of spin-1/2 rather than spinless tachyonic theories.) Furthermore, around Eq. (5.7) of Ref. [17], it is argued that the vacuum state in a tachyonic theory cannot be Lorentz invariant, but is filled with those (real) anti-fermions whose energies are "pushed down" to energies below zero, from initially positive-energy states, under the Lorentz transformation. Our Fig. 2 illustrates how the decay and energy loss rates, under a Lorentz transformation, turn into neutrino-antineutrino collision rates (leading to decay and energy loss) with the "downshifted" real antiparticle states which are the result of the Lorentz transformation, finally restoring the Lorentz invariance of the results for the decay and energy loss rates, given in Eqs. (III.55) and (IV.11).

D. Superluminal Signal Propagation
A very important question regarding the conceivable existence of tachyonic neutrinos concerns the possibility of superluminal signal propagation. We thus follow Appendix A of Ref. [47] and ask how difficult it is to reliably "stamp" any information onto the superluminal neutrinos. When assuming the dispersion relation E = ( k 2 − m 2 ν ) 1/2 with its classical equivalent E = m ν / v 2 ν − 1 = 0, the dilemma is that high-energy tachyonic neutrinos approach the light cone and travel only infinitesimally faster than light itself. In the high-energy limit, their interaction cross sections may be sufficiently large to allow for good detection efficiency but this is achieved at the cost of sacrificing the "speed advantage" in comparison to the speed of light. Low-energy tachyonic neutrinos may a substantially faster than light, but their interaction cross sections are small and the information sent via them may be lost. The smallness of the cross sections sets important boundaries for the possibility to transmit information, as follows. In Appendix A of Ref. [47], it has been shown that, by postulating that superluminal particles should not have the capacity to transport any "imprinted" information into the past, one is naturally led to the assumption that any conceivable superluminal particles have to be very light, and weakly interacting.
[58], we now supplement these considerations with a numerical estimate. Neutrino-electron cross sections for 1 GeV < E ν < 1 PeV can be estimated to good accuracy using the formula with A 0 ∼ 0.0095 fb and E 0 = 1 GeV. By order-of-magnitude, Eq. (V.8) remains valid for neutrino scattering off electrons, for all three neutrino flavors, even if additional charged-current interactions exist for electron neutrinos, due to exchange graphs with virtual W bosons (for muon and tau neutrinos, only the Z boson contributes at tree level). A particle typically cannot be localized to better than an area equal to the square of its (reduced) Compton wavelength (we temporarily restore factors of and c), The detection probability P for a perfectly focused particle therefore cannot exceed If we are to send information reliably, then the detection probability should be of order unity. Setting P = 1 leads to When traveling at a speed c + δc for a path length s, the neutrino acquires a path length difference of δs, which compares to its Compton wavelengthλ = /(m ν c) as follows, The (V.13) Even at a (larger-than-realistic) mass square m 2 ν = 1 eV 2 for the tachyonic neutrino flavor eigenstate, the value of s 0 ∼ 10 74 m far exceeds the commonly assumed size of the Universe of 10 26 m by many orders of magnitude. The permissibility of slightly superluminal propagation on small length and distance scales has been discussed in the literature previously (see, e.g., Ref. [59]). Furthermore, we refer to the experiments in the group of Nimtz [60][61][62], which also use a compact apparatus and rely on the quantum mechanical tunneling effect, which lies outside the regime of classical mechanics. Thus, a very slightly superluminal neutrino flavor eigenstate with a light mass does not necessarily lead to a detectable violation of causality.

VI. CONCLUSIONS
In the current article, we have considered the tachyonic neutrino decay width against lepton-pair and neutrinopair Cerenkov radiation (LPCR and NPCR, see Fig. 1), via the exchange of a virtual Z 0 boson. This process is kinematically allowed for a fast-than-light, oncoming neutrino. We use the hypothesis of tachyonic neutrinos described by the tachyonic Dirac (not Majorana) equation [1][2][3][4][5]. Various kinematic considerations are summarized in Sec. II. The tachyonic threshold is found according to Eq. (II.5), E th ≈ 2 m 2 e /m ν , in the limit m e ≫ m ν . Specificities of the tachyonic decay are studied in Sec. II A (threshold calculation), Sec. II B (absence of threshold for NPCR), Sec. II C (maximum q 2 of the Z 0 boson and validity of Fermi theory), and Sec. II D (rest frame of the tachyon). In Sec. II E, it is shown that, because tachyonic particle states may transform into antiparticle states upon a Lorentz transformation, it is indispensable to carry out the calculation directly in the lab frame [17,18].
We continue with a discussion of the interaction Lagrangians relevant for our studies, from the GWS (Glashow-Weinberg-Salam) model in Sec. III A. After a brief digression on the degrees of freedom of three-particle decay processes in Sec. III B and a discussion on the general rationale of the investigation in Sec. III C, the calculation of the tachyonic decay width is approached in two steps. In Sec. III D, we first demonstrate that it is possible to carry out standard decay rate calculations of the electroweak theory, directly in the lab frame, using the muon decay width as an example. We are finally in the position (see Sec. III E) to carry out the integration of the decay rate, for the tachyonic dispersion relation, in the lab frame. We find an explicit dependence of the formulas for the decay rate, Γ, and the energy loss rate, dE/dx, on the energy E 1 of the incoming neutrino (all decay processes are studied within the conventions from Fig. 4). The main results of our investigations are summarized in Eqs. (III.55), (III.56) and (III.57); these formulas describe the decay width of a tachyonic neutrino against LPCR, and the energy loss per distance of an incoming tachyonic neutrino beam. This investigation is supplemented, in Sec. IV, by a calculation of NPCR, culminating in the results given in Eq. (IV.11) for the decay and energy loss rates.
In Sec. V, we find that the neutrino pair Cerenkov radiation (NPCR) process, even if threshold-less, has such a low probability due to the weak-interaction physics involved, that it cannot constrain the tachyonic models, even for large tachyonic neutrino mass parameters of the order of 10 −2 eV. The lifetime of a tachyonic neutrino against LPCR and NPCR, assuming a realistic magnitude of the mass parameter, far exceeds the age of the Universe. Even a "Big Bird" neutrino of energy of E ν ≈ 2 PeV, would easily survive the travel from the blazar PKS B1424-418 (see Ref. [63]. In contrast to Lorentz-violating models, neutrino pair Cerenkov radiation does not pressure the tachyonic neutrino hypothesis. According to Sec. V, we should take the opportunity to clarify that in contrast to Ref. [25], it is actually impossible to relate a hypothetical cutoff of the cosmic neutrino spectrum at the "Big Bird" energy of 2 PeV to the threshold energy for (charged) lepton pair Cerenkov radiation, and thus, to a neutrino mass parameter. The reasons are twofold: First, a further decay process exists for tachyonic neutrinos which is not subject to a threshold condition, namely neutrino pair Cerenkov radiation. Second, the decay and energy loss rates for both (charged) lepton as well as neutrino pair Cerenkov radiation simply are too small to lead to any appreciable energy loss for an oncoming tachyonic neutrino flavor eigenstate, over cosmic distances and time scales. Formulated differently, we can say that that neither lepton nor neutrino pair Cerenkov radiation processes pressure the tachyonic model in any way.
Finally, we hope that the detailed outline of the calculation of the decay processes given in Sec. III and IVk could be of interest in a wider context, regarding decay processes and cross sections involving tachyonic spin-1/2 particles. It is indispensable to introduce further helicity projectors in the calculation of the bispinor matrix elements relevant to the process, and the calculations become a little more complex than for ordinary Dirac spinors [see Eqs. (III.33)-(III. 35)]. Our approach relies on a consistent formalism developed for the fundamental tachyonic bispinor solutions, as reported in various recent investigations [6][7][8][9][10].

ACKNOWLEDGMENTS
Helpful conversation with J. H. Noble and B. J. Wundt are gratefully acknowledged. This research was supported by the National Science Foundation (Grant PHY-1403973 and PHY-1710856). This work was also supported by a János Bolyai Research Scholarship of the Hungarian Academy of Sciences.
This result also is in agreement with Eq. (5.57) on p. 153 of Ref. [64], up to an overall minus sign which is fixed by the conventions. According to p. 107 of Ref. [65], the effective Weinberg angle reads as The W and Z masses are connected by virtue of the Weinberg angle, according to Eq. (III.5).

B. COVARIANT PAIR PRODUCTION INTEGRALS
In our evaluation, for the outgoing electron-positron pair in the decay of the tachyonic neutrino, we shall need a few integrals. In the conventions of Fig. 1, the outgoing momenta p 2 and p 4 are on the mass shell, E 2 = k 2 2 + m 2 e and E 4 = k 2 4 + m 2 e . Let us anticipate the results for the integrals I, J, and K, which are defined as follows, By symmetry, one immediately has J ρλ (q) = J λρ (q). The evaluation of these integrals is essentially simplified because of the Lorentz invariance of the integration measures, which entails the possibility to choose a coordinate system where q = (q 0 , q = 0), and then, identify the occurrences of (q 0 ) 2 with q 2 . The derivation of the results is discussed below.
We observe that because E 2 = k 2 2 + m 2 e and E 4 = k 2 4 + m 2 e , we have We shall go through the calculation if the integral I(q) in great detail, For J λρ (q), we write Projection onto the tensors g λρ and q λ q ρ leads to g λρ J λρ (q) = q 2 (4A + B) , q λ q ρ J λρ (q) = q 4 (A + B) .  Yet,