Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Statistical Analyses of a Class of Random Cyclooctatetraene Chain Networks with Respect to Several Topological Properties
Next Article in Special Issue
Analogies between Lattice QCD and the Truncated Nambu–Jona-Lasinio Model
Previous Article in Journal
Symmetric Adaptive Higher-Order Energy-Preserving Methods for a Charged Particle System and Guiding Center System
Previous Article in Special Issue
Dynamical Chiral Symmetry Breaking in Quantum Chromo Dynamics: Delicate and Intricate
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Neutrino Oscillations in Finite Time Path Out-of-Equilibrium Thermal Field Theory

1
Rudjer Bošković Institute, Bijenička cesta 54, 10000 Zagreb, Croatia
2
Physics Department, Faculty of Science, University of Zagreb, Bijenička cesta 32, 10000 Zagreb, Croatia
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(11), 1970; https://doi.org/10.3390/sym15111970
Submission received: 3 July 2023 / Revised: 11 September 2023 / Accepted: 21 September 2023 / Published: 24 October 2023

Abstract

:
We demonstrate that the Finite-Time-Path Field Theory is an adequate tool for calculating neutrino oscillations. We apply this theory using a mass-mixing Lagrangian which involves the correct Dirac spin and chirality structure and a Pontecorvo–Maki–Nakagawa–Sakata (PMNS)-like mixing matrix. The model is exactly solvable. The Dyson–Schwinger equations transform propagators of the input free (massless) flavor neutrinos into a linear combination of oscillating (massive) neutrinos. The results are consistent with the predictions of the PMNS matrix while allowing for extrapolation to early times.

1.  Introduction

The evidence for neutrino oscillations, which solve the solar neutrino problem [1], started with the study of atmospheric oscillations [2]. Evidence has been collected from various sources and for various distances from the sources over a wide range of neutrino energies and detectors [3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43].
Neutrino oscillations, as well as similar problems such as kaon oscillations and decays, suffer from the absence of a proper formalism for the features in which a finite-time description is essential.
The most important aspect of this paper is to demonstrate that the Finite Time Path Field Theory (FTPFT) is an adequate tool [44,45] to treat such problems. In particular, we start with the free Lagrangian L 0 of massless neutrinos with three flavors; however, we mix them dynamically through an interaction Lagrangian L I , as in addition to the standard weak interaction term L W it has the term L M i x , which contains a Pontecorvo–Maki–Nakagawa–Sakata (PMNS)-like matrix involving neutrino masses. The term L W is of course crucial for the processes of creation and detection of flavor neutrinos. Nevertheless, we adopt the approximation where the neutrino mixing is fully due to the term L M i x , that is, we neglect the residual influence of L W on the neutrino masses and mixing. Of course, we must keep in mind the caveat that not calculating these corrections misses the contributions from the vacuum condensate, which Blasone et al. [46,47,48] found to be important corrections to the PMNS result. Nevertheless, they vanish in the relativistic limit, while the PMNS result is recovered [46,47,48]. Because the case of nonrelativistic neutrino kinematics seems to be far from present experimental capabilities, we proceed to study neutrino propagators in our aforedescribed framework, which provides us with an exactly solvable model suitable for studying how FTPFT can be applied to neutrino oscillations.
Through the exact solution of the Dyson–Schwinger equations (DSE), we then obtain propagators (with the retarded S R and advanced S A components) of oscillating (massive) neutrinos which couple to the weak interactions. The DSE equation for the Keldysh component ( S K ) of neutrino propagators leads to two possible histories distinguished by two of three masses (say, m i and m j ) for which the interference leads to oscillations. An equal time limit of S K provides the number of neutrinos at time t.
The origin of neutrino masses has not yet been answered conclusively. The variety of extant approaches includes the Majorana mass term, right-handed neutrinos with a see-saw mechanism, supersymmetric extensions of Standard model, and more. One intriguing possibility is the generation of neutrino masses and mixing via the aforementioned flavor vacuum condensate [46,47,48]. However, in the present paper we do not address the issue of the origin of neutrino masses, and focus on demonstrating the calculation of neutrino oscillations in the FTPFT framework.
Whatever the mechanism of neutrino mass generation, in this paper we assume that it manifests itself through an additional term in the Lagrangian: the mass mixing term L M i x . Here, the interaction Lagrangian contains the mass mixing matrix in a close analogy to the PMNS matrix except with an explicit Dirac spin and chirality structure. The model is exactly solvable.
To calculate neutrino oscillations, we use FTPFT; [49,50] in addition, see [51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68], often called the closed time path (CTP) thermal field theory (TFT). We find this to be an appropriate tool for calculating a wide variety of problems, including equilibrium and nonequilibrium TFT problems, for which it was originally constructed [49,50], as well as scattering processes, decays, and oscillations.
How can field theory describe oscillations? It depends on which version of field theory is employed (for example, see [69,70] and references therein). If we look at the predominant field theory approach, namely, the S-matrix approach, it exhibits excellent properties: adiabatic switching of interaction, Gell-Mann–Low (G-L) theorem, Fermi’s golden rule, reduction technique, renormalization, etc. However, for time-dependent processes, notably oscillations, it has nothing to say, in the sense that it includes adiabatic switching (on and off) of the interaction and acts between the infinite initial time t i = and final time t f = .
FTPFT, on the other hand, uses a finite time path, no adiabatic switching, and no G-L theorem. The particle states are not exactly the eigenstates of the full Hamiltonian. Renormalization, as we have indicated in [67,68], is somewhat more involved; nevertheless, it can be performed in finite time with essentially the same technique of dimensional regularization.
The drawback of this theory is the lack of a Gell-Mann–Low-type theorem. Nevertheless, FTPFT does not contain disconnected subdiagrams owing to the absence of the maximal-time vertices [61,62,64]. Thus, this drawback may be cured after carrying out all of the calculations in a proper way by considering the large time behavior of the contributions. Then, in this limit one can expect the G-L theorem to again be valid.
The particle number is obtained as the equal time limit (ETL) of the re-summed S K propagator (or better) as the average of ETL (AETL) for the contributions distinguished by t 1 f t 2 f (see Appendix A, especially Appendix A.7). These processes contribute inclusively to the particle number. Fermi’s golden rule need not be put in by hand, as it emerges naturally during calculation of the convolution products (∗-products).
In certain applications, the solution of DSE can be written in a closed form. The solution of DSE is illustrative: conserving energy and momentum along the chain, the chain of retarded functions extends from the final time t towards the earliest time t m i n 0 . There, it meets a S K propagator containing the initial distribution function. At that time, there is a (single) convolution product, which does not conserve energy. Instead, it creates just enough oscillating deviations from the energy-conserving value to satisfy the uncertainty relations. From the lowest time, another chain of advanced functions progresses back to the final time. It conserves energy and momentum as well; however, the energy of the advanced chain differs from the energy of the retarded chain. With the help of this energy difference, the convolution product creates oscillations. It is easy to see that with increasing time the uncertainty in energy becomes smaller, which is in the full agreement with the uncertainty relations.
In certain special cases, energy is conserved at the remaining ∗-product. For instance, this is the case when Fermi’s golden rule emerges.
The remainder of this paper is organized as follows: first, we provide a short overview of the PMNS-approach; then, we define the dynamics of our approach through the mass mixing interaction Lagrangian L M i x (9).
Our calculation is based on the DSE (12)–(15) and their closed form solution. This solution defines the oscillating neutrino which possesses flavor as massive and weakly interacting.
The equal time limit of the S K propagator is directly connected to the neutrino number, as it is measured at time t. In the process of calculation, we show how the number of complicated ∗-products appearing in the formal solution simplifies to a single ∗-product, which is necessary for the time evolution and energy–time uncertainty relations.
As a result of this calculation, we obtain the neutrino number at time t. Projection to the flavor degree of freedom exhibits oscillating behavior. This result is consistent with Heisenberg’s uncertainty relations between time and energy. We reproduce the PMNS result exactly within the limit of ultrarelativistic neutrinos and large time.
Finally, we discuss prospects for further research.
Further details and calculation can be found in Appendix A.

2. PMNS Theory of Neutrino Oscillation

In the PMNS theory of neutrino oscillation [71,72,73,74,75], three flavor neutrino states that interact weakly are mixed to three different superpositions of the neutrino states of definite mass (for mixing in gauge theories, see [6]). In the flavor states, neutrinos are emitted and absorbed through weak processes; however, they travel as mass eigenstates.
This is mathematically expressed as
ν α = i U α i ν i , ν i = α U α i ν α ,
where α = e , μ , τ are the flavor indices, respectively, e-electron neutrinos, μ -muon neutrinos, and τ -tauon neutrinos, that label neutrinos with definite flavors, while m i , i = 1 , 2 , 3 , are the indices of the neutrino mass states. To describe antineutrinos, it is necessary to use complex-conjugate matrices U α i U α i . The matrix U α i is the PMNS-matrix, introduced in a way analogous to the CKM matrix describing the mixing of quarks.
The size of the matrix depends on theory. For the standard three-neutrino theory, the matrix is 3 × 3 . For more neutrinos (including sterile ones), the matrix could be larger. There are anomalies [4,5,76,77,78,79,80,81,82,83,84,85,86], suggesting that the model requires further refinement.
The 3 × 3 PMNS matrix (2), for example, is provided in [22]:
U = U e 1 U e 2 U e 3 U μ 1 U μ 2 U μ 3 U τ 1 U τ 2 U τ 3 = 1 0 0 0 c 23 s 23 0 s 23 c 23 c 13 0 s 13 e i δ 0 1 0 s 13 e i δ 0 c 13 c 12 s 12 0 s 12 c 12 0 0 0 1 e i α 1 / 2 0 0 0 e i α 2 / 2 0 0 0 1 ,
where c i j = cos θ i j and s i j = sin θ i j . The phase factors α 1 and α 2 are relevant only if the neutrinos are Majorana particles, i.e., if the neutrino is identical to its antineutrino, otherwise they can be ignored. The phase factor δ measures the degree of violation of CP symmetry, which has not yet been observed experimentally.
Note that Relation (1) is valid only approximately, as particles with the same energy and momentum but different masses cannot simultaneously be on the mass shell. The same feature appears in the usual treatment of neutrino oscillations in field theory, where it is known as mass-shell approximation [69,70,87,88].
A neutrino travels through space as a massive neutrino. Its time evolution is described (in units c = 1 , = 1 ) by ν i ( t ) = e i E i t p i · x ν i ( 0 ) for each mass- and energy-eigenstate ν i with mass m i and eigenenergy E i = p i 2 + m i 2 . However, weak processes actually produce flavor states ν α , which are, in accordance with Equation (1), superpositions of these mass- and energy-eigenstates ν i . Due to their energy differences E i E j , such superpositions have an energy uncertainty Δ E (5).
In the quantum mechanical treatment of neutrino oscillations, it is standard to assume (for example, see [89,90]) that all mass eigenstates ν i have the same momentum, i.e., their energies E i differ by their masses m i and the flavor neutrino state ν α is determined by the momentum.
In the ultrarelativistic limit, p i = p i m i . Thus, p i E , the neutrino energy in the limit m i 0 , such that i :
E i p i + m i 2 2 p i E + m i 2 2 E , and t L ,
where t is the time from the beginning of evolution and L is the distance traveled.
In the process of measurement, the neutrino is projected back to the flavor states. The probability that the initial neutrino with flavor α will be detected later as having flavor β is defined as
P α β = ν β ( L ) ν α 2 = i U α i U β i e i m i 2 L 2 E 2 .
This can be written as
P α β = δ α β 4 i > j Re U α i U β i U α j U β j sin 2 Δ m i j 2 L 4 E + 2 i > j Im U α i U β i U α j U β j sin Δ m i j 2 L 2 E ,
where Δ m i j 2 m i 2 m j 2 .
The second term is related to CP asymmetry:
A C P ( α β ) = P ( ν α ν β ) P ( ν ¯ α ν ¯ β ) = 4 i > j I m U α i U β i U α j U β j sin Δ m i j 2 L 2 E .
With the help of the Jarlskog invariant I m U α i U β i U α j U β j = J γ , k ε α β γ ε i j k , the CP asymmetry is expressed as
A C P ( α β ) = 16 J γ ε α β γ sin Δ m 21 2 L 4 E sin Δ m 32 2 L 4 E sin Δ m 31 2 L 4 E .
Note that CP asymmetry in neutrino oscillations has not yet been observed.
Because
Δ E max i , j = 1 , 2 , 3 ( E i E j ) m i 2 m j 2 2 E ,
Equation (4) satisfies the time–energy uncertainty relation Δ E Δ t 1 of Heisenberg (and of Mandelstam and Tamm, as discussed in [90]) for
t 2 E Δ m i j 2 .
For times shorter than (6), the uncertainty relations are fulfilled with the help of neutrino production processes. Specifically, they produce the flavor states ν α , which are superpositions (1) of ν i , the eigenstates of neutrino masses m i , which are so close that they cannot be resolved in a time shorter than (6) [90].
The ultrarelativistic limit applies to all currently observed neutrinos, as it is known that the differences of their squared masses are of the order 10 4 eV 2 and their energies are at least 1 MeV. Measured oscillation distances L are on the order of kilometers.

3. Neutrino Oscillation as a Dynamical Process

3.1. The Mass Mixing Term in the Interaction Lagrangian

In order to treat neutrino oscillation as a dynamical process, we start with the neutrino mixing interaction L M i x of the Dirac spinor Lagrangian, defined by:
L 0 ( x ) = α ν ¯ α ( x ) i / ν α ( x )
L I = L W + L M i x
L M i x ( x ) = α , i ν ¯ α ( x ) U α , i M i j U β , j ν β ( x ) + a n t i n e u t r i n o s
M = m 1 0 0 0 m 2 0 0 0 m 3 ,
where U α , i is a 3 × 3 matrix analogous to the PMNS-matrix. Here, ν α are neutrino spinor wave functions for flavor α and m i represents the propagating neutrino masses. The initial conditions are fixed through the neutrino distribution function (including either the chirality or a handedness projection operator) which are built in the S f , K component of the flavor neutrino propagator. The massless chiral “flavor neutrinos” (with propagators S f ( p ) and weak interaction) enter the DSE. The self-energy is identified from the next-to-lowest order to the DSE through the use of (9) and (10), and is simply the constant matrix
Σ α , β = Σ α , β , R = Σ α , β , A = U α , i M i j U β , j , Σ α , β , K = 0
The DSE for fermions and their formal solutions are provided in Appendix A, specifically, (A15)–(A19).
The DSEs for oscillating neutrinos are
S ˜ β , α , R = S β , R δ β , α + i S β , R Σ β , η , R S ˜ η , α , R
S ˜ β , α , A = S β , A δ β , α + i S β , η , R Σ η , ζ , A S ˜ ζ , α , A
S ˜ β , α , K = S β , K δ β , α + i [ S β , η , R Σ η , ζ , R S ˜ ζ , α , K
+ S β , η , K Σ η , ζ , A S ˜ ζ , α , A + S β , η , R Σ η , ζ , K S ˜ ζ , α , A ] ,
where S denotes the lowest-order Green function and the same symbol with a tilde ( S ˜ ) denotes the corresponding re-summed Green function, as explained in Appendix A.6.

3.2. Solution of the Dyson–Schwinger Equations for Oscillating Neutrinos

The formal solution for (12)–(15) is
S ˜ β , α , R = [ 1 i S R Σ R ] β , η 1 S η , α , R = S β , η , R [ 1 i Σ R S R ] η , α 1 , S ˜ β , α , A = [ 1 i S A Σ A ] β , η 1 S η , α , A = S β , η , A [ 1 i Σ A S A ] η , α 1 , S ˜ β , α , K = S β , η , K , A ( 1 i Σ A S A ) η , α 1 + ( 1 i S R Σ R ) β , η 1 S η , α , K , R + i ( 1 i S R Σ R ) β , η 1 [ S η , ζ , R Σ ζ , ω , K S ω , θ , A S η , ζ , R Σ ζ , ω , R S ω , θ , K , A + S η , ζ , K , R Σ ζ , ω , A S ω , θ , A ] ( 1 i Σ A S A ) θ , α 1 .
The solution simplifies owing to the following three reasons:
  • The self-energy is a simple matrix, not a retarded or advanced function;
  • The ∗-products among the bare propagators turn to algebraic products except for the case where one factor is retarded (R or K,R) and the other is advanced (A or K,A);
  • The matrix U is unitary.
Next, we calculate ( 1 i Σ R ( A ) S R ( A ) ) θ , α 1 ; the matrix inversion is simple, and we obtain
( 1 i Σ S R ) β , α 1 ( p ) = i U β , i p 2 + m i p p 2 m i 2 + i p 0 ϵ U α , i .
Then, we obtain the re-summed retarded propagator component:
S ˜ β , α , R ( p ) = ( 1 i Σ S R ) β , η 1 S η , α , R ( p ) = i U β , i i ( m i + p ) p 2 m i 2 + i p 0 ϵ U α , i
and the resummed advanced propagator component
S ˜ β , α , A ( p ) = ( 1 i Σ S A ) β , η 1 S η , α , A ( p ) = i U β , i i ( m i + p ) p 2 m i 2 i p 0 ϵ U α , i .
We now have an important conclusion. Relations (18) and (19) express the re-summed flavor propagator through the linear combination of “propagating” neutrino propagators. This should be contrasted to the “similar” (mathematically ill-defined) relations between wave functions of flavor and propagating neutrinos.
Now, we turn to the S ˜ β , α , K propagator. Insertion of (17)–(19) into (16) provides us with
S ˜ β , α , K = δ f , α i U β , i U α , i i p 2 m i 2 i p 0 ϵ × { [ 1 2 n f ( p ) ] p 0 + | p | 2 | p | ( γ 0 | p | γ p ) [ 1 2 n f ¯ ( p ) ] p 0 | p | 2 | p | ( γ 0 | p | + γ p ) } 1 γ 5 2 + δ β , f i U β , i U α , i i p 2 m i 2 + i p 0 ϵ × { [ 1 2 n f ( p ) ] p 0 + | p | 2 | p | ( γ 0 | p | γ p ) [ 1 2 n f ¯ ( p ) ] p 0 | p | 2 | p | ( γ 0 | p | + γ p ) } 1 γ 5 2 + i δ f , α i , j U β , i p 2 + m i p ( p 2 m i 2 + i p 0 ϵ ) ( p 2 + i p 0 ϵ ) U α , i [ i m i ( i p ) { [ 1 2 n f ( p ) ] p 0 + | p | 2 | p | ( γ 0 | p | γ p ) [ 1 2 n f ¯ ( p ) ] p 0 | p | 2 | p | ( γ 0 | p | + γ p ) } 1 γ 5 2 + { [ 1 2 n f ( p ) ] p 0 + | p | 2 | p | ( γ 0 | p | γ p ) [ 1 2 n f ¯ ( p ) ] p 0 | p | 2 | p | ( γ 0 | p | + γ p ) } 1 γ 5 2 ] i m j ( i p ) ] × U α , j p 2 + m j p ( p 2 m j 2 i p 0 ϵ ) ( p 2 i p 0 ϵ ) U β , j ,
where f is the flavor of the initial neutrino beam.
The resummed Keldysh component S ˜ β , α , K of the neutrino propagator consists of two contributions: (a) terms without ∗-products (which we call algebraic) and (b) the terms with a single ∗-product (which we call convolutional).
This propagator carries the information which, after the equal time limit, provides the number and momentum distribution of all types of flavor neutrinos measured at time t.

3.3. The ∗-Products and the Average of Equal Time Limits

The remaining ∗-product in (20) between the retarded and advanced functions is expressed (see Equation (A2) for details) as follows:
C X 0 ( p 0 , p ) = d p 01 d p 02 P X 0 ( p 0 , p 01 + p 02 2 ) i 2 π e i X 0 ( p 01 p 02 + i ϵ ) 1 p 01 p 02 + i ϵ A , R ( p 01 , p ) B , A ( p 02 , p ) , P X 0 ( p 0 , p 0 ) = 1 π Θ ( X 0 ) sin 2 X 0 ( p 0 p 0 ) p 0 p 0 ,
where the label means that the values of A and B should be taken at infinite time. After time t, the number of oscillating neutrinos is expressed through the average equal-time limit of the re-summed Keldysh component S ˜ β , α , K of the neutrino propagator (see Appendix A):
1 N β , p ( t ) = 1 2 π [ lim 0 < Δ 0 + lim 0 > Δ 0 ] d p 0 e i p 0 Δ T r [ γ 0 4 S ˜ β , β , K , t ( p ) 1 γ 5 2 ] ,
where Δ = s 01 s 02 and X 0 = ( s 01 + s 02 ) / 2 = t . In this expression, we have taken into account that the initial condition (A13) contains only flavor neutrinos of type f, not antineutrinos, and we calculate number of flavor neutrinos of type β . The equal time limit removes projector P from the above ∗-product:
1 π lim Δ 0 d p 0 e i Δ p 0 C s 02 = t ( p 0 , p ) = d p 01 d p 02 i 2 π e i t ( p 01 p 02 ) 1 p 01 p 02 A , R ( p 01 , p ) B , A ( p 02 , p ) .

3.4. Contributions to Neutrino Oscillation

To calculate neutrino oscillation, we start with chiral fermion number (22) and insert the solution of the Dyson–Schwinger Equation (20).
Upon managing all ∗-products appearing in (20), we end up with two types of contributions:
  • Contributions without any ∗-product. These contributions (the first pair in (20)) are independent in time. They should reproduce the initial (input) density of neutrinos. By calculating geometrical series term by term, one would obtain the lowest order providing the input neutrino density. All the higher term would vanish, as they are equal time limit of the product of two or more retarded functions (or two or more advanced functions). The re-summed propagator is nonperturbative, and consequently the result renormalizes the initial density. Nevertheless, we obtain the input neutrino density.
  • Contributions containing S α , K (the second pair in (20)). These refer to the initial input of flavor neutrinos of type α . They contribute to oscillation.

3.5. The Algebraic Term

Note here that we have to decide whether we approach equal time from above or from below. The final result does not depend on our choice; however, for our choice Δ = s 01 s 02 > 0 the advanced contribution should vanish.
We calculate the contribution of [ 1 i S R Σ R ] 1 S α , K , R (at n f ¯ ( p ) = 0
I α , a l g = 1 N α , a l g ( p ) ( t ) = i π lim s 01 s 02 = t d p 0 e i p 0 ( s 01 s 02 ) 1 4 T r γ 0 δ α , f i U α , i U α , i i p 2 m i 2 + i p 0 ϵ [ 1 2 n f ( p ) ] p 0 + | p | 2 | p | ( γ 0 | p | γ p ) 1 γ 5 2 1 4 T r γ 0 ( θ + η p ) ( γ 0 | p | + γ p ) [ 1 γ 5 ] = θ | p | I α , a l g = lim s 01 s 02 = t d p 0 e i p 0 ( s 01 s 02 ) i | U α , i | 2 i p 2 m i 2 + i p 0 ϵ [ 1 2 n f ( p ) ] p 0 + | p | 2 π
Finally, we find
I α , a l g = δ α , f [ 1 n f ( p ) ] .
This contribution is the only time-independent contribution to neutrino yield. As we will see, it confirms the conservation of the total number of neutrinos in our model. It indicates that the eventual (finite) wave function renormalization is not necessary.

3.6. ∗-Product, Term Containing S K , R

In the terms with a convolution product, the chirality projector ( 1 γ 5 2 ) appears twice: once to select chiral neutrinos in initial distribution function, and a second time to be selected by the weak interaction measurement device:
I S c o n v , S K , R = 1 4 T r γ 0 i 2 π 2 d p 01 d p 02 e i ( p 01 p 02 ) t 1 p 01 p 02 i [ 1 i Σ S R ] 1 S K , R } ( p 01 ) { Σ S A [ 1 i Σ S A ] 1 } ( p 02 ) [ 1 γ 5 ] / 2 = δ α , f i , j U β , i U α , i U α , j U β , j T r γ 0 i 16 π 2 d p 01 d p 02 e i ( p 01 p 02 ) t 1 p 01 p 02 i p 1 2 + m i p 1 ( p 1 2 m i 2 + i p 01 ϵ ) ( p 1 2 + i p 01 ϵ ) { [ 1 2 n f ( p ) ] p 01 + | p | 2 | p | ( γ 0 | p | γ p ) [ 1 2 n f ¯ ( p ) ] p 01 | p | 2 | p | ( γ 0 | p | + γ p ) } [ 1 γ 5 ] m j ( i p 2 ) ] × p 2 2 + m j p 2 ( p 2 2 m j 2 i p 02 ϵ ) ( p 2 2 i p 02 ϵ ) 1 γ 5 2 p 1 = ( p 01 , p ) , p 2 = ( p 02 , p ) .
The trace is
1 4 T r γ 0 ( p 1 2 + m i p 1 ) [ γ 0 | p | γ p ] [ 1 γ 5 ] ( m j + p 2 ) [ 1 γ 5 ] = 2 m j p 1 2 | p | .
Now, we have
I S c o n v , S K , R = i δ α , f i , j U β , i U α , i U α , j U β , j i m j 2 4 π 2 d p 01 d p 02 e i ( p 01 p 02 ) t 1 p 01 p 02 { [ 1 2 n f ( p ) ] ( p 01 + | p | ) [ 1 2 n f ¯ ( p ) ] ( p 01 | p | ) } ( p 1 2 m i 2 + i p 01 ϵ ) ( p 2 2 m j 2 i p 02 ϵ ) .
The integrals are performed by closing the integration path from below for d p 01 and from above for d p 02 .
d p 01 e i ( p 01 p 02 ) t 1 p 01 p 02 p 01 ± | p | p 1 2 m i 2 + i ϵ p 01 = i π λ i = ± 1 ( e i ( λ i ω i p 02 ) t 1 ) ( ω i ± λ i | p | ) ω i ( λ i ω i p 02 ) d p 02 e i ( λ i ω i p 02 ) t 1 ( λ i ω i p 02 ) ( p 2 2 m j 2 i p 02 ϵ ) = i π λ j = ± 1 λ j e i ( λ i ω i λ j ω j ) t 1 ω j ( λ i ω i λ j ω j )
Here, ω i = [ | p | 2 + m i 2 ] 1 / 2 and ω j = [ | p | 2 + m j 2 ] 1 / 2
Thus, the contribution consists of contributions from all singularities (four terms):
I S c o n v , S K , R = m j 2 4 δ α , f i j U β , i U α , i U α , j U β , j λ i , λ j = ± 1 λ j e i ( λ i ω i λ j ω j ) t 1 ω i ω j ( λ i ω i λ j ω j ) [ λ i | p | n f ( p ) ( ω i + λ i | p | ) + n f ¯ ( p ) ] ( ω i λ i | p | ) ] .
After setting n α ( p ) = 0 , the integral is
I S c o n v , S K , R = m j 2 4 δ α , f i j U β , i U α , i U α , j U β , j λ i , λ j = ± 1 λ j e i ( λ i ω i λ j ω j ) t 1 ω i ω j ( λ i ω i λ j ω j ) [ λ i | p | n f ( p ) ( ω i + λ i | p | ) ω i = [ | p | 2 + m i 2 ] 1 / 2 , ω j = [ | p | 2 + m j 2 ] 1 / 2 .

3.7. ∗-Product, Term Containing S K , A and i j Terms

This integral is similar to the previous one:
I S c o n v , S K , A = m i 2 4 δ α , f i j U β , i U α , i U α , j U β , j λ i , λ j = ± 1 λ i e i ( λ i ω i λ j ω j ) t 1 ω i ω j ( λ i ω i λ j ω j ) [ λ j | p | n f ( p ) ( ω j + λ j | p | ) ] .
With i and j interchanged, the contributions are
I S c o n v , S K , R j i = m i 2 4 δ α , f i j U β , j U α , j U α , i U β , i λ i , λ j = ± 1 λ i e i ( λ j ω j λ i ω i ) t 1 ω i ω j ( λ j ω j λ i ω i ) [ λ j | p | n f ( p ) ( ω k i + λ i | p | ) ]
and
I S c o n v , S K , A j i = m j 2 4 δ α , f i j U β , j U α , j U α , i U β , i λ i , λ j = ± 1 λ j e i ( λ j ω j λ i ω i ) t 1 ω i ω j ( λ j ω j λ i ω i ) [ λ i | p | n f ( p ) ( ω k i + λ i | p | ) ] .

3.8. The Dominant Contribution

The results of the previous section are consistent with Heisenberg’s uncertainty condition between energy and time. In our further calculations, we assume that the time is large compared to the energy differences | ω i ω j | . Together with the fact that the measured neutrinos ( m i < < | p | ) are mostly ultrarelativistic, with | p | ranging from 1 MeV to 10 MeV, the following approximations are justified:
ω i ω j m i 2 m j 2 2 | p | , ω i + ω j 2 | p | , ω i + | p | 2 | p | and ω i | p | m i 2 2 | p | .
Owing to the above, all contributions proportional to [ 1 2 n f ( p ) ] are dominated by λ i = λ j = + 1 , while the contributions proportional to [ 1 2 n α , ( p ) ] are dominated by λ i = λ j = 1 . The contribution from the constant 1 in [ 1 2 n f ] is killed by this procedure.
Notice, however, that for short times ( t h 4 π | ω i ω j | , i j ) the approximations are not justified, and it is necessary to deal with the full expression. in order to satisfy Heisenberg’s uncertainty relations among time and energy.
The above contributions then become
I S c o n v , S K , R d = m j 2 δ α , f i j U β , i U α , i U α , j U β , j e i m i 2 m j 2 2 | p | t 1 m i 2 m j 2 n f ( p )
I S c o n v , S K , A d = m i 2 δ α , f i j U β , i U α , i U α , j U β , j e i m i 2 m j 2 2 | p | t 1 m i 2 m j 2 n f ( p )
I S c o n v , S K , R j i d = m i 2 δ α , f i j U β , j U α , j U α , i U β , i e i m i 2 m j 2 2 | p | t 1 m i 2 m j 2 n f ( p )
and
I S c o n v , S K , A j i d = m j 2 δ α , f i j U β , j U α , j U α , i U β , i e i m i 2 m j 2 2 | p | t 1 m i 2 m j 2 n f ( p ) .
By adding the contributions, we obtain
I S c o n v , S K = I S c o n v , S K , R d I S c o n v , S K , A d + I S c o n v , S K , R j i d I S c o n v , S K , A j i d = δ α , f i j { U β , i U α , i U α , j U β , j [ e i m i 2 m j 2 2 | p | t 1 ] [ U β , j U α , j U α , i U β , i e i m i 2 m j 2 2 | p | t 1 ] } n f ( p ) ,
which can be further written as
I S c o n v , S K = n f ( p ) i j δ α , f n f ( p ) [ 4 R e ( U β , i U α , i U α , j U β , j ) sin 2 m i 2 m j 2 4 | p | t + 2 I m ( U β , j U α , j U α , i U β , i ) sin m i 2 m j 2 2 | p | t ] .

4. Final Result

From the preceding section, notably Equation (22), and from Appendix A, notably Appendix A.8, it can be seen that at time t the total number of particles of the flavor β stemming from the initial flavor α is
N β , f , p ( t ) = δ α , f n α , + ( p ) + n f ( p ) δ α , f i j [ 4 R e ( U β , i U α , i U α , j U β , j ) sin 2 m k i 2 m l j 2 4 | p | t + 2 I m ( U β , j U α , j U α , i U β , i ) sin m i 2 m j 2 2 | p | t ] .
Let us point out a few features of this expression:
  • The result (42) is identical to the standard PMNS expression (4). The ultrarelativistic relation (3) reveals the equality of the arguments of the sines, while division by the initial distribution of the number of particles n f recasts (42) in terms of probability, as in (4). Thus, with the same presently available inputs, our result (42) would provide the same numerical results as, for example, [86].
  • If we sum over β , the oscillating contribution vanishes. This reflects the fact that the total neutrino number is conserved within the realm of chiral neutrinos; notably, the sterile neutrinos are not involved! Notice that our conclusion is valid for low energy neutrino beams as well; this is easily verified by looking at (31)–(34).
  • The results are valid for moderate energies, although it is necessary to skip the simplifications in (3) and (35) and include all contributions.

Conclusions

In this paper, we apply the Finite Time Path Field Theory (FTPFT), originally designed to deal with out-of-equilibrium many-body statistical ensembles, to the problem in particle physics. We demonstrate that FTPFT is an appropriate tool for the treatment of neutrino oscillations.
We calculate neutrino oscillation within a simple model, with the interaction Lagrangian containing the term L M i x built as mass mixing through the PMNS-matrix with a built-in Dirac spinor and chirality structure. The model is exactly solvable. This is an extension of standard PMNS-case, as it involves sterile neutrinos, at least those with the same mass as flavor and propagating neutrinos. Of course, more sophisticated sterile neutrinos would require additional model building.
The result is consistent with Heisenberg’s uncertainty relation between energy and time. The flavor neutrinos are chiral spin 1 / 2 particles, though a massive neutrino propagator should additionally contain the right-handed component along with the re-summed oscillating propagator. Nevertheless, the result is chirally invariant, as it is dictated from the weak interaction Lagrangian ( L W ) taken into account through the factors ( 1 γ 5 ) appearing at the beginning (production) and end (detection) stages. The sterile neutrinos do not contribute to oscillations, even for low energy beams. Our result coincides with the PMNS-formula for large times and for the ultrarelativistic case.
In the derivation of our result, the PMNS relation (1) is not used. Instead of using the relation which mixes flavor and propagating neutrino states (and in which at least one state cannot be on mass shell), the mass mixing Lagrangian through DSE provides an equally significant relation among the oscillating neutrino propagators obtained from DSE (18) and (19). They enable the calculation of the Keldysh component of the propagator S ˜ β , α , K (20), which by application of AETL provides measurable particle numbers (42).
The FTPFT approach successfully passes the test of neutrino oscillations, in the sense that we demonstrate how to perform calculations of these oscillations using the presented approach. We reproduce the standard PMNS results in the present model, which introduces the neutrino masses precisely through the PMNS-matrix in the mixing part of the model Lagrangian. This means that while the present approach does not address the issue of the origin of neutrino masses, and in general does not provide answers about the dynamics and physics of processes, the existing knowledge can be used as an input in the form of masses, matrices, self-energies, etc.
The application of the FTPFT approach to neutrino oscillations shows that it is an interesting candidate for a complementary tool to the S-matrix, which is formulated for infinite times and involves switching interactions on and off adiabatically. In the case of phenomena where finite times are essential, such as the presently pertinent neutrino oscillations as well as the oscillation of kaons and B and D mesons, decays, and symmetry violations, it has been necessary to use largely heuristic methods, such as elementary quantum mechanics in the PMNS approach and the Gell-Mann–Pais approach for kaons. In such cases, the FTPFT approach is obviously a candidate for a more rigorous description.
Further tests of the FTPFT approach could proceed along the following main lines:
  • Improving the model by taking into account the eventually-confirmed anomalies.
  • If suitable for considering decays of heavier neutrinos, the model could be easily adapted to build these features in. In this case, the self-energies should again be provided through an adequate calculation.
  • Applying the formalism to other oscillating and decay processes (e,g., decays of K 0 , D 0 , and B 0 mesons, positronium, the Cabbibo angle, etc.). This work is in progress and almost completed. It roughly confirms the Gell-Mann–Pais results. The factor which limits the predictive power is the rudimentary knowledge of the self-energies in the existing literature. Authors have mostly been concerned with obtaining the imaginary parts (decay rates), while the real parts (mass shifts) often involve renormalization.
  • In a classical out-of-equilibrium problem, the damping rates are the first thing to address. Braaten–Pisarski re-summation has provided a good start. Even for this case, two-loop self-energy diagrams contain minimal time vertices, and possible “upgrades” could be very tricky. This is another area where work is in progress.

Author Contributions

Conceptualization, I.D. and D.K.; Methodology, I.D.; Software, I.D.; Validation, I.D. and D.K.; Formal analysis, D.K.; Investigation, D.K.; Data curation, D.K.; Writing—original draft, D.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

No new data were created in this study.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Appendix A.1. Finite Time Path Field Theory

To calculate neutrino oscillations, we use Finite Time Path Field Theory (FTPFT), often called the closed time path thermal field theory (CTP-TFT). “Thermal” is used here for historical reasons, as the theory was first constructed to treat ensembles at the thermal equilibrium or very close to it. Now, we find it to be an appropriate tool for calculation of a wide variety of problems, including equilibrium and nonequilibrium TFT problems (within the linear response approximation) as well as scattering processes, decays, and oscillations. There are a number of specific features distinguishing it from S-matrix theory:
  • The time path C is closed and finite: C = ( 0 + i ϵ , t + i ϵ ) U ( t + i ϵ , t i ϵ ) U ( t i ϵ , 0 i ϵ ) .
  • The subject of the S-matrix is amplitude (“wave function”), while in FTPFT it is a two-point function.
  • The product of two point functions is not algebraic, instead being a convolution product (see Appendix A.2); only under special conditions does it become an algebraic product.
  • Instead of Feynman propagators, matrix propagators S i j , i = 1 , 2 , j = 1 , 2 are obtained. These are linearly transformed into S R (retarded), S A (advanced), and S K (Keldysh) propagators. Our method further separates S K into its retarded and advanced pieces ( S K = S K , R S K , A ); S K contains single particle distribution functions of the unperturbed system (i.e., as they are determined at t = 0 ).
  • A measured quantity is obtained as an equal time limit of S K . Compared to a scattering matrix, these measured quantities are more inclusive: one particle is separated (and measured), while the others are integrated over. This is equivalent to the exclusive S-matrix approach. In addition, the calculated quantities correspond to a yield, i.e., the number of particles found at time t, while the equivalent in the S-matrix approach is the cross-section, i.e., related to the time derivative of the yield.
  • While primarily developed for thermal equilibrium and out-of-equilibrium (particularly “almost equilibrated”) ensembles, nothing prevents it from being applied to decays and oscillations (as in this paper), or to scattering processes.
  • For application to scattering processes, it is necessary to choose initial single particle distributions (i.e., for incoming particles) as two plane waves of extremely low intensities. After calculation it is then necessary to carry out the t limit. In our experience, the results are physically equivalent to S-matrix calculation. Adiabatic switching (on and off) of the interaction is not possible in finite time. With an infinite time limit ( t ), the lack of adiabatic switching does not matter.

Appendix A.2. Convolution Product of Two Two-Point Functions

The convolution product [61] of two Green functions is defined as
C = A B C ( x , y ) = d z A ( x , z ) B ( z , y ) .
In terms of Wigner transforms (see [61] for more details), it becomes
C X 0 ( p 0 , p ) = d p 01 d p 02 P X 0 ( p 0 , p 01 + p 02 2 ) i 2 π e i X 0 ( p 01 p 02 + i ϵ ) p 01 p 02 + i ϵ A ( p 01 , p ) B ( p 02 , p ) ,
where the projector P X 0 is
P X 0 ( p 0 , p 0 ) = 1 2 π Θ ( X 0 ) 2 X 0 2 X 0 d s 0 e i s 0 ( p 0 p 0 ) = 1 π Θ ( X 0 ) sin 2 X 0 ( p 0 p 0 ) p 0 p 0
and
e i s 0 p 0 Θ ( X 0 ) Θ ( 2 X 0 + s 0 ) Θ ( 2 X 0 s 0 ) = d p 0 e i s 0 p 0 P X 0 ( p 0 , p 0 ) .
It is important to note that
lim X 0 P X 0 ( p 0 , p 0 ) = lim X 0 1 π sin 2 X 0 ( p 0 p 0 ) p 0 p 0 = δ ( p 0 p 0 ) .
The retarded (advanced) function is supposed to satisfy the following properties [61]:
(1)
The function of p 0 is analytic above (below) the real axis;
(2)
The function vanishes as | p 0 | approaches infinity in the upper (lower) semiplane.
If A is an advanced projected operator, we can integrate the expression (A2) even further. After closing the p 01 integration contour in the lower semi-plane, we obtain
C X 0 ( p 0 , p ) = d p 01 P X 0 ( p 0 , p 01 ) A ( p 01 , p ) B ( p 01 , p ) .
If B is a retarded projected operator, we can achieve the same result by closing the p 02 integration contour in the upper semi-plane.
In the X 0 limit (A5), Equation (A6) becomes the simple algebraic product
C ( p 0 , p ) = A ( p 0 , p ) B ( p 0 , p ) .
Note that Equations (A6) and (A7) are valid for combinations A A B A , A A B R , A R B R , but not for A R B A .
For A R (retarded) and B A (advanced), the Equal Time Limit (ETL) of the convolution product is obtained as follows:
1 π lim Δ 0 d p 0 e i Δ p 0 C X 0 ( p 0 , p ) = d p 01 d p 02 i 2 π e i X 0 ( p 01 p 02 ) 1 p 01 p 02 A , R ( p 01 , p ) B , A ( p 02 , p )
where Δ = s 01 s 02 and X 0 = ( s 01 + s 02 ) / 2 . The ETL is finally changed to AETL (the average of ETL), where lim Δ 0 is replaced by 1 2 [ lim 0 < Δ 0 + lim 0 > Δ 0 ] . From the numerator, the constant term ( 1 ) could have been subtracted, because for this term it is possible to close the integration path d p 01 from below (or to close the integration path d p 02 from above), in order to find that the contribution vanishes. The resulting “kernel” is not singular, and we omit the i ϵ prescription.
From now on, we skip the index as self-understandable.

Appendix A.3. Massive Neutrino Propagator

The massive neutrino propagators with masses m i 0 , i = 1 , 2 , 3 for the case when the fermion and antifermion distribution are equal:
n + ( ω p , p ) = n ( ω p , p ) is simple,
D R ( p ) = ( p + m ) G R ( p , m ) , D A ( p ) = ( p + m ) G A ( p , m ) , G R ( A ) ( p , m ) = i p 2 m 2 ± 2 i p 0 ϵ , D K ( p ) = D K , R ( p ) D K , A ( p ) , D K , R ( p ) = [ 1 2 n ( ω p ) ] ( P + m p 0 ω p ) G R ( p , m ) , D K , A ( p ) = [ 1 2 n ( ω p ) ] ( P + m p 0 ω p ) G A ( p , m ) , p = γ μ p μ , p = ( p 0 , p ) , P = γ μ P μ , P = ( ω p , p 0 ω p p ) , ω p = p 2 + m 2 ,
where n ( ω p ) is the initial fermion distribution function.

Appendix A.4. “Flavor Neutrino” Propagator

Flavor neutrinos ν e , ν μ , and ν τ are, by assumption, massless and chiral; only left-handed flavor neutrinos and right-handed antineutrinos exist:
S f , R ( p ) = p G R ( p , 0 ) , S f , A ( p ) = p G A ( p , 0 ) , G R ( A ) ( p , 0 ) = i p 2 ± 2 i p 0 ϵ , p = γ μ p μ , p = ( p 0 , p ) , P = γ μ P μ , P = ( | p | , p 0 | p | p ) .
For unequal flavor neutrino and flavor antineutrino distributions, we have
n ( p 0 , p ) = 1 γ 5 2 Θ ( p 0 ) n + ( p ) 1 γ 5 2 Θ ( p 0 ) n ( p ) ,
where n ( p 0 , p ) is now a 4 × 4 matrix.
Now, we decompose Keldysh propagator into its retarded and advanced parts:
S f , K ( p ) = S f , K , R ( p ) S f , K , A ( p ) , S f , K , R ( p ) = G R ( p , 0 ) L f ( p 0 , p ) , S f , K , A ( p ) = G A ( p , 0 ) L f ( p 0 , p ) ,
where
L f ( p 0 , p ) = [ 1 2 n f ( p ) ] p 0 + | p | 2 | p | ( γ 0 | p | γ p ) 1 γ 5 2 [ 1 2 n f ¯ ( p ) ] p 0 | p | 2 | p | ( γ 0 | p | + γ p ) 1 γ 5 2 .
We define the propagator for flavor antineutrinos by S ¯ , which is obtained from S via the replacement of n f n f ¯ . These propagators satisfy the following properties under inversion:
S f , R ( p ) = S ¯ f , A ( p ) , S f , K , R ( p ) = S ¯ f , K , A ( p ) .

Appendix A.5. Oscillating Neutrino Propagator

In addition to flavor and propagating neutrinos, we have defined “oscillating neutrinos” above by solving Equations (12)–(15) through partial re-summation of DSE over the powers of the mixing mass interaction self-energies.

Appendix A.6. Dyson–Schwinger Equation for Fermions

Here, S is the lowest order Green function and S ˜ is the re-summed Green function:
[ S 1 i Σ ] S ˜ = 1 ,
where
S = S R S K 0 S A , Σ = Σ R Σ K 0 Σ A , S ˜ = S ˜ R S ˜ K 0 S ˜ A
or in components,
S ˜ R = S R + i S R Σ R S ˜ R S ˜ A = S A + i S A Σ A S ˜ A S ˜ K = S K + i [ S R Σ K S ˜ A + S K Σ A S ˜ A + S R Σ R S ˜ K ] .
The formal solution is
S ˜ R = [ 1 i S R Σ R ] 1 S R = S R [ 1 i Σ R S R ] 1 , S ˜ A = [ 1 i S A Σ A ] 1 S A = S A [ 1 i Σ A S A ] 1 , S ˜ K = + i S ˜ R Σ K S ˜ A + S ˜ R S R 1 S K S A 1 S ˜ A
S ˜ K = S K , A ( 1 i Σ A S A ) 1 + ( 1 i S R Σ R ) 1 S K , R + i ( 1 i S R Σ R ) 1 [ S R Σ K S A S R Σ R S K , A + S K , R Σ A S A ] ( 1 i Σ A S A ) 1 ,
where we have used
S K = S K , R S K , A , Σ K = Σ K , R Σ K , A .

Appendix A.7. Fermion Particle Number

The number of fermions of momentum p at the time t is obtained from the time evolution of the of number operator, which is the Average of the Equal Time Limit (AETL) of the propagator S K , s 01 , s 02 . The average is necessary because we must distinguish the limit from above s 01 s 02 = Δ > 0 and the limit from below Δ < 0 . If the propagator S ˜ K , s 01 , s 02 ( p ) has a pole which is apart from the real axis for the finite imaginary part, these two limits will be different, in which case we take the average. Using AETL, we are able to understand the limit:
A E T L = 1 2 [ lim 0 < Δ 0 + lim 0 > Δ 0 ] .
The particle number is defined as
N p ( t ) = ( 2 π ) 3 d N / ( d 3 x d 3 p )
and
1 N f , p ( t ) = 1 2 π [ lim 0 < Δ 0 + lim 0 > Δ 0 ] d p 0 e i p 0 Δ T r [ γ 0 4 S ˜ K , t ( p ) ] ,
where Δ = s 01 s 02 and X 0 = ( s 01 + s 02 ) / 2 = t .
The contribution (A22) corresponds to a single polarization. To sum over the polarizations, an additional factor of 2 is necessary.
When acting on the momentum and spin eigenstates of fermions | + , p , s and antifermions | , p , s (both normalized to ± , ± p , s | ± , ± p , s = 1 / m ), the projectors Λ ± satisfy
Λ + ( ω p , p ) | + , p , s = | + , p , s and Λ ( ω p , p ) | + , p , s = 0 , while
Λ + ( ω p , p ) | , p , s = 0 and Λ ( ω p , p ) | , p , s = | , p , s .
In the rest of this paper, the distribution functions are assumed to depend on | p | and not on the direction of p . The lowest order contribution ( N f + f ¯ , p 0 ( t ) ) is identical to the initial distribution (for n f ¯ ( ω p ) = 0 ):
1 N f , p 0 ( t ) = 1 n + f ( ω p ) .

Appendix A.8. Massless Chiral Fermion Particle Number

For a massless chiral particle there is one more helicity projector ( 1 γ 5 ) / 2 with respect to the trace in Equation (A22). The particle number becomes (22):
1 N β , p ( t ) = 1 2 π [ lim 0 < Δ 0 + lim 0 > Δ 0 ] d p 0 e i p 0 Δ T r [ γ 0 4 S ˜ β , β , K , t ( p ) 1 γ 5 2 ] .
Here, N β , p ( t ) is the number of neutrinos of flavor f and momentum p detected at time t.
Note that one more helicity projector 1 γ 5 2 is necessary to ensure the proper helicity (chirality) of the detected neutrinos. The lowest-order contribution for a massless chiral particle is identical to the initial distribution:
1 N f p 0 ( t ) = 1 n f ( | p | ) .

References

  1. Davis, R., Jr.; Harmer, D.S.; Hoffman, K.C. Search for neutrinos from the sun. Phys. Rev. Lett. 1968, 20, 1205–1209. [Google Scholar] [CrossRef]
  2. Fukuda, Y.; Hayakawa, T.; Ichihara, E.; Inoue, K.; Ishihara, K.; Ishino, H.; Itow, Y.; Kajita, T.; Kameda, J.; Kasuga, S.; et al. Evidence for oscillation of atmospheric neutrinos. Phys. Rev. Lett. 1998, 81, 1562. [Google Scholar] [CrossRef]
  3. Apollonio, M.; Baldini, A.; Bemporad, C.; Caffau, E.; Cei, F.; D’eclaise, Y.; de Kerret, H.; Dieterle, B.; Etenko, A.; George, J.; et al. Limits on neutrino oscillations from the CHOOZ experiment. Phys. Lett. B 1999, 466, 415–430. [Google Scholar] [CrossRef]
  4. Abazajian, K.N.; Acero, M.A.; Agarwalla, S.K.; Aguilar-Arevalo, A.A.; Albright, C.H.; Antusch, S.; Arguelles, C.A.; Balantekin, A.B.; Barenboim, G.; Barger, V.; et al. Light Sterile Neutrinos: A White Paper. arXiv 2012, arXiv:1204.5379. [Google Scholar]
  5. Aguilar, A.; Auerbach, L.B.; Burman, R.L.; Caldwell, D.O.; Church, E.D.; Cochran, E.D.; Donahue, J.B.; Fazely, A.; Garvey, G.T.; Gunasingha, R.M.; et al. Evidence for neutrino oscillations from the observation of νe appearance in a νμ beam. Phys. Rev. D 2001, 64, 112007. [Google Scholar] [CrossRef]
  6. Schechter, J.; Valle, J.W.F. Neutrino Masses in SU(2) x U(1) Theories. Phys. Rev. D 1980, 22, 2227. [Google Scholar] [CrossRef]
  7. Eguchi, K.; Enomoto, S.; Furuno, K.; Goldman, J.; Hanada, H.; Ikeda, H.; Ikeda, K.; Inoue, K.; Ishihara, K.; Itoh, W.; et al. First results from KamLAND: Evidence for reactor anti-neutrino disappearance. Phys. Rev. Lett. 2003, 90, 021802. [Google Scholar] [CrossRef]
  8. Hosaka, J.; Ishihara, K.; Kameda, J.; Koshio, Y.; Minamino, A.; Mitsuda, C.; Miura, M.; Moriyama, S.; Nakahata, M.; Namba, T.; et al. Solar neutrino measurements in super-Kamiokande-I. Phys. Rev. D 2006, 73, 112001. [Google Scholar] [CrossRef]
  9. Cervera, A.; Donini, A.; Gavela, M.B.; Cadenas, J.J.G.; Hernandez, P.; Mena, O.; Rigolin, S. Golden measurements at a neutrino factory. Nucl. Phys. B 2001, 579, 17–55. [Google Scholar] [CrossRef]
  10. Barger, V.; Marfatia, D.; Whisnant, K. Breaking eight fold degeneracies in neutrino CP violation, mixing, and mass hierarchy. Phys. Rev. D 2002, 65, 073023. [Google Scholar] [CrossRef]
  11. Yasuda, O. New plots and parameter degeneracies in neutrino oscillations. New J. Phys. 2004, 6, 83. [Google Scholar] [CrossRef]
  12. Burguet-Castell, J.; Gavela, M.B.; Gomez-Cadenas, J.J.; Hernandez, P.; Mena, O. On the Measurement of leptonic CP violation. Nucl. Phys. B 2001, 608, 301. [Google Scholar] [CrossRef]
  13. Minakata, H.; Nunokawa, H. Exploring neutrino mixing with low-energy superbeams. J. High Energy Phys. 2001, 001. [Google Scholar] [CrossRef]
  14. Fogli, G.L.; Lisi, E. Tests of three flavor mixing in long baseline neutrino oscillation experiments. Phys. Rev. D 1996, 54, 3667. [Google Scholar] [CrossRef] [PubMed]
  15. Abe, K.; Abgrall, N.; Ajima, Y.; Aihara, H.; Albert, J.B.; Andreopoulos, C.; Andrieu, B.; Aoki, S.; Araoka, O.; Argyriades, J.; et al. Indication of Electron Neutrino Appearance from an Accelerator-produced Off-axis Muon Neutrino Beam. Phys. Rev. Lett. 2011, 107, 041801. [Google Scholar] [CrossRef] [PubMed]
  16. Adamson, P.; Auty, D.J.; Ayres, D.S.; Backhouse, C.; Barr, G.; Betancourt, M.; Bishai, M.; Blake, A.; Bock, G.J.; Boehnlein, D.J.; et al. Active to sterile neutrino mixing limits from neutral-current interactions in MINOS. Phys. Rev. Lett. 2011, 107, 181802. [Google Scholar] [CrossRef] [PubMed]
  17. Abe, Y.; Aberle, C.; Akiri, T.; dos Anjos, J.C.; Ardellier, A.F.; Barbosa, F.; Baxter, A.; Bergevin, M.; Bernstein, A.; Bezerra, T.J.C.; et al. Indication of Reactor νe Disappearance in the Double Chooz Experiment. Phys. Rev. Lett. 2012, 108, 131801. [Google Scholar] [CrossRef]
  18. An, F.P.; Bai1, J.Z.; Balantekin, A.B.; Band, H.R.; Beavis, D.; Beriguete, W.; Bishai, M.; Blyth, S.; Boddy, K.; Brown, R.L.; et al. Observation of electron-antineutrino disappearance at Daya Bay. Phys. Rev. Lett. 2012, 108, 171803. [Google Scholar] [CrossRef]
  19. Ahn, J.K.; Chebotaryov, S.; Choi, J.H.; Choi, S.; Choi, W.; Choi, Y.; Jang, H.I.; Jang, J.S.; Jeon, E.J.; Jeong, I.S.; et al. Observation of Reactor Electron Antineutrino Disappearance in the RENO Experiment. Phys. Rev. Lett. 2012, 108, 191802. [Google Scholar] [CrossRef]
  20. Abe, K.; Akhlaq, N.; Akutsu, R.; Ali, A.; Alt, C.; Andreopoulos, C.; Antonova, M.; Aoki, S.; Arihara, T.; Asada, Y.; et al. Neutron-antineutron oscillation search using a 0.37 megaton-years exposure of Super-Kamiokande. Phys. Rev. D 2021, 103, 112008. [Google Scholar] [CrossRef]
  21. Acero, M.A.; Adamson, P.; Aliaga, L.; Anfimov, N.; Antoshkin, A.; Arrieta-Diaz, E.; Asquith, L.; Aurisano, A.; Back, A.; Backhouse, C.; et al. Improved measurement of neutrino oscillation parameters by the NOvA experiment. Phys. Rev. D 2022, 106, 032004. [Google Scholar] [CrossRef]
  22. Gonzalez-Garcia, M.C.; Maltoni, M.; Schwetz, T. NuFIT: Three-Flavour Global Analyses of Neutrino Oscillation Experiments. Universe 2021, 7, 459. [Google Scholar] [CrossRef]
  23. Capozzi, F.; Valentino, E.D.; Lisi, E.; Marrone, A.; Melchiorri, A.; Palazzo, A. Unfinished fabric of the three neutrino paradigm. Phys. Rev. D 2021, 104, 083031. [Google Scholar] [CrossRef]
  24. de Salas, P.F.; Forero, D.V.; Gariazzo, S.; Martínez-Miravé, P.; Mena, O.; Ternes, C.A.; Tórtola, M.; Valle, J.W.F. 2020 global reassessment of the neutrino oscillation picture. J. High Energy Phys. 2021, 2, 71. [Google Scholar] [CrossRef]
  25. Cleveland, B.T.; Daily, T.; Davis, R., Jr.; Distel, J.R.; Lande, K.; Lee, C.K.; Wildenhain, P.S.; Ullman, J. Measurement of the solar electron neutrino flux with the Homestake chlorine detector. Astrophys. J. 1998, 496, 505–526. [Google Scholar] [CrossRef]
  26. Kaether, F.; Hampel, W.; Heusser, G.; Kiko, J.; Kirsten, T. Reanalysis of the GALLEX solar neutrino flux and source experiments. Phys. Lett. B 2010, 685, 47–54. [Google Scholar] [CrossRef]
  27. Abdurashitov, J.N.; Gavrin, V.N.; Gorbachev, V.V.; Gurkina, P.P.; Ibragimova, T.V.; Kalikhov, A.V.; Khairnasov, N.G.; Knodel, T.V.; Mirmov, I.N.; Shikhin, A.A.; et al. Measurement of the solar neutrino capture rate with gallium metal. III: Results for the 2002–2007 data-taking period. Phys. Rev. C 2009, 80, 015807. [Google Scholar] [CrossRef]
  28. Cravens, J.P.; Abe, K.; Iida, T.; Ishihara, K.; Kameda, J.; Koshio, Y.; Minamino, A.; Mitsuda, C.; Miura, M.; Moriyama, S.; et al. Solar neutrino measurements in Super-Kamiokande-II. Phys. Rev. D 2008, 78, 032002. [Google Scholar] [CrossRef]
  29. Abe, K.; Hayato, Y.; Iida, T.; Ikeda, M.; Ishihara, C.; Iyogi, K.; Kameda, J.; Kobayashi, K.; Koshio, Y.; Kozuma, Y.; et al. Solar neutrino results in Super-Kamiokande-III. Phys. Rev. D 2011, 83, 052010. [Google Scholar] [CrossRef]
  30. Nakajima, Y. Recent results and future prospects from Super-Kamiokande. In Proceedings of the XXIX International Conference on Neutrino Physics and Astrophysics, Chicago, IL, USA, 22 June–2 July 2020. [Google Scholar] [CrossRef]
  31. Aharmim, B.; Ahmed, S.N.; Anthony, A.E.; Barros, N.; Beier, E.W.; Bellerive, A.; Beltran, B.; Bergevin, M.; Biller, S.D.; Boudjemline, K.; et al. Combined Analysis of all Three Phases of Solar Neutrino Data from the Sudbury Neutrino Observatory. Phys. Rev. C 2013, 88, 025501. [Google Scholar] [CrossRef]
  32. Bellini, G.; Benziger, J.; Bick, D.; Bonetti, S.; Bonfini, G.; Avanzini, M.B.; Caccianiga, B.; Cadonati, L.; Calaprice, F.; Carraro, C.; et al. Precision measurement of the 7Be solar neutrino interaction rate in Borexino. Phys. Rev. Lett. 2011, 107, 141302. [Google Scholar] [CrossRef] [PubMed]
  33. Bellini, G.; Benziger, J.; Bonetti, S.; Buizza Avanzini, M.; Caccianiga, B.; Cadonati, L.; Calaprice, F.; Carraro, C.; Chavarria, A.; Chepurnov, A.; et al. Measurement of the solar 8B neutrino rate with a liquid scintillator target and 3 MeV energy threshold in the Borexino detector. Phys. Rev. D 2010, 82, 033006. [Google Scholar] [CrossRef]
  34. Borexino Collaboration; Bellini, G.; Benziger, J.; Bick, D.; Bonfini, G.; Bravo, D.; Caccianiga, B.; Cadonati, L.; Calaprice, F.; Caminata, A.; et al. Neutrinos from the primary proton–proton fusion process in the Sun. Nature 2014, 512, 383–386. [Google Scholar]
  35. Gando, A.; Gando, Y.; Hanakago, H.; Ikeda, H.; Inoue, K.; Ishidoshiro, K.; Ishikawa, H.; Koga, M.; Matsuda, R.; Matsuda, S.; et al. Reactor On-Off Antineutrino Measurement with KamLAND. Phys. Rev. D 2013, 88, 033001. [Google Scholar] [CrossRef]
  36. Abe, K.; Bronner, C.; Haga, Y.; Hayato, Y.; Ikeda, M.; Iyogi, K.; Kameda, J.; Kato, Y.; Kishimoto, Y.; Marti, L.; et al. Atmospheric neutrino oscillation analysis with external constraints in Super-Kamiokande I-IV. Phys. Rev. D 2018, 97, 072001. [Google Scholar] [CrossRef]
  37. Aartsen, M.G.; Ackermann, M.; Adams, J.; Aguilar, J.A.; Ahlers, M.; Ahrens, M.; Altmann, D.; Anderson, T.; Arguelles, C.; Arlen, T.C.; et al. Measurement of the Atmospheric νe Spectrum with IceCube. Phys. Rev. D 2015, 91, 072004. [Google Scholar] [CrossRef]
  38. Bezerra, T. New Results from the Double Chooz Experiment. In Proceedings of the XXIX International Conference on Neutrino Physics and Astrophysics, Chicago, IL, USA, 22 June–2 July 2020. [Google Scholar] [CrossRef]
  39. Adey, D.; An, F.P.; Balantekin, A.B.; Band, H.R.; Bishai, M.; Blyth, S.; Cao, D.; Cao, G.F.; Cao, J.; Chan, Y.L. Measurement of the Electron Antineutrino Oscillation with 1958 Days of Operation at Daya Bay. Phys. Rev. Lett. 2018, 121, 241805. [Google Scholar] [CrossRef]
  40. Yoo, J. Recent Results from RENO Experiment. In Proceedings of the XXIX International Conference on Neutrino Physics and Astrophysics, Chicago, IL, USA, 22 June–2 July 2020. [Google Scholar] [CrossRef]
  41. Adamson, P.; Anghel, I.; Backhouse, C.; Barr, G.; Bishai, M.; Blake, A.; Bock, G.J.; Bogert, D.; Cao, S.V.; Castromonte, C.M. Measurement of neutrino and antineutrino oscillations using beam and atmospheric data in MINOS. Phys. Rev. Lett. 2013, 110, 251801. [Google Scholar] [CrossRef]
  42. Dunne, P. Latest Neutrino Oscillation Results from T2K. In Proceedings of the XXIX International Conference on Neutrino Physics and Astrophysics, Chicago, IL, USA, 22 June–2 July 2020. [Google Scholar] [CrossRef]
  43. Himmel, A. New Oscillation Results from the NOvA Experiment. In Proceedings of the XXIX International Conference on Neutrino Physics and Astrophysics, Chicago, IL, USA, 22 June–2 July 2020. [Google Scholar] [CrossRef]
  44. Altherr, T.; Salati, P. The Electric charge of neutrinos and plasmon decay. Nucl. Phys. B 1994, 421, 662–682. [Google Scholar] [CrossRef]
  45. Ho, C.M.; Boyanovsky, D.; de Vega, H.J. Neutrino oscillations in the early universe: A Real time formulation. Phys. Rev. D 2005, 72, 085016. [Google Scholar] [CrossRef]
  46. Blasone, M.; Jizba, P.; Mavromatos, N.E.; Smaldone, L. Dynamical generation of field mixing via flavor vacuum condensate. Phys. Rev. D 2019, 100, 045027. [Google Scholar] [CrossRef]
  47. Blasone, M.; Smaldone, L.; Vitiello, G. Neutrino mixing, entanglement and the gauge paradigm in quantum field theory. PoSCORFU2021 2022, 049. [Google Scholar] [CrossRef]
  48. Blasone, M.; Smaldone, L.; Vitiello, G. Non-Abelian gauge structure and flavor mixing in quantum field theory. J. Phys. Conf. Ser. 2023, 2533, 012010. [Google Scholar] [CrossRef]
  49. Schwinger, J. Brownian motion of a quantum oscillator. J. Math. Phys. 1961, 2, 407. [Google Scholar] [CrossRef]
  50. Keldysh, L.V. Diagram technique for nonequilibrium processes. Sov. Phys. JETP 1965, 20, 1018–1026. [Google Scholar]
  51. Kadanoff, L.P.; Baym, G. Quantum Statistical Mechanics; Benjamin: New York, NY, USA, 1962. [Google Scholar]
  52. Danielewicz, P. Quantum Theory of Nonequilibrium Processes. 1. Ann. Phys. 1984, 152, 239. [Google Scholar] [CrossRef]
  53. Chou, K.-C.; Su, Z.-B.; Hao, B.-L.; Yu, L. Equilibrium and Nonequilibrium Formalisms Made Unified. Phys. Rep. 1985, 118, 1. [Google Scholar] [CrossRef]
  54. Rammer, J.; Smith, H. Quantum field-theoretical methods in transport theory of metals. Rev. Mod. Phys. 1986, 58, 323. [Google Scholar] [CrossRef]
  55. Landsman, N.P.; van Weert, C.G. Real and Imaginary Time Field Theory at Finite Temperature and Density. Phys. Rep. 1987, 145, 141. [Google Scholar] [CrossRef]
  56. Calzetta, E.; Hu, B.L. Nonequilibrium Quantum Fields: Closed Time Path Effective Action, Wigner Function and Boltzmann Equation. Phys. Rev. D 1988, 37, 2878. [Google Scholar] [CrossRef]
  57. Bellac, M.L. Thermal Field Theory; Cambridge University Press: Cambridge, UK, 1996. [Google Scholar]
  58. Brown, D.A.; Danielewicz, P. Partons in phase space. Phys. Rev. D 1998, 58, 094003. [Google Scholar] [CrossRef]
  59. Dadić, I. Two mechanisms for elimination of pinch singularities in/out of equilibrium thermal field theories. Phys. Rev. D 1999, 59, 125012. [Google Scholar] [CrossRef]
  60. Blaizot, J.-P.; Iancu, E. The Quark gluon plasma: Collective dynamics and hard thermal loops. Phys. Rep. 2002, 359, 355–528. [Google Scholar] [CrossRef]
  61. Dadić, I. Out of equilibrium thermal field theories: Finite time after switching on the interaction and Wigner transforms of projected functions. Phys. Rev. D 2000, 63, 025011. [Google Scholar] [CrossRef]
  62. Dadić, I. Out-of-equilibrium TFT-energy nonconservation at vertices. Nucl. Phys. A 2002, 702, 356. [Google Scholar] [CrossRef]
  63. Arleo, F.; Aurenche, P.; Bopp, F.W.; Dadić, I.; David, G.; Delagrange, H.; d’Enterria, D.G.; Eskola, K.J.; Gelis, F.; Guillet, J.P. Hard probes in heavy-ion collisions at the LHC: Photon physics in heavy ion collisions at the LHC. CERN Yellow Book CERN-2004-009-D. arXiv 2004, arXiv:hep-ph/0311131. [Google Scholar]
  64. Dadić, I. Retarded propagator representation of out-of-equilibrium thermal field theories. Nucl. Phys. A 2009, 820, 267C. [Google Scholar] [CrossRef]
  65. Millington, P.; Pilaftsis, A. Perturbative nonequilibrium thermal field theory. Phys. Rev. D 2013, 88, 085009. [Google Scholar] [CrossRef]
  66. Millington, P.; Pilaftsis, A. Thermal field theory to all orders in gradient expansion. J. Phys. Conf. Ser. 2013, 447, 012071. [Google Scholar] [CrossRef]
  67. Dadić, I.; Klabučar, D. Causality and Renormalization in Finite-Time-Path Out-of-Equilibrium ϕ3 QFT. Particles 2019, 2, 92–102. [Google Scholar] [CrossRef]
  68. Dadić, I.; Klabučar, D.; Kuić, D. Direct photons from Quark Matter in Renormalized Finite-Time-Path QED. Particles 2020, 3, 676–692. [Google Scholar] [CrossRef]
  69. Beuthe, M. Oscillations of neutrinos and mesons in quantum field theory. Phys. Rept. 2003, 375, 105–218. [Google Scholar] [CrossRef]
  70. Akhmedov, E.K.; Kopp, J. Neutrino Oscillations: Quantum Mechanics vs. Quantum Field Theory. J. High Energy Phys. 2013, 4, 8, Erratum in J. High Energy Phys. 2013, 10, 52. [Google Scholar] [CrossRef]
  71. Pontecorvo, B. Mesonium and anti-mesonium. Sov. Phys. JETP 1957, 6, 429. [Google Scholar]
  72. Maki, Z.; Nakagawa, M.; Sakata, S. Remarks on the unified model of elementary particles. Prog. Theor. Phys. 1962, 28, 870. [Google Scholar] [CrossRef]
  73. Wolfenstein, L. Neutrino Oscillations in Matter. Phys. Rev. D 1978, 17, 2369. [Google Scholar] [CrossRef]
  74. Gribov, V.N.; Pontecorvo, B. Neutrino astronomy and lepton charge. Phys. Lett. B 1969, 28, 493–496. [Google Scholar] [CrossRef]
  75. Mikheyev, S.P.; Smirnov, A.Y. Resonance Amplification of Oscillations in Matter and Spectroscopy of Solar Neutrinos. Sov. J. Nucl. Phys. 1985, 42, 913–917. [Google Scholar]
  76. Aguilar-Arevalo, A.A.; Brown, B.C.; Conrad, J.M.; Dharmapalan, R.; Diaz, A.; Djurcic, Z.; Finley, D.A.; Ford, R.; Garvey, G.T.; Gollapinni, S. Updated MiniBooNE neutrino oscillation results with increased data and new background studies. Phys. Rev. D 2021, 103, 052002. [Google Scholar] [CrossRef]
  77. Mueller, T.A.; Lhuillier, D.; Fallot, M.; Letourneau, A.; Cormon, S.; Fechner, M.; Giot, L.; Lasserre, T.; Martino, J.; Mention, G.; et al. Improved Predictions of Reactor Antineutrino Spectra. Phys. Rev. C 2011, 83, 054615. [Google Scholar] [CrossRef]
  78. Huber, P. On the determination of anti-neutrino spectra from nuclear reactors. Phys. Rev. C 2011, 84, 024617. [Google Scholar] [CrossRef]
  79. Giunti, C.; Laveder, M. Neutron oscillations. Phys. Rev. D 2008, 77, 093002. [Google Scholar] [CrossRef]
  80. Aartsen, M.G.; Abbasi, R.; Ackermann, M.; Adams, J.; Aguilar, J.A.; Ahlers, M.; Ahrens, M.; Alispach, C.; Amin, N.M.; Andeen, K. eV-Scale Sterile Neutrino Search Using Eight Years of Atmospheric Muon Neutrino Data from the IceCube Neutrino Observatory. et al. Phys. Rev. Lett. 2020, 125, 141801. [Google Scholar] [CrossRef] [PubMed]
  81. Abe, K.; Haga, Y.; Hayato, Y.; Ikeda, M.; Iyogi, K.; Kameda, J.; Kishimoto, Y.; Marti, L.; Miura, M.; Moriyama, S. Solar Neutrino Measurements in Super-Kamiokande-IV. Phys. Rev. D 2016, 94, 052010. [Google Scholar] [CrossRef]
  82. Gonzalez-Garcia, M.C.; Maltoni, M. Determination of matter potential from global analysis of neutrino oscillation data. J. High Energy Phys. 2013, 1309, 152. [Google Scholar] [CrossRef]
  83. Maltoni, M.; Smirnov, A.Y. Solar neutrinos and neutrino physics. Eur. Phys. J. A 2016, 52, 87. [Google Scholar] [CrossRef]
  84. Guzzo, M.M.; Masiero, A.; Petcov, S.T. On the MSW effect with massless neutrinos and no mixing in the vacuum. Phys. Lett. B 1991, 260, 154. [Google Scholar] [CrossRef]
  85. Roulet, E. MSW effect with flavor changing neutrino interactions. Phys. Rev. D 1991, 44, 935. [Google Scholar] [CrossRef]
  86. Yasuda, O. Neutrino oscillation phenomenology and the impact of Professor Masatoshi Koshiba. Prog. Theor. Exp. Phys. 2022, 12, 12B105. [Google Scholar] [CrossRef]
  87. Alves, G.F.S.; Bertuzzo, E.; Salla, G.M. On-shell approach to neutrino oscillations. Phys. Rev. D 2022, 106, 036028. [Google Scholar] [CrossRef]
  88. Banks, H.; Kelly, K.J.; McCullough, M. How broad is a neutrino. J. High Energy Phys. 2023, 02, 136. [Google Scholar] [CrossRef]
  89. Bilenky, S.M. Flavor neutrino oscillations and time-energy uncertainty relation. Phys. Scr. T 2006, 127, 8. [Google Scholar] [CrossRef]
  90. Bilenky, S.M.; von Feilitzsch, F.; Potzel, W. Neutrino Oscillations and Uncertainty Relations. J. Phys. G 2011, 38, 115002. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Dadić, I.; Klabučar, D. Neutrino Oscillations in Finite Time Path Out-of-Equilibrium Thermal Field Theory. Symmetry 2023, 15, 1970. https://doi.org/10.3390/sym15111970

AMA Style

Dadić I, Klabučar D. Neutrino Oscillations in Finite Time Path Out-of-Equilibrium Thermal Field Theory. Symmetry. 2023; 15(11):1970. https://doi.org/10.3390/sym15111970

Chicago/Turabian Style

Dadić, Ivan, and Dubravko Klabučar. 2023. "Neutrino Oscillations in Finite Time Path Out-of-Equilibrium Thermal Field Theory" Symmetry 15, no. 11: 1970. https://doi.org/10.3390/sym15111970

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop