Nonextensive Tsallis statistics in Unruh effect for Dirac neutrinos

Flavor mixing of quantum fields was found to be responsible for the breakdown of the thermal Unruh effect. Recently, this result was revisited in the context of nonextensive Tsallis thermostatistics, showing that the emergent vacuum condensate can still be featured as a thermal-like bath, provided that the underlying statistics is assumed to obey Tsallis prescription. This was analyzed explicitly for bosons. Here we extend this study to Dirac fermions and in particular to neutrinos. Working in the relativistic approximation, we provide an effective description of the modified Unruh spectrum in terms of the q-generalized Tsallis statistics, the q-entropic index being dependent on the mixing parameters sin {\theta} and \Delta m. As opposed to bosons, we find q>1, which is indicative of the subadditivity regime of Tsallis entropy. An intuitive understanding of this result is discussed in relation to the nontrivial entangled structure exhibited by the quantum vacuum for mixed fields, combined with the Pauli exclusion principle.


I. INTRODUCTION
The Unruh effect states that an observer moving through the inertial vacuum with uniform acceleration a (Rindler observer) perceives a thermal radiation at temperature [1].
Recently, the growing interest in the Unruh effect has led to unveil unconventional features of this phenomenon. For instance, in [28] it was argued that the Casimir-Polder forces between two relativistic uniformly accelerated atoms exhibit a transition from the short distance thermal-like behavior predicted by the Unruh effect to a long distance non-thermal behavior. Similarly, in [18] the Unruh temperature T U was claimed to be sensitive to the degree of departure from linearity of the dispersion relation for the case of massive fields. Exotic behaviors of Unruh radiation also appear when considering the superposition of fields with different masses, as stated in [29,30]. Specifically, in [29] it was shown that the vacuum density of mixed particles detected by the Rindler observer deviates from the pure Planckian spectrum, the correction being quantified by the mass difference and the mixing angle. In turn, this feature can be ascribed to the peculiar nature of the vacuum for mixed fields (flavor vacuum), which appears as a condensate of entangled particle/antiparticle pairs due to the nontrivial structure of the mixing transformations at level of the ladder operators [31,32].
The mixing-induced departure of Unruh condensate from the Planckian profile was originally interpreted as a breakdown of the thermality of Unruh effect [29]. In a recent work [33], this result has been revisited in the context of the q-generalized Tsallis statistics based on the nonadditive Tsallis entropy [34][35][36][37]. As well-known, Tsallis entropy is a generalization of the standard Boltzmann-Gibbs definition parameterized by the q-index (in particular, q = 1 corresponds to the nonextensive Tsallis statistics, while q → 1 gives back Boltzmann-Gibbs framework). The concept was introduced as a basis for extending the usual statistical mechanics and is formally identical to Havrda-Charvát structural α-entropy considered in the information theory [38]. The rationale behind the study of the modified Unruh effect within Tsallis statistics is that the q-generalized entropy well describes systems exhibiting long-range interactions and/or spacetime entanglement, as is the case for mixed fields [31,32,39,40]. In this vein, in [33] it has been shown that the Unruh condensate for mixed particles can still be described as a thermal-like distribution, provided that the underlying statistics is assumed to obey Tsallis's prescription. This picture allows to relate the q-entropic index and the characteristic mixing parameters sin θ and ∆m in a nontrivial way. Clearly, in the absence of mixing (i.e. for vanishing θ and/or ∆m), the Boltzmann-Gibbs theory is recovered. In passing, we mention that applications of Tsallis entropy have been widely considered in literature. Among the various systems that clarify the physical conditions under which Tsallis entropy and the associated statistics apply, we quote self-gravitating stellar systems [41,42], black holes [36], the cosmic background radiation [43,44], low-dimensional dissipative systems [37], solar neutrinos [45], polymer chains [46] and cosmological models [47,48].
The description of Unruh effect in Tsallis's theory has been developed for mixing of boson fields. In that case, it has been found that q < 1, indicating a superadditive behavior of Tsallis entropy. In what follows, we extend this study to the case of Dirac fermions and in particular to neutrinos, which are the paradigmatic example of mixed particles. Again, we relate the q-index to the mixing angle and the mass difference, in such a way that the standard Fermi-Dirac distribution based on the extensive Boltzmann-Gibbs entropy is restored for vanishing mixing. Similarly to [33], we find a running behavior of q as a function of the energy scale, which is typical for QFT systems [47]. Nevertheless, as opposed to bosons, we now obtain q > 1, that is to say, the entropy function turns out to be subadditive. Following [33], we explore our result in connection with the entangled condensate structure exhibited by the quantum vacuum for mixed fields.
The remainder of the work is organized as follows: in Sec. II we review the quantization of Dirac field in Rindler spacetime and the origin of the Unruh effect. Section III is devoted to the study of the QFT formalism of field mixing in Rindler spacetime. This prepares the ground for the analysis of Sec. IV, where we introduce the framework of Tsallis thermodynamics and investigate the connection with the phenomenon of field mixing. To avoid technicalities, we make use of some approximations, such as those of relativistic neutrinos and small deviations of q from unity. Discussion of the results and outlook are summarized in Sec. V. Two Appendices containing details of calculations conclude the work. Throughout all the manuscript, we use the 4-dimensional metric with the mostly negative signature Furthermore we adopt natural units = c = k B = 1.

II. FIELD QUANTIZATION IN RINDLER SPACETIME: THE UNRUH EFFECT
Let us consider an observer moving in Minkowski spacetime along the x-axis with constant proper acceleration a > 0. From [49], it is known that his trajectory is given by parameterized by τ . It is easy to see that this is an hyperbola in the (t, x)-plane of asymptotes x = ±t, which act as event horizons for the accelerated observer. By varying a, we obtain different hyperbolas with the same geometric (and, thus, physical) properties. Starting from the Minkowski metric ds 2 = dt 2 − dx 2 − dy 2 − dz 2 , we now introduce the following set of coordinates with the other coordinates left unchanged. The metric then becomes which describes a stationary spacetime of timelike Killing vector ∂/∂η. By comparing Eqs. (3) and (4), the wordlines ρ = const, y = const, z = const correspond to the trajectories of uniformly accelerated observers with a = ρ −1 and proper time τ = aη. The collection of these worldlines forms the so-called Rindler spacetime (see the right R-wedge in Fig. 1). Notice that the hypersurfaces of η = const describe events which are simultaneous from the point of view of the Rindler (uniformly accelerated) observer.
In spite of the minimalistic setting, the geometric structure of Rindler manifold differs significantly from that of Minkowski spacetime. Indeed, an observer moving along a fixed hyperbola in R is causally disconnected from the L-sector, which means that he can neither send nor receive any signal from a source in that region. Of course, this does not apply to inertial (Minkowski) observers, which are indeed connected to the entire spacetime.
Causal structure of Rindler metric in 1 + 1 dimensions. We are assuming that the acceleration is along the x-axis. The right branch of hyperbola represents the worldine of the Rindler observer. The arrow indicates the direction of motion.
Strictly speaking, the set of coordinates (4)-(5) covers only the right sector of Minkowski spacetime. To describe the left L-wedge, one has to require ρ = − √ x 2 − t 2 < 0, with η running in the opposite direction. On the other hand, the future F and past P regions are covered by the charts where the ± signs refer to F and P , respectively. We emphasize that these charts are obtained by inverting the rôles of time and space coordinates in Eq. (4). For our next purposes, it is enough to consider the right and left wedges only.

A. Minkowski spacetime
Let us now turn to the problem of quantizing the Dirac field in Rindler spacetime. Toward this end, it is useful to review the canonical planewave expansion in Minkowski space. This is given by where x ≡ (x, y, z) denotes the set of spatial Minkowski coordinates. Here, ψ r + k = N u r k e −iω k t+ik·x and ψ r − k = N v r k e iω k t−ik·x are the positive and negative frequency modes with three-momentum k ≡ (k x , k y , k z ) and frequency ω k = m 2 + |k| 2 , respectively. The normalization coefficient is fixed to N = (2π) −3/2 . The following explicit form of the spinors u r k and v r k is adopted where we have used the usual notation / p = p µ γ µ , with γ µ being the gamma matrices in the Dirac representation and u r (0), v r (0) the unit spinors in the particle rest frame of k = (m, 0) (see [50]). They are normalized according to and similarly for v r k . In the second quantization formalism the operators a r k and b r k act as annihilators of the ordinary Minkowski vacuum |0 M , i.e.
while a r † k and b r † k create a particle/antiparticle with quantum numbers k and r. Such operators are assumed to satisfy the canonical anticommutation relations a r k , a s † with all other anticommutators vanishing. In turn, this implies the following relation for the equal-time anticommutator between the field and its hermitian conjugate where α and β are the spinorial indices.

B. Rindler spacetime
We now properly aim at performing the second quantization of the spinor field in Rindler spacetime. To fix the notation, we closely follow [51] and consider the solutions of the Dirac equation in Rindler coordinates in the R-wedge. We emphasize that these modes can be extended to the remaining sectors F , L and P via analytic continuation across the event horizons, which represent branch points for these functions.
By using Eq. (4), the Dirac equation in R takes the form This admits as solutions with and with In the above equations we have used the shorthand notation x ≡ (y, z) and k ≡ (k y , k z ) for the transverse positon and momentum coordinates, respectively. The Hamiltonian eigenvalue has been denoted by Ω > 0, while κ ≡ {Ω, k}.
For later purposes, it is also useful to consider the solutions (15) and (17) expressed in Minkowski coordinates. By taking into account the spinor transformations under coordinate changes, these solutions take the form [51] We stress that these functions are still defined in the R-wedge only.
Exploiting the above tools, we now expand the Dirac field in R as whereκ ≡ (−Ω, k). As for the Minkowski quantization, the operators a j,κ and b j,κ are assumed to be canonical, i.e.
with all other anticommutators being trivial. They now annihilate the Rindler vacuum state while a † j,κ and b † j,κ create a particle/antiparticle for the Rindler observer in R. To make it easier the comparison between the field expansions (9) and (20), let us introduce an alternative field quantization in Minkowski spacetime, namely a quantization in terms of the Lorentz momentum eigenfunctions. The physical reason behind this formalism is that the Rindler Hamiltonian is exactly the Lorentz boost generator written in Minkowski coordinates. This means that the time evolution of a Rindler observer can be properly described as an infinite succession of infinitesimal boost transformations (in our case along the x-axis).
We consider the integral representations of the Bessel function K ν (ρ) given by From the last representation, by using the coordinate transformation (4) one can show that [51] K ν (µ k ρ(t, x)) e −νη = where P − θ (t, x) is the two-dimensional positive-frequency planewave with ω = µ k cosh θ and k x = µ k sinh θ. On the other hand, Eq. (23) leads to where P + θ (t, x) is the two-dimensional negative-frequency planewave. Remarkably, Eqs. (25) and (26) show that the function K ν (µ k ρ(t, x)) e −νη can be expressed equivalently as superposition of either positive-or negative-frequency planewaves.
If we now plug the above representations of the Bessel function into Eq. (19), we obtain (27) As opposed to Eq. (19), these functions are globally defined (and analytical) on the entire Minkowski manifold, except for the origin. Furthermore, consistently with our previous considerations, they are eigenfunctions of the boost generator, with eigenvalue Ω. The coefficient N ∓ κ = (2π √ µ k ) −1 e ±πΩ/2 is now fixed by requiring the orthonormality with respect to the ordinary inner product in Minkowski spacetime. We also notice that the existence of two different integral representations of the Bessel function is reflected into the existence of two different global representations of the Rindler modes in Minkowski coordinates 1 . In turn, this is related to the possibility of extending such modes across the event horizons along two different paths (see [51] for more details).
At this stage, we can exploit the fact that the functions ψ + j,κ (ψ − j,κ ) are linear combinations of only positive-(negative-) frequency planewaves to introduce a field expansion equivalent to Eq. (9) as regards the associated Fock space. This expansion reads Since a field excitation in this quantization is still a positive-frequency mode with respect to the time t, it is easy to understand that the vacuum for c k,κ and d k,κ is the same as that for a r k and b r k defined in Eq. (12). Mathematically, this can be seen by the following relations which show that the c j,κ -annihilator is a linear superposition of the a r k 's only. Similarly, for d † j,κ ψ + j,κ we have which show that the d j,κ -annihilator is a linear superposition of the b r k 's only. To simplify the notation, in the above formulas we have omitted the dependence of F j,r and G j,r by the transverse momentum k, since it does not play any significant rôle in our next calculations. Furthermore, the field expansion (9) has been used as right entry for the computation of the Minkowski inner products (29) and (31).
By using Eq. (12), it goes without saying that The proof of the equivalence between the two quantum constructions is completed by noticing that the transformations (29) and (31) are canonical, which implies that the new ladder operators c j,κ and d j,κ still obey the canonical anticommutation relations.

C. Unruh effect
We can now derive the Unruh effect for the Dirac field as originally proposed in [1]. Such computation requires the formulation of a quantum construction which is valid for the Minkowski spacetime and gives back the quantization (20) when restricted to the R-wedge. For this task, we make use of the global modes ψ ∓ j,κ and define the two combinations These functions have all the required properties. Indeed, they are analytical and orthonormal in the whole Minkowski space. Furthermore R j,κ is such that it is defined everywhere but in the L-sector and reduces to Eqs. (15) and (17) in R (depending on whether j = 1, 2). The opposite behavior holds for L j,κ . By inverting the above relations with respect to ψ ∓ j,κ and replacing into Eq. (28), we finally get where the operators r j,κ and l j,κ are connected to c j,κ and d j,κ via the Bogoliubov transformations If we now identify r j,κ with the creation operator for the Rindler observer in R, we can evaluate the spectrum of particles detected by such an observer in the Minkowski vacuum quite straightforwardly, obtaining where in the last step we have introduced the Rindler proper energy E = aΩ and the inverse Unruh temperature β = 1/T U . Clearly, the same result would be obtained by considering the vacuum expectation value of l † j,κ l i,κ , which represents the particle number operator for an accelerated observer in L. Thus, from Eq. (37) we conclude that the Minkowski vacuum appears as a thermal state for the Rindler observer. The spectrum of particle is distributed according to the Fermi-Dirac statistics with temperature proportional to the magnitude of the acceleration. This is the well-known Unruh effect [1].
A comment is in order here: the spectrum (37) diverges for κ = κ . This is due to the fact the states created by the action of r † j,κ and a r, † k on the vacuum are ill-defined, in the sense that they are not properly normalizable. To avoid unphysical divergences, a wavepacket based approach should be performed. This has been done explicitly for bosons in [33,53], where it has been shown that the ensuing spectrum retains the same profile as in the planewave formalism, with the Dirac delta function in energy-momentum being now replaced by the Kronecker delta. However, since wavepackets do not affect the spinorial structure of the field, a similar result is expected to be obtained for fermions as well. Then, in the wavepacket approach we can formally write without any pathological behavior for fixed κ = κ . In the next Section, we shall go on working with planewaves and only at the end we take care of this subtlety.

III. FLAVOR MIXING IN RINDLER SPACETIME
Quantum mixing is among the most challenging and important topics in Particle Physics. In the well-established quantum mechanical version, it consists in the superposition of two states of particles with different masses, rotated each other of an angle θ to give rise to the so-called flavor states. Since its prediction by Pontecorvo [54], the theoretical basis of flavor mixing, and in particular of neutrino mixing, have been largely investigated [55], leading to the discovery of the phenomenon of neutrino oscillations. Experimental developments have later provided convincing evidences supporting Pontecorvo's proposal [56,57]. This has inevitably opened up new scenarios in the physics beyond the Standard Model, which account for the neutrino having a non-vanishing mass.
Two decades later Pontecorvo's quantum mechanical approach, a field theoretical formalism for mixed fields has been developed [31]. This has shown the shortcomings of the original description by pointing out the ortoghonality between the vacuum for fields with definite flavors and that for fields with definite masses. Remarkably, this analysis has revealed a rich structure of the interacting-field vacuum as a su(2) coherent state, which in turn underlies the unitary inequivalence between the flavor and mass Fock spaces and the consequent alteration of the oscillation formula to include the antiparticle degrees of freedom (see Appendix A for more details).
All of the above field theoretical considerations on flavor mixing have been carried out in Minkowski background. Recently, they have been generalized to Rindler metric for the case of boson mixing [29]. To figure out how the flavor/mass inequivalence appears to the Rindler observer, let us consider the mixing transformations for Dirac fields in a simplified two-flavor model where Ψ e , Ψ µ are the interacting fields with definite flavors (electron and muon), while Ψ 1 , Ψ 2 denote the free fields with definite masses m 1 and m 2 (with m 2 > m 1 ), satisfying the canonical anticommutator (13).
By assuming the standard planewave expansion (9) for Ψ 1 , Ψ 2 , Eqs. (39) take the form [31] where (σ) is the flavor (mass) index and we have rewritten the expansion (9) in a formally different (but equivalent) way by embedding the time-dependent exponential factor in the spinors u r k, The flavor annihilators a r k, (θ, t) and b r k, (θ, t) have been computed explicitly in Appendix A. Furthermore, the algebraic nature of Eqs. (39) has been explored in Eq. (A1), showing that each of them exhibits the structure of a rotation nested into a Bogoliubov transformation when considered at the level of ladder operators. As remarked above, the most striking implication of this peculiar structure turns out to be the inequivalence between the mass vacuum |0 1,2 = |0 1 ⊗ |0 2 (such that a r k,i |0 1,2 = b r k,i |0 1,2 = 0, i = 1, 2) and the flavor vacuum |0(θ, t) e,µ (such that a r k, (θ, t)|0(θ, t) e,µ = b r k, (θ, t)|0(θ, t) e,µ = 0, ∀t, = e, µ), which becomes a condensate of entangled massive particle/antiparticle pairs of density given by Eq. (A21). Now, in the previous Section we have seen that a useful tool to quantize the Dirac field for an accelerated observer is the boost-mode expansion (28). Thus, to extend the QFT formalism of mixing to the Rindler metric, we first consider the analogue of Eq. (28) for flavor fields. By using the same approach as that adopted for the derivation of Eq. (40), we infer the following relation where the flavor ladder operators c j,κ, and d j,κ, are related to the flavor operators a r k, and b r k, in the planewave expansion (40) by with F j,r,σ being given in Eq. (30), and with G j,r,σ being given in Eq. (32) (see Eqs. (29) and (31) for the corresponding relations in the mass basis). By relying on Eq. (33) and the related discussion, we stress that Eqs. (42) and (43) are canonical transformations which leave the flavor vacuum |0(θ, t) e,µ unchanged. The expression (41) provides us with the springboard we were looking for. Indeed, we can closely follow the Rindler quantization scheme developed for massive fields at the end of Sec. (II) and apply to flavor fields. In doing so, we arrive at the following expansion where ( , σ) = {(e, 1), (µ, 2)}, as usual. The flavor ladder operators for the Rindler observer read which have indeed the same structure as the transformations (36), the mass index being now replaced by the flavor one. Two comments deserve attention here: first, as opposed to the Minkowski frequencies ω k,i , i = 1, 2, the Rindler frequency Ω does not carry any index since it is independent of the field masses. Furthermore, we emphasize that the flavor ladder operators for the Rindler observer are given by the combinations of two Bogoliubov transformations, the thermal one responsible for the Unruh effect and encoded by the Ω-dependent Bogoliubov coefficients, the second one arising from mixing and implicitly contained in c j,κ, and d j,κ, trough Eqs. (42) and (43). This is the same structure found for the case of boson field mixing in Rindler space (see [29,58]). By using Eqs. (A11) and (A12), we can now explicitly write down r j,κ, as and similarly for l j,κ, . As in Eq. (37), we can now derive the vacuum distribution of mixed particles for the Rindler observer by computing the expectation value 1,2 0|r † j,κ, (θ, t) r j,κ, (θ, t)|0 1,2 . However, this faces us with the computation of nontrivial integrals, which cannot be solved analytically. Despite such technicalities, interesting implications can still be derived by setting t = η = 0 and in the realistic approximation of relativistic neutrinos, i.e. ≡ m 2 2 − m 2 1 /ω k,1 1. Under these assumptions, it is a matter of direct calculations to show that where the U k and V k coefficients have been defined in Eqs. (A13)-(A17). We stress that all the time-dependent quantities in the above equation must be intended as evaluated for t = η = 0. Therefore, the particle spectrum at t = η = 0 becomes 1,2 0|r † j,κ, (θ) r j,κ, (θ)|0 1,2 = 1 e 2πΩ + 1 + sin 2 θ 2 cosh(πΩ) e πΩ N F F j,κ,σ − e −πΩ N GG j,κ,σ + N F G j,κ,σ + c.c. , where the first contribution is the standard Fermi-Dirac distribution, while the terms multiplied by sin 2 θ are the mixing-induced corrections. Here, we have used the relation (A20) for the Bogoliubov coefficients U k and V k . Furthermore, we have implicitly taken into account the advantages of the wavepacket approach, which allows us to get rid of any pathological divergency in the evaluation of the spectrum for fixed momentum κ (see the discussion at the end of Sec. II). The explicit expressions of N F F j,κ,σ , N GG j,κ,σ and N F G j,κ,σ are rather awkward to exhibit. Here we provide their integral expressions For concreteness, we consider the case of spin j = 1 and flavor = e in Eq. (48) (the same considerations hold true for j = 2 and/or = µ). In turn, this amounts to setting the mass index σ equal to one. Then, by replacing Eqs. (30), (32) and (A19) into (49) and (50) and rewriting one can show that On the other hand, from Eq. (51) it follows that where the function I κ (k x ) is explicitly given in Appendix B. By further manipulation, we finally arrive at where we have factorized out the dependence on both |∆m 2 | = |m 2 2 − m 2 1 | and Ω, and H(µ k,1 ) denotes the (dimensionless) k x -integral of the residual function. Although such integral cannot be evaluated analytically, numerical estimations show that it assumes finite values and its real part is positive, at least in the regime of interest (see the discussion below Eq. (B4)).
Within the framework of Boltzmann-Gibbs thermodynamics, the result (48) (or equivalently (56)) is a signature of the broken thermality of Unruh radiation for mixed fields. Such an effect has been discussed in more detail in [29] for the case of boson mixing. Regardless of the specific nature of the fields, it can essentially be traced back to the complex nature of the quantum vacuum for mixed fields. Indeed, whilst in the absence of mixing this state is made up of one type of pairwise correlated particles/antiparticles, for mixed fields it contains pairs of both equal and different types (see the discussion in Appendix A). This gives rise to a nontrivial modification of the Unruh distribution, which however is recovered for θ = 0 and/or m 2 = m 1 , consistently with the vanishing of mixing. We further elaborate on the physical meaning of Eq. (56) in the next Section, where we revisit the above effect in the context of nonextensive Tsallis theory.

IV. Q-GENERALIZED TSALLIS STATISTICS FOR NEUTRINO MIXING
In Boltzmann-Gibbs statistical mechanics, entropy is as a measure of the number of possible microscopic states of a system in thermodynamic equilibrium, consistent with its macroscopic thermodynamic properties. This definition is enclosed in the well-known formula for a set of W discrete microstates, where p i is the probability of the i-th microscopic configuration with the condition W i=1 p i = 1. Clearly, if probabilities are all equal, this takes the form S BG = log W . It is immediate to check that the entropy as defined above is additive, which means that, given two probabilistically independent systems A and B with entropies S BG (A) and S BG (B), the total entropy is S BG (A + B) = S BG (A) + S BG (B). Intuitively, this reflects the fact that, in the absence of interactions, the number of total states of a composite systems is simply given by the sum of the states of its constituent parts.
Although Boltzmann-Gibbs statistics is the correct way of approaching thermodynamic ergodic systems, it has been argued that strong correlations may be sometime responsible for the emergence of exotic features that lye outside the domain of description of the standard theory. The most emblematic example of systems falling within this category are the gravitational systems, as already pointed out by Gibbs at the beginning of the last century. Hence, in such cases a generalization of Boltzmann-Gibbs statistics and the related additive entropy is unavoidable.
Along this line, in [34][35][36][37] it has been shown that a proper generalization of Eq. (57) for complex systems exhibiting long-range interactions and/or spacetime entanglement is given by the following nonadditive entropy where This is the so-called q-Tsallis entropy and the framework built upon it is the nonextensive (or Tsallis) statistical mechanics. We notice that the Boltzmann-Gibbs definition (57) is a special case of S q for q → 1. Moreover, by considering again two independent systems such that p A+B which shows that S q is superadditive or subadditive, depending on whether q < 1 or q > 1. Since the q-entropic index quantifies the departure from Boltzmann-Gibbs, it is typically named nonextensive Tsallis parameter. As remarked above, the definition (58) works quite well for strongly gravitating systems, and in particular for black holes [36], although in recent years encouraging results have been obtained in the study of the cosmic background radiation [43,44], low-dimensional dissipative systems [37], solar neutrinos [45] and polymer chains [46], among others. The q-generalized entropy (58) naturally leads to nontrivial modifications of the standard tools of Boltzmann-Gibbs theory. For instance, in [59] it has been discussed how the removal of the assumption of extensivity of entropy from statistical mechanics affects the core thermodynamic relations, including the zero-th law of thermodynamics and the concepts of temperature and pressure. Similarly, in [60][61][62][63][64] Tsallis entropy has been used to derive the following generalized average occupational numbers which are obtained through the usual procedure of maximizing Tsallis entropy under the constraints of keeping constant the average internal energy and the average number of particles. Here the upper (lower) sign refers to the generalized Fermi-Dirac (Bose-Einstein) distribution 2 . Furthermore β = 1/T , where T is the temperature of the system. We remark that the usual Fermi-Dirac and Bose-Einstein distributions are recovered in the q → 1 limit.
To make physical sense, the distributions (61) must be non-negative and real-valued. This gives rise to the following constraint Following [66], we also point out that Eq. (61) only provides a preliminary but still useful approximation. Indeed, one cannot derive the exact analytical expression of the generalized distribution for arbitrary values of q, due to the nonseparability of the partition function. Nevertheless, the error committed turns out to be fairly negligible at very low temperatures for systems with large numbers of particles, as it is for quantum fields (see [66] for more detailed numerical estimations). Accordingly, since the Unruh temperature is extremely small even for huge accelerations 3 , the use of Eq. (61) is well justified.
In the previous Section we have emphasized that the Unruh spectrum for mixed fields loses its typical thermal profile owing to the entangled structure acquired by the quantum vacuum, which becomes a (time-dependent) su(2) coherent state (see Eqs. (48) or (56)). Based on the observation that the generalized Tsallis framework well describes strongly correlated systems, either on quantum or classical grounds, in [33] the question arose of how such a result could appear when rephrased in the language of nonextensive Tsallis statistics. Remarkably, for mixing of bosons it has been shown that the modified Unruh distribution can be mapped into the q-generalized Bose-Einstein distribution (61), provided that the q-entropic index is properly related to the characteristic mixing parameters sin θ and ∆m.
Following this analysis, here we extend the Tsallis paradigm to mixing of fermions as well. Clearly, since the correction in Eq. (56) slightly affects the Fermi-Dirac spectrum at both high and low energy regimes, it is reasonable to expand the generalized distribution (61) (with the + sign in the denominator) for tiny departures of q from unity.
To the leading order, we then obtain to be compared with the corresponding Tsallis-induced correction for bosons (see Eq. (39) of [33]) In passing, we note that a similar perturbative analysis in nonextensive statistics has been carried out in the context of high-energy heavy-ion collisions [67], starting from the experimental evidence that high-energy nucleusnucleus collisions cannot be described in terms of superpositions of elementary nucleon-nucleon interactions. In that case, the use of Tsallis statistics is motivated by the hypothesis that collective effects in the hadronic medium can be traced back to some kind of memory effects and long-range forces that may occur during high-energy ion collisions. Now, Eq. (63) can be compared to the modified spectrum (56) by setting E = aΩ (the proper energy of Rindler modes) and β = 1/T U = 2π/a. Straightforward calculations lead to Although we are only retaining the first term in the expansion, we notice that the energy dependence of the qgeneralized spectrum gets nontrivially altered.
Let us now investigate the correspondence between the modification of Unruh distribution arising in the context of flavor mixing and that obtained in the nonextensive Tsallis theory, respectively. By equating Eqs. (56) and (65), we infer where in the last step we have approximated apart from an irrelevant factor, and we have defined The function H(µ k,1 ) is defined in Eqs. (B2)-(B4) of Appendix B. Now, since Re{H(µ k,1 )} > 0 (at least in the regime we are considering, see Appendix B), from the above relations we infer that q > 1, indicating the subadditive regime of Tsallis entropy, according to Eq. (60). Furthermore, from Eq. (62) we see that there is no upper limit on the energy scale. By contrast, for the case of boson mixing one has q < 1, which corresponds to a flavor vacuum entropy higher than the sum of the entropies of the two mass vacua taken separately (Tsallis superadditive regime) [33].
Therefore, in the above effective picture, the properties of Unruh condensate for mixed particles can be described in terms of the nonextensive Tsallis statistics, with the q-entropic index satisfying the condition (66). As expected, the departure of q from unity depends on the mixing angle and the mass difference in such a way that the usual Boltzmann-Gibbs theory with q = 1 is recovered for θ = 0 and/or ∆m = 0, consistently with the vanishing of mixing in both cases. The same behavior occurs in the ultrarelativistic limit ∆m/µ k,1 → 0, since Re{H(µ k,1 )} is finite even for large momenta. This can be understood by observing that in such approximation the field theoretical description of flavor mixing reproduces the standard quantum mechanical framework, in which the flavor and mass representations are connected by a simple rotation, and thus the related Fock spaces are unitarily equivalent to each other (see the discussion in Appendix A).
Following [33], we now exploit Eq. (66) to provide an alternative interpretation for the result (56). Indeed, instead of regarding the mixing-induced distortion of the spectrum as a nonthermal effect within Boltzmann-Gibbs framework, we can trace its origin back to a modification of the vacuum distribution of particles/antiparticle pairs at a more fundamental statistical level. We remark that a similar connection between the Unruh effect and Tsallis statistics has been exhibited in [68] in the context of deformed uncertainty relations accounting for a minimal length at Planck scale (Generalized Uncertainty Principle).
As for the case of bosons, the identification in Eq. (66) shows that the q-entropic index exhibits a running behavior as a function of the Rindler frequency Ω. This is quite a common behavior for quantum field theoretical and quantum gravity systems. An analogous study in Tsallis statistics with a varying nonextensive parameter has been recently developed in cosmology in [47]. Interestingly, we notice that, while for mixing of bosons the deviation of q from unity depends on Ω −1 [33], in the present case it scales as Ω −2 . Of course, this does not entail any dimensional problem, being the Rindler frequency Ω dimensionless.
In this connection, one may argue that Eq. (66) exhibits a pathological behavior for arbitrarily small frequencies, as it diverges in the Ω → 0 limit. However, as discussed in [33], this must be interpreted as a warning concerning the domain of validity of our approximation. Indeed, for Ω below a certain threshold Ω min , it happens that the qentropic index might largely deviate from unity, which a posteriori would invalidate the expansion (65). To keep the self-consistency of our analysis, we have to require that |q − 1| 1. In turn, this gives rise to the infrared frequency cutoff which means that the more accurate the approximation of relativistic neutrinos ∆m/µ k,1 1, the higher the number of Rindler modes that fit with the q-generalized Fermi-Dirac distribution (61). Clearly, the entire frequency spectrum is spanned for F θ → 0. To give some numerical estimations, here we notice that, by setting the sampling values θ = π/4, ∆m 2 ∼ 10 −5 eV 2 and k y = k z 10 GeV, we get Ω min 10 −8 . On the other hand, the departure from extensivity becomes increasingly negligible for large Ω. In particular the Boltzmann-Gibbs statistics with q = 1 is recovered for Ω → ∞. As explained in [33], this is somehow consistent with the fact that the higher the energy of the level we are considering, the lower the expected number of particles it can accommodate, with both the Fermi-Dirac and q-generalized distributions approaching zero as Ω → ∞. Therefore, it is reasonable to expect that the difference between the two spectra narrows for relatively large Ω. This is exactly the same behavior found in [33] for mixing of bosons.
The connection established between the modification (56) induced by flavor mixing and the q-generalized Fermi-Dirac distribution (63) can be explained in terms of the nontrivial structure acquired by the vacuum for mixed fields (see Appendix A). Indeed, due to the peculiar features of mixing transformations in quantum field theory, this state turns out to be a condensate of entangled particle/antiparticle pairs of both equal and different types. As a consequence, while the standard Unruh effect may be seen as due to single-type fluctuations spontaneously appearing close to the Rindler horizon (one component falling back into the horizon, the other escaping as Unruh radiation), in the presence of field mixing it can be generated by different types of virtual pairs, resulting in a deformed spectrum of Unruh radiation. This has been pictorially represented in Fig. 2 by associating different (online) colors to different types of particles. We have shown above that such effect can be framed in the context of Tsallis statistics, with a modified distribution given by the q-generalized distribution (63) based on nonadditive Tsallis entropy (58). In particular, the departure q − 1 from extensivity is quantified by the characteristic mixing parameters and the energy scale (see Eq. (66)). As a remark, from Eq. (60) we also emphasize that q > 1 implies that mixed particles are packaged in the Unruh condensate in such a way that the total entropy of the vacuum for mixed fields is lower than the sum of the vacuum entropies for the corresponding free fields. Heuristically speaking, this can be understood as follows: suppose we dispose of a detector of red-type particles only (analogous considerations hold for a detector of blue-type particles). In the absence of mixing, the detection of a particle definitely corresponds to the same type of antiparticle fallen back into the horizon. By contrast, for mixed fields, if we see a red particle, it may be associated to a red-antiparticle, as well as to a blue-antiparticle crossing inward the horizon, with the sum of probabilities being normalized to unity (see Fig. 2). However, a blue-antiparticle that falls into the horizon can also lead to the escape of a blue-particle, which will remain undetected under our hypothesis. Since, for the Pauli exclusion principle, two identical fermions cannot occupy the same quantum state simultaneously, we infer that the existence of blue-antiparticle/blue particle pairs in the vacuum affects the microscopic configuration of the system by reducing the available phase space volume for the blue-antiparticle/red particle pairs. Indeed, it may happen that the blue antiparticle of a blue-blue pair with given quantum numbers crosses inward the horizon before the blue antiparticle of a blue-red pair with the same quantum numbers does. Then, for its part the blue antiparticle of the blue-red pair will be prevented from entering the horizon and its companion from escaping as Unruh radiation. This results in a lower number of possible microstates for the observed system and a related decrease of the total entropy with respect to the case of unmixed fields. Of course, such a behavior does not occur for bosons, for which the condensate structure of vacuum made up For free (i.e unmixed) fields, the vacuum state is made up of single-type pairs (either blue-blue or red-red). By contrast, in the presence of field mixing its structure becomes far richer due to the simultaneous appearance of hybrid pairs (red-blue and blue-red). The ensuing spectrum can still be featured as a thermal-like bath, provided that vacuum constituents are assumed to be distributed according to nonextensive Tsallis statistics.
of hybrid types of pairs, combined with the fact that bosons are not subject to the Pauli exclusion principle, leads to the opposite superadditive regime for Tsallis entropy [33].
It is interesting to observe that a similar comparative analysis of the generalized thermostatistic properties of bosons and fermions has been carried out in [69] in the context of the quantum q-deformed algebra. This differs from Tsallis framework in that it relies on the modification of the quantum algebra of the creation and annihilation operators based on the q-calculus [70], rather than the modification of the partitioning of the microstates of the many body system. In that case, it has been shown that q-deformed bosons and fermions have an enhancement of the quantum statistical effects compared to standard behavior. In addition, in [71] the thermodynamic geometry of an ideal q-deformed boson and fermion gas has been constructed, investigating some properties such as the stability and statistical interaction.
The above considerations allow us to state that the correlations induced by mixing do spoil the macroscopic properties of Unruh condensate by affecting the statistical behavior of its microscopic configurations. In this connection, we notice that nonextensive statistics based on Tsallis entropies, such as the relative entropy and the Peres criterion, have been widely used in the study of entanglement in the last years [72].

V. CONCLUSIONS AND OUTLOOK
The Unruh effect is a notorious prediction built in quantum field theory, which states that empty space appears as a thermal bath of particles from the point of view of an accelerating observer [1]. Recently, it has been shown that this effect is nontrivially spoilt when dealing with mixing of fields with different masses [29,30], leading to a nonthermal distortion of the vacuum condensate. Following the analysis of [33] for boson mixing, here we have revisited this phenomenon in the context of nonextensive Tsallis thermodynamics. By considering the superposition of two Dirac neutrino fields in the relativistic approximation, we have provided an effective description of the modified Unruh spectrum in terms of the q-generalized Fermi-Dirac distribution based on the nonadditive Tsallis entropy. In this framework, the departure from the standard Boltzmann-Gibbs theory turns out to be quantified by the mass difference of the mixed fields and the mixing angle, in such a way that the ordinary extensive mechanical statistics is recovered for vanishing mixing. Furthermore, the condition q > 1 indicates the subadditive feature of Tsallis entropy. In this sense, the present result differs from that of [33], where the opposite superadditive regime has been obtained.
As remarked in [33], the above analysis provides us with a more conservative interpretation for the result of [29,30], since it allows us to extend the usual thermal behavior of Unruh effect to the case of mixed fields. In other terms, one can maintain a thermal-like picture for the vacuum condensate, provided that the distribution of particle/antiparticle pairs is assumed to obey Tsallis (rather than Boltzmann) prescription. In turn, the origin of this statistical effect can be ascribed to the nontrivial structure acquired by the vacuum for mixed fields, which becomes a hybrid condensate of entangled particle-antiparticle pairs. Therefore, in this scenario it is still possible to associate an effective temperature to the vacuum condensate, which would however be different from the standard definition (1), since it would depend on the entropic index (see [59] for more details). Clearly, this aspect deserves careful attention and will be investigated in more detail elsewhere, along with the possibility of having a non-stationary Unruh condensate in the presence of field mixing. In this respect, it should be emphasized that a q-dependent expression for the Unruh temperature has been derived in the context of the Generalized Uncertainty Principle (GUP) in [68]. In that case, the departure of Unruh effect from Boltzmann statistics is induced by the emergence of a minimal length of the order of Planck scale, which affects the phase-space structure by changing the elementary cell volume. In light of this similarity, it would be interesting to investigate whether there exists any kind of connection between the two frameworks. For instance, a suggestive idea could be to rethink the phenomenon of mixing as an effect induced by background geometry, just like the GUP arises from a modified geometry of spacetime at quantum gravity scale.
Apart from these long-term perspectives, the present analysis can be refined by looking at the following issues: to avoid technicalities, in our computations we have considered a simplified model involving only two fields. However, we expect that the generalization to the more realistic three-flavor scenario does not affect our results at the conceptual level. Furthermore, we have made use of the assumption of relativistic neutrinos and, at the same time, we have expanded the q-generalized Fermi-Dirac distribution for small deviations of q from unity. To strengthen our outcome, the exact calculation would be required. Then, the question arises as to how Eq. (66) appears when relaxing any approximation. We also developed computations for t = 0, neglecting the time-dependence of the vacuum for mixed fields. Clearly, a more complete analysis should take into account the time evolution of the vacuum (and any possible deviation from thermal equilibrium) as well. A further direction would be to analyze the above framework by considering neutrinos as Majorana fields. As well-known, indeed, the very nature of neutrinos -Dirac or Majorana -is still unclear, the only practical way to discriminate being the search for neutrinoless double β-decays. Based on [73], where the Unruh effect was shown to be insensitive to the Majorana/Dirac dichotomy, also in this case we envisage no significant deviations from the current result.
Finally, it is easy to understand that experimental tests of our prediction are extremely challenging at the present, the reason being that Unruh radiation itself has not yet been detected. Nevertheless, promising proposals have been put forward in the context of analog models. Among these, we mention experiments on water waves [22,23], graphene [25] and other condensed matter systems [24,26]. Therefore, such models represent the only source to address experimentally the Unruh effect -and any possible deviation from the standard behavior -to date.
Again, the density (A21) vanishes for θ = 0 and/or m 2 = m 1 , thus giving back the expected equality between the flavor and mass representations when there is no mixing. The same occurs in the limit of large momenta with respect to the masses m 1 and m 2 , consistently with the recovery of the standard quantum mechanical picture.
Of course, by exploiting the symmetric structure of Eq. (A10), one can reverse the above reasoning and study the inherent structure of the mass vacuum in terms of the definite flavor quanta. What is obtained in that case is that such state would appear as a condensate of particle/antiparticle pairs of both equal and different flavors. For more details on the QFT formalism of field mixing and the related physical implications, we refer the reader to [31,39,40,74]. Specifically, in [39] it has been shown that the flavor Fock space cannot be obtained by the direct product of the Fock spaces for massive fields. Therefore, the nontrivial nature of QFT mixing is a feature boiling down to the nonfactorizability of the flavor states in terms of those with definite massses, including the vacuum. The entanglement content of the flavor vacuum has been quantified in [40] for the case of boson mixing and in the limit of relatively small difference of masses.
where we have separated out the dependence on both the mass difference |∆m 2 | and the Rindler energy Ω. H(µ k,1 ) is a shorthand notation for the (dimensionless) k x -integral of the residual function 5 .
The integral (B2) can be estimated numerically. In the limit of relativistic neutrinos and for values of Ω such that the unmodified Fermi-Dirac spectrum in Eq. (48) differs appreciably from zero, N F G j,κ,σ is finite and its real part assumes positive values. As Ω increases, the rapidly oscillating behavior of H(µ k,1 ) makes it difficult even numerical evaluations. However, this regime turns out to be less relevant, since the larger the energy Ω, the lower the expected number of particles in the vacuum condensate. Further investigation of these technicalities is under active consideration and will be carried out elsewhere.