Scientific Research

An Academic Publisher

Fluid State of Dirac Quantum Particles ()

^{}

Keywords

Share and Cite:

*Journal of Modern Physics*,

**9**, 2402-2419. doi: 10.4236/jmp.2018.914154.

1. Introductory Summary

In our previous works on spacetime structures of quantum particles, we suggest that quantum particles should be endowed with geometric and topological structures of differentiable manifolds and their motion should be described as isometric embeddings in higher Euclidean space. We also suggest that all quantum particles are formed from mass points which are joined together by contact forces as a consequence of viewing quantum particles as CW-complexes [1] [2] [3]. Fundamentally, we show that the three main dynamical descriptions of physical events in classical physics, namely Newton mechanics, Maxwell electromagnetism and Einstein gravitation, can be formulated in the same general covariant form and they can be represented by the general equation

${\nabla}_{\beta}M=kJ$ (1)

where M is a mathematical object that represents the corresponding physical system and ${\nabla}_{\beta}$ is a covariant derivative. For Newton mechanics,

$M=\frac{1}{2}m{\displaystyle {\sum}_{\mu =1}^{3}{\left(\text{d}{x}^{\mu}/\text{d}t\right)}^{2}}+V$ and $J=0$. For Maxwell electromagnetism,

$M={F}^{\alpha \beta}={\partial}^{\mu}{A}^{\nu}-{\partial}^{\nu}{A}^{\mu}$, with the four-vector potential ${A}^{\mu}\equiv \left(V,A\right)$ and J can be identified with the electric and magnetic currents. And for Einstein gravitation, $M={R}^{\alpha \beta}$ and J can be defined in terms of a metric ${g}_{\alpha \beta}$ and the Ricci

scalar curvature using the Bianchi identities ${\nabla}_{\beta}{R}^{\alpha \beta}=\frac{1}{2}{g}^{\alpha \beta}{\nabla}_{\beta}R$, that is,

$J=\frac{1}{2}{g}^{\alpha \beta}{\nabla}_{\beta}R$. If we use the Bianchi identities as field equations for the gravitational field then Einstein field equations ${T}_{\mu \nu}=k\left({R}_{\mu \nu}-\frac{1}{2}R{g}_{\mu \nu}+\text{\Lambda}{g}_{\mu \nu}\right)$, as in

the case of the electromagnetic field, should be regarded as a definition for the energy-momentum tensor ${T}_{\mu \nu}$ for the gravitational field [4]. An interesting feature that emerges from Equation (1) for the gravitational field is that we can derive the Ricci flow $\partial {g}_{\alpha \beta}/\partial t=\kappa {R}_{\alpha \beta}$ for a vacuum field $J=0$. Mathematically, the Ricci flow is a geometric process that can be employed to smooth out irregularities of a Riemannian manifold [5]. From the definition of the four-current ${j}^{\alpha}=\left(\rho ,{j}_{i}\right)=\frac{1}{2}{g}^{\alpha \beta}{\nabla}_{\beta}R$ for the gravitational field, by comparing with the Poisson equation for a potential V in classical physics, ${\nabla}^{2}V=4\text{\pi}\rho $, we can identify the scalar potential V with the Ricci scalar curvature R and then obtain a diffusion equation ${\partial}_{t}R=k{\nabla}^{2}R$ whose solutions can be found to take the form $R\left(x,y,z,t\right)=\left(M/{\left(\sqrt{4\text{\pi}kt}\right)}^{3}\right){\text{e}}^{-\left({x}^{2}+{y}^{2}+{z}^{2}\right)/4kt}$, which determines the probabilistic distribution of an amount of geometrical substance M which is defined via the Ricci scalar curvature R and manifests as observable matter [6]. It is worth mentioning that in fact a similar diffusion equation can also be derived from the Ricci flow $\partial {g}_{\alpha \beta}/\partial t=\kappa {R}_{\alpha \beta}$ of the form $\partial R/\partial t=\Delta R+2{\left|\text{Ric}\right|}^{2}$, where $\Delta $ is the Laplacian defined as $\Delta ={g}^{\alpha \beta}{\nabla}_{\alpha}{\nabla}_{\beta}$ and $\left|\text{Ric}\right|$ is a shorthand for a mathematical expression that we will not be concerned with in this work [7]. Therefore, the Bianchi field equations of general relativity in the covariant form given in Equation (1) can be used to formulate quantum particles as differentiable manifolds. For example, we showed that the Ricci scalar curvature R associated with a differentiable manifold that represents a quantum system, such as the hydrogen atom, can be expressed in terms of the Schrödinger wavefunction $\psi $ in quantum mechanics as $R=k\left({\displaystyle {\sum}_{\mu =1}^{3}{\left(\text{d}{x}^{\mu}/\text{d}t\right)}^{2}}-\left(\hslash /m\right)\left({\partial}_{t}\psi +{\displaystyle {\sum}_{\mu =1}^{3}{\partial}_{\mu}\psi \left(\text{d}{x}^{\mu}/\text{d}t\right)/\psi}\right)\right)$.

On the other hand, we have also shown that Maxwell field equations of electromagnetism and Dirac relativistic equation of quantum mechanics can be formulated covariantly from a general system of linear first order partial differential equations [8] [9]. An explicit form of a system of linear first order partial differential equations can be written as follows [10] [11]

$\underset{i=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\underset{j=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{a}_{ij}^{r}\frac{\partial {\psi}_{i}}{\partial {x}_{j}}={k}_{1}\underset{l=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}{b}_{l}^{r}{\psi}_{l}+{k}_{2}{c}^{r},r=1,2,\cdots ,n$ (2)

The system of equations given in Equation (2) can be rewritten in a matrix form as

$\left(\underset{i=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{A}_{i}\frac{\partial}{\partial {x}_{i}}\right)\psi ={k}_{1}\sigma \psi +{k}_{2}J$ (3)

where $\psi ={\left({\psi}_{1},{\psi}_{2},\cdots ,{\psi}_{n}\right)}^{\text{T}}$, $\partial \psi /\partial {x}_{i}={\left(\partial {\psi}_{1}/\partial {x}_{i},\partial {\psi}_{2}/\partial {x}_{i},\cdots ,\partial {\psi}_{n}/\partial {x}_{i}\right)}^{\text{T}}$, ${A}_{i}$, $\sigma $ and J are matrices representing the quantities ${a}_{ij}^{k}$, ${b}_{l}^{r}$ and ${c}^{r}$, and ${k}_{1}$ and ${k}_{2}$ are undetermined constants. Now, if we apply the operator $\sum}_{i=1}^{n}{A}_{i}\partial /\partial {x}_{i$ on the left on both sides of Equation (3) then we obtain

$\begin{array}{l}\left(\underset{i=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{A}_{i}^{2}\frac{{\partial}^{2}}{\partial {x}_{i}^{2}}+\underset{i=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\underset{j>i}{\overset{n}{{\displaystyle \sum}}}\left({A}_{i}{A}_{j}+{A}_{j}{A}_{i}\right)\frac{{\partial}^{2}}{\partial {x}_{i}\partial {x}_{j}}\right)\psi \\ ={k}_{1}^{2}{\sigma}^{2}\psi +{k}_{1}{k}_{2}\sigma J+{k}_{2}\underset{i=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{A}_{i}\frac{\partial J}{\partial {x}_{i}}\end{array}$ (4)

In order for the above systems of partial differential equations to be used to describe physical phenomena, the matrices ${A}_{i}$ must be determined. We have shown that for both Dirac and Maxwell field equations, the matrices ${A}_{i}$ must take a form so that Equation (4) reduces to the following equation

$\left(\underset{i=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{A}_{i}^{2}\frac{{\partial}^{2}}{\partial {x}_{i}^{2}}\right)\psi ={k}_{1}^{2}{\sigma}^{2}\psi +{k}_{1}{k}_{2}\sigma J+{k}_{2}\underset{i=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{A}_{i}\frac{\partial J}{\partial {x}_{i}}$ (5)

To obtain Dirac equation we simply set
${A}_{i}{A}_{j}+{A}_{j}{A}_{i}=0$ with
${A}_{i}^{2}=\pm 1$, and in this case the matrices
${A}_{i}$ are the matrices

$\begin{array}{l}{\gamma}_{1}=\left(\begin{array}{cccc}1& 0& 0& 0\\ 0& 1& 0& 0\\ 0& 0& -1& 0\\ 0& 0& 0& -1\end{array}\right),{\gamma}_{2}=\left(\begin{array}{cccc}0& 0& 0& 1\\ 0& 0& 1& 0\\ 0& -1& 0& 0\\ -1& 0& 0& 0\end{array}\right),\\ {\gamma}_{3}=\left(\begin{array}{cccc}0& 0& 0& -i\\ 0& 0& i& 0\\ 0& i& 0& 0\\ -i& 0& 0& 0\end{array}\right),{\gamma}_{4}=\left(\begin{array}{cccc}0& 0& 1& 0\\ 0& 0& 0& -1\\ -1& 0& 0& 0\\ 0& 1& 0& 0\end{array}\right)\end{array}$ (6)

For Maxwell field equations, in order to specify the matrices ${A}_{i}$ we need to use the form of Maxwell field equations established in classical electrodynamics [13] [14]. And the matrices ${A}_{i}$ take the forms

$\begin{array}{l}{A}_{1}=\left(\begin{array}{cccccc}-1& 0& 0& 0& 0& 0\\ 0& -1& 0& 0& 0& 0\\ 0& 0& -1& 0& 0& 0\\ 0& 0& 0& 1& 0& 0\\ 0& 0& 0& 0& 1& 0\\ 0& 0& 0& 0& 0& 1\end{array}\right),{A}_{2}=\left(\begin{array}{cccccc}0& 0& 0& 0& 0& 0\\ 0& 0& 0& 0& 0& -1\\ 0& 0& 0& 0& 1& 0\\ 0& 0& 0& 0& 0& 0\\ 0& 0& -1& 0& 0& 0\\ 0& 1& 0& 0& 0& 0\end{array}\right),\\ {A}_{3}=\left(\begin{array}{cccccc}0& 0& 0& 0& 0& 1\\ 0& 0& 0& 0& 0& 0\\ 0& 0& 0& -1& 0& 0\\ 0& 0& 1& 0& 0& 0\\ 0& 0& 0& 0& 0& 0\\ -1& 0& 0& 0& 0& 0\end{array}\right),{A}_{4}=\left(\begin{array}{cccccc}0& 0& 0& 0& -1& 0\\ 0& 0& 0& 1& 0& 0\\ 0& 0& 0& 0& 0& 0\\ 0& -1& 0& 0& 0& 0\\ 1& 0& 0& 0& 0& 0\\ 0& 0& 0& 0& 0& 0\end{array}\right),\\ {A}_{5}=\left(\begin{array}{cccccc}\mu & 0& 0& 0& 0& 0\\ 0& \mu & 0& 0& 0& 0\\ 0& 0& \mu & 0& 0& 0\\ 0& 0& 0& 1& 0& 0\\ 0& 0& 0& 0& 1& 0\\ 0& 0& 0& 0& 0& 1\end{array}\right)\end{array}$ (7)

Besides the covariant formulations of classical and quantum physics as described above, we have also discussed the topological transformation of quantum dynamics by showing the wave dynamics of a quantum particle on different types of topological structures in various dimensions from the fundamental polygons of the corresponding universal covering spaces [15]. We presented our discussions in the form of Bohr model in one, two and three dimensions using linear wave equations. For the clarity of our presentation in terms of Bohr model, we want to mention here that in order to successfully construct a model for the hydrogen atom which predicts correctly the spectrum of the energy radiated from the atom, Bohr proposed three postulates which state that the centripetal force required for the electron to orbit the nucleus in a stable circle is the Coulomb force $m{v}^{2}/r=k{q}^{2}/{r}^{2}$, the permissible orbits are those that satisfy the condition that the angular momentum of the electron equals $n\hslash $, that is $mvr=n\hslash $, and when the electron moves in one of the stable orbits it does not radiate, however, it will radiate when it makes a transition between the stable orbits [16]. Furthermore, in his work on the concept of matter wave, de Broglie proposed that an electron has both a wave and a particle nature by regarding the electron as a standing wave around the circumference of an orbit, as shown in Figure 1 [17].

It is seen that de Broglie’s requirement leads to the wave condition $2\text{\pi}r=n\lambda $. This is equivalent to assuming that the standing wave around a circle, which is a 1-sphere, is similar to a standing wave on the fundamental interval of a straight line which is the universal covering space of the circle ${S}^{1}$, where the translations taking the interval to the next images will generate the holonomy group [18]. If we apply de Broglie wavelength $\lambda $ defined in terms of the momentum

Figure 1. de Broglie waves around a circle.

of a quantum particle $p=mv$ as $\lambda =h/mv$, then using the wavelength also given by $\lambda =2\text{\pi}R/n$ we obtain $h/mv=2\text{\pi}R/n$, that leads to the Bohr’s postulate of the quantisation of angular momentum $mvR=n\hslash $. Then the energy spectrum ${E}_{n}$ can be calculated from Coulomb’s law $m{v}^{2}/R=k{q}^{2}/{R}^{2}$ as ${E}_{n}=-m{k}^{2}{q}^{4}/2{\hslash}^{2}{n}^{2}$, where R now is the radius of the nth stationary orbit.

In mathematics, the fundamental polygon in one dimension is an interval and the universal covering space is the straight line and in this case the standing wave on a finite string is transformed into the standing wave on a circle which can be applied into the Bohr model of the hydrogen atom. In two dimensions, the fundamental polygon is a square and the universal covering space is the plane and in this case the standing wave on the square is transformed into the standing wave on different surfaces that can be formed by gluing opposite sides of the square, which include a 2-sphere, a 2-torus, a Klein bottle and a projective plane. This may be seen as an extension of the Bohr model of the hydrogen atom from one-dimensional manifolds of the 1-sphere and 1-torus embedded in the ambient two-dimensional Euclidean space R^{2} into two-dimensional manifolds embedded or immersed in the ambient three-dimensional Euclidean space R^{3}. In three dimensions, the fundamental polygon is a cube and the universal covering space is the three-dimensional Euclidean space. It is shown that a 3-torus and the manifold
$K\times {S}^{1}$ defined as the product of a Klein bottle and a circle can be constructed by gluing opposite faces of a cube therefore in three-dimensions the standing wave on a cube is transformed into the standing wave on a 3-torus or on the manifold
$K\times {S}^{1}$. We also discuss a transformation of a stationary wave on the fundamental cube into a stationary wave on a 3-sphere despite it still remains unknown whether a 3-sphere can be constructed directly from a cube by gluing its opposite faces. In spite of this uncertainty, however, we speculate that mathematical degeneracy in which an element of a class of objects degenerates into an element of a different but simpler class may play an important role in quantum dynamics. For example, a 2-sphere is a degenerate 2-torus when the axis of revolution passes through the centre of the generating circle. Therefore, it seems reasonable to assume that if an n-torus degenerates into an n-sphere then wavefunctions on an n-torus may also be degenerated into wavefunctions on an n-sphere. Furthermore, since an n-sphere can degenerate itself into a single point, therefore the mathematical degeneracy may be related to the concept of wavefunction collapse in quantum mechanics where the classical observables such as position and momentum can only be obtained from the collapse of the associated wavefunctions for physical measurements. This consideration suggests that quantum particles associated with differentiable manifolds may possess the more stable mathematical structures of an n-torus rather than those of an n-sphere.

The above formulation of quantum particles in terms of differentiable manifolds and the consideration of their intrinsic geometric and topological characteristics raise the question of how the standing waves that represent quantum particles could be established physically. The aim of this work is to answer this question by showing that Dirac equation in fact describes quantum particles as composite structures that are in a fluid state in which the components of the wavefunction can be identified with the stream function and the velocity potential of a potential flow formulated in the theory of classical fluids. In this case Dirac quantum particles can manifest as standing waves which are the result of the superposition of two fluid flows moving in opposite directions. For example, if two opposite waves are represented by the function ${\psi}_{1}\left(r,t\right)=a\mathrm{sin}\left(kr-vt\right)$ and ${\psi}_{2}\left(r,t\right)=a\mathrm{sin}\left(kr+vt\right)$ then the resultant standing wave can be obtained as ${\psi}_{1}\left(r,t\right)+{\psi}_{2}\left(r,t\right)=2a\mathrm{sin}\left(kr\right)\mathrm{cos}\left(vt\right)$. We also show that even though a Dirac quantum particle has the potential to transfer energy in opposite directions to establish a standing wave within its physical structure, for a steady motion without an external disturbance a Dirac quantum particle can be considered as a classical particle defined in classical physics.

2. Two-Dimensional Hydrogen-Like Physical System

In this work we will show that quantum particles can be described as physical systems in a state of fluids that can be formulated in terms of the fluid dynamics in two dimensions. Since such systems can be seen to have the physical structure of a two-dimensional hydrogen-like atom therefore in this section we will examine further how they can be described in terms of quantum mechanics. First we need to extend our formulation of Maxwell field equations of electromagnetism and Dirac relativistic equation from a general system of linear first order partial differential equations to that with an external field. Such system of equations is given as follows [19]

$\underset{i=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\underset{j=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{a}_{ij}^{r}\frac{\partial {\psi}_{i}}{\partial {x}_{j}}=\underset{i=1}{\overset{n}{{\displaystyle \sum}}}\left(\underset{j=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{b}_{ij}^{r}{V}_{j}+{c}_{i}^{r}\right){\psi}_{i}+{d}^{r},r=1,2,\cdots ,n$ (8)

The system of equations given in Equation (8) can be rewritten in a matrix form as

$\left(\underset{i=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{A}_{i}\frac{\partial}{\partial {x}_{i}}\right)\psi =-i\left(\underset{i=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}q{B}_{i}{V}_{i}+m\sigma \right)\psi +J$ (9)

where $\psi ={\left({\psi}_{1},{\psi}_{2},\cdots ,{\psi}_{n}\right)}^{\text{T}}$, $\partial \psi /\partial {x}_{i}={\left(\partial {\psi}_{1}/\partial {x}_{i},\partial {\psi}_{2}/\partial {x}_{i},\cdots ,\partial {\psi}_{n}/\partial {x}_{i}\right)}^{\text{T}}$ with ${A}_{i}$, ${B}_{i}$, σ and J are matrices representing the quantities ${a}_{ij}^{r}$, ${b}_{ij}^{r}$, ${c}_{j}^{r}$ and ${d}^{r}$, which are assumed to be constant in this work. While the quantities q, m and J represent physical entities related directly to the physical properties of the particle, and the quantities ${V}_{i}$ represent an external field, such as the potentials of an electromagnetic field. By applying the operator $\sum}_{i=1}^{n}{A}_{i}\partial /\partial {x}_{i$ on the left on both sides of Equation (9), with the assumption that the coefficients ${a}_{ij}^{k}$, ${b}_{i}^{r}$ and ${c}^{r}$ are constants, then we obtain

$\left(\underset{i=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{A}_{i}\frac{\partial}{\partial {x}_{i}}\right)\left(\underset{j=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{A}_{j}\frac{\partial}{\partial {x}_{j}}\right)\psi =\left(\underset{i=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{A}_{i}\frac{\partial}{\partial {x}_{i}}\right)\left(-i\left(\underset{j=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}q{B}_{j}{V}_{j}+m\sigma \right)\psi +J\right)$ (10)

Since the quantities ${A}_{i}$, ${B}_{i}$, σ, q, m and J are assumed to be constant, Equation (10) becomes

$\begin{array}{l}\left(\underset{i=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{A}_{i}{}^{2}\frac{{\partial}^{2}}{\partial {x}_{i}^{2}}+\underset{i=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\underset{j>i}{\overset{n}{{\displaystyle \sum}}}\left({A}_{i}{A}_{j}+{A}_{j}{A}_{i}\right)\frac{{\partial}^{2}}{\partial {x}_{i}\partial {x}_{j}}\right)\psi \\ =\left(-i\left(\underset{i=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{A}_{i}\frac{\partial}{\partial {x}_{i}}\right)\left(\underset{j=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}q{B}_{j}{V}_{j}+m\sigma \right)\right)\psi \\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}-i\left(\underset{i=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}q{B}_{i}{V}_{i}+m\sigma \right)\left(\left(\underset{j=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{A}_{j}\frac{\partial}{\partial {x}_{j}}\right)\psi \right)+\underset{i=1}{\overset{n}{{\displaystyle \sum}}}{A}_{i}\frac{\partial J}{\partial {x}_{i}}\\ =-i\left(\underset{i=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\underset{j=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}q{A}_{i}{B}_{j}\frac{\partial {V}_{j}}{\partial {x}_{i}}\right)\psi -(\underset{i=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\underset{j>i}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{q}^{2}\left({B}_{i}{B}_{j}+{B}_{j}{B}_{i}\right){V}_{i}{V}_{j}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}-2i\underset{i=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}qm{B}_{i}{V}_{i}\sigma -{m}^{2}{\sigma}^{2})\psi -i\left(\underset{i=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}q{B}_{i}{V}_{i}+m\sigma \right)J+\underset{i=1}{\overset{n}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{A}_{i}\frac{\partial J}{\partial {x}_{i}}\end{array}$ (11)

Dirac equation for an arbitrary field can be formulated from the system of linear first order partial differential equations given in Equation (9) by setting ${B}_{i}={A}_{i}={\gamma}_{i}$, $\sigma =1$, $J=0$ and ${A}_{i}{A}_{j}+{A}_{j}{A}_{i}=0$. In this case, in terms of the operator ${\gamma}_{i}$, Equation (9) becomes

$\left(\underset{i=1}{\overset{4}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\gamma}_{i}\frac{\partial}{\partial {x}_{i}}\right)\psi =-i\left(\underset{i=1}{\overset{4}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}q{\gamma}_{i}{V}_{i}+m\right)\psi $ (12)

Equation (12) can be written in a covariant form as Dirac equation for an arbitrary field as

$\left({\gamma}^{\mu}\left(i{\partial}_{\mu}-q{V}_{\mu}\right)-m\right)\psi =0$ (13)

Equation (11) also reduces to the following equation

$\left(\underset{i=1}{\overset{4}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\gamma}_{i}{}^{2}\frac{{\partial}^{2}}{\partial {x}_{i}^{2}}\right)\psi =\left(-i\underset{i=1}{\overset{4}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\underset{j>i}{\overset{4}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}q{\gamma}_{i}{\gamma}_{j}\left(\frac{\partial {V}_{j}}{\partial {x}_{i}}-\frac{\partial {V}_{i}}{\partial {x}_{j}}\right)+2i\underset{i=1}{\overset{4}{{\displaystyle \sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}qm{\gamma}_{i}{V}_{i}-{m}^{2}\right)\psi $ (14)

Even though in the following we will examine only physical states of Dirac quantum particles in which ${V}_{i}=0$ where the physical quantity ${V}_{i}$ is assumed to be associated with an external field, however, if we consider quantum particles as differentiable manifolds which are formed by mass points joined together by contact forces then we may suggest that they are endowed with intrinsic geometric and topological structures and in this case the quantity ${V}_{i}$ may be considered as an internal field that is responsible for the stability of the physical structure of a quantum particle. As we will show below this is in fact the case when at least part of a quantum particle exists as a two-dimensional structure in which the intrinsic angular momentum can take half-integer values. The problem that we considered can actually be started with Dirac equation given in Equation (13). It can be shown that in the non-relativistic limit, Dirac equation reduces to the Pauli equation for stationary system as [20]

$\left(\frac{1}{2\mu}{\left(-i\hslash \nabla -qA\right)}^{2}-\frac{q\hslash}{2\mu}\left(\sigma \cdot B\right)+q\varphi \right)\psi \left(r\right)=E\psi \left(r\right)$ (15)

In the case when $A=0$, $B=0$, and $q\varphi =k/r$ then we have

$-\frac{{\hslash}^{2}}{2\mu}{\nabla}^{2}\psi \left(r\right)-\frac{k}{r}\psi \left(r\right)=E\psi \left(r\right)$ (16)

Now let us examine a physical system that is described by the Schrödinger wave equation given in Equation (16) from the viewpoint of an observer who sees it as a planar system [21] [22] [23]. It is shown that if we consider physical systems whose configuration space is multiply connected, such as the physical system of a hydrogen-like atom in two rather than three dimensions, then multivalued wavefunctions can be used to describe the system [24]. In two-dimensional space, the Schrödinger equation in the planar polar coordinates takes the form

$-\frac{{\hslash}^{2}}{2\mu}\left[\frac{1}{r}\frac{\partial}{\partial r}\left(r\frac{\partial}{\partial r}\right)+\frac{1}{{r}^{2}}\frac{{\partial}^{2}}{\partial {\varphi}^{2}}\right]\psi \left(r,\varphi \right)-\frac{k}{r}\psi \left(r,\varphi \right)=E\psi \left(r,\varphi \right)$ (17)

Solutions of the form $\psi \left(r,\varphi \right)=R\left(r\right)\Phi \left(\varphi \right)$ then reduce the above equation to two separate equations for the functions $\Phi $ and R

$\frac{{\text{d}}^{2}\Phi}{\text{d}{\varphi}^{2}}+{m}^{2}\Phi =0$ (18)

$\frac{{\text{d}}^{2}R}{\text{d}{r}^{2}}+\frac{1}{r}\frac{\text{d}R}{\text{d}r}-\frac{{m}^{2}}{{r}^{2}}R+\frac{2\mu}{{\hslash}^{2}}\left(\frac{k}{r}-E\right)R=0$ (19)

where m is identified as the angular momentum of the system. From the system of ordinary differential equations given in Equations (18) and (19), the energy spectrum can be found as

$E=-\frac{{k}^{2}\mu}{2{\hslash}^{2}{\left(n+m+1/2\right)}^{2}}$ (20)

It is seen that if the physical system is the Bohr model of two-dimensional hydrogen-like atom then the angular momentum m must take half-integral values. Hence, the topological structure of a configuration space of a physical system can determine the quantum nature of an observable of the system. This result should be expected in quantum mechanics since we know that the quantum behaviour of a particle depends almost entirely on the configuration of an experiment. If, in a particular experiment, the electron of a hydrogen-like atom is constrained to move in a plane, then the orbital angular momentum of the electron must take half-integral values if we use the Schrödinger equation to study the dynamics of the electron and want to retain the same energy spectrum as the Bohr model. As a consequence, it might seem possible to invoke the result to explain the Stern-Gerlach experiment without the necessity of introduction of spin into the quantum theory.

3. Dirac Real Equation

In this section we show that Dirac equation for a free particle can be used to describe the state of a fluid of the quantum particle formulated in the theory of classical fluids. For free Dirac quantum particles, Equation (13) reduces to

$\left(i{\gamma}^{\mu}{\partial}_{\mu}-m\right)\psi =0$ (21)

By expanding Equation (21) using the matrices ${\gamma}_{i}$ given in Equation (6), we obtain

$-\frac{\partial {\psi}_{1}}{\partial t}-im{\psi}_{1}=\left(\frac{\partial}{\partial x}-i\frac{\partial}{\partial y}\right){\psi}_{4}+\frac{\partial {\psi}_{3}}{\partial z}$ (22)

$-\frac{\partial {\psi}_{2}}{\partial t}-im{\psi}_{2}=\left(\frac{\partial}{\partial x}+i\frac{\partial}{\partial y}\right){\psi}_{3}-\frac{\partial {\psi}_{4}}{\partial z}$ (23)

$\frac{\partial {\psi}_{3}}{\partial t}-im{\psi}_{3}=\left(-\frac{\partial}{\partial x}+i\frac{\partial}{\partial y}\right){\psi}_{2}-\frac{\partial {\psi}_{1}}{\partial z}$ (24)

$\frac{\partial {\psi}_{4}}{\partial t}-im{\psi}_{4}=\left(-\frac{\partial}{\partial x}-i\frac{\partial}{\partial y}\right){\psi}_{1}+\frac{\partial {\psi}_{2}}{\partial z}$ (25)

First, it is observed that with the form of the field equations given in Equations (22)-(25), we may interpret that the change of the field $\left({\psi}_{1},{\psi}_{2}\right)$ with respect to time generates the field $\left({\psi}_{3},{\psi}_{4}\right)$, similar to the case of Maxwell field equations in which the change of the electric field generates the magnetic field. With this observation it may be suggested that, like the Maxwell electromagnetic field which is composed of two essentially different physical fields, the Dirac field of massive particles may also be viewed as being composed of two different physical fields, namely the field $\left({\psi}_{1},{\psi}_{2}\right)$, which plays the role of the electric field in Maxwell field equations, and the field $\left({\psi}_{3},{\psi}_{4}\right)$, which plays the role of the magnetic field. The similarity between Maxwell field equations and Dirac field equations can be carried further by showing that it is possible to reformulate Dirac equation as a system of real equations. When we formulate Maxwell field equations from a system of linear first order partial differential equations we rewrite the original Maxwell field equations from a vector form to a system of first order partial differential equations by equating the corresponding terms of the vectorial equations [9]. Now, since, in principle, a complex quantity is equivalent to a vector quantity therefore in order to form a system of real equations from Dirac complex field equations we equate the real parts with the real parts and the imaginary parts with the imaginary parts. In this case Dirac equation given in Equations (22)-(25) can be rewritten as a system of real equations as follows

$-\frac{\partial {\psi}_{1}}{\partial t}=\frac{\partial {\psi}_{4}}{\partial x}+\frac{\partial {\psi}_{3}}{\partial z}$ (26)

$-\frac{\partial {\psi}_{2}}{\partial t}=\frac{\partial {\psi}_{3}}{\partial x}-\frac{\partial {\psi}_{4}}{\partial z}$ (27)

$-\frac{\partial {\psi}_{3}}{\partial t}=\frac{\partial {\psi}_{2}}{\partial x}+\frac{\partial {\psi}_{1}}{\partial z}$ (28)

$-\frac{\partial {\psi}_{4}}{\partial t}=\frac{\partial {\psi}_{1}}{\partial x}-\frac{\partial {\psi}_{2}}{\partial z}$ (29)

$\frac{\partial {\psi}_{4}}{\partial y}=m{\psi}_{1}$ (30)

$\frac{\partial {\psi}_{3}}{\partial y}=-m{\psi}_{2}$ (31)

$\frac{\partial {\psi}_{2}}{\partial y}=-m{\psi}_{3}$ (32)

$\frac{\partial {\psi}_{1}}{\partial y}=m{\psi}_{4}$ (33)

If the wavefunction $\psi $ satisfies Dirac field equations given in Equations (26)-(33) then we can derive the following system of equations for all components

$\frac{{\partial}^{2}{\psi}_{i}}{{\partial}^{2}y}-{m}^{2}{\psi}_{i}=0$ (34)

$\frac{{\partial}^{2}{\psi}_{i}}{\partial {t}^{2}}-\frac{{\partial}^{2}{\psi}_{i}}{\partial {x}^{2}}-\frac{{\partial}^{2}{\psi}_{i}}{\partial {z}^{2}}=0$ (35)

Solutions to Equation (34) are

${\psi}_{i}={c}_{1i}\left(x,z\right){\text{e}}^{my}+{c}_{2i}\left(x,z\right){\text{e}}^{-my}$ (36)

where ${c}_{1i}$ and ${c}_{2i}$ are undetermined functions of $\left(x,z\right)$. The solutions given in Equation (36) give a distribution of a physical quantity along the y-axis. On the other hand, Equation (35) can be used to describe the dynamics, for example, of a vibrating membrane in the $\left(x,z\right)$ -plane. If the membrane is a circular membrane of radius a then the domain D is given as $D=\left\{{x}^{2}+{z}^{2}<{a}^{2}\right\}$. In the polar coordinates given in terms of the Cartesian coordinates $\left(x,z\right)$ as $x=r\mathrm{cos}\theta $, $z=r\mathrm{sin}\theta $, the two-dimensional wave equation given in Equation (35) becomes

$\frac{1}{{c}^{2}}\frac{{\partial}^{2}\psi}{\partial {t}^{2}}-\frac{{\partial}^{2}\psi}{\partial {r}^{2}}-\frac{1}{r}\frac{\partial \psi}{\partial r}-\frac{1}{{r}^{2}}\frac{{\partial}^{2}\psi}{\partial {\theta}^{2}}=0$ (37)

The general solution to Equation (37) for the vibrating circular membrane with the condition $\psi =0$ on the boundary of D can be found as [6] [11]

$\begin{array}{c}\psi \left(r,\theta ,t\right)=\underset{m=1}{\overset{\infty}{{\displaystyle \sum}}}{J}_{0}\left(\sqrt{{\lambda}_{0m}}r\right)\left({C}_{0m}\mathrm{cos}\sqrt{{\lambda}_{0m}}ct+{D}_{0m}\mathrm{sin}\sqrt{{\lambda}_{0m}}ct\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\underset{m,n=1}{\overset{\infty}{{\displaystyle \sum}}}{J}_{n}\left(\sqrt{{\lambda}_{nm}}r\right)\left({A}_{nm}\mathrm{cos}n\theta +{B}_{nm}\mathrm{sin}n\theta \right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\times \left(\left({C}_{nm}\mathrm{cos}\sqrt{{\lambda}_{nm}}ct+{D}_{nm}\mathrm{sin}\sqrt{{\lambda}_{nm}}ct\right)\right)\end{array}$ (38)

where ${J}_{n}\left(\sqrt{{\lambda}_{nm}}r\right)$ is the Bessel function of order n and the quantities ${A}_{nm}$, ${B}_{nm}$, ${C}_{nm}$ and ${D}_{nm}$ can be specified by the initial and boundary conditions. It is also observed that at each moment of time the vibrating membrane appears as a 2D differentiable manifold which is a geometric object whose geometric structure can be constructed using the wavefunction given in Equation (38) as

$R=\left(2({\psi}_{11}{\psi}_{22}-{\left({\psi}_{12}\right)}^{2}\right)/{\left(1+{\psi}_{1}^{2}+{\psi}_{2}^{2}\right)}^{2}$

where ${\psi}_{\mu}=\partial \psi /\partial {x}^{\mu}$ and ${\psi}_{\mu \nu}={\partial}^{2}\psi /\partial {x}^{\mu}\partial {x}^{\nu}$ [25]. Even though elementary particles may have the geometric and topological structures of a 3D differentiable manifold, it is seen from the above descriptions via the Schrödinger wave equation and Dirac equation that they appear as 3D physical objects that embedded in three-dimensional Euclidean space. Interestingly, in the following we will show that the solution given in Equation (38) can be used to describe a standing wave in a fluid due to the motion of two waves in opposite directions. For a steady state in which the system is time-independent, the system of equations given in Equations (26)-(29) reduces to the following system of equations

$\frac{\partial {\psi}_{2}}{\partial x}+\frac{\partial {\psi}_{1}}{\partial z}=0,\frac{\partial {\psi}_{1}}{\partial x}-\frac{\partial {\psi}_{2}}{\partial z}=0$ (39)

$\frac{\partial {\psi}_{4}}{\partial x}+\frac{\partial {\psi}_{3}}{\partial z}=0,\frac{\partial {\psi}_{3}}{\partial x}-\frac{\partial {\psi}_{4}}{\partial z}=0$ (40)

It is observed from Dirac equation for steady states that the field $\left({\psi}_{1},{\psi}_{2}\right)$ and the field $\left({\psi}_{3},{\psi}_{4}\right)$ satisfy the Cauchy-Riemann equations in the $\left(x,z\right)$ -plane. We will now discuss whether it is possible to consider Dirac quantum particles as physical systems which exist in a fluid state as defined in the classical fluid dynamics as substances that retain a definite volume, have the ability to flow and deform continually, hence they can exhibit a wave motion. For references in the next section we will outline the main features in the theory of classical fluids, especially, in two dimensions.

4. Fluid Dynamics in Two-Dimensions

In fluid dynamics, quantities that satisfy the Cauchy-Riemann equations can be identified with the velocity potential and the stream function of an incompressible and irrotational flow [26] [27] [28] [29]. In two-dimensional fluid dynamics, a streamline is a theoretical line that is assumed to be tangential to the instantaneous velocity, therefore there is no flow that can cross the streamline. For a continuous stream of fluid, the streamlines can form continuous lines or closed curves. As an illustration, in the following we will consider a free or potential vortex flow whose streamlines are concentric circles in the $\left(x,z\right)$ -plane as shown in Figure 2.

In two-dimensional incompressible flow, the stream function $\Psi $ is defined as a volume flux through a curve given by

$\Psi ={{\displaystyle \int}}^{\text{}}\left({v}_{x}\text{d}z-{v}_{z}\text{d}x\right)$ (41)

From the definition given in Equation (41), we have $\text{d}\Psi ={v}_{x}\text{d}z-{v}_{z}\text{d}x$. On the other hand, the total derivative of the stream function $\Psi $ is given by $\text{d}\Psi =\left(\partial \Psi /\partial x\right)\text{d}x+\left(\partial \Psi /\partial z\right)\text{d}z$, therefore we obtain the following relationships between the velocity components $\left({v}_{x},{v}_{z}\right)$ and the stream function $\Psi $

${v}_{x}=\frac{\partial \Psi}{\partial z},{v}_{z}=-\frac{\partial \Psi}{\partial x}$ (42)

If the stream function $\Psi $ is defined in terms of polar coordinates $\left(r,\theta \right)$ as $\Psi =\Psi \left(r,\theta \right)$, then we have $\text{d}\Psi =\left(\partial \Psi /\partial r\right)\text{d}r+\left(\partial \Psi /\partial \theta \right)\text{d}\theta $. As shown in Figure 3, we also have $\text{d}\Psi ={v}_{r}\left(r\text{d}\theta \right)-{v}_{\theta}\text{d}r$.

Figure 2. Circular streamlines.

Figure 3. Fluid flow in terms of stream function in polar coordinates.

Therefore we also obtain the following relationships between the velocity components $\left({v}_{r},{v}_{\theta}\right)$ and the stream function $\Psi $

${v}_{r}=\frac{1}{r}\frac{\partial \Psi}{\partial \theta},{v}_{\theta}=-\frac{\partial \Psi}{\partial r}$ (43)

From the definition of the stream function we see that the radial component of the velocity of a vortex flow is equals to zero, ${v}_{r}=0$, since there is no flow that can cross the streamlines. In fluid dynamics, the circulation $\Gamma $ around a closed curve is a line integral of velocity $v$ defined as

$\Gamma ={\displaystyle \oint v\cdot \text{d}s}={\displaystyle \oint {v}_{s}\text{d}s}$ (44)

where ${v}_{s}$ is the tangential velocity. By Stokes’ theorem, the circulation $\Gamma $ is related to the vorticity $\omega =\nabla \times v$ as

$\text{\Gamma}={\displaystyle {\oint}_{\partial S}v\cdot \text{d}s}={\displaystyle {\iint}_{S}\omega \cdot \text{d}S}$ (45)

It is seen from the above equation that the flux of vorticity is the circulation. In particular, for a two-dimensional flow in the $\left(x,z\right)$ -plane, the circulation becomes

$\text{\Gamma}={\displaystyle {\int}_{\text{A}}\left(\frac{\partial {v}_{z}}{\partial x}-\frac{\partial {v}_{x}}{\partial z}\right)\text{d}x\text{d}z}$ (46)

A flow for which the circulation is equal to zero, $\omega =\nabla \times v=0$, is called a potential or irrotational flow. In two dimensions it is seen from Equation (46) that the condition for potential flow is

$\frac{\partial {v}_{z}}{\partial x}-\frac{\partial {v}_{x}}{\partial z}=0$ (47)

It should be mentioned that in potential flow we have $\text{\Gamma}={\displaystyle {\oint}_{\partial S}v\cdot \text{d}s}=0$ therefore closed streamlines cannot exist in such flow. In general this result prevents us from identifying the components of Dirac equation ${\psi}_{1}$ as closed stream function at the same time identifying ${\psi}_{2}$ as velocity potential. However, it is seen that the result of $\text{\Gamma}={\displaystyle \oint v\cdot \text{d}s}=0$ may not be valid if the region of space is multiply-connected since the velocity circulation may not be zero if the closed contour cannot be contracted to a point. This important feature can be discussed further as follows. For a two-dimensional irrotational flow given in polar coordinates, it can be shown that the flow velocity ${v}_{\theta}$ and the radius r satisfy the following relationship $r{v}_{\theta}=C$, where C is a constant. The constant C can be established by using the singularity in the irrotational vortex flow where to velocity becomes infinite at the centre of the vortex with the vorticity is given by the relation $\omega =-\left(\partial {v}_{\theta}/\partial r+{v}_{\theta}/r\right)$. In this case the circulation around a circular streamline can be found as

$\text{\Gamma}={\displaystyle \oint v\cdot \text{d}s}={\displaystyle \oint {v}_{s}\text{d}s}=2\text{\pi}r{v}_{\theta}$ (48)

It is also interesting to note that in the Bohr model of the hydrogen atom in which the electron is assumed to move around the nucleus in stationary circular orbits with ${v}_{\theta}=v$ then the angular momentum is quantised as $mrv=n\hslash $. If we now also assume that Dirac quantum particles are in fluid states whose circulation $\Gamma $ is also quantised as the angular momentum then we obtain the following quantisation for the circulation

$\text{\Gamma}=2\text{\pi}rv=\frac{2\text{\pi}n\hslash}{m}=\frac{nh}{m}$ (49)

Since ${v}_{r}=0$, the stream function $\text{\Psi}=\text{\Psi}\left(r,\theta \right)$ can be obtained as follows

$\Psi =-{{\displaystyle \int}}^{\text{}}\text{\hspace{0.05em}}{v}_{\theta}\text{d}r=-{{\displaystyle \int}}^{\text{}}\frac{\Gamma}{2\text{\pi}r}\text{d}r=-\frac{\Gamma}{\text{2\pi}}\mathrm{ln}\left({k}_{1}r\right)=-\frac{n\hslash}{m}\mathrm{ln}\left({k}_{1}r\right)$ (50)

where ${k}_{1}$ is an undetermined constant.

In fluid dynamics, another important concept that is connected with an irrotational flow is the concept of the velocity potential $\Phi $ which is defined in the $\left(x,z\right)$ -plane as

$\Phi ={{\displaystyle \int}}^{\text{}}\text{\hspace{0.05em}}{v}_{s}\text{d}s={{\displaystyle \int}}^{\text{}}\text{\hspace{0.05em}}{v}_{x}\text{d}x+{v}_{z}\text{d}z$ (51)

It is seen from Equation (51) that the velocity components can be expressed in terms of the velocity potential as follows

${v}_{x}=\frac{\partial \Phi}{\partial x},{v}_{z}=\frac{\partial \Phi}{\partial z}$ (52)

In polar coordinates $\left(r,\theta \right)$, the velocity potential and its relationship with the velocity components are given as

$\Phi ={{\displaystyle \int}}^{\text{}}\text{\hspace{0.05em}}{v}_{s}\text{d}s={{\displaystyle \int}}^{\text{}}\text{\hspace{0.05em}}{v}_{r}\text{d}r+{v}_{\theta}r\text{d}\theta ,$ (53)

${v}_{r}=\frac{\partial \text{\Phi}}{\partial r},{v}_{\theta}=\frac{1}{r}\frac{\partial \text{\Phi}}{\partial \theta}$ (54)

Similarly, the velocity potential $\Phi $ can also be obtained using the relation $\text{d}\Phi =r{v}_{\theta}\text{d}\theta +{v}_{r}\text{d}r$ as

$\Phi ={{\displaystyle \int}}^{\text{}}\text{\hspace{0.05em}}r{v}_{\theta}\text{d}\theta ={{\displaystyle \int}}^{\text{}}\text{\hspace{0.05em}}r\left(\frac{\Gamma}{2\text{\pi}r}\right)\text{d}\theta ={{\displaystyle \int}}^{\text{}}\text{\hspace{0.05em}}r\left(\frac{n\hslash}{mr}\right)\text{d}\theta =\frac{n\hslash}{m}\theta +{k}_{2}$ (55)

where ${k}_{2}$ is an undetermined constant. From the relationships given in Equations (42) and (52) we then obtain the Cauchy-Riemann equations

$\frac{\partial \Phi}{\partial x}+\frac{\partial \Psi}{\partial z}=0,\frac{\partial \Psi}{\partial x}-\frac{\partial \Phi}{\partial z}=0$ (56)

5. Fluid state of Dirac Quantum Particles

By comparing Equation (56) to Dirac equations given in Equations (39) and (40), the field $\left({\psi}_{1},{\psi}_{2}\right)$ may be identified as the stream function and the velocity potential of one fluid flow and the field $\left({\psi}_{3},{\psi}_{4}\right)$ with another fluid flow. However, the main problem that we want to deal with now is whether the two fields $\left({\psi}_{1},{\psi}_{2}\right)$ and $\left({\psi}_{3},{\psi}_{4}\right)$ are connected and, most importantly, how such connection would lead to the prospect of using them to describe a Dirac quantum particle as a standing wave. In the following we will show that in fact this is the case by using the relationships between the components of these two fields given in Equations (30)-(33). For convenience we rewrite these equations as follows

$\frac{\partial {\psi}_{4}}{\partial y}=m{\psi}_{1},\frac{\partial {\psi}_{1}}{\partial y}=m{\psi}_{4}$ (57)

$\frac{\partial {\psi}_{3}}{\partial y}=-m{\psi}_{2},\frac{\partial {\psi}_{2}}{\partial y}=-m{\psi}_{3}$ (58)

$\frac{{\partial}^{2}{\psi}_{i}}{{\partial}^{2}y}-{m}^{2}{\psi}_{i}=0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}i=1,2,3,4$ (59)

If the physical quantity m, which is identified with the inertial mass of a quantum particle, is assumed to be positive, $m>0$, then it is observed that it is possible to describe the physical structure of a Dirac quantum particle as a spinning top if we consider solutions to Equation (59) as hybrid functions of the form

${\psi}_{i}=\{\begin{array}{l}{c}_{1i}\left(x,z\right){\text{e}}^{my}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}y<0\\ {c}_{2i}\left(x,z\right){\text{e}}^{-my}\text{\hspace{0.05em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}y\ge 0\end{array}$ (60)

For simplicity, instead of the hybrid form given in Equation (60), in the following we will show only for the case in which $y\ge 0$ since similar results can be obtained for which $y<0$. The solutions given in Equation (60) can be rewritten in the following forms

${\psi}_{1}={c}_{21}\left(x,z\right){\text{e}}^{-my},{\psi}_{2}={c}_{22}\left(x,z\right){\text{e}}^{-my}$ (61)

${\psi}_{3}={c}_{23}\left(x,z\right){\text{e}}^{-my},{\psi}_{4}={c}_{24}\left(x,z\right){\text{e}}^{-my}$ (62)

Using the equations given in Equations (57) and (58), we further obtain the conditions ${c}_{24}=-{c}_{21}$ and ${c}_{23}={c}_{22}$. If we write ${c}_{21}=f\left(x,z\right)$ and ${c}_{22}=g\left(x,z\right)$ then we have

${\psi}_{1}=f\left(x,z\right){\text{e}}^{-my},{\psi}_{2}=g\left(x,z\right){\text{e}}^{-my}$ (63)

${\psi}_{3}=g\left(x,z\right){\text{e}}^{-my},{\psi}_{4}=-f\left(x,z\right){\text{e}}^{-my}$ (64)

From the above forms of solutions given to the components ${\psi}_{i}$ of the wavefunction $\psi $ we can show how a standing wave can be established from the superposition of a wave associated with the field $\left({\psi}_{1},{\psi}_{2}\right)$ and a wave associated with the field $\left({\psi}_{3},{\psi}_{4}\right)$. Let ${\psi}_{1}=f\left(x,z\right){\text{e}}^{-my}$ be identified with the velocity potential and ${\psi}_{2}=g\left(x,z\right){\text{e}}^{-my}$ with the stream function of one fluid flow. Now we have two different descriptions that can be given to the field $\left({\psi}_{3},{\psi}_{4}\right)$. If we identify the component ${\psi}_{3}=g\left(x,z\right){\text{e}}^{-my}$ with the velocity potential and ${\psi}_{4}=-f\left(x,z\right){\text{e}}^{-my}$ with the stream function of another fluid flow then we have the stream function of the first flow equals the velocity potential of the second flow, and the stream function of the second flow is a reflection of the velocity of the first flow. Even though this kind of identification may be used to describe a particular type of fluid flow of Dirac quantum particles, it does not give rise to the physical structure that we are looking for, that is a standing wave. However, if we now identify the component ${\psi}_{3}=g\left(x,z\right){\text{e}}^{-my}$ with the stream function and ${\psi}_{4}=-f\left(x,z\right){\text{e}}^{-my}$ with the velocity potential of the second flow then the two flows are identical except for their flow directions, which are opposite to each other, and in fact this is what we want to obtain because they can form a required standing wave. It is also observed that for a steady motion a Dirac quantum particle does not exhibit a wave motion even though it has the potential to establish a wave within its physical structure. Therefore, without an external disturbance a Dirac quantum particle may be considered as a classical particle defined in classical physics. Furthermore, we may also speculate that the two opposite fluid flows associated with the physical structure of a Dirac quantum particle may be related to the concept of spin-half that is introduced into quantum mechanics.

6. Conclusion

In our previous works on spacetime structures of quantum particles, we suggest that all quantum particles are formed from mass points which are joined together by contact forces, which is a consequence of viewing quantum particles as CW-complexes. Being identified with differentiable manifolds, quantum particles therefore should be endowed with geometric and topological structures of differentiable manifolds and their motion should be described as isometric embeddings in higher Euclidean space. In particular, we show that quantum particles may have the geometric and topological structures of a 3D differentiable manifold which can be described as standing waves which are solutions to the Schrödinger wave equation and Dirac equation. In this work we have extended our previous discussions by showing that Dirac equation can be used to describe quantum particles as composite structures that are in a fluid state in which the components of the wavefunction can be identified with the stream function and the velocity potential of a potential flow formulated in the theory of classical fluids. With this fluid composition, physically, Dirac quantum particles can manifest as standing waves which are the result of the superposition of two fluid flows moving in opposite directions. However, for a steady motion, a quantum particle whose physical structure is constructed in terms of Dirac equation does not exhibit a wave motion even though it has the potential to establish a wave within its physical structure. Therefore, if there are no external fields acting on it, a Dirac quantum particle may be considered as a classical particle defined in classical physics. It is also noted from the fact that there are two identical fluid flows in opposite directions within their physical structures, the fluid state model of Dirac quantum particles can be invoked to explain why fermions are spin-half particles as discussed in Section 2.

Acknowledgements

We would like to thank Vanessa Cheng and Jane Gao of the administration of JMP for their editorial advice during the preparation of this work.

Conflicts of Interest

The author declares no conflicts of interest regarding the publication of this paper.

[1] | Ho, V.B. (2018) International Journal of Physics, 6, 105-115. |

[2] | Ho, V.B. (2018) GJSFR-A, 18, 37-58. |

[3] | Ho, V.B. (2018) International Journal of Physics, 6, 139-146. |

[4] | Einstein, A. (1952) The Principle of Relativity. Dover Publications, New York. |

[5] |
Hamilton, R.S. (1982) Journal of Differential Geometry, 17, 255-306. https://doi.org/10.4310/jdg/1214436922 |

[6] | Strauss, W.A. (1992) Partial Differential Equation. John Wiley & Sons, Inc., New York. |

[7] |
Cao, H.D. and Zhu, X.P. (2006) Asian Journal of Mathematics, 10, 165-492. https://doi.org/10.4310/AJM.2006.v10.n2.a2 |

[8] | Ho, V.B. (2017) A Derivation of Dirac Equation from a General System of Linear First Order Partial Differential Equations. Preprint, ResearchGate, viXra 1712.0404v1. |

[9] | Ho, V.B. (2018) Formulation of Maxwell Field Equations from a General System of Linear First Order Partial Differential Equations. Preprint, ResearchGate, viXra 1802.0055v1. |

[10] | Melshko, S.V. (2005) Methods for Constructing Exact Solutions of Partial Differential Equations. Springer Science & Business Media, Inc. |

[11] | Sobolev, S.L. (1964) Partial Differential Equations of Mathematical Physics. Dover Publications, Inc, New York. |

[12] |
Dirac, P.A.M. (1928) Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 117, 610-624. https://doi.org/10.1098/rspa.1928.0023 |

[13] | Jackson, J.D. (1975) Classical Electrodynamics. John Wiley & Sons, New York. |

[14] | Landau, L.D. and Lifshitz, E.M. (1987) The Classical Theory of Fields. Pergamon Press, Sydney. |

[15] | Ho, V.B. (2018) Topological Transformation of Quantum Dynamics. Preprint, ResearchGate, viXra 1810.0300v1. |

[16] | Bohr, N. (1913) Philosophical Magazine, 26, 1-25. |

[17] | De Broglie, L. (1925) Annales de Physique, 10, 22-128. |

[18] | Lachieze-Rey, M. and Luminet, J.P. (2003) Cosmic Topology. |

[19] | Ho, V.B. (2018) Formulation of Dirac Equation for an Arbitrary Field from a System of Linear First Order Partial Differential Equations. |

[20] | Bransden, B.H. and Joachain, C.J. (1986) Physics of Atoms and Molecules. Longman Scientific & Technical, New York. |

[21] |
Ho, V.B. (1994) Journal of Physics A: Mathematical and General, 27, 6237-6241. https://doi.org/10.1088/0305-4470/27/18/031 |

[22] |
Ho, V.B. and Morgan, M.J. (1996) Journal of Physics A: Mathematical and General, 29, 1497. https://doi.org/10.1088/0305-4470/29/7/019 |

[23] | Ho, V.B. (1996) Geometrical and Topological Methods in Classical and Quantum Physics. PhD Thesis, Monash University, Clayton. |

[24] | Balachandran, A.P., Marmo, G., Skagerstam, B.S. and Stern, A. (1991) Classical Topology and Quantum States. World Scientific, Singapore. |

[25] | Kreyszig, E. (1975) Introduction to Differential Geometry and Riemannian Geometry. University of Toronto Press, Toronto. |

[26] | Landau, L.D. and Lifshitz, E.M. (1989) Fluid Mechanics. Pergamon Press, Sydney. |

[27] | Douglas, J.F., Gasiorek, J.M. and Swaffield, J.A. (1990) Fluid Mechanics. Longman Scientific & Technical, New York. |

[28] | Batchelor, G.K. (1988) An Introduction to Fluid Dynamics. Cambridge University Press, Melbourne. |

[29] | White, F.M. (2011) Fluid Mechanics. McGraw Hill, New York. |

Copyright © 2020 by authors and Scientific Research Publishing Inc.

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.