Next Article in Journal
Semi-Hyers–Ulam–Rassias Stability of a Volterra Integro-Differential Equation of Order I with a Convolution Type Kernel via Laplace Transform
Previous Article in Journal
An Eco-Epidemiological Model Incorporating Harvesting Factors
Previous Article in Special Issue
Hidden Dynamical Symmetry and Quantum Thermodynamics from the First Principles: Quantized Small Environment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Neutrino Oscillations in the Model of Interaction of Spinor Fields with Zero-Range Potential Concentrated on a Plane

1
Department of High Energy and Elementary Particle Physics, Faculty of Physics, St. Petersburg State University, Ulyanovskaya Str. 1-3, Peterhof, 198504 St. Petersburg, Russia
2
Military Educational and Scientific Center of the Navy, Naval Polytechnic Institute, Naval Academy, Kadetsky Boulevard, 1, Pushkin, 196604 St. Petersburg, Russia
*
Author to whom correspondence should be addressed.
Symmetry 2021, 13(11), 2180; https://doi.org/10.3390/sym13112180
Submission received: 23 July 2021 / Revised: 19 October 2021 / Accepted: 4 November 2021 / Published: 15 November 2021
(This article belongs to the Special Issue Symmetry in Particle Physics II)

Abstract

:
Symanzik’s approach to the description of quantum field systems in an inhomogeneous space-time is used to construct a model for the interaction of neutrino fields with matter. In this way, the problem of the influence of strong inhomogeneities of the medium on the processes of oscillations is considered. As a simple example, a model of neutrino scattering on a material plane is investigated. Within this model, in the collisions of particles with planes, a special filtration mechanism can be formed. It has a significant impact on the dynamics of subsequent neutrino oscillations which are analogous to the Mikheev-Smirnov-Wolfenstein effect at propagation of these particles in an adiabatic medium. Taking into account the possibility of the filtration process in a highly inhomogeneous environment can be useful in planning and carrying out experimental studies of neutrino physics. It can also be considered by investigations of the role of neutrino in astrophysical processes by means of numerical simulations methods.

1. Introduction

Recently, much attention has been paid to research in the field of neutrino physics. Here, the effects of neutrino oscillations predicted theoretically more than 60 years ago [1,2,3,4,5] are finding new applications [6,7,8,9], and after the experimental confirmation of their existence, new interest has arisen to such processes [10,11,12,13]. It should be noted that the current stage of theoretical studies of neutrino oscillations is characterized by the desire to develop for them a quantum-field formalism that is more adequate for their physical essence instead of the phenomenological approach that has been used so far, based on the general principles of quantum mechanics [11,14,15].
The features of the interaction of neutrinos with the material environment are also being intensively investigated [16,17,18,19,20,21]. In order to describe such phenomena, we propose in this paper to apply the Symanzik approach for modeling interactions of quantized fields with space-time inhomogeneities [22]. It was applied to study the effects of interactions of quantum electrodynamics (QED) fields with two-dimensional materials [23,24]. Employing basic physical requirements of locality, renormalizability, and gauge invariance, which are assumed in this approach, we impose significant restrictions on the possible form of the model [25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43].
The interaction of photon fields with a two-dimensional (2D) object is carried out by the Chern-Simons potential, including only one dimensionless parameter and concentrated on a 2D-subspace of full Minkowski space-time [23,24,27,28,29,30,31,32]. To describe the interaction of homogeneous and isotropic material plane with spinor field, one needs, in the most general case, no more than eight dimensionless constants [33,34,35,36,37,38,39]. Within the framework of models constructed in this way, a number of unusual effects have been found. In particular, it is shown that the Casimir force between thin plates of an uncharged capacitor can not only be attractive, but also repulsive [23,24]. Using the Symanzik approach, we constructed models of interaction of Dirac particles with 2D materials, studied their bound states, scattering on a material plane [33,34,35,36,37,38,39,40,41], and calculated the Casimir force generated by vacuum fluctuations of quantum massless Dirac fields between two parallel thin plates [42].
Taking into account this experience, we apply the Symanzik approach for modeling the interaction of neutrino fields with strongly inhomogeneous material. As an example of this situation, we consider the interaction of neutrinos with a material plane and analyse the influence of collisions with it on the oscillation processes.

2. Neutrino Dynamics in the Free Field Approximation

Over the past 20 years, we have seen an intensive development in the physics of materials. One of the theoretical problems in this area is the construction of models of the interaction between 2D objects and fields of QED. In order to find possible methods for its solution, it was suggested to employ the Symansik approach [22] for construction of the QED models with space-time inhomogeneities interpreted as a description of material environments [23,24]. On this basis, modifications of QED were developed for modeling the interaction between QED fields and 2D materials. Some effects of this interaction were investigated and are presented in [23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43].
To describe the processes called neutrino oscillations, one uses a model of a system with three pairs of four-component spinor fields ψ ¯ λ j ( x ) , ψ λ j ( x ) , λ = 1 , 2 , 3 , j = 1 , 2 , 3 , 4 which a free-action functional reads as [6,11,14,15]
S 0 ( Ψ ¯ , Ψ , M ) = Ψ ¯ ( x ) ( i ^ + M ) Ψ ( x ) d x = λ , λ = 1 3 j , j = 1 4 ψ ¯ λ j ( x ) ( i ^ j j δ λ λ + M λ λ δ j j ) ψ λ j ( x ) d x
where δ λ λ , δ j j are the Kronecker delta-symbols, M is a Hermitean ( 3 × 3 ) - matrix with three eigenvalues m μ and corresponding normalized three-component eigenvectors e μ :
δ λ λ = 0 , b y λ λ , δ λ λ = 1 , λ = 1 3 M λ λ e μ λ = m μ e μ λ , ( e μ * e μ ) = λ = 1 3 e μ λ * e μ λ = δ μ μ , ( e λ * e λ ) = μ = 1 3 e μ λ * e μ λ = δ λ λ , M λ λ = μ = 1 3 m μ e μ λ e μ λ * , λ , λ ; μ , μ = 1 , 2 , 3 .
We assume that 0 m 1 m 2 m 3 . In (1), we used the notation Ψ ( x ) = { ψ 1 ( x ) , ψ 2 ( x ) , ψ 3 ( x ) } , Ψ ¯ ( x ) = { ψ ¯ 1 ( x ) , ψ ¯ 2 ( x ) , ψ ¯ 3 ( x ) } , and ^ is a Lorentz invariant scalar product of four-differential vector with a four-Dirac matrix vector
^ = ν = 0 3 γ ν ν , ν = x ν , γ 0 = τ 0 0 0 τ 0 , γ 1 = 0 τ 1 τ 1 0 , γ 2 = 0 τ 2 τ 2 0 , γ 3 = 0 τ 3 τ 3 0 ,
where τ 0 is the unique ( 2 × 2 ) matrix, and τ 1 , τ 2 , τ 3 are the Pauli matrices
τ 0 = 1 0 0 1 , τ 1 = 0 1 1 0 , τ 2 = 0 i i 0 , τ 3 = 1 0 0 1 .
It is supposed that the spinor fields ψ ¯ λ ( x ) , ψ λ ( x ) and the matrix M λ λ are convenient for direct description of neutrino physics and its experimentally observed features.
Using notations
ϕ μ ( x ) = λ = 0 3 e μ λ * ψ λ ( x ) , ϕ ¯ μ ( x ) = λ = 0 3 e μ λ ψ ¯ λ ( x ) , M μ μ = λ , λ = 1 3 e μ λ * M λ λ e μ λ = m μ δ μ μ ,
we can write the free action (1) of the model in terms of the fields Φ ¯ ( x ) = { ϕ ¯ 1 ( x ) , ϕ ¯ 2 ( x ) , ϕ ¯ 3 ( x ) } , Φ ( x ) = { ϕ 1 ( x ) , ϕ 2 ( x ) , ϕ 3 ( x ) } as
S 0 ( Ψ ¯ , Ψ , M ) = S 0 ( Φ ¯ , Φ , M ) = Φ ¯ ( i ^ M ) Φ = μ = 1 3 j j = 1 4 ϕ ¯ μ j ( x ) ( i ^ j j m μ δ j j ) ϕ μ j ( x ) .
One says that the system is considered in a lepton (also-called flavor) representation, if fields Ψ ¯ ( x ) , Ψ ( x ) and the non-diagonal mass matrix M are used for its description. In the so-called mass representation, the system states are characterized by the fields Φ ¯ ( x ) , Φ ( x ) and diagonal mass matrix M (i.e., by the masses m μ , μ = 1 , 2 , 3 ). For writing indices, we will use the letter λ in the lepton representation and the letter μ in the mass one.
The considered system of spinor fields can be characterised by the local, independent from representation, bilinear function G Γ ( μ ) ( x ) = G Γ ( λ ) ( x ) defined by a ( 4 × 4 ) —matrix Γ as follows
G Γ ( μ ) ( x ) = Φ ¯ ( x ) Γ Φ ( x ) = μ = 1 3 j , j = 1 4 ϕ ¯ μ j ( x ) Γ j j ϕ μ j ( x ) = λ , λ , μ = 1 3 e μ λ * e μ λ ( ψ ¯ λ ( x ) Γ ψ λ ( x ) ) = λ , λ = 1 3 δ λ λ ( ψ ¯ λ ( x ) Γ ψ λ ( x ) ) = λ 3 ψ ¯ λ ( x ) Γ ψ λ ( x ) = Ψ ¯ ( x ) Γ Ψ ( x ) = G Γ ( λ ) ( x ) .
However, properties of the components ϕ ¯ μ ( x ) Γ ϕ μ ( x ) , ψ ¯ λ ( x ) Γ ψ λ ( x ) of G Γ ( μ ) ( x ) and G Γ ( λ ) ( x ) appear to be essentially different.
In the stationarity point of S 0 ( Φ ¯ , Φ , M ) the fields ϕ ¯ μ ( x ) , ϕ μ ( x ) satisfy the Dirac equations
( i ^ m μ ) ϕ μ ( x ) = 0 , i ν ϕ ¯ μ ( x ) γ ν + m μ ϕ ¯ μ ( x ) = 0 , μ = 1 , 2 , 3 .
If one chooses ϕ ¯ μ ( x ) , ϕ μ ( x ) as their plane wave solutions
ϕ μ ( x ) = e i p μ x χ μ ( p μ ) , ϕ ¯ μ ( x ) = e i p μ x χ ¯ μ ( p μ ) , ( p ^ μ m μ ) χ μ ( p μ ) = 0 , χ ¯ μ ( p μ ) ( p ^ μ m μ ) = 0 , p μ = { p μ 0 , p μ 1 , p μ 2 , p μ 3 } , p μ 2 = p μ 0 2 p μ 1 2 p μ 2 2 p μ 3 2 = m μ 2 ,
then ϕ ¯ μ ( x ) Γ ϕ μ ( x ) = χ ¯ μ ( p μ ) Γ χ μ ( p μ ) does not depend on the spase-time point x. For a similar quantity of flavor representation, one obtains
ψ ¯ λ ( x ) Γ ψ λ ( x ) = μ , μ 3 e μ λ * e μ λ ϕ ¯ μ ( x ) Γ ϕ μ ( x ) = μ = 1 3 e μ λ * e μ λ χ ¯ μ ( p μ ) Γ χ μ ( p μ ) + μ μ = 1 3 e μ λ * e μ λ e i ( p μ p μ ) x χ ¯ μ ( p μ ) Γ χ μ ( p μ ) .
The dependence on the point x of this expression is determined by the factors e i ( p μ p μ ) x . If space parts p μ = { p μ 1 , p μ 2 , p μ 3 } of the four moments p μ = { p μ 0 , p μ 1 , p μ 2 , p μ 3 } coincide by μ = 1 , 2 , 3 : p 1 = p 2 = p 3 = { p 1 , p 2 , p 3 } = p , then exp { i ( p μ p μ ) x } = exp { i ( p μ 0 p μ 0 ) x 0 } and ψ ¯ λ ( x ) Γ ψ λ ( x ) do not depend on the space coordinates of x = { x 0 , x } = { x 0 , x 1 , x 2 , x 3 } . For given p , m μ , the moment component p μ 0 is defined as p μ 0 = m μ 2 + p 2 , and exp { i ( p μ p μ ) x } is a periodic function of the time coordinate x 0 with period T μ μ = 2 π / | m μ 2 + p 2 m μ 2 + p 2 | . Thus, the function ψ ¯ λ ( x ) Γ ψ λ ( x ) describes an evolution of the system which is characterized by three periods, T 12 , T 13 , T 23 . It is an example of a typical process called neutrino oscillations within the flavor description of the system.
If m 1 m 2 m 3 , and m μ m μ , then
p μ 0 p μ 0 = m μ 2 + p 2 m μ 2 + p 2 = ( m μ m μ ) 1 | p | 2 2 m μ m μ + O | p | 4 m μ 2 m μ 2 ,
for small | p | 2 / ( m μ m μ ) and
p μ 0 p μ 0 = m μ 2 m μ 2 2 | p | 1 m μ 2 + m μ 2 4 | p | 2 + O m μ 2 m μ 2 | p | 4
for small ( m μ m μ ) / | p | 2 . The free field approximation of the action functional enables to describe the propagation of neutrino in vacuum.
For processes in which the influence of the material environment is significant, it was proposed to represent this in the model by an additional potential in the Hamiltonian. In this way, models with constant and adiabatically varying density of the matter were constructed and studied by Mikheev, Smirnov, and Wolfenstein [44,45,46,47,48,49,50,51]. It was shown that the effective masses of neutrino are changed by their interaction with material media. This can cause resonance effects in the processes of neutrino oscillations (MSW resonance), which significantly change their characteristics.
The problem of modeling the interaction of neutrinos with external media attracts the attention of many researchers. It remains actual at the present time. In developing the methods used in [44,45,46,47,48,49,50], many models describing the interactions of neutrino and matter with constant and adiabatically distributed density have been constructed [14,16,17,18,19,20,21].
However, little attention has been paid to the study of boundary effects and phenomena generated by the strong inhomogeneous medium, for modeling of which it is necessary to take into account the interaction of neutrinos with singular density distribution concentrated in a d < 4 -dimensional subspace of the Minkowski space-time. In this paper, we will demonstrate the possibility of applying the methods of quantum field theory to such problems.

3. Interaction of Neutrinos with Matter

The main idea of Symanzik’s approach in constructing renormalizable models of quantum field theory in a non-uniform space-time is to use the possibility of modifying the action functional of a usual renormalizable quantum field model which is invariant in respect to the space-time translations and Lorentz transformations by appending an additional so-called defect action functional (DAF) obeying some general requirements [22].
The most important of these is that the modified model must remain renormalizable. It is a formal mathematical requirement that imposes strong restrictions on the possible form of the DAF. It should naturally also be assumed that the basic physical principles of interaction laws in the original model also remain non-broken in the modified one. In the gauge theory models, these could be the basic postulate about locality and local gauge invariance. In addition to that, some common physical requirements can be taken into account. For example, the DAF does not break the unitarity of the scattering matrix.
In the framework of Symansik’s approach, one constructed the model describing the interaction of the QED fields with two-dimensional material, which form is defined by the solution of equation f ( x ) = 0 . The full action functional of that reads as
S ( ψ ¯ , ψ , A , f ) = S ( ψ ¯ , ψ , A ) + S d e f ( ψ ¯ , ψ , A , f )
where S ( ψ ¯ , ψ , A ) is the usual action functional of QED
S ( ψ ¯ , ψ , A ) = 1 4 F μ ν F μ ν + ψ ¯ ( i ^ + i e A ^ m ) ψ , F μ ν = μ A ν ν A μ ,
with an electromagnetic field A and spinor fields ψ ¯ , ψ , electron charge e, and mass m. The DAF is the sum of two terms: S d e f ( ψ ¯ , ψ , A , f ) = S f ( A , a ) + S f ( ψ ¯ , ψ , Q ) written as
S f ( A , a ) = a 2 ε α β μ ν α A β μ f ( x ) A ν δ ( f ( x ) ) d x , S f ( ψ ¯ , ψ , Q ) = ψ ¯ ( x ) Q ψ ( x ) δ ( f ( x ) ) d x .
We used here the notation ε α β μ ν for the totally antisymmetric Levi-Civita tensor ( ε 0123 = 1 ), a and the elements of the matrix Q are dimensionless parameters. The matrix Q satisfies the condition γ 0 Q γ 0 = Q . The parameter a is a real number. The delta-function δ ( f ( x ) ) describes a subspace f ( x ) = 0 of ( 3 + 1 ) -space-time filled with 2D material [23,24,33,34]. Any ( 4 × 4 ) matrix can be represented as a linear combination of 16 linearly independent matrices with complex coefficients. As such basic elements, we will use the matrices Γ k , 1 k 16 of the following form
I , γ j , γ 5 = i γ 0 γ 1 γ 2 γ 3 , γ j γ 5 , σ j k = i γ j γ k = i 2 [ γ j , γ k ] , j < k , j , k = 0 , 1 , 2 , 3
where I is the ( 4 × 4 ) identity matrix. These Γ i can be considered as matrices that form a basis for a linear (reducible) representation of the Lorentz group. The Dirac matrices γ 0 , γ 1 , γ 2 , γ 3 are transformed as components of a Lorentz contravariant vector, I is the scalar, and γ 5 is the pseudoscalar. The matrices γ ν γ 5 , σ α β are represented as contravariant components of the pseudovector and antisymmetric tensor of the second rank, respectively.
Thus, the differential matrix operator ^ = γ ν ν and the QED action functional S ( ψ ¯ , ψ , A ) are invariant in respect to Lorentz transformations, and the S d e f ( ψ ¯ , ψ , A , f ) describing the interaction of the QED fields with the extended object is breaking this symmetry. The remaining symmetry properties of the system are defined by the form of the surface f ( x ) = 0 and the choice of parameters r k of the matrix Q = k = 1 16 r k Γ k .
The action functional (3) was proposed in [23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43] as a realization of the opportunity to construct a model of the interaction of QED fields with two-dimensional materials within the framework of the Symanzik approach, unless the electron mass, S ( ψ ¯ , ψ , A , f ) does not contain other dimensional parameters. This model is local, gauge-invariant, and renormalizable. For a material with a given shape (function f ( x ) ) and with given material properties (parameters r k of matrix Q and a), it is possible to investigate theoretically within the framework of the model a large class of various problems. For example, it can be used for calculating the characteristics of scattering processes and bound states of particles. In this case, using conventional methods of QED and corresponding to their specific problem modifications, one can obtain quantitative results with a high degree of accuracy that are suitable for experimental verification and various predictions.
In our work, we propose to generalize this approach to the case of interaction of neutrino fields with matter whose distribution of density would be a local function concentrated in a subspace with dimension d < 4 of the four-dimensional Minkowski space-time.
Since, by definition, the action functionals are dimensionless and the dimension of the product of two spinor fields is equal to three, for the DAF, it is necessary by d < 3 to have parameters with negative dimensions. Adding such a DAF to the basic action of the model will violate the renormalizability of that. Therefore, the only valid value for d < 4 is d = 3 . Therefore, a possible generalization of the QED functional S d e f ( ψ ¯ , ψ , A , f ) for a description of interaction of neutrino fields with a singularly distributed medium could be proposed for mass representation, as
S d e f ( Φ ¯ , Φ , L , Q , f ) = Φ ¯ ( x ) L Q Φ ( x ) δ ( f ( x ) ) d x = μ , μ = 1 3 L μ μ ϕ ¯ μ ( x ) Q ϕ μ ( x ) δ ( f ( x ) ) d x
Here, the elements of a hermitian ( 3 × 3 ) matrix L and ( 4 × 4 ) matrix Q are constant dimensionless parameters. The matrix Q is supposed to be presented as Q = j = 1 16 r j Γ j with 16 complex numbers r j and linear independent matrices Γ j of the form (4). The solution of equation f ( x ) = 0 describes a region of Minkowski space filled with the matter that interact with neutrinos. Its properties are presented by the parameters r k . In this paper, we consider as the extended material object the plane x 3 = 0 . It corresponds to choosing f ( x ) = x 3 .
We put on the matrix Q the restriction γ 0 Q γ 0 = Q , which is necessary for the scattering matrix unitarity. It follows from γ 0 γ 0 = I , γ 0 γ i γ 0 = γ i = γ i by i = 1 , 2 , 3 that γ 0 γ 5 γ 0 = γ 5 = γ 5 , γ 0 σ j k γ 0 = σ j k , γ 0 γ 5 γ j γ 0 = ( γ 5 γ j ) . Therefore, the coefficient r i by the matrix Γ i = γ 5 in the representation Q = j = 1 16 r j Γ j is imaginary, and all other coefficients r j , j i are real.
The matrix Q is simplified if there is a symmetry in the interaction of plane x 3 = 0 and spinor fields. If it is assumed that the material plane is isotropic and homogeneous, that is, the DAF (5) is invariant with respect to rotation about the x 3 —axis and to translations along x 1 , x 2 —directions, then Q has the form [34]
Q = r 1 I + i r 2 γ 5 + r 3 γ 3 + r 4 γ 5 γ 3 + r 5 γ 0 + r 6 γ 5 γ 0 + r 7 σ 03 + r 8 σ 12
where r k , k = 1 , . . . , 8 , are real numbers.
The free action functional (1) has the form
S 0 ( Φ ¯ , Φ , M ) = k = 1 n j , j = 1 4 μ , μ = 1 3 ϕ ¯ μ j ( x ) L μ μ k Q j j k ϕ μ j ( x ) d x , n = 2 .
Here, L μ μ k is the Hermitian ( 3 × 3 ) matrix and Q j j k is the Dirac ( 4 × 4 ) operator matrix. The proposed contribution in DAF of spinor fields Φ ¯ , Φ (5) is the quadratic form of the same structure with n = 1 . From this point of view, (5) can be considered as a minimal model of interaction of neutrinos with material planes. We will use it in this paper.
This model contains 17 real parameters: 9 of them define the matrix L , and there are 8 in Q. The unitary flavor transformation of fields Φ ¯ , Φ enables one to diagonalize the matrix L in (5). In this presentation, the eigenvalues of L appear to be the coupling constants of three neutrino mixes with the plane x 3 = 0 .
In order to obtain an analogue of the Chern-Simons potential describing the coupling of electromagnetic fields with 2D material for the theory of weak interactions, it is sufficient to use the fact that the strength tensor of the non-abelian Yang-Mills field A ^ μ has the form F ^ μ ν = μ A ^ ν ν A ^ μ i g [ A ^ μ , A ^ ν ] , and changes as F ^ μ ν ( x ) Θ ( x ) F ^ μ ν ( x ) Θ ( x ) 1 under the local gauge transformation with the matrix Θ ( x ) . It means that
S d e f ( A ^ , f , λ , g ) = λ 4 ε μ ν ρ σ T r ( F ^ μ ν ( x ) F ^ ρ σ ( x ) ) θ ( f ( x ) ) d x = λ ε μ ν ρ σ θ ( f ( x ) ) μ T r A ^ ν ( x ) ρ A ^ σ ( x ) 2 i g 3 A ^ ν ( x ) A ^ ρ ( x ) A ^ σ ( x ) d x ,
where θ ( f ( x ) ) is the Heaviside step function of f ( x ) and λ , g are constant parameters, is a gauge invariant functional.
If the field A ^ μ ( x ) disappears at large x, then S d e f ( A ^ , f , λ , g ) = S d e f ( A ^ , f , λ , g ) , where
S d e f ( A ^ , f , λ , g ) = λ ε μ ν ρ σ T r A ^ ν ( x ) ρ A ^ σ ( x ) 2 i g 3 A ^ ν ( x ) A ^ ρ ( x ) A ^ σ ( x ) μ f ( x ) δ ( f ( x ) ) d x .
In the framework of the perturbation theory the functional S d e f ( A ^ , f , λ , g ) is also equivalent to S d e f ( A ^ , f , λ , g ) , and the last one can be used as DAF in modeling the interaction of the Yang-Mills field A ^ μ ( x ) with 2D material concentrated on the surface f ( x ) = 0 . It can be considered as a possible generalization of the abelian Chern-Simons action functional for the non-abelian gauge vector field. If it disagrees with non-perturbative results obtained by using S d e f ( A ^ , f , λ , g ) as an alternative versions of DAFs, this situation will impose a special investigation.
It is important to note that in the theory of gauge interactions of bosonic vector and fermionic spinor fields, their interaction with 2D materials is described within the framework of the proposed approach by the sum of functionals, each of which contains only bosonic or fermionic fields. Therefore, the influence of fields of one type on the effects of interaction of fields of another type with extended 2D objects is, in the main approximation, insignificant.
We assume that this is also true for the processes of interaction of neutrinos with a strongly inhomogeneous medium, and to study their features, we will use a model with DAF (5) which contains only neutrino fields.
The invariant in respect to all not affecting the axis x 3 transformations of the Lorentz group interaction of plane x 3 = 0 with a Dirac field was considered in [33]. For symmetry of such a kind, one needs to put r 5 = r 6 = r 7 = r 8 = 0 in (6) and the matrix Q obtain the form
Q = r 1 I + i r 2 γ 5 + r 3 γ 3 + r 4 γ 5 γ 3 .
If one takes into account only the properties of the plane material which are invariant in respect to all rotations and busts, then one can put r 3 = r 4 = 0 and obtain
Q = r 1 I + i r 2 γ 5 .
This matrix depends on two real parameters r 1 , r 2 . If the parity symmetry is supposed to not be broken by the DAF, then r 2 = 0 and
Q = r 1 I .
Thus, the full action functional describing the interaction of the material plane x 3 = 0 with the system of Dirac fields Φ ¯ , Φ in the mass representation reads as
S ( Φ ¯ , Φ , M , L , Q ) = S 0 ( Φ ¯ , Φ , M ) + S d e f ( Φ ¯ , Φ , L , Q ) = Φ ¯ ( x ) ( i ^ + M + L Q δ ( x 3 ) ) Φ ( x ) d x .
Here, the fields Φ ¯ ( x ) , Φ ( x ) have three mass components ϕ ¯ μ ( x ) , ϕ μ ( x ) , and μ = 1 , 2 , 3 , and each of them has four spinor components. For notations of spinor and flavor indices, we will use the Latin and Greek letters, respectively. The matrices M , L , and Q do not depend on x coordinates. The M one is diagonal on the spinor and mass indices: M k μ k μ = m μ δ k k δ μ μ , but for Q it is so only in the special case (9). The matrix L is supposed to be Hermitian of general form.
The action functional (10) describes three free Dirac particles with masses 0 m 1 m 2 m 3 interacting on the plane x 3 = 0 . The matrix Q represents the properties that are material of this plane which are essential for its interaction with spinor fields. The diagonal part of the matrix L defines, for each particle, its interaction constant with the plane. The non-diagonal elements of L can be considered as induced by the plane coupling constants between different μ -components of the fields Φ ¯ , Φ .
In the flavor representation, the action functional (10) is written as
S ( Ψ ¯ , Ψ , M , L , Q ) = Ψ ¯ ( x ) ( i ^ + M + L Q δ ( x 3 ) ) Ψ ( x ) d x .
For Q = 0 , it is presented in (1) and
S d e f ( Ψ ¯ , Ψ , M , L , Q ) = Ψ ¯ ( x ) L Q δ ( x 3 ) Ψ ( x ) d x = λ λ = 1 3 L λ λ ψ ¯ λ ( x ) Q ψ λ ( x ) δ ( x 3 ) d x .
The matrix elements of M, L and M , L are connected by relations
M λ λ = μ , μ = 1 3 e μ λ M μ μ e μ λ * = μ = 1 3 e μ λ m μ e μ λ * , L λ λ = μ , μ = 1 3 e μ λ L μ μ e μ λ * .
Thus, we constructed the model of interaction of neutrino fields with strong inhomogeneous matter based on the Symanzik approach in quantum field theory. The DAF (5) is supposed to be used as the addition term to the action functional of renormalized models describing neutrino physics. The constructed model is an analog of the model of interaction of the QED fields with 2D matter, of which the investigation in the Gaussian approximation enabled one to obtain non-trivial theoretical results about Casimir and Casimir-Polder effects [23,24,27,31,42], scattering processes [32,40,41], and the bound state of photons and Dirac particles [37,39,43]. Within the Gaussian approximation of the proposed model, we consider the scattering of neutrinos on the material plane and analyze the influence of collisions with it on their oscillations.

4. Statement of the Problem

Although the action functional (10) is Gaussian, the processes, which it describes, are nontrivial. We will study the scattering on the plane x 3 = 0 of particles, which are presented by the fields Φ ¯ , Φ , by using the modified Dirac equations
δ δ Φ ¯ S ( Φ ¯ , Φ , M , L , Q ) = ( i ^ M + L Q δ ( x 3 ) ) Φ ( x ) = 0 ,
δ δ Φ S ( Φ ¯ , Φ , M , L , Q ) = i μ Φ ¯ ( x ) γ μ + Φ ¯ ( M + L Q δ ( x 3 ) ) = 0 ,
characterizing the point of stationarity of the functional S ( Φ ¯ , Φ , M , L , Q ) . The ordinary way to do it is to find the solution of (11) and (12) and applying that to construct the currents of incident, reflected, and transmitted particles. It enables one to calculate the characteristics of the scattering process. Such a problem was solved for the interaction of one Dirac field with the plane x 3 = 0 defined by the DAF with matrices Q of the form (6), (7). Our task is to obtain such results for the model with an action functional (10).
If Φ ¯ + ( x ) , Φ + ( x ) , and Φ ¯ ( x ) , Φ ( x ) denote solutions of (11) and (12) by x 3 > 0 and x 3 < 0 , respectively; then they must satisfy the free Dirac equations
( i ^ M ) Φ ± ( x ) = 0 , i ν Φ ¯ ± ( x ) γ ν + M Φ ¯ ± ( x ) = 0 ,
and conditions on the plane x 3 = 0 :
lim x 3 + 0 Φ + ( x ) = Λ S lim x 3 0 Φ ( x ) ,
with matrix S corresponding to the symmetry of considered interaction defined by the matrix Q, and a ( 3 × 3 ) flavor matrix Λ .
We suppose that in the scattering process, the incident and reflected particles are in the subspace x 3 < 0 , and the transmitted ones are in the region x 3 > 0 . The incident and transmitted particles move in the positive direction of the x 3 axes, and we denote by Φ ( x ) and Φ + ( x ) the describing them spinors. Reflected particles moving in the opposite direction will be represented by the spinor Φ ( x ) . Thus, the fields Φ ± ( x ) in (14) have the form: Φ ( x ) = Φ ( x ) + Φ ( x ) , Φ + ( x ) = Φ + ( x ) .
If the reflection of particles exists, then Φ ( x ) 0 , the functions Φ ± ( x ) are not continuous by x 3 = 0 , and Λ S 1 in (14). It seems to be in contradiction with (11), since (11) is correctly defined if Φ ( x ) is a continuous by x 3 = 0 function.
This problem is solved by an auxiliary regularization of the δ -function in the interaction action [33,39]. It enables one to construct a regularized version of the conditions (14), and it is possible after removing regularization in this expression, to obtain a finite limit for S in terms of the coupling constants of the plane with a Dirac field. For the matrix Q defining this interaction in the model with one spinor field, one received S = exp { i γ 3 Q } = S 1 [39]. However, in the framework of other regularization schemes, it appears to be S = ( I + i γ 3 Q / 2 ) ( I i γ 3 Q / 2 ) 1 = S 2 [33].
Thus, the matrix S which must be expressed in terms of Q and used in such an approach depends on choosing the regularization, but it is essential that both S = S 1 and S = S 2 obey the requirement
S γ 0 γ 3 S = γ 0 γ 3 .
If (14) is fulfilled, and Λ Λ = 1 , then the equality (15) ensures that
lim x 3 + 0 Φ ¯ + ( x ) γ 3 Φ + ( x ) = lim x 3 0 Φ ¯ ( x ) γ 3 Φ ( x ) ,
that is, no additional current is created on the plane x 3 = 0 along the x 3 axis.
In constructing a solution to the problem proposed by us, we assume that its determining parameters are the elements of the matrices S and Λ . In this case, it is supposed that there is a regularization procedure for the delta function in the action of the model, which makes it possible to establish a one-to-one relationship between elements of matrices S , Λ and Q , L . In this respect, S and Λ can be considered to be directly related to the observables and independent from the choice of the regularization procedure. Calculations based on the use of boundary condition (14) do not require any additional regularization scheme. Therefore, values of elements S and Λ can be expressed in terms of experimental data. Values of the elements of the matrices Q , L may depend on the choice of regularization. It can be compared with renormalization in models of quantum field theory, where the observable values are expressed in terms or renormalized parameters which are considered as independent from bare parameters of Lagrangian and used regularization.
We suppose to calculate the characteristic of the scattering process by using the boundary condition (14) and to obtain results in terms of matrices S , Λ . The problem here is that for our model, we are given matrix Q, but we do not know the matrix S independent of the choice of regularization. In this situation, it is natural to try to construct this matrix by analyzing the properties of the matrices S 1 and S 2 .
To reveal their structural features, which may also be the same for S, we introduce convenient notation that was used in [37,39]. If M is a ( 2 × 2 ) matrix with elements M i j , i , j = 1 , 2 , then M ( ± ) are the ( 4 × 4 ) matrices
M ( + ) = M 11 0 M 12 0 0 0 0 0 M 21 0 M 22 0 0 0 0 0 , M ( ) = 0 0 0 0 0 M 11 0 M 12 0 0 0 0 0 M 21 0 M 22 .
We will use the notations τ ν ( + ) , τ ν ( ) , ν = 0 , 1 , 2 , 3 for the ( 4 × 4 ) matrices corresponding to the unit ( 2 × 2 ) matrix τ 0 and Pauli matrices τ j , j = 1 , 2 , 3 . The matrix Q (6) can be written as
Q = Q + ( + ) + Q ( ) , Q ± = j = 0 3 q ± j τ j , Q ± ( ± ) = j = 0 3 q ± j τ ( ± )
where q ± 0 = r ± 18 , q ± 1 = i r ± 27 , q ± 2 = ± i r 36 , q ± 3 = r 45 , and r ± i j = r i ± r j .
In virtue of γ 3 = i τ 2 ( + ) i τ 2 ( ) , we receive
Q = i γ 3 Q = τ 2 ( + ) Q + ( + ) τ 2 ( ) Q ( ) = Q + ( + ) + Q ( ) , Q ± ( ± ) = j = 0 3 q ± j τ j ( ± ) , Q ± = ± j = 0 3 q ± j τ 2 τ j = j = 0 3 q ± j τ j , q ± 0 = ± q ± 2 , q ± 1 = i q ± 3 , q ± 2 = ± q ± 0 , q ± 3 = ± i q ± 1 .
As γ 0 = τ 3 ( + ) + τ 3 ( ) , the condition γ 0 Q γ 0 = Q is written for the ( 2 × 2 ) -matrices Q ± as τ 3 Q ± τ 3 = Q ± . This means that q ± 1 , q ± 2 are imaginary and q ± 0 , q ± 3 are real numbers, since τ j = τ j , and τ j τ k + τ k τ j = 2 δ j k τ 0 with j , k = 1 , 2 , 3 . It is fulfilled for real parameters r k , k = 1 , . . . , 8 in (6). It follows from (16) that the parameters q ± 0 , q ± 1 are imaginary and q ± 2 , q ± 3 are real.
If M and N are ( 2 × 2 ) -matrices, and f ( x ) is an analytical function at x = 0 , then the following relations are fulfilled for the matrices M ( ± ) , N ( ± )
M ( ± ) N ( ) = 0 , M ( ± ) N ( ± ) = ( M N ) ( ± ) , f ( M ( + ) + N ( ) ) = f ( M ) ( + ) + f ( N ) ( ) .
The matrix S from the boundary condition (14) can be obtained in the regularised model. It appears to be dependent on the regularisation scheme function of the matrix Q . Examples of such could be S 1 = exp { Q } and S 2 = ( I + Q / 2 ) ( I Q / 2 ) 1 .
For the description of neutrino scattering on the plane x 3 = 0 , we will use the matrix S of the form
S = S + ( + ) + S ( ) ,
S ± ( ± ) = e i η ± ς 0 ± τ 0 ( ± ) + i ς 1 ± τ 1 ( ± ) + ς 2 ± τ 2 ( ± ) + ς 3 ± τ 3 ( ± ) , ς 0 ± 2 + ς 1 ± 2 ς 2 ± 2 ς 3 ± 2 = 1 ,
where 0 η ± π and ς k ± , k = 0 , 1 , 2 , 3 are real numbers. Employing this parameterization does not generate difficulties for constructing a complete solution to the problem posed by us. On the other hand, there is no reason to expect that the appearance of results obtained in this way cannot be received within the framework of the approach using regularization (see the Appendix A). To present formulas in compact form, it is convenient to also use the following notations
u 1 + = ς 0 + + ς 3 + , u 2 + = ς 0 + ς 3 + , u 1 = ς 0 + ς 3 , u 2 = ς 0 ς 3 ,
v 1 + = ς 1 + + ς 2 + , v 2 + = ς 1 + ς 2 + , v 1 = ς 1 + ς 2 , v 2 = ς 1 ς 2 .
The condition ς 0 ± 2 + ς 1 ± 2 ς 2 ± 2 ς 3 ± 2 = 1 is written for these parameters as u 1 ± u 2 ± + v 1 ± v 2 ± = 1 .
It can also be useful to present ς j ± , j = 0 , 1 , 2 , 3 as
ς 0 ± = cosh ( r ± ) cos ( α ± ) , ς 1 ± = cosh ( r ± ) sin ( α ± ) , ς 2 ± = sinh ( r ± ) cos ( β ± ) , ς 3 ± = sinh ( r ± ) sin ( β ± ) .
Here, r ± 0 , 0 α ± 2 π , 0 β ± 2 π and there are not conditions connecting r ± , α ± , β ± with each other.

5. Scattering of Plane Waves

The solution of free Dirac equations (13) can be presented in the form
ϕ μ ( x ) = e i p x ϕ μ ( p ¯ ) d p ¯ , μ = 1 , 2 , 3 .
Here, we used the notation p ¯ = ( p 0 , p 1 , p 2 ) and p 3 to obey the condition p 2 = p 0 2 p 1 2 p 2 2 p 3 2 = m μ 2 . If p 0 m μ , then the spinor ϕ μ ( p ¯ , x 3 ) = exp ( i p 3 x 3 ) ϕ μ ( p ¯ ) in the integrand of (21) describes by p 3 > 0 the particle moving in the positive direction of x 3 -axes, and p 3 < 0 corresponds to movement in the opposite direction. The spinor ϕ μ ( p ¯ ) fulfills the Dirac equation ( p ^ m μ ) ϕ μ ( p ¯ ) = 0 . For the scattering process described by (13) and (14), the most general plane wave presentation of spinors ϕ μ ± ( p ¯ , x 3 ) can be chosen as
ϕ μ + ( p ¯ , x 3 ) = ϕ μ + ( p ¯ , x 3 ) , f o r x 3 0 , ϕ μ ( p ¯ , x 3 ) = ϕ μ ( p ¯ , x 3 ) + ϕ μ ( p ¯ , x 3 ) , f o r x 3 0 .
Here,
ϕ μ + ( p ¯ , x 3 ) = ( e i ω μ ϕ μ 1 ( p ¯ ) c μ 1 + ϕ μ 2 ( p ¯ ) c μ 2 ) e i κ μ x 3 , ϕ μ ( p ¯ , x 3 ) = ( e i ω μ ϕ μ 1 ( p ¯ ) a μ 1 + ϕ μ 2 ( p ¯ ) a μ 2 ) e i κ μ x 3 , ϕ μ ( p ¯ , x 3 ) = ( e i ω μ ϕ μ 1 ( p ¯ ) b μ 1 + ϕ μ 2 ( p ¯ ) b μ 2 ) e i κ μ x 3 ,
where κ μ = p 0 2 p 1 2 p 2 2 m μ 2 0 ,
ϕ μ 1 ( p ¯ ) = 1 0 k μ f μ e i ω μ , ϕ μ 2 = 0 1 f μ e i ω μ k μ , ϕ μ 1 ( p ¯ ) = 1 0 k μ f μ e i ω μ , ϕ μ 2 = 0 1 f μ e i ω μ k μ , k μ = κ μ p 0 + m μ , p 1 + i p 2 p 0 + m μ = f μ e i ω μ , f μ = | f μ | , 0 ω μ 2 π .
and a 1 , a 2 , b 1 , b 2 , c 1 , c 2 are arbitrary complex parameters. Functions ϕ ± ( p ¯ , x 3 ) describe the incidend and transmitted particles moving in the positive direction of x 3 -axes, and ϕ ( p ¯ , x 3 ) corresponds to reflected particles moving from the plane x 3 = 0 in the negative direction of x 3 -axes. The boundary condition (14) is written for ϕ ± ( p ¯ ) as
ϕ μ + ( p ¯ μ , 0 ) = μ = 1 3 Λ μ μ S ϕ μ ( p ¯ μ , 0 )
with matrix S presented in (17) and (18) and
μ = 1 3 Λ μ μ Λ μ μ = δ μ μ .
Thus, (22) is a system of 12 linear equations which enables one to express the amplitudes b 1 μ , b 2 μ , c 1 μ , c 2 μ , μ = 1 , 2 , 3 of reflected and transmitted particles in terms of amplitudes a 1 μ , a 2 μ of incident ones. Substituting the spinors φ μ ± ( p ¯ , 0 ) in (22) and using the notations
a μ = a 1 μ a 2 μ , b μ = b 1 μ b 2 μ , c μ = c 1 μ c 2 μ , m μ + = 1 0 k μ f μ , n μ + = 1 0 k μ f μ , , m μ = 0 1 f μ k μ , n μ = 0 1 f μ k μ
for μ = 1 , 2 , 3 , we obtain six equations of the form
m μ ± c μ = μ = 1 3 Λ μ μ s ± ( m μ ± a μ + n μ ± b μ ) .
It follows from (24) that
c μ = μ = 1 3 Λ μ μ m ¯ μ ± s ± ( m μ ± a μ + n μ ± b μ ) , μ = 1 3 Λ μ μ n ¯ μ ± s ¯ ± m μ ± c μ = n ¯ μ ± m μ ± a μ + b μ ,
where we used the notations s ± = e i η ± ς 0 ± τ 0 + i ς 1 ± τ 1 + ς 2 τ 2 + ς 3 τ 3 ,
m ¯ μ ± = m μ ± 1 , n ¯ μ ± = n μ ± 1 , s ¯ ± = s ± 1 = e i η ± ς 0 ± τ 0 i ς 1 ± τ 1 ς 2 τ 2 ς 3 τ 3 .
Excluding b μ , c μ in (25), one receives the equations
μ = 1 3 Ω μ μ c μ = Θ μ a μ , μ = 1 3 Ω μ μ b μ = μ = 1 3 Ω μ μ a μ
with ( 2 × 2 )-matrices
Ω μ μ = Λ μ μ ( n ¯ μ + s ¯ + m μ + n ¯ μ s ¯ m μ ) , Ω μ μ = Λ μ μ ( m ¯ μ + s + n μ + m ¯ μ s n μ ) Ω μ μ = Λ μ μ ( m ¯ μ s n μ m ¯ μ + s + n μ + ) , Θ μ = n ¯ μ + m μ + n ¯ μ m μ .
The solution for b and c can be written as
c μ = μ = 1 3 Ω ˜ μ μ Θ μ a μ , b μ = μ , μ = 1 3 Ω ˜ μ μ Ω μ μ a μ ,
and the problem is to construct the matrices Ω ˜ μ μ , Ω ˜ μ μ in an explicit form.
It can be solved in the more general formulation. Let K k k be ( n × n )-matrices, V k , U k - n-component vectors, k , k = 1 , 2 , 3 , and
k = 1 3 K k k V k = U k .
One needs to construct the solution of the system of equations (27)
V k = k = 1 3 K ˜ k k U k
with an explicit form of the ( n × n ) -matrices K ˜ k k , k , k = 1 , 2 , 3 .
Finding the components V 1 , V 2 , V 3 of the vector V directly from Equation (27), one obtains them in the form
V i = H i 1 ( U i ( K i j K i l K l l 1 K l j ) ( K j j K j l K l l 1 K l j ) 1 U j ( K i l K i j K j j 1 K j l ) ( K l l K l j K j j 1 K j l ) 1 U l ) , H i = K i i K i j ( K j j K j l K l l 1 K l j ) 1 ( K j i K j l K l l 1 K l i ) K i l ( K l l K l j K j j 1 K j l ) 1 ( K l i K l j K j j 1 K j i ) .
Here, the sets { i j l } of the indexes are assumed to be chosen as { i j l } = { 123 } , { 231 } , { 312 } . Comparing (28) and (29), we receive the following expressions for the matrices K ˜ j , j
K ˜ 11 = H 1 1 , K ˜ 12 = H 1 1 ( K 12 K 13 K 33 1 K 32 ) ( K 22 K 23 K 33 1 K 32 ) 1 , K ˜ 13 = H 1 1 ( K 13 K 12 K 22 1 K 23 ) ( K 33 K 32 K 22 1 K 23 ) 1 , K ˜ 22 = H 2 1 , K ˜ 21 = H 2 1 ( K 21 K 23 K 33 1 K 31 ) ( K 11 K 13 K 33 1 K 31 ) 1 , K ˜ 23 = H 2 1 ( K 23 K 31 K 11 1 K 13 ) ( K 33 K 31 K 11 1 K 13 ) 1 , K ˜ 33 = H 3 1 , K ˜ 31 = H 3 1 ( K 31 K 32 K 22 1 K 21 ) ( K 11 K 12 K 22 1 K 21 ) 1 , K ˜ 32 = H 3 1 ( K 32 K 31 K 11 1 K 12 ) ( K 22 K 21 K 11 1 K 12 ) 1 .
Using Ω μ μ , Ω μ μ and (26)–(29), one can calculate the matrices Ω ˜ μ μ , Ω ˜ μ μ and obtain the right-hand sides in the representations (26) of b μ , c μ in an evident form. However, for the unitary ( 3 × 3 ) matrix Λ of general form, they turn out to be rather cumbersome and inconvenient for analyzing their properties. Therefore, to study the most simple effects of the neutrino interaction with planes, we restrict ourselves in this paper to the case of a diagonal matrix, Λ .

6. Explicit Results in a Simplified Model

In virtue of Λ Λ = 1 , a diagonal part of Λ has the form Λ = diag { exp ( i ρ 1 ) , exp ( i ρ 2 ) , exp ( i ρ 3 ) } with real parameters ρ 1 , ρ 2 , ρ 3 . For Λ of such a kind, the matrices Ω μ μ , Ω μ μ , Ω μ μ are diagonal, and (26) is written as
b μ = e i ρ μ ( m ¯ μ s n μ m ¯ μ + s + n μ + ) 1 ( m ¯ μ + s + m μ + m ¯ μ s m μ ) a μ ,
c μ = e i ρ μ ( n ¯ μ + s ¯ + m μ + n ¯ μ s ¯ m μ ) 1 ( n ¯ μ + m μ + n ¯ μ m μ ) a μ .
Using the notations (19) and (20), one can present (30) and (31) as
b μ = R ( μ ) a μ , c μ = T ( μ ) a μ , R ( μ ) = e i ρ μ d μ r 11 ( μ ) r 12 ( μ ) r 21 ( μ ) r 22 ( μ ) , T ( μ ) = e i ( ρ μ + η + ) d μ t 11 ( μ ) t 12 ( μ ) t 21 ( μ ) t 22 ( μ )
with
d μ = f μ 2 ( 2 cos ( ξ ) u 2 + u 2 u 1 + u 1 + v 1 v 2 + + v 1 + v 2 ) ( ( k μ 2 + f μ 2 ) v 2 i k μ ( u 2 + u 1 ) + v 1 ) ( ( k μ 2 + f μ 2 ) v 2 + + i k μ ( u 2 + + u 1 + ) + v 1 + ) , r 11 ( μ ) = f μ 4 v 2 + v 2 ( k μ 2 v 2 + i k μ ( u 1 + u 2 + ) v 1 + ) ( k μ 2 v 2 i k μ ( u 1 + u 2 ) + v 1 ) + + f μ 2 ( u 1 + u 1 + u 2 + u 2 + i k μ ( ( u 1 + + u 2 + ) v 2 + ( u 1 u 2 ) v 2 + ) 2 cos ( ξ ) , r 12 ( μ ) = 2 f μ k μ ( e i ξ u 2 + u 2 ( f μ 2 + k μ 2 ) v 2 + v 2 i k μ ( u 2 + v 2 v 2 + u 2 ) ) , r 21 ( μ ) = 2 f μ k μ ( u 2 u 2 + + ( f μ 2 + k μ 2 ) v 2 + v 2 + i k μ ( u 2 + v 2 v 2 + u 2 ) e i ξ ) , r 22 ( μ ) = f μ 4 v 2 + v 2 ( k μ 2 v 2 + + i k μ ( u 2 + + u 1 + ) + v 1 + ) ( k μ 2 v 2 + i k μ ( u 1 u 2 ) v 1 ) + + f μ 2 ( u 1 + u 1 + u 2 + u 2 i k μ ( ( u 1 + u 2 + ) v 2 + ( u 1 + u 2 ) v 2 + ) 2 cos ( ξ ) ) ,
and
t 11 ( μ ) = 2 i k μ ( f μ 2 v 2 + e i ξ k μ 2 v 2 + i k μ ( u 1 + u 2 ) v 1 ) , t 12 ( μ ) = 2 f μ k μ ( ( u 1 + i k μ v 2 + ) e i ξ u 2 i k μ v 2 ) , t 21 ( μ ) = 2 f μ k μ ( ( u 2 + i k μ v 2 + ) e i ξ u 1 i k μ v 2 ) , t 22 ( μ ) = 2 i k μ ( ( k μ 2 v 2 + + i k μ ( u 1 + + u 2 + ) + v 1 + ) e i ξ f μ 2 v 2 )
where ξ = η + η .
The currents J i n ν , J r ν , J t ν of incident, reflected, and transmitted particles are of the form
J i n μ ν = ϕ ¯ i n μ γ ν ϕ i n μ = p i n μ ν 2 a μ a μ * m μ + p μ 0 , J r μ ν = ϕ ¯ r μ γ ν ϕ r μ = p r μ ν 2 b μ b μ * m μ + p μ 0 , J t μ ν = ϕ ¯ t μ γ ν ϕ t μ = p t μ ν 2 c μ c μ * m μ + p μ 0 .
Here, p i n μ ν , p r μ ν , p t μ ν are components of the momentum vectors p i n μ = p t μ = { p μ 0 , p μ 1 , p μ 2 , p μ 3 } , p r μ = { p μ 0 , p μ 1 , p μ 2 , p μ 3 } , p r μ 2 = p r μ 2 = p μ 0 2 p μ 1 2 p μ 2 2 p μ 3 2 = m μ 2 . In the used parametrization, p μ 1 = ( m μ + p μ 0 ) f μ cos ( ω μ ) , p μ 2 = ( m μ + p μ 0 ) f μ sin ( ω μ ) , p μ 3 = ( m μ + p μ 0 ) k μ .
In virtue of (23) and ϕ ¯ i n μ γ 3 ϕ r μ + ϕ μ ¯ r μ γ 3 ϕ i n μ = 0 ,
J t μ 3 = lim x 3 + 0 ϕ ¯ μ + ( x ) γ 3 ϕ μ + ( x ) = lim x 3 0 ϕ ¯ μ ( x ) γ 3 ϕ μ ( x ) = J i n μ 3 + J r μ 3 .
Thus, the continuity of the current component J μ 3 at the plane x 3 = 0 implies that J i n μ 3 = J t μ 3 J r μ 3 . Taking into account the direction of the current J r μ 3 , we come to the conclusion that the equality
a μ * a μ = b μ b μ * + c μ c μ * = a μ * R ( μ ) R ( μ ) a μ + a μ * T ( μ ) T ( μ ) a μ
must be fulfilled. To verify (32), it is sufficient to use the above given expressions for the matrices R ( μ ) , T ( μ ) and to take into account the relations u 1 ± u 2 ± + v 1 ± v 2 ± = 1 for the parameters on which they depend. It follows from (32) that if one denotes the reflection and transmission coefficients for the considered scattering process of the μ -th particle as
K r μ = b μ b μ * a μ a μ * = a μ * R ( μ ) R ( μ ) a μ a μ a μ * , K t μ = c μ c μ * a μ a μ * = a μ * T ( μ ) T ( μ ) a a μ a μ * ,
then K r μ + K t μ = 1 . The matrices R ( μ ) R ( μ ) , T ( μ ) T ( μ ) are Hermitian. The elements T μ i j = { T ( μ ) T ( μ ) } i j , i , j = 1 , 2 , i j of the matrix T ( μ ) T ( μ ) can be presented as
T μ 11 = 4 k μ 2 ( k μ 2 ( 2 + u 2 2 u 2 + 2 ) + v 1 2 + ( f μ 2 + k μ 2 ) ( u 2 + 2 + u 1 2 + k μ 2 v 2 2 + f μ 2 v 2 + 2 ) 2 f μ 2 ( ( u 2 + u 1 + v 2 + v 1 ) cos ( ξ ) + k μ ( u 2 + v 2 v 2 + u 2 ) sin ( ξ ) ) , T μ 12 = 4 i f μ k μ 2 ( u 2 + ( i k μ u 2 + + v 1 + ) u 2 ( i k μ u 2 + v 1 ) + + ( ( u 1 + v 1 u 1 v 1 + ) i k μ ( u 2 + u 1 u 2 u 1 + + v 2 + v 1 v 2 v 1 + ) ) e i ξ ) + + ( f μ 2 + k μ 2 ) ( v 2 ( u 1 i k μ v 2 ) v 2 + ( u 1 + i k μ v 2 + ) ) + ( u 2 v 2 + u 2 + v 2 ) ( f μ 2 e i ξ + k μ 2 e i ξ ) , T μ 22 = 4 k μ 2 ( k μ 2 ( 2 + u 2 + 2 u 2 2 ) + v 1 + 2 + ( f μ 2 + k μ 2 ) ( u 1 + 2 + u 2 2 + f μ 2 v 2 2 + k μ 2 v 2 + 2 ) 2 f μ 2 ( ( u 1 + u 2 + v 1 + v 2 ) cos ( ξ ) + k μ ( u 2 + v 2 u 2 v 2 + ) sin ( ξ ) ) ) ,
where ξ = η + η . The elements R μ i j = R μ j i * of the matrix R ( μ ) R ( μ ) for i j , i , j = 1 , 2 are the following
R μ 11 = d μ * d μ T μ 11 , R μ 12 = T μ 12 , R μ 22 = d μ * d μ T μ 22 .
This is a consequence of the relation (32).
We see that in the mass representation, the obtained results for characteristics of the scattering process do not depend on parameters ρ i of diagonal matrix Λ . However, they influence the oscillation events of transmitted and reflected neutrinos in the flavor description of the system, taking the contribution to the oscillating part of (2) in the form
μ μ = 1 3 c μ λ * c μ λ e i ( p μ p μ ) x χ ¯ μ ( p μ ) Γ χ μ ( p μ ) .
It can essentially change the characteristics of neutrino oscillations after their collision with the material plane.
In the model with diagonal matrix Λ and matrix S of the form (17), we constructed a plane wave solution of Equations (13) and (14), describing the motion of particles with an arbitrary angle of incidence to the plane x 3 = 0 . Now, we consider a special process of neutrino collision with a plane with an angle of incidence equal to zero.

7. Moving of Particles along the x 3 -Axes

If the particles move orthogonally to the plane’s x 3 = 0 direction, then p i n μ = p t μ = { p μ 0 , 0 , 0 , p 3 } , p r μ = { p μ 0 , 0 , 0 , p 3 } . This scattering process is described by reflection and transition matrices
R ( μ ) = R ( μ ) | f = 0 = e i ρ μ r μ + 0 0 r μ , T ( μ ) = T ( μ ) | f = 0 = e i ( ρ μ + η + ) t μ + 0 0 t μ
where
r μ + = r + ( k μ ) = k μ 2 v 2 + i k μ ( u 1 + u 2 + ) v 1 + k μ 2 v 2 + + i k μ ( u 1 + + u 2 + ) + v 1 + , r μ = r ( k μ ) = k μ 2 v 2 + i k μ ( u 1 u 2 ) v 1 k μ 2 v 2 i k μ ( u 1 + u 2 ) + v 1 , t μ + = t + ( k μ ) = 2 i k μ k μ 2 v 2 + + i k μ ( u 1 + + u 2 + ) + v 1 + , t μ = t ( k μ ) = 2 i k μ e i ( η η + ) i k μ ( u 1 + u 2 ) k μ 2 v 2 v 1 , k μ = p μ 0 m μ p μ 0 + m μ = p 3 m μ + p 3 2 + m μ 2 = ρ 1 + ρ 2 + 1 | ρ = p 3 m μ .
The matrices R ( μ ) R ( μ ) , T ( μ ) T ( μ ) are also diagonal. Their elements
| r μ ± | 2 = ( v 1 ± k μ 2 v 2 ± ) 2 + k μ 2 ( u 1 ± u 2 ± ) 2 ( v 1 ± + k μ 2 v 2 ± ) 2 + k μ 2 ( u 1 ± + u 2 ± ) 2 , | t μ ± | 2 = 4 k μ 2 ( v 1 ± + k μ 2 v 2 ± ) 2 + k μ 2 ( u 1 ± + u 2 ± ) 2
fulfill the relations | r μ ± | 2 + | t μ ± | 2 = 1 . The reflection and transition coefficient could be written as
K r μ = K r ( k μ ) = | r μ + | 2 cos ( θ μ ) 2 + | r μ | 2 sin ( θ μ ) 2 , K t μ = K t ( k μ ) = | t μ + | 2 cos ( θ μ ) 2 + | t μ | 2 sin ( θ μ ) 2 .
Here, the angle θ μ is defined by relation cos ( θ μ ) 2 = a μ 1 a μ 1 * / ( a μ 1 a μ 1 * + a μ 2 a μ 2 * ) .
Let us denote
f ( q ) = f ( u , w 1 , w 2 ; q ) = 4 u + w 1 q + q w 2 ,
where u > 0 , w 1 0 , w 2 0 are constant parameters determining this function, which will be used to analyze the scattering process under consideration. It follows directly from definition that f ( q ) > 0 . For q 0 , this function reaches its maximal value f m a x by q = q m a x = w 1 / w 2 :
f m a x = f ( q m a x ) = f w 1 w 2 = 4 u + 2 w 1 w 2 .
The equation f ( q ) = f m a x / 2 has, by w 1 w 2 0 and 0 < q , two solutions, q = q , q + , 0 < q < q m a x < q + :
q ± = 8 f m a x u ± ( 8 f m a x u ) 2 4 f m a x 2 w 1 w 2 f m a x w 2 = u + 4 w 1 w 2 ± ( u + 4 w 1 w 2 ) 2 4 w 1 w 2 2 w 2 .
The neighbourhood q < q < q + of the point q max can be considered as characterizing the function f ( q ) region. We denote an estimation of its extension as f e x t = q + q . Then,
f e x t 2 = u + 2 w 1 w 2 w 2 .
The parameters u , w 1 , w 2 can be expressed in terms of q m a x , f m a x , f e x t as
u = 4 ( f e x t 2 2 q m a x ) f e x t 2 f m a x , w 1 = 4 q m a x 2 f e x t 2 f m a x , w 2 = 4 f e x t 2 f m a x .
Substituting them in f ( u , w 1 , w 2 ; q ) , we receive the following presentation of this function
f ( q ) = f m a x f e x t 2 q f e x t 2 q + ( q m a x q ) 2 = f m a x g q m a x f m a x 2 ; q q m a x , g ( c ; x ) = 1 1 + c ( x x 1 ) 2 .
Let us put u = 2 + u 1 2 + u 2 2 , w 1 = v 1 2 , w 2 = v 2 2 , where u 1 , u 2 , v 1 , v 2 are real constants which can be parameterized by 0 r , 0 α < 2 π , 0 β < 2 π as follows:
u 1 = cosh ( r ) sin ( α ) + sinh ( r ) cos ( β ) , u 2 = cosh ( r ) sin ( α ) sinh ( r ) cos ( β ) , v 1 = cosh ( r ) cos ( α ) + sinh ( r ) sin ( β ) , v 2 = cosh ( r ) cos ( α ) sinh ( r ) sin ( β ) .
Then, u 1 u 2 + v 1 v 2 = 1 and
f ( u , w 1 , w 2 ; k 2 ) = 4 2 + u 1 2 + u 2 2 + v 1 2 k 2 + k 2 v 2 2 = 4 ( u 1 + u 2 ) 2 + ( v 1 k + k v 2 ) 2 = 4 4 + ( u 1 u 2 ) 2 + ( v 1 k k v 2 ) 2 = h ( u 1 , u 2 , v 1 , v 2 ; k ) = h ( k ) = h ˜ ( r , α , β ; k ) = h ˜ ( k ) .
The function h ˜ ( k ) has the form h ˜ ( k ) = f m a x g ( k m a x 2 / f e x t 2 ; k / k m a x ) , where
f m a x = 2 2 c a 2 b 2 ( a 2 b 2 ) 2 , f e x t = 2 ( 2 c a 2 b 2 ( a 2 b 2 ) 2 ) ( a b ) 2 , k m a x 2 = ( a + b ) 2 ( a 2 b 2 ) 2
and a = cosh ( r ) sin ( α ) , b = sinh ( r ) cos ( β ) , c = cosh r 2 . Since this relation between parameters f m a x , f e x t , k m a x , and r , α , β is symmetric in respect to replacement a b , there are two possibilities of inverting it, as
cosh ( r ) 2 = f e x t 2 + ( k m a x 2 1 ) 2 f e x t 2 f m a x , sin ( α ) cos ( β ) = k m a x 4 1 f e x t 2 + ( k m a x 2 1 ) 2 f e x t 2 ( 1 f m a x ) + ( k m a x 1 ) 2 and ( 1 ) sin ( α ) 2 = ( k m a x 2 1 ) 2 f e x t 2 + ( k m a x 2 1 ) 2 , cos ( β ) 2 = ( k m a x 2 + 1 ) 2 f e x t 2 ( 1 f m a x ) + ( k m a x 2 1 ) 2 ; ( 2 ) sin ( α ) 2 = ( k m a x 2 + 1 ) 2 f e x t 2 + ( k m a x 2 1 ) 2 , cos ( β ) 2 = ( k m a x 2 1 ) 2 f e x t 2 ( 1 f m a x ) + ( k m a x 2 1 ) 2 .
Since sin ( α ) 2 1 , cos ( β ) 2 1 , the relation 2) imposes the restriction f e x t 2 k m a x on f e x t , and in case 1), the inequality f e s t 2 k m a x / 1 f m a x must be fulfilled. Thus, the function h ˜ ( k ) can be presented in the form h ˜ ( k ) = c 1 g ( c 2 ; k / c 3 ) with 0 < c 1 < 1 , 0 < c 2 1 / 4 , 0 < c 3 . It is even: h ˜ ( k ) = h ˜ ( k ) has the maximum value h ˜ m a x = c 1 by | k | = c 3 = k m a x , and the extension of the neighborhood of k m a x is h e x t = c 3 / c 2 . By given c 3 , the function h ˜ ( k ) has the minimal value of h e x t by c 2 = 1 / 4 . If c 2 = 1 / 4 , then h ˜ ( k ) = h m a x g ( 1 / 4 , k / c 3 ) .
The graphs of the functions g ( 1 / 4 , k / c ) by different c = k m a x are shown on Figure 1. With an increase of c, the vicinity of the maxima becomes increasingly more flat, and the positions of the maxima become increasingly less noticeable. Using the notation
h ( u 1 ± , u 2 ± , v 1 ± , v 2 ± ; k ) = h ± ( k ) = h ˜ ( r ± , α ± , β ± ; k ) = h ˜ ± ( k ) = c 1 ± g ( c 2 ± ; k / c 3 ± ) ,
one can present the functions | t μ ± | 2 of k μ in the form | t μ ± | 2 = h ± ( k μ ) = g ± ( k μ ) = c 1 ± g ( c 2 ± ; k μ / c 3 ± ) and obtain the following expression for the transmission coefficient:
K t ( k ) = c 1 + g ( c 2 + ; k / c 3 + ) cos ( θ ) 2 + c 1 g ( c 2 ; k / c 3 ) sin ( θ ) 2 .
Let us consider, as an example, a scattering process for which θ = 0 and u 1 + = 0.5 , u 2 + = 0.3 , v 1 + = 0.01 , v 2 + = 85 . In this case, α + = 0.00941038 , β + = 1.56844 , r + = 4.442675 and
K t ( k ) = h + ( k ) = 4 4 + 0 . 2 2 + ( 0.01 k 85 k ) 2 = c 1 + 1 + c 2 + c 3 + k k c 3 + 2 = c 1 + g ( c 2 + ; k / c 3 + ) ,
c 1 + = 4 4.04 0.990099 , c 2 + = 0.85 4.04 0.210396 , c 3 + = 1 10 85 0.0108465 .
The maximal value of K t ( k ) is h m a x + = c 1 + = h + ( k m ) , k m a x = c 3 + , h e x t + = c 3 + c 2 + = 4.04 85 0.0236468 . The parameters c 1 , c 2 , c 3 (34) are very close to ones defining the function presented by Graphic 1 on the Figure 1. The difference between it and the graph of the function K t ( k ) (33) is inessential. Note that c 2 + (34) fulfill the restrictions c 2 + 1 / 4 and h e x t = 4.04 / 85 0.0236 > 2 k m a x 0.0216 . However, 2 k m a x / 1 f m a k 0.218 and the inequality h e x t 1 f m a k k m a x is not satisfied. Therefore, there is only one { r , α , β } -parametrization h ˜ + ( r , α , β ; k ) of the function h + ( u 1 , u 2 , v 1 , v 2 ; k ) in this case.
Outside the neighborhood, 0.00162 < k < 0.0726 ( 0.00324 < ρ < 0.146 ) of the point k m a x = 0.0108 ( ρ m a x = 0.0217 ) , h + ( k ) (33) is 10 times less than its maximum value. For 0.1 < k < 1 , it decreases monotonically from h + ( 0.1 ) = 0.0536193 to h + ( 1 ) = 0.000553454 and has the maximal value h + ( m a x ) at the point 0 < c 3 + < 1 (34).
The graph of h + ( k ) (33) differs significantly from the (graphs (2)–(6)) shown in Figure 2 for K t ( k ) by cos ( θ ) 2 1 , 0 . These ones can have two maxima (graphs (2) and (3)), one smoothed maximum (graph (6)), or change very little on the most part of the interval 0 k 1 (graphs (4) and (5)).
In the general case, the process of neutrino scattering on the plane is described by one universal function g ( k ) = g ( α , β , k ) parameterized by constants 0 α < 0.25 , β > 0 : K t ( k ) = g ( α , β , k ) γ + g ( α + , β + , k ) γ + with factors γ ± > 0 . In the terms of parameters c 1 ± , c 2 ± , c 3 ± , θ used above, α ± = c 2 ± , β ± = c 3 ± , γ + = c 1 + cos ( θ ) 2 , γ = c 1 sin ( θ ) 2 . The maximal value of g ( k ) is 1 = g ( α , β ; β ) . The argument k of this function is k = p 0 m / p 0 + m , where m is the mass of the particle and p 0 its energy. Hence, 0 k < 1 by the physical values m p 0 < of p 0 . If 0 < β 1 , then g ( α , β ; k ) has the maximum on the interval 0 < k 1 . If β > 1 , then g ( α , β ; k ) grows from g ( α , β ; 0 ) = 0 to g ( α , β ; β ) = α with an increase of k by 0 < k 1 .
Since, in the considered scattering process m, p 3 and k fulfill the relations k = p 0 m / p 0 + m , p 0 2 p 3 2 = m 2 , one can present p 3 , p 0 as p 0 = ( 1 + k 2 ) m / ( 1 k 2 ) , p 3 = 2 k m / ( 1 k 2 ) . Hence, if g ( α , β ; k ) = 1 , then k = β and p 0 = ( 1 + β 2 ) m / ( 1 β 2 ) , p 3 = 2 β m / ( 1 β 2 ) .
Corresponding to (33) function K t ( p ) , the transmission coefficient as functions of p 3 has the maximum by p 3 = 2 c 3 + m / ( 1 c 3 + 2 ) = 20 85 m / 8499 0.0216956 m, and as the function of p 0 , it is maximal by p 0 = ( 1 + c 3 + 2 ) m / ( 1 c 3 + 2 ) = 8501 m / 8499 1.00024 m .
Now, we analyze the features of system dynamics in the framework of the flavor representation. As the matrices Γ used in (2) when discussing the oscillation process of neutrinos, we consider Γ = I , γ 0 , γ 3 ; γ 5 , γ 0 γ 5 , γ 3 γ 5 . The function χ ¯ μ ( p μ ) χ μ ( p μ ) ( χ ¯ μ ( p μ ) γ 5 χ μ ( p μ ) ) is interpreted as the scalar (pseudo-scalar) density of the μ - th particles with momentum p μ = { p μ 0 , 0 , 0 , p 3 } , χ ¯ μ ( p μ ) { γ 0 , 0 , 0 , γ 3 } χ μ ( p μ ) ( χ ¯ μ ( p μ ) { γ 0 γ 5 , 0 , 0 , γ 3 γ 5 } χ μ ( p μ ) ) is the four-current vector (axial four-current pseudo-vector) of the μ - th particles moving along the x 3 -axis. We will use the plane wave presentation for spinors χ ¯ μ ( p μ ) , χ μ ( p μ ) describing the particles with momentum p 3 moving orthogonal to the plane x 3 = 0 :
χ μ = e i ρ μ ( σ 1 μ α 1 μ + σ 2 μ α 2 μ ) , χ ¯ μ = e i ρ μ ( σ 1 μ α 1 μ * + σ 2 μ α 2 μ * ) γ 0 , σ 1 μ = 1 0 k μ 0 , σ 2 μ = 0 1 0 k μ , k μ = p 3 p μ 0 + m μ , p μ 0 = ( p 3 ) 2 + m μ 2 , k μ = 1 2 p 3 m μ 1 8 p 3 m μ 3 + O p 3 m μ 4 , k μ = 1 m μ p 3 + 1 2 m μ p 3 2 + O m μ p 3 4 .
By means of notation a μ μ ± = α 1 μ * α 1 μ ± α 2 μ * α 2 μ , the results obtained for the χ ¯ μ ( p μ ) Γ χ μ ( p μ ) are written in the form
χ ¯ μ χ μ = 1 k μ k μ e i ( ρ μ ρ μ ) a μ μ + , χ ¯ μ γ 0 χ μ = 1 + k μ k μ e i ( ρ μ ρ μ ) a μ μ + , χ ¯ μ γ 3 χ μ = k μ + k μ e i ( ρ μ ρ μ ) a μ μ +
and
χ ¯ μ γ 5 χ μ = k μ k μ e i ( ρ μ ρ μ ) a μ μ , χ ¯ μ γ 0 γ 5 χ μ = k μ + k μ e i ( ρ μ ρ μ ) a μ μ , χ ¯ μ γ 3 γ 5 χ μ = 1 + k μ k μ e i ( ρ μ ρ μ ) a μ μ .
Oscillations of incident waves are presented by ψ ¯ λ ( x ) Γ ψ λ ( x ) = μ , μ 3 e μ λ * e μ λ ϕ ¯ μ ( x ) Γ ϕ μ ( x ) as
ψ ¯ λ ( x ) ψ λ ( x ) = μ = 1 3 e μ λ * e μ λ 1 k μ 2 a μ μ + + μ μ = 1 3 e i ( p μ 0 p μ 0 ) x 0 e μ λ * e μ λ 1 k μ k μ e i ( ρ μ ρ μ ) a μ μ + , ψ ¯ λ ( x ) γ 0 ψ λ ( x ) = μ = 1 3 e μ λ * e μ λ 1 + k μ 2 a μ μ + + μ μ = 1 3 e i ( p μ 0 p μ 0 ) x 0 e μ λ * e μ λ 1 + k μ k μ e i ( ρ μ ρ μ ) a μ μ + , ψ ¯ λ ( x ) γ 3 ψ λ ( x ) = 2 μ = 1 3 e μ λ * e μ λ k μ a μ μ + + μ μ = 1 3 e i ( p μ 0 p μ 0 ) x 0 e μ λ * e μ λ k μ + k μ e i ( ρ μ ρ μ ) a μ μ + .
and
ψ ¯ λ ( x ) γ 5 ψ λ ( x ) = μ μ = 1 3 e i ( p μ 0 p μ 0 ) x 0 e μ λ * e μ λ k μ k μ e i ( ρ μ ρ μ ) a μ μ , ψ ¯ λ ( x ) γ 0 γ 5 ψ λ ( x ) = 2 μ = 1 3 e μ λ * e μ λ k μ a μ μ + μ μ = 1 3 e i ( p μ 0 p μ 0 ) x 0 e μ λ * e μ λ k μ + k μ e i ( ρ μ ρ μ ) a μ μ , ψ ¯ λ ( x ) γ 3 γ 5 ψ λ ( x ) = μ = 1 3 e μ λ * e μ λ 1 + k μ 2 a μ μ + μ μ = 1 3 e i ( p μ 0 p μ 0 ) x 0 e μ λ * e μ λ 1 + k μ k μ e i ( ρ μ ρ μ ) a μ μ .
For transmitted waves, one needs to replace a 1 μ c 1 μ = t μ + a 1 μ , a 2 μ c 2 μ = t μ a 2 μ . Thus, changing in (35) and (36)
a μ μ ± c μ μ ± = c 1 μ * c 1 μ ± c 2 μ * c 2 μ = ( t μ + * t μ + ) a 1 μ * a 1 μ ± ( t μ * t μ ) a 2 μ * a 2 μ ,
one receives the corresponding results for transmitted neutrino oscillations. If m 1 < m 2 < m 3 and momentum p 3 is the same for all particles, then p 3 / m 1 > p 3 / m 2 > p 3 / m 3 . For given p 3 , we denote k μ = p 3 / ( m μ + ( p 3 ) 2 + m μ 2 ) , μ = 1 , 2 , 3 . If | t ± | 2 ( k ) monotonically increases from zero to | t ± | 2 ( 1 ) by 0 k 1 , then | t ± | 2 ( k 1 ) > | t ± | 2 ( k 2 ) > | t ± | 2 ( k 3 ) , since k 1 > k 2 > k 3 . In virtue of
| t ± | 2 ( k 1 ) | t ± | 2 ( k 2 ) = ( k 1 k 2 ) k | t ± | 2 ( k ) | k = k , k 1 > k > k 2 , | t ± | 2 ( k 2 ) | t ± | 2 ( k 3 ) = ( k 2 k 3 ) k | t ± | 2 ( k ) | k = k , k 2 > k > k 3 ,
the greater the derivatives of the function | t ± | 2 ( k ) are at points k and k for given k 1 , k 2 , k 3 , that is, the faster | t ± | 2 ( k ) grows on the interval k 1 > k > k 3 , the more there will be an influence of the neutrino collisions with the plane on the process of their oscillations.
If there is a maximum of | t ± | 2 ( k ) by 0 < k < 1 , then each | t ± | 2 ( k μ ) , μ = 1 , 2 , 3 can be maximal, and the following situations are possible:
(1)
| t ± | 2 ( k 1 ) is maximal, | t ± | 2 ( k 3 ) < | t ± | 2 ( k 2 ) < | t ± | 2 ( k 1 ) ;
(2)
| t ± | 2 ( k 2 ) is maximal, | t ± | 2 ( k 3 ) < | t ± | 2 ( k 1 ) < | t ± | 2 ( k 2 )
or | t ± | 2 ( k 1 ) < | t ± | 2 ( k 3 ) < | t ± | 2 ( k 2 ) ;
(3)
| t ± | 2 ( k 3 ) is maximal, | t ± | 2 ( k 1 ) < | t ± | 2 ( k 2 ) < | t ± | 2 ( k 3 ) .
The more the function | t ± | 2 ( k ) resembles an approximation of the delta-function and the nearer the maximal | t ± | 2 ( k μ ) is to the maximum of | t ± | 2 ( k ) , the greater is the difference of maximal | t ± | 2 ( k μ ) from non-maximal | t ± | 2 ( k μ ) , and the more significant will be the changes in neutrino oscillations as a result of their collisions with the plane.
We considered the model with 8 real parameters in the matrix Q (6). It could be interesting to compare its properties with features of the system with matrices Q of the form (7)–(9). In the case (7), there are four real parameters defining Q and
S ± ( ± ) = e i η ( ς 0 τ 0 ( ± ) + i ς 1 τ 1 ( ± ) ± ς 2 τ 2 ( ± ) ± ς 3 τ 3 ( ± ) ) ,
u 1 + = u 2 = ς 0 + ς 3 = u 1 , u 2 + = u 1 = ς 0 ς 3 = u 2 ,
v 1 + = v 2 = ς 1 + ς 2 = v 1 , v 2 + = v 1 = ς 1 ς 2 = v 2 ,
ς 0 = cosh ( r ) cos ( α ) , ς 1 = cosh ( r ) sin ( α ) , ς 2 = sinh ( r ) cos ( β ) , ς 3 = sinh ( r ) sin ( β ) .
Here, ς 0 , ς 1 , ς 2 , ς 3 , u 1 , u 2 , v 1 , v 2 are real parameters fulfilled the conditions ς 0 2 + ς 1 2 ς 2 2 ς 3 2 = 1 , u 1 u 2 + v 1 v 2 = 1 and r 0 , 0 α 2 π , 0 β 2 π , 0 η π . Hence, for such materials,
t μ + = t μ * = 2 i k μ k μ 2 v 2 + i k μ ( u 1 + u 2 ) + v 1 = t μ , | t μ | 2 = 4 4 + ( u 1 u 2 ) 2 + ( v 1 k μ k μ v 2 ) 2 = K t μ , | r μ + | 2 = | r μ | 2 = | r μ | 2 = K r μ = 1 K t μ , c μ μ ± = ( t μ * t μ ) a 1 μ * a 1 μ ± ( t μ t μ * ) a 2 μ * a 2 μ .
Putting η = 0 , ς 1 = 0 , α = 0 in (37)–(40), one obtains for the model given by matrix Q (8) with two real parameters
S ± ( ± ) = ς 0 τ 0 ( ± ) ± ς 2 τ 2 ( ± ) ± ς 3 τ 3 ( ± ) , ς 0 = cosh ( r ) , ς 2 = sinh ( r ) cos ( β ) , ς 3 = sinh ( r ) sin ( β )
u 1 = ς 0 + ς 3 , u 2 = ς 0 ς 3 , v 1 = ς 2 = v , v 2 = ς 2 = v ,
with real ς 0 , ς 2 , ς 3 , u 1 , u 2 , v fulfilling the conditions ς 0 2 ς 2 2 ς 3 2 = 1 , u 1 u 2 v 2 = 1 , and r 0 , 0 β 2 π . In this model,
t μ + = t μ * = 2 i k μ i k μ ( u 1 + u 2 ) + v ( 1 k μ 2 ) = t μ , | t μ | 2 = 4 4 + ( u 1 u 2 ) 2 + v 2 ( 1 k μ + k μ ) 2 = K t μ .
If one sets in (41) and (42) ς 3 = 0 , β = 0 , then one receives, for the model with a matrix Q (9) containing one real parameter,
S ± ( ± ) = ς 0 τ 0 ( ± ) ± ς 2 τ 2 ( ± ) , ς 0 = cosh ( r ) , ς 2 = sinh ( r ) , u 1 = u 2 = ς 0 = u , v 1 = v 2 = ς 2 = v .
Here, ς 0 , ς 2 , u , v are real parameters, ς 0 2 ς 2 2 = 1 , u 2 v 2 = 1 and r 0 . In this case,
t μ = 2 i k μ 2 i u k μ + v ( 1 k μ 2 ) , | t μ | 2 = 4 4 + v 2 ( 1 k μ + k μ ) 2 = K t μ .
The essential difference between the models with matrix Q of the form (6)–(9) is that the functions | t μ ± | 2 ( k ) , | t μ | 2 ( k ) , in the models (6), (7) can have the maximum by 0 < k < 1 , but | t μ | 2 ( k ) in models (8), (9) are, by 0 k 1 , the monotonously growing functions.

8. Conclusions

In our work, we considered the problem of neutrino interaction with matter. Using our experience in constructing models of QED in singular background fields, we have proposed a quantum-field approach, which may be useful for the theoretical description of neutrino propagation in a highly inhomogeneous medium. It assumes taking into account the basic symmetry principles of modern physics of fundamental interactions that underlie the Standard Model and can, in principle, be generalized to describe the interaction of all lepton fields with the external environment. Mainly, attention was paid to the problem of neutrino scattering on a material plane, considered as a simplest example of process in the space with a strongly inhomogeneous distribution of substance.
In a general form, for a model with an off-diagonal unitary matrix Λ mixing Dirac fields in the mass representation, expressions for the reflected and transmitted waves were obtained. For them, in the model with a diagonal Λ , an explicit solution was obtained, which was used to analyze the oscillations of neutrinos in the case of their motion, orthogonally to the plane x 3 = 0 . It was shown that the parameters that determine the material properties of the plane can be chosen so that its effect on the neutrino flux is similar to a filter that transmits particles in a narrow interval of low energies and almost completely reflects all other ones. As a result of the neutrino collision with the plane, the parallel component of the momentum does not change, and the orthogonal one does not change in absolute value. Only the amplitudes of fields can change significantly.
The example we have considered with parameters of the model (33) and (34) shows that the interaction of neutrinos with a plane can effect their filtration. A characteristic feature of the filtration process of particles upon collision with a plane is the possibility of essentially different transmission coefficients for neutrinos of different masses. In this case, due to filtration of their flux, the regimes of the neutrino oscillations before the plane and behind it can be strongly different. This phenomenon can be used to estimate the masses of neutrinos of various types in carrying out analyses of experimental data.
Although the filtering and MSW-resonance results are similar in many ways, their mechanics are not the same. The MSW effect is formed non-locally in space and time. This requires a certain volume of matter and a certain period of time, generally speaking, that are different for various substances. In order for the filtration of the neutrino flux to occur, their collision with the plane is sufficient, which (in the framework of the considered model) occurs instantaneously and locally at x 3 = 0 . From the point of view of the possibility of verifying the adequacy of the proposed model, it would be interesting to reveal in the dynamics of neutrino oscillations an effect, which cannot arise as a consequence of MSW resonance and is a manifestation of the filtration process.
One of the current theoretical problems in neutrino astro-physics is constructing numerical models of dynamics of supernova explosions. Many research teams have been working on this issue [52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67]. Perhaps, taking into account the filtering mechanism in such models will be useful to achieve a better understanding of the features of the process of collapse of the super-heavy star core.
In general terms, a possible scenario of its evolution can be presented as follows. If, in the core of the star, its shell filters neutrinos by energies, they can be divided into two classes. Particles with energies from the narrow range of the low-energy region leave the core unhindered. These neutrinos can be called free. For all the others, which we will call bound, the core shell is impermeable.
In the process of the star’s evolution, its core is subjected to the pressure of the increasing gravitational forces. In it, neutrinos are born, the free ones are emitted, and the bound neutrinos are accumulated in the core. This can go on until the main features of the interaction between the core shell and neutrinos changes, the class of free neutrinos expands, and the star will emit them with further contraction without a significant change in its structure. In our model, this can be described by changing the function K t ( k ) . For instance, if the plots of possible K t ( k ) , is shown in Figure 2, then within the change ( 1 ) ( 4 ) , a large fraction of the bound high-energy neutrinos become free, they will leave the star, essentially changing the intensity and spectrum of its neutrino emission.
If the main interaction features of neutrinos with the core shell are not changed and at least some of the bound neutrinos do not become free, then the enormous energy accumulated by them in the core destroys its shell sooner or later. After that, the core and the star are exploded.
There is great interest among experimentalists to determine neutrino masses directly [68,69,70,71,72]. One can assume that the employment of 2D materials and special surface treatment techniques in the construction of neutrino detectors would enable one to efficiently use the filtering mechanism in experiments of such a kind.
The maximal number of parameters in the model, the properties of which we have studied in detail, was eight for the matrix Q and three for the diagonal matrix Λ . Not all of them are included in our results, and the question arises whether it is possible to reduce the number of parameters in the model without limiting its area of applicability. We considered versions of the model with four, two, and one parameters in the Q matrix, simplified for symmetry reasons, and found a difference in the properties of their predicted transmission coefficients. This raises the question of whether one can confine oneself to using the simplest models to describe real neutrino oscillation data.
We suppose that the proposed method for modeling the processes of interaction of neutrinos with matter can be useful for theoretical studies and analysis of the obtained experimental data.

Author Contributions

All authors contributed equally to this paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Regularization and Parameters of Model

If the matrix S in (14) is a function of Q = Q + ( + ) + Q ( ) , then S ( Q ) = S + ( + ) + S ( ) where S ± ( ± ) = S ( Q ± ( ± ) ) , and the corresponding ( 2 × 2 ) -matrices S ± can be presented as
S ± = e i η ± ( ς 0 ± τ 0 + i ς 1 ± τ 1 + ς 2 ± τ 2 + ς 3 ± τ 3 )
with real parameters ς k ± , k = 0 , 1 , 2 , 3 , 0 η ± π and ς 0 ± 2 + ς 1 ± 2 ς 2 ± 2 ς 3 ± 2 = 1 .
Indeed, for the ( 2 × 2 ) matrices Q ± we can use the parametrisation of the form S ± = k = 0 3 s ± k τ k . In virtue of γ 0 γ 3 = τ 1 ( + ) τ 1 ( ) , the condition (15) is written for S ± as S ± τ 1 S ± = τ 1 and S ± τ 1 = τ 1 S ± 1 . Since S ± 1 = ( s ± 0 τ 0 s ± τ ) / det S ± , one obtaines
S ± τ 1 = k = 0 3 s ± * k τ k τ 1 = τ 1 ( s ± 0 τ 0 s ± τ ) det S ± , det S ± = s ± 0 2 s ± 1 2 s ± 2 2 s ± 3 2 .
The solution of (A2) for s ± * k has the form
s ± * k = s ± k det S ± , k = 0 , 2 , 3 , s ± * 1 = s ± 1 det S ± .
In virtue of | s ± * k | = | s ± k | , it follows from (A2) and (A3) that | det S ± | = 1 , det S ± = e 2 i η ± and s ± k = e i η ± ς k ± , k = 0 , 2 , 3 , s ± 1 = e i η ± i ς 1 ± with 0 η ± π , real ς k ± , k = 0 , 1 , 2 , 3 and ς 0 ± 2 + ς 1 ± 2 ς 2 ± 2 ς 3 ± 2 = 1 .
Thus, the representation (A1) for S ± is proven, and the matrices S ± ( ± ) are written as
S ± ( ± ) = e i η ± ς 0 ± τ 0 ( ± ) + i ς 1 ± τ 1 ( ± ) + ς 2 ± τ 2 ( ± ) + ς 3 ± τ 3 ( ± ) , ς 0 ± 2 + ς 1 ± 2 ς 2 ± 2 ς 3 ± 2 = 1 .
Here, 0 η ± π and ς μ ± , μ = 0 , 1 , 2 , 3 are real numbers which can be expressed in terms of parameters r j , 1 j 8 , of the matrix Q (6) in accordance with chosen regularisation.
If the parameters h 0 = h 0 * , h 1 = h 1 * , h 2 = h 2 * , h 3 = h 3 * define the matrix h = k = 0 3 h i τ i = h 0 τ 0 + h τ , then
h τ 1 = ( h 0 τ 0 h 1 τ 1 + h 2 τ 2 + h 3 τ 3 ) τ 1 = τ 1 ( h 0 τ 0 + h 1 τ 1 + h 2 τ 2 + h 3 τ 3 ) = τ 1 h ,
and for the function F ( h ) one obtains F ( h ) τ 1 = τ 1 F ( h ) . Therefore F ( h ) τ 1 F ( h ) = τ 1 , if F ( h ) = F ( h ) and F ( h ) F ( h ) = τ 0 .
The parameters q ± i of the matrices Q ± (16) are imaginary by i = 0 , 1 and real by i = 2 , 3 . For the functions F ( Q ± ) = ( τ 0 Q ± / 2 ) ( τ 0 + Q ± / 2 ) 1 and F ( Q ± ) = exp ( Q ± ) , the conditions F ( Q ± ) = F ( Q ± ) and F ( Q ± ) F ( Q ± ) = τ 0 are fulfilled. Hence, the expressions
S = S + ( + ) + S ( ) = ( I + i γ 3 Q / 2 ) ( I i γ 3 Q / 2 ) 1 , S = S + ( + ) + S ( ) = exp ( i γ 3 Q )
obtained for S ± ( ± ) by means of the regularization procedures [33,39] are in agrement with (15) and (A4).
The solutions of equations exp { Q ± } = S ± with S ± given in the form (A1) can be written as
Q ± = i ( η ± + 2 n π ) τ 0 + 2 ( i ς 1 ± τ 1 + ς 2 ± τ 2 + ς 3 ± τ 3 ) ln 1 + z ± 2 + z ± z ± , z ± = ς 2 ± 2 + ς 3 ± 2 ς 1 ± 2
where n is an arbitrary integer number.
The parameters z ± fulfill the restriction 1 + z ± 2 = 1 ς 1 ± 2 + ς 2 ± 2 + ς 3 ± 2 = ς 0 ± 2 0 . Hense, if z ± is real then 1 + z ± 2 + z ± > 0 and ln 1 + z ± 2 + z ± / z ± is real. If z ± is imaginary, then
| 1 + z ± 2 + z ± | 2 = 1 | z ± | 2 + | z ± | 2 = 1 .
Therefore, in this case, ln 1 + z ± 2 + z ± is imaginary and ln 1 + z ± 2 + z ± / z ± is real.
The solutions of equations ( τ 0 Q ± / 2 ) ( τ 0 + Q ± / 2 ) 1 = S ± in respect to Q ± are of the form
Q ± = 2 ( S ± τ 0 ) ( τ 0 + S ± ) 1 = 2 ( i sin ( η ± ) τ 0 + i ς 1 ± τ 1 + ς 2 ± τ 2 + ς 3 ± τ 3 ) cos ( η ± ) + ς 0 ±
For the matrix Q ± written as (16), the coefficients q 0 ± , q 1 ± are imaginary and q 2 ± , q 3 ± are real. The matrices Q ± of the form (A5), (A6) fulfill this condition. The calculation of the characteristics of the neutrino scattering on the plane x 3 = 0 can be carried out using the regularization of the modified Dirac equations. In this case, the expression of the results in terms of the parameters of the Lagrangian turns out to depend on the choice of the regularization scheme. The same results are obtained if we use the solutions of the Dirac equations for free particles and the boundary condition of the form (14) with the appropriately, in accordance with the used regularization scheme chosen matrix S.
Since regularization is not required to solve the boundary problem (14) with arbitrary S defined as (A4), the results turn out to be expressed directly in terms of the parameters of the matrix S. They can be considered as the only solutions of the posed problem that have a physical meaning. This is similar to what happens in quantum field theory, when the results of calculations are expressed in terms of the independent of the regularization scheme renormalized parameters.

References

  1. Pontecorvo, B. Mesonium and antimesonium. Sov. Phys. JETP 1957, 6, 429–431. [Google Scholar]
  2. Pontecorvo, B. Inverse β processes and nonconservation of lepton charge. Sov. Phys. JETP 1958, 7, 172–173. [Google Scholar]
  3. Gribov, V.; Pontecorvo, B. Neutrino astronomy and lepton charge. Phys. Lett. B 1969, 28, 493–496. [Google Scholar] [CrossRef]
  4. Bilenky, S.M.; Pontecorvo, B. Again on neutrino oscillations. Lett. Nuovo C. 1976, 17, 569–574. [Google Scholar] [CrossRef]
  5. Bilenky, S.M.; Pontecorvo, B. Lepton mixing and neutrino oscillations. Phys. Rep. 1978, 41, 225–261. [Google Scholar] [CrossRef]
  6. Bilenky, S.M.; Petcov, S.T. Massive Neutrinos and Neutrino Oscillations. Rev. Mod. Phys. 1987, 59, 671–754. [Google Scholar] [CrossRef]
  7. Bilenky, S.M.; Giunti, C.; Grimus, W. Phenomenology of neutrino oscillations. Prog. Part Nucl. Phys. 1999, 43, 1–86. [Google Scholar] [CrossRef] [Green Version]
  8. Fisher, P.; Kayser, B.; McFarland, K.S. Neutrino mass and oscillation. Ann. Rev. Nucl. Part Sci. 1999, 49, 481–528. [Google Scholar] [CrossRef] [Green Version]
  9. Bilenky, S.M. Neutrino masses, mixing and oscillations. Proc. Roy. Soc. London A 2004, 460, 403–444. [Google Scholar] [CrossRef]
  10. Bilenky, S. Neutrino oscillations: From a historical perspective to the present status. Nucl. Phys. B 2016, 908, 2–13. [Google Scholar] [CrossRef] [Green Version]
  11. Lobanov, A.E. Oscillations of particles in the Standard Model. Theor. Math. Phys. 2017, 192, 1000–1015. [Google Scholar] [CrossRef]
  12. Bilenky, S. Introduction to the physics of massive and mixed neutrinos. Lect. Notes Phys. 2018, 947, 1–273. [Google Scholar]
  13. Kudenko, Y.G. Neutrino oscillations: Recent results and future prospects. Phys.-Uspekhi 2018, 61, 739–747. [Google Scholar] [CrossRef]
  14. Lobanov, A.E. Particle quantum states with indefinite mass and neutrino oscillations. Ann. Phys. 2019, 403, 82–105. [Google Scholar] [CrossRef] [Green Version]
  15. Naumov, D.V.; Naumov, V.A. Quantum-field theoretical approach to neutrino oscillations. Phys. Part Nuclei 2020, 51, 5–209. [Google Scholar] [CrossRef]
  16. Giunti, C.; Studenikin, A. Neutrino electromagnetic interactions: A window to new physics. Rev. Mod. Phys 2015, 87, 531–591. [Google Scholar] [CrossRef] [Green Version]
  17. Kartavtsev, A.; Raffelt, G.; Vogel, H. Neutrino propagation in media: Flavor, helicity, and pair correlations. Phys. Rev. D 2015, 91, 125020. [Google Scholar] [CrossRef] [Green Version]
  18. Chukhnova, A.V.; Lobanov, A.E. Neutrino evolution in dense matter and electromagnetic field. EPJ Web Conf. 2019, 222, 03001. [Google Scholar] [CrossRef]
  19. Smirnov, A.Y.; Xu, X.J. Wolfenstein potentials for neutrinos induced by ultra-light mediators. J. High Energy Phys. 2019, 1912, 046. [Google Scholar] [CrossRef] [Green Version]
  20. Chukhnova, A.V.; Lobanov, A.E. Neutrino flavor oscillations and spin rotation in matter and electromagnetic field. Phys. Rev. D 2020, 101, 013003. [Google Scholar] [CrossRef] [Green Version]
  21. Agarwalla, S.K.; Das, S.; Masud, M.; Swaina, P. Evolution of Neutrino Mass-Mixing Parameters in Matter with Non-Standard Interactions. arXiv 2021, arXiv:2103.13431. [Google Scholar]
  22. Symanzik, K. Schrödinger representation and casimir effect in renormalizable quantum field theory. Nucl. Phys. B 1981, 190, 1–44. [Google Scholar] [CrossRef]
  23. Markov, V.N.; Pis’mak, Y.M. Casimir effect for thin films from imperfect materials. arXiv 2005, arXiv:hep-th/0505218v3. [Google Scholar]
  24. Markov, V.N.; Pis’mak, Y.M. Casimir effect for thin films in QED. J. Phys. A Math. Gen. 2006, 39, 6525. [Google Scholar] [CrossRef]
  25. Fialkovsky, I.V.; Markov, V.N.; Pis’mak, Y.M. Renormalizable mean field calculation in QED with fermion background. J. Phys. A Math. Gen. 2006, 39, 6357. [Google Scholar] [CrossRef]
  26. Fialkovsky, I.V.; Markov, V.N.; Pis’mak, Y.M. Field of homogeneous plane in quantum electrodynamics. Int. J. Mod. Phys. A 2006, 21, 2601–2616. [Google Scholar] [CrossRef]
  27. Marachevsky, V.N.; Pis’mak, Y.M. Casimir-Polder effect for a plane with Chern-Simons interaction. Phys. Rev. D 2010, 81, 065005. [Google Scholar] [CrossRef] [Green Version]
  28. Pismak, D.Y.; Pismak, Y.M. Scattering of electromagnetic waves on a planar surface in a model with the Chern–Simons potential. Theor.Math. Phys. 2011, 169, 1423–1431. [Google Scholar] [CrossRef]
  29. Pismak, D.Y.; Pismak, Y.M. Modeling the interaction of an electromagnetic field with a two-dimensional surface in the framework of Symanzik’s approach. Theor. Math. Phys. 2013, 175, 442–454. [Google Scholar] [CrossRef]
  30. Pis’mak, D.Y.; Pis’mak, Y.M. Chern-Simon potential in models of interaction of electromagnetic field with thin films. Phys. Part Nuclei. 2013, 4, 450–461. [Google Scholar] [CrossRef]
  31. Pis’mak, D.Y.; Pis’mak, Y.M. Chern-Simons potential in models of Casimir effect. AIP Conf. Proc. 2014, 1606, 337. [Google Scholar]
  32. Pis’mak, D.Y.; Pis’mak, Y.M.; Wegner, F.J. Electromagnetic waves in a model with Chern-Simons potential. Phys. Rev. E 2015, 92, 013204. [Google Scholar] [CrossRef] [Green Version]
  33. Pismak, D.Y.; Pismak, Y.M. Modeling the interaction of a material plane with a spinor field in the framework of Symanzik’s approach. Theor. Math. Phys. 2015, 184, 1329–1341. [Google Scholar] [CrossRef]
  34. Pismak, Y.M.; Shukhobodskaia, D.Y. Bound states in a model of interaction of Dirac field with material plane. EPJ Web Conf. 2016, 125, 05022. [Google Scholar] [CrossRef]
  35. Pismak, Y.M.; Shukhobodskaia, D.Y. Model of Dirac field interacting with material plane within Symanzik’s approach. EPJ Web Conf. 2016, 126, 05012. [Google Scholar] [CrossRef] [Green Version]
  36. Pismak, Y.M.; Shukhobodskaia, D.Y. Symanzik approach in modeling of bound states of Dirac particle in singular background. EPJ Web Conf. 2017, 158, 07005. [Google Scholar] [CrossRef] [Green Version]
  37. Pismak, Y.M.; Wegner, F.J. Dispersion relations and dynamic characteristics of bound states in the model of a Dirac field interacting with a material plane. EPJ Web Conf. 2018, 191, 06015. [Google Scholar] [CrossRef]
  38. Pismak, Y.M. Modelling of bound states of Dirac particles in singular background in framework of Symanzik approach. Phys. Part Nucl. Lett. 2018, 15, 380–383. [Google Scholar] [CrossRef]
  39. Wegner, F.J.; Pismak, Y.M. Wave packets of bound states of a Dirac field on a material plane. Theor. Math. Phys. 2019, 200, 553–565. [Google Scholar] [CrossRef]
  40. Pismak, Y.M.; Shakhova, O.Y. Symanzik approach in modeling the interaction of quantum fields with extended objects: Scattering of Dirac particles on material plane. Phys. Part Nucl. Lett. 2019, 16, 441–444. [Google Scholar] [CrossRef]
  41. Pismak, Y.M. Zero-range potential in the model of the spinor field interacting with a material plane: Conditions for transparency of the plane to Dirac particles. EPJ Web Conf. 2019, 222, 03026. [Google Scholar] [CrossRef]
  42. Pismak, Y.M.; Shakhova, O.Y. Interaction of Massless Dirac Field with 2D-Material in a Model Constructed in the Framework of Symanzik Approach. Phys. Part Nucl. 2020, 51, 578–582. [Google Scholar] [CrossRef]
  43. Pismak, Y.M.; Shukhobodskaia, D.Y. Bound state in a model of electromagnetic field interacting with a material plane. Mod. Phys. Lett. A 2020, 35, 2050170. [Google Scholar] [CrossRef]
  44. Wolfenstein, L. Neutrino oscillations in matter. Phys. Rev. D 1978, 17, 2369–2374. [Google Scholar] [CrossRef]
  45. Wolfenstein, L. Neutrino oscillations and stellar collapse. Phys. Rev. D 1979, 20, 2634–2635. [Google Scholar] [CrossRef]
  46. Mikheev, 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]
  47. Mikheev, S.P.; Smirnov, A.Y. Neutrino oscillations in a variable density medium and neutrino bursts due to the gravitational collapse of stars. Sov. Phys. JETP 1986, 64, 4–7. [Google Scholar]
  48. Mikheev, S.P.; Smirnov, A.Y. Resonant amplification of neutrino oscillations in matter and solar neutrino spectroscopy. Nuovo Cim. C 1986, 9, 17–26. [Google Scholar] [CrossRef]
  49. Mikheev, S.P.; Smirnov, A.Y. Neutrino Oscillations in an Inhomogeneous Medium: Adiabatic Regime. Sov. Phys. JETP 1987, 65, 230–236. [Google Scholar]
  50. Mikheyev, S.P.; Smirnov, A.Y. Resonant neutrino oscillations in matter. Prog. Part Nucl. Phys. 1989, 23, 41–136. [Google Scholar] [CrossRef]
  51. Ioannisian, A.N.; Kazarian, N.A.; Smirnov, A.Y.; Wyler, D. A precise analytical description of the Earth matter effect on oscillations of low energy neutrinos. Phys. Rev. D 2002, 65, 064023. [Google Scholar] [CrossRef] [Green Version]
  52. Burbidge, E.M.; Burbidge, G.R.; Fowler, W.A.; Hoyle, F. Synthesis of the Elements in Stars. Rev. Mod. Phys. 1957, 29, 547–650. [Google Scholar] [CrossRef] [Green Version]
  53. Colgate, S.A.; White, R.H. The hydrodynamic behavior of supernovae explosions. Astrophys. J. 1966, 143, 626–681. [Google Scholar] [CrossRef]
  54. Colgate, S.A. Hot bubbles drive explosions. Nature 1989, 341, 489–490. [Google Scholar] [CrossRef]
  55. Burrows, A.; Livne, E.; Dessart, L.; Ott, C.D.; Murphy, J. Features of the Acoustic Mechanism of Core-Collapse Supernova Explosions. Astrophys. J. 2007, 655, 416–433. [Google Scholar] [CrossRef]
  56. Strumia, A.; Vissani, F. Neutrino masses and mixings and…. arXiv 2010, arXiv:hep-ph/0606054v3. [Google Scholar]
  57. Melson, T.; Janka, H.-T.; Marek, A. Neutrino-driven supernova of a low-mass iron-core progenitor boosted by threedimensional turbulent convection. Astrophys. J. Lett. 2015, 801, L24. [Google Scholar] [CrossRef] [Green Version]
  58. Lentz, E.J.; Bruenn, S.W.; Hix, W.R.; Mezzacappa, A.; Bronson Messer, O.E.; Endeve, E.; Blondin, J.M.; Harris, J.A.; Marronetti, P.; Yakunin, K.N. Three-dimensional core-collapse supernova simulated using a 15 m progenitor. Astrophys. J. Lett. 2015, 807, L31. [Google Scholar] [CrossRef] [Green Version]
  59. Gibney, E. How to blow up a star. Nature 2018, 556, 287–289. [Google Scholar] [CrossRef]
  60. Stapleford, C.J.; Fröhlich, C.; Kneller, J.P. Coupling neutrino oscillations and simulations of core-collapse supernovae. Phys. Rev. D 2020, 102, 081301. [Google Scholar] [CrossRef]
  61. Johns, J.L.; Nagakura, H.; Fuller, G.M.; Burrows, A. Neutrino oscillations in supernovae: Angular moments and fast instabilities. Phys. Rev. D 2020, 101, 043009. [Google Scholar] [CrossRef] [Green Version]
  62. Azari, M.D.; Yamada, S.; Morinaga, T.; Nagakura, H.; Furusawa, S.; Harada, A.; Okawa, H.; Iwakami, W.; Sumiyoshi, K. Fast collective neutrino oscillations inside the neutrino sphere in core-collapse supernovae. Phys. Rev. D 2020, 101, 023018. [Google Scholar] [CrossRef] [Green Version]
  63. Nagakura, H.; Johns, L. Constructing angular distributions of neutrinos in core collapse supernova from zero-th and first moments calibrated by full Boltzmann neutrino transport. Phys. Rev. D 2021, 103, 123025. [Google Scholar] [CrossRef]
  64. Abbar, S.; Capozzi, F.; Glas, R.; Janka, H.-T.; Tamborra, I. On the characteristics of fast neutrino flavor instabilities in three-dimensional core-collapse supernova models. Phys. Rev. D 2021, 103, 063033. [Google Scholar] [CrossRef]
  65. Bizouard, M.-A.; Maturana-Russel, P.; Torres-Forné, A.; Obergaulinger, M.; Cerdá-Durán, P.; Christensen, N.; Font, J.A.; Meyer, R. Inference of proto-neutron star properties from gravitational-wave data in core-collapse supernovae. Phys. Rev. D 2021, 103, 063006. [Google Scholar] [CrossRef]
  66. Román-Garza, J.; Bavera, S.S.; Fragos, T.; Zapartas, E.; Misra, D.; Andrews, J.; Coughlin, S.; Dotter, A.; Kovlakas, K.; Serra, J.G.; et al. The role of core-collapse physics in the observability of black-hole neutron-star mergers as multi-messenger sources. Astrophys. J. Lett. 2021, 912, L23. [Google Scholar] [CrossRef]
  67. Portilla, M.L.; Di Palma, I.; Drago, M.; Cerda-Duran, P.; Ricci, F. Deep learning for multimessenger core-collapse supernova detection. Phys. Rev. D 2021, 103, 063011. [Google Scholar]
  68. Robertson, R.G.H.; Knapp, D.A. Direct measurements of neutrino mass. Annu. Rev. Nucl. Part Sci. 1988, 38, 185. [Google Scholar] [CrossRef]
  69. Holzschuh, E. Measurement of the neutrino mass from tritium beta decay. Rept. Prog. Phys. 1992, 55, 1035. [Google Scholar] [CrossRef]
  70. Otten, E.; Weinheimer, C. Neutrino mass limit from tritium beta decay. Rept. Prog. Phys. 2008, 71, 086201. [Google Scholar] [CrossRef] [Green Version]
  71. Drexlin, G.; Hannen, S.; Mertens, V.; Weinheimer, C. Current direct neutrino mass experiments. Adv. High Energy Phys. 2013, 2013, 293986. [Google Scholar] [CrossRef] [Green Version]
  72. Formaggio, J.A.; de Gouvêa, A.L.C.; Robertson, R.G.H. Direct measurements of neutrino mass. Phys. Rep. 2021, 914, 1–54. [Google Scholar] [CrossRef]
Figure 1. The graphs of function g ( 1 / 4 , k / c ) by 0 k 1 and different values of parameter c:
(1) c = 0.01 ; (2) c = 0.05 ; (3) c = 0.25 ; (4) c = 0.5 ; (5) c = 0.9 ; (6) c = 2.5 .
Figure 1. The graphs of function g ( 1 / 4 , k / c ) by 0 k 1 and different values of parameter c:
(1) c = 0.01 ; (2) c = 0.05 ; (3) c = 0.25 ; (4) c = 0.5 ; (5) c = 0.9 ; (6) c = 2.5 .
Symmetry 13 02180 g001
Figure 2. Transmission coefficient K t ( k ) = c 1 + g ( c 2 + ; k / c 3 + ) cos ( θ ) 2 + c 1 g ( c 2 ; k / c 3 ) sin ( θ ) 2 by different values of c 1 ± , c 2 ± , c 3 ± , θ : (1) c 1 + = 0.99 , c 2 + = 0.225 , c 3 + = 0.01 , c 1 = 0.8 , c 2 = 0.025 , c 3 = 0.1 , cos ( θ ) 2 = 0.95 ; (2) c 1 + = 0.95 , c 2 + = 0.225 , c 3 + = 0.07 , c 1 = 0.9 , c 2 = 0.25 , c 3 = 0.5 , cos ( θ ) 2 = 0.55 ; (3) c 1 + = 0.8 , c 2 + = 0.25 , c 3 + = 0.02 , c 1 = 0.9 , c 2 = 0.2 , c 3 = 0.6 , cos ( θ ) 2 = 0.35 ; (4) c 1 + = 0.9 , c 2 + = 0.025 , c 3 + = 0.5 , c 1 = 0.8 , c 2 = 0.0025 , c 3 = 0.7 , cos ( θ ) 2 = 0.9 ; (5) c 1 + = 0.8 , c 2 + = 0.025 , c 3 + = 0.7 , c 1 = 0.7 , c 2 = 0.0025 , c 3 = 0.9 , cos ( θ ) 2 = 0.9 ; (6) c 1 + = 0.6 , c 2 + = 0.25 , c 3 + = 0.7 , c 1 = 0.8 , c 2 = 0.25 , c 3 = 0.8 , cos ( θ ) 2 = 0.7 .
Figure 2. Transmission coefficient K t ( k ) = c 1 + g ( c 2 + ; k / c 3 + ) cos ( θ ) 2 + c 1 g ( c 2 ; k / c 3 ) sin ( θ ) 2 by different values of c 1 ± , c 2 ± , c 3 ± , θ : (1) c 1 + = 0.99 , c 2 + = 0.225 , c 3 + = 0.01 , c 1 = 0.8 , c 2 = 0.025 , c 3 = 0.1 , cos ( θ ) 2 = 0.95 ; (2) c 1 + = 0.95 , c 2 + = 0.225 , c 3 + = 0.07 , c 1 = 0.9 , c 2 = 0.25 , c 3 = 0.5 , cos ( θ ) 2 = 0.55 ; (3) c 1 + = 0.8 , c 2 + = 0.25 , c 3 + = 0.02 , c 1 = 0.9 , c 2 = 0.2 , c 3 = 0.6 , cos ( θ ) 2 = 0.35 ; (4) c 1 + = 0.9 , c 2 + = 0.025 , c 3 + = 0.5 , c 1 = 0.8 , c 2 = 0.0025 , c 3 = 0.7 , cos ( θ ) 2 = 0.9 ; (5) c 1 + = 0.8 , c 2 + = 0.025 , c 3 + = 0.7 , c 1 = 0.7 , c 2 = 0.0025 , c 3 = 0.9 , cos ( θ ) 2 = 0.9 ; (6) c 1 + = 0.6 , c 2 + = 0.25 , c 3 + = 0.7 , c 1 = 0.8 , c 2 = 0.25 , c 3 = 0.8 , cos ( θ ) 2 = 0.7 .
Symmetry 13 02180 g002
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pismak, Y.; Shakhova, O. Neutrino Oscillations in the Model of Interaction of Spinor Fields with Zero-Range Potential Concentrated on a Plane. Symmetry 2021, 13, 2180. https://doi.org/10.3390/sym13112180

AMA Style

Pismak Y, Shakhova O. Neutrino Oscillations in the Model of Interaction of Spinor Fields with Zero-Range Potential Concentrated on a Plane. Symmetry. 2021; 13(11):2180. https://doi.org/10.3390/sym13112180

Chicago/Turabian Style

Pismak, Yury, and Olga Shakhova. 2021. "Neutrino Oscillations in the Model of Interaction of Spinor Fields with Zero-Range Potential Concentrated on a Plane" Symmetry 13, no. 11: 2180. https://doi.org/10.3390/sym13112180

APA Style

Pismak, Y., & Shakhova, O. (2021). Neutrino Oscillations in the Model of Interaction of Spinor Fields with Zero-Range Potential Concentrated on a Plane. Symmetry, 13(11), 2180. https://doi.org/10.3390/sym13112180

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