Relativistic scattering of fermions in quaternionic quantum mechanics

In this article, we propose a quaternionic version of the Dirac equation in the presence of scalar and vector potentials. It has been shown that in complex limit of such an equation, the complex version of this equation can be covered. After setting a quaternionic form for the Dirac delta potential, scattering due to the considered interaction has been studied. Wave functions and discontinuity conditions of the problem considered have been derived in detail. Using the continuity equation, we have found a constraint implying the conservation law of the probability current.


Introduction
The mathematical structure of quantum mechanics, one of the fundamental pillars of modern physics of the description of new aspects of nature, consists of Hilbert spaces defined over the field of complex numbers. Physicists have witnessed that the underlying foundations of the theory are enormously successful in describing different quantal phenomena [1][2][3][4][5][6]. For the first time quantum theory was extended to include the field of quaternions by Hamilton [7,8].
Mathematically the quaternions can be expressed by extending complex numbers in the following form: where i, j and k are imaginary units and q l , l = 1, 2, 3, 4 are real. The imaginary units are interrelated by i j = − ji = k, jk = −k j = i, ki = −ik = j, (1.2) which signals the non-commutative character of the field of quaternions with respect to the multiplication. a e-mail: hadisobhani8637@gmail.com Quaternions have a representation similar to complex numbers as follows: and the conjugation of quaternion is defined as Having appropriate mathematical tools using quaternions [9,10], a lot of efforts have been made to reconstruct quantum mechanics in terms of quaternion functions during the past few decades such as those by Kaneno [11], Finkelstein and Jauch [12,13], Emch [14], Horwitz and Biedenharn [15], De Leo [16][17][18][19][20], to name a few; maybe the best-known person in this research field is Adler [21][22][23][24][25][26]. In recent times, scattering for non-relativistic and spinless quantum particles has been studied in [27][28][29] in the presence of a quaternionic Dirac delta potential in the direction proposed by De Leo [20]. Motivated by the above progress, in this article, we want to investigate the relativistic scattering of fermions in the presence of scalar and vector potentials. In order to do this, we investigate the quaternionic version of the Dirac equation in Sect. 2. In Sect. 3 we introduce the scattering potential and present related calculations. We derive the reflection and transmission coefficients to establish the conservation of the current probability in Sect. 4, and finally in Sect. 5 we present the conclusion of our work.

Quaternionic version of the Dirac equation
According to the rules of quaternionic quantum mechanics we should write the Dirac equation in terms of an antihermitian operator; this version of the Dirac equation [20] in the presence of a quaternionic vector and scalar potential can be written as where we have usedh = c = 1 and supposed that S a (r), V a (r) are real functions and S b (r), V b (r) are complex functions in Eq. (2.1). Let us test Eq. (2.1) carefully.
In the complex limit this equation must recover the famous Dirac equation. In order to check this important point, we set 2) is the complex version of the Dirac equation [30,31]. In Eqs. (2.1) and (2.2) α and β are defined by the Pauli and unity matrices as follows: Setting the wave function (r, t) = (r)e −i Et , we get to the eigenvalue form of Eq. (2.1): (r) is a spinor with two components, up and down, that we denote + (r) and − (r), respectively. By some calculation and setting i S a (r) , we can find a connection between these two components: where the interaction direction has been considered in the x-direction. The quaternionic wave function components are in the form of Considering such a form of the wave function components in Eq. (2.6) yields Substitution of Eqs. (2.7) and (2.8) into Eq. (2.5) provides a system of coupled differential equations, where p 2 = E 2 − m 2 and * means the complex conjugation. Now, we are in a position to investigate scattering states due to the quaternionic Dirac delta potential.

Investigation of quaternionic Dirac delta scattering
The interaction which we want to study consists of a Dirac delta function with a quaternion coefficient where V a and V b are real constants. Inserting Eq. (3.1) into Eqs. (2.9) and (2.10) and taking an integral in the neighborhood of x = 0 result in the discontinuity condition of the wave function components, We would like to assume that the particles are coming from region I (x < 0). After experiencing the scattering potential at the origin, some of them will be reflected into region I and the others will be transmitted into region II (x > 0). For the free particles we have where the coefficients are complex constants in general. Therefore, according to our assumption as regards the particles, the physical wave functions can be written as in which the coefficients will be determined using continuity and discontinuity conditions of the wave functions at x = 0. Continuity of the wave functions necessitates Using discontinuity of derivatives of the wave functions results in The coefficients can be derived by solving the system of equations (3.8), (3.9), (3.10) and (3.11). After solving the system of equations we have (3.14) In the next section, the conservation law of probability will be studied.

Reflection and transmission coefficients
In quaternionic quantum mechanics, similar to the complex version, we have the continuity equation .

(4.5)
It is straightforward to prove that by using the definition J x = † α x , we derive the constraint Details of the derivation of Eq. (4.6) are presented in the appendix. We plotted treatments of the reflection and transmission coefficients in Fig. 1, considering m = 1, V a = V b = 0.5. Furthermore, in Figs. 2 and 3, different numerical values for elements of the considered potential have been assumed and effects of these elements on reflection |r | 2 , transmission |t| 2 coefficients have been plotted.
As can be understood from Figs. 2 and 3, the effects of V a on the reflection and transmission coefficients are stronger than V b .

Conclusion
In this article the quaternionic form of the Dirac equation was studied. Because of the non-commutative nature of the quaternions, the derivation of the wave function in the formalism considered is different from the usual one. In this study a quaternionic Dirac delta potential was considered. We have shown how such a potential provided a discontinuity condition for the derivative of the wave functions. Using boundary conditions, we were able to find explicit form of reflection and transmission coefficients. The conservation law of probability was checked. The correctness of the results was checked and proved as we expected from ordinary quantum mechanics. Furthermore, in some plots,  . (6.1) It should be noted that the spinor of a wave function is a quaternion and since the coefficients are complex constants, their order and j, the imaginary unit, is important. Hence in the daggered form of the spinor we are faced with a reversed order in Eq. (6.1). Proceeding in the matrix multiplication form we have