Dirac field, van der Waals gas, Weyssenhoff fluid, Newton particle

1. Introduction

Both in geometric construction, and for its far-reaching applications, the Dirac field is among the most important fields in mathematics and physics. Still, when confronted to possible interpretations, there appears to be the spread consensus that no-one really knows what a spinor is.

To help understanding, visualizing spinor fields would be essential. But in this respect, difficulties arise because of the intrinsically complex nature they have.

On the other hand, all complex quantities may always be written in polar form, in which complex functions are re-expressed as product of modules times unitary phases, with modules and phases being real. Relativistic spinors, having two chiral components, need extra care in under-going the polar decomposition since, under Lorentz transformations, the various components would mix. Still, the relativistic polar formulation is doable just as well [1,2].

The advantage on the polar decomposition of relativistic spinor fields is that it converts the entire Dirac theory into a form that is genuinely hydrodynamic [3]. This does not only mean that all variables are real, but also that in themselves all variables are perfectly visualizable.

In fact, of the four sets of variables in terms of which the spinor field can be decomposed, two are the density and velocity, exactly as those we have in standard hydrodynamics, while another is the spin, which is also a very well-known concept of quantum mechanics. Interpreting the chiral angle is instead not yet easy, although we hope that such an object will be better understood in the light of the investigation that we intend to do in this work.

2. Dirac Field in Polar Form

2.1. Dirac Spinors

We start with a brief summary of the Dirac spinors to set our convention, and to establish the relations that we are going to need later on. So to begin, let ${\gamma }^i$ be matrices belonging to the Clifford algebra $\{{\gamma }^i,{\gamma }^j\}=2\mathbb{I}{\eta }^{ij}$ with ${\eta }^{ij}$ the Minkowski matrix. Then ${\sigma }_{ik}=[{\gamma }_i,{\gamma }_k]/4$ are defined as generators of the Lorentz group. In $2i{\sigma }_{ab}={\epsilon }_{abcd}\pi {\sigma }^{cd}$ we find the implicit definition of the parity-odd matrix $\pi$ and whose existence stipulates that the Lorentz group is reducible.1 The exponentiation of the generators gives an element of the Lorentz group $\Lambda$ and therefore $\mathbf{S}=\Lambda e^{iq\alpha }$ is an element of the spinor group accounting also for gauge transformations for which q is the charge. A spinor field is an object that transforms like $\psi \to \mathit{S}\psi$ and $\overline{\psi }\to \overline{\psi }{\mathit{S}}^{-1}$ where $\overline{\psi }={\psi }^{†}{\gamma }^0$ is the adjoint operation. With a pair of adjoint spinors we can form the spinorial bi-linears

$$\begin{array}{cc}& {\Sigma }^{ab}=2\overline{\psi }{\sigma }^{ab}\pi \psi M^{ab}=2i\overline{\psi }{\sigma }^{ab}\psi \end{array}$$

(1)

$$\begin{array}{cc}& S^a=\overline{\psi }{\gamma }^a\pi \psi U^a=\overline{\psi }{\gamma }^a\psi \end{array}$$

(2)

$$\begin{array}{cc}& \Theta =i\overline{\psi }\pi \psi \Phi =\overline{\psi }\psi \end{array}$$

(3)

which are all real tensors. They verify the Hodge duality

$$\begin{array}{cc}& {\Sigma }^{ab}=-\frac{1}{2}{\epsilon }^{abij}{M}_{ij}{M}^{ab}=\frac{1}{2}{\epsilon }^{abij}{\Sigma }_{ij}\end{array}$$

(4)

beside the constitutive relations

$$\begin{array}{cc}& M_{ik}{U}^i=\Theta S_k{\Sigma }_{ik}{U}^i=\Phi S_k\end{array}$$

(5)

$$\begin{array}{cc}& M_{ik}{S}^i=\Theta U_k{\Sigma }_{ik}{S}^i=\Phi U_k\end{array}$$

(6)

as well as

$$\begin{array}{cc}& M_{ab}\Phi -{\Sigma }_{ab}\Theta =U^j{S}^k{\epsilon }_{jkab}\end{array}$$

(7)

$$\begin{array}{cc}& M_{ab}\Theta +{\Sigma }_{ab}\Phi =U_{[a}{S}_{b]}\end{array}$$

(8)

together with

$$\begin{array}{cc}& \frac{1}{2}{M}_{ab}{M}^{ab}=-\frac{1}{2}{\Sigma }_{ab}{\Sigma }^{ab}={\Phi }^2-{\Theta }^2\end{array}$$

(9)

$$\begin{array}{cc}& \frac{1}{2}{M}_{ab}{\Sigma }^{ab}=-2\Theta \Phi \end{array}$$

(10)

and

$$\begin{array}{cc}& U_a{U}^a=-S_a{S}^a={\Theta }^2+{\Phi }^2\end{array}$$

(11)

$$\begin{array}{cc}& U_a{S}^a=0\end{array}$$

(12)

called Fierz re-arrangements. They show that not all the bi-linears are independent, and in fact if ${\Phi }^2+{\Theta }^2\ne 0$ then both antisymmetric tensors $M_{ab}$ and ${\Sigma }_{ab}$ can be dropped in favour of the two vectors and the two scalars. In turn, under the same condition, the axial-vector and the vector $S_a$ and $U_a$ are space-like and time-like, showing that they can be recognized as spin and velocity, respectively [3].

The spinorial covariant derivative is defined as

$$\begin{array}{cc}& {\nabla }_{\mu }\psi ={\partial }_{\mu }\psi +\frac{1}{2}{C}_{ab\mu }{\sigma }^{ab}\psi +iq{A}_{\mu }\psi \end{array}$$

(13)

in which $A_{\mu }$ is the gauge potential and $C_{\phantom{ab}\mu }^{ab}$ is the space-time spin connection. We are here in the torsionless case although full generality will be recovered by introducing torsion as an axial-vector field in the dynamics.

As usual, the commutator

$$\begin{array}{cc}& [{\nabla }_{\mu },{\nabla }_{\nu }]\psi =\frac{1}{2}{R}_{ab\mu \nu }{\sigma }^{ab}\psi +iq{F}_{\mu \nu }\psi \end{array}$$

(14)

defines the Riemann curvature and the Maxwell strength.

The dynamics is assigned by the torsion field equations

$$\begin{array}{cc}& {\nabla }_{\rho }{(\partial W)}^{\rho \mu }+M^2{W}^{\mu }=X{S}^{\mu }\end{array}$$

(15)

together with the gravitational field equations

$$\begin{array}{cc}\hfill & R^{\rho \sigma }-\frac{1}{2}R{g}^{\rho \sigma }-\Lambda g^{\rho \sigma }=\frac{1}{2}[\frac{1}{4}{F}^2{g}^{\rho \sigma }-F^{\rho \alpha }{F}_{\phantom{\sigma }\alpha }^{\sigma }+\\ \hfill & +\frac{1}{4}{(\partial W)}^2{g}^{\rho \sigma }-{(\partial W)}^{\sigma \alpha }{(\partial W)}_{\phantom{\rho }\alpha }^{\rho }+M^2(W^{\rho }{W}^{\sigma }-\frac{1}{2}{W}^2{g}^{\rho \sigma })+\\ & +\frac{i}{4}(\overline{\psi }{\gamma }^{\rho }{\nabla }^{\sigma }\psi -{\nabla }^{\sigma }\overline{\psi }{\gamma }^{\rho }\psi +\overline{\psi }{\gamma }^{\sigma }{\nabla }^{\rho }\psi -{\nabla }^{\rho }\overline{\psi }{\gamma }^{\sigma }\psi )-\frac{1}{2}X(W^{\sigma }{S}^{\rho }+W^{\rho }{S}^{\sigma })]\end{array}$$

(16)

and the electrodynamic field equations

$$\begin{array}{cc}& {\nabla }_{\sigma }{F}^{\sigma \mu }=q{U}^{\mu }\end{array}$$

(17)

where ${(\partial W)}_{\alpha \nu }={\nabla }_{\alpha }{W}_{\nu }-{\nabla }_{\nu }{W}_{\alpha }$ and M the torsion mass, and where we define $R_{\phantom{\alpha }\rho \alpha \sigma }^{\alpha }=R^{\rho \sigma }$ and $R^{\rho \sigma }{g}_{\rho \sigma }=R$ as the Ricci tensor and scalar and $\Lambda$ the cosmological constant.

As for matter, the dynamics is assigned in terms of the Dirac spinor field equation given by

$$\begin{array}{cc}& i{\gamma }^{\mu }{\nabla }_{\mu }\psi -X{W}_{\sigma }{\gamma }^{\sigma }\pi \psi -m\psi =0\end{array}$$

(18)

where $W_{\sigma }$ is the Hodge dual of the torsion tensor and X the torsion-spin coupling constant, added to recover full generality as we have already anticipated [4].

The set of field equations (15-16) with (17) is conceived in this way so to give rise to conservation laws that turn out to be automatically satisfied when the Dirac spinorial field equations (18) are valid, and so it is consistent.

2.2. Polar Decomposition

In the aforementioned case in which ${\Phi }^2+{\Theta }^2\ne 0$ we can perform what is called polar decomposition of the spinor field. Specifically, it is possible to demonstrate [1,2] that under the above condition any spinor field can always be written, in chiral representation, in the form

$$\begin{array}{cc}& \psi =\varphi e^{-\frac{i}{2}\beta \pi }{\mathit{L}}^{-1}\left(\begin{array}{c}\hfill 1\\ \hfill 0\\ \hfill 1\\ \hfill 0\end{array}\right)\end{array}$$

(19)

for a pair of functions $\varphi$ and $\beta$ and for some $\mathit{L}$ with the structure of a spinor transformation. As anticipated, the two antisymmetric tensors are expressed by means of the two vectors and the two scalars, and these are given by

$$\begin{array}{cc}& S^a=2{\varphi }^2{s}^a{U}^a=2{\varphi }^2{u}^a\end{array}$$

(20)

and

$$\begin{array}{cc}& \Theta =2{\varphi }^{2}sin\beta \Phi =2{\varphi }^{2}cos\beta \end{array}$$

(21)

when the polar form is implemented. The last two show that $\varphi$ and $\beta$ are a scalar and a pseudo-scalar, known as module and chiral angle. Then (11-12) reduce to

$$\begin{array}{cc}& u_a{u}^a=-s_a{s}^a=1\end{array}$$

(22)

$$\begin{array}{cc}& u_a{s}^a=0\end{array}$$

(23)

showing that the velocity has only 3 independent components, the 3 spatial rapidities, whereas the spin has only 2 independent components, the 2 angles that, in the rest-frame, its spatial part forms with the third axis. As for $\mathit{L}$ we can read its meaning as that of the specific transformation that takes a given spinor to its rest-frame with spin aligned along the third axis. For the spinorial fields in polar form, the 8 real components are re-configured in such a way that the 2 scalars $\varphi$ and $\beta$ are isolated from the 6 parameters of $\mathit{L}$ that can always be transferred into the frame and which are thus the Goldstone fields.

Because in general one can prove that

$$\begin{array}{cc}& {\mathit{L}}^{-1}{\partial }_{\mu }\mathit{L}=iq{\partial }_{\mu }\zeta \mathbb{I}+\frac{1}{2}{\partial }_{\mu }{\zeta }_{ij}{\sigma }^{ij}\end{array}$$

(24)

for some $\zeta$ and ${\zeta }_{ij}$ then we can define

$$\begin{array}{cc}& R_{ij\mu }:={\partial }_{\mu }{\zeta }_{ij}-C_{ij\mu }\end{array}$$

(25)

$$\begin{array}{cc}& P_{\mu }:=q({\partial }_{\mu }\zeta -A_{\mu })\end{array}$$

(26)

which are real tensors. By reading these expressions one can see that after the Goldstone fields are transferred into the frame, they combine with spin connection and gauge potential to become the longitudinal components of the $P_{\mu }$ and $R_{ij\mu }$ tensors, hence called gauge and space-time tensorial connections. From (19) with (26-25) we get

$$\begin{array}{cc}& {\nabla }_{\mu }\psi =({\nabla }_{\mu }ln\varphi \mathbb{I}-\frac{i}{2}{\nabla }_{\mu }\beta \pi -\frac{1}{2}{R}_{\alpha \nu \mu }{\sigma }^{\alpha \nu }-i{P}_{\mu }\mathbb{I})\psi \end{array}$$

(27)

as the polar form of the covariant derivative. Notice that

$$\begin{array}{cc}& {\nabla }_{\mu }{s}_{\nu }=s^{\alpha }{R}_{\alpha \nu \mu }{\nabla }_{\mu }{u}_{\nu }=u^{\alpha }{R}_{\alpha \nu \mu }\end{array}$$

(28)

as general identities. The covariant derivative of the velocity is the object with which one builds the strain-rate tensor in continuum mechanics. Expressions (28) are the extension to both velocity and spin of relationships that make $R_{ab\mu }$ interpretable as the strain-rate tensor.

The tensorial connections are such that

$$\begin{array}{cc}& -R_{\phantom{i}j\mu \nu }^i={\nabla }_{\mu }{R}_{\phantom{i}j\nu }^i-{\nabla }_{\nu }{R}_{\phantom{i}j\mu }^i+R_{\phantom{i}k\mu }^i{R}_{\phantom{k}j\nu }^k-R_{\phantom{i}k\nu }^i{R}_{\phantom{k}j\mu }^k\end{array}$$

(29)

$$\begin{array}{cc}& -q{F}_{\mu \nu }={\nabla }_{\mu }{P}_{\nu }-{\nabla }_{\nu }{P}_{\mu }\end{array}$$

(30)

therefore being the covariant potentials of these tensors.

In the gravitational field equations, the right-hand side aside for the factor $1/2$ is the energy density tensor, and it is expressed in polar variables according to

$$\begin{array}{cc}\hfill & T^{\rho \sigma }=\frac{1}{4}{F}^2{g}^{\rho \sigma }-F^{\rho \alpha }{F}_{\phantom{\sigma }\alpha }^{\sigma }+\frac{1}{4}{(\partial W)}^2{g}^{\rho \sigma }-{(\partial W)}^{\sigma \alpha }{(\partial W)}_{\phantom{\rho }\alpha }^{\rho }+M^2(W^{\rho }{W}^{\sigma }-\frac{1}{2}{W}^2{g}^{\rho \sigma })+\\ \hfill & +{\varphi }^2[P^{\rho }{u}^{\sigma }+P^{\sigma }{u}^{\rho }+({\nabla }^{\rho }\beta /2-X{W}^{\rho })s^{\sigma }+({\nabla }^{\sigma }\beta /2-X{W}^{\sigma })s^{\rho }-\\ & -\frac{1}{4}{R}_{\alpha \nu }^{\phantom{\alpha \nu }\sigma }{s}_{\kappa }{\epsilon }^{\rho \alpha \nu \kappa }-\frac{1}{4}{R}_{\alpha \nu }^{\phantom{\alpha \nu }\rho }{s}_{\kappa }{\epsilon }^{\sigma \alpha \nu \kappa }]\end{array}$$

(31)

where as explicit presence of the space-time tensorial connection and its coupling to the spin can be seen.

The Dirac spinor field equations in polar form are

$$\begin{array}{cc}& {\nabla }_{\mu }\beta -2X{W}_{\mu }+B_{\mu }-2P^{\iota }{u}_{[\iota }{s}_{\mu ]}+2m{s}_{\mu }cos\beta =0\end{array}$$

(32)

$$\begin{array}{cc}& {\nabla }_{\mu }ln{\varphi }^2+R_{\mu }-2P^{\rho }{u}^{\nu }{s}^{\alpha }{\epsilon }_{\mu \rho \nu \alpha }+2m{s}_{\mu }sin\beta =0\end{array}$$

(33)

in which $R_{\mu \nu }^{\phantom{\mu \nu }\nu }=R_{\mu }$ and $\frac{1}{2}{\epsilon }_{\mu \alpha \nu \iota }{R}^{\alpha \nu \iota }=B_{\mu }$ were defined.

Upon the introduction of the potentials

$$\begin{array}{cc}& 2Y_{\mu }={\nabla }_{\mu }\beta -2X{W}_{\mu }+B_{\mu }\end{array}$$

(34)

$$\begin{array}{cc}& 2Z_{\mu }={\nabla }_{\mu }ln{\varphi }^2+R_{\mu }\end{array}$$

(35)

it becomes easier to work the polar spinor field equations (32-33) in order to isolate the gauge tensorial connection

$$\begin{array}{cc}& P^{\eta }=mcos\beta u^{\eta }+Y_{\mu }{u}^{[\mu }{s}^{\eta ]}+Z_{\mu }{u}_{\pi }{s}_{\tau }{\epsilon }^{\mu \pi \tau \eta }\end{array}$$

(36)

which is recognized to be the momentum of the field and with which the energy (31) acquires the form

$$\begin{array}{cc}\hfill & T^{\rho \sigma }=\frac{1}{4}{F}^2{g}^{\rho \sigma }-F^{\rho \alpha }{F}_{\phantom{\sigma }\alpha }^{\sigma }+\frac{1}{4}{(\partial W)}^2{g}^{\rho \sigma }-{(\partial W)}^{\sigma \alpha }{(\partial W)}_{\phantom{\rho }\alpha }^{\rho }+M^2(W^{\rho }{W}^{\sigma }-\frac{1}{2}{W}^2{g}^{\rho \sigma })+\\ \hfill & +{\varphi }^2[2mcos\beta u^{\rho }{u}^{\sigma }-2Y_{\mu }{s}^{\mu }{u}^{\rho }{u}^{\sigma }+Y_{\mu }{u}^{\mu }(s^{\rho }{u}^{\sigma }+s^{\sigma }{u}^{\rho })+Y^{\rho }{s}^{\sigma }+Y^{\sigma }{s}^{\rho }+\\ \hfill & +Z_{\mu }{u}_{\pi }{s}_{\tau }({\epsilon }^{\mu \pi \tau \sigma }{u}^{\rho }+{\epsilon }^{\mu \pi \tau \rho }{u}^{\sigma })-\frac{1}{4}(R_{\alpha \nu \pi }{\epsilon }^{\rho \alpha \nu \pi }{g}^{\sigma \kappa }+R_{\alpha \nu \pi }{\epsilon }^{\sigma \alpha \nu \pi }{g}^{\rho \kappa }+\\ & +R_{\alpha \nu }^{\phantom{\alpha \nu }\sigma }{\epsilon }^{\rho \alpha \nu \kappa }+R_{\alpha \nu }^{\phantom{\alpha \nu }\rho }{\epsilon }^{\sigma \alpha \nu \kappa }\left)s_{\kappa }\right]\end{array}$$

(37)

in terms of the $R_{ab\mu }$ tensor and the $Y_{\mu }$ and $Z_{\mu }$ potentials.

3. Torsion Effective Approximation and van der Waals gas

3.1. General Thermodynamic Variables

When combining the two principles of thermodynamics into the single relation $dU=TdS-pdV$ and considering that $dS$ is an exact differential form, one can extract

$$\begin{array}{c}\hfill {\left(\frac{\partial U}{\partial V}\right)}_T={\left(\frac{\partial p}{\partial T}\right)}_{V}T-p\end{array}$$

(38)

in case V and T are the independent variables. With this equation, called internal energy equation, one can deduce the internal energy once the equation of state is assigned.

For example, take the simplest non-perfect gas, that is the van der Waals gas, whose equation of state is

$$\begin{array}{c}\hfill \left(p+\frac{a}{V^2}\right)(V-b)=RT\end{array}$$

(39)

in which a is a constant related to the effective pressure due to forces between the molecules, positive in the case of attraction, and b is the volume that is occupied by the molecules. By means of (38) one can deduce that

$$\begin{array}{c}\hfill U=\omega +C_{V}T-\frac{a}{V}\end{array}$$

(40)

where $\omega$ is a generic constant.

With the equation of the internal energy we are giving an axiomatic definition of thermodynamical variables in the sense that we are assigning a meaning to the different terms entering (38) according to the role they play in such an equation. For example, if we knew that U had a given dependence on V then the right-hand side of (38) would be known, and any pair of variables satisfying the right-hand side of (38) in exactly the way p and T are would respectively be interpreted as pressure and temperature.

With this in mind, we are now going to investigate the thermodynamic structure of the Dirac spinor field theory.

3.2. Massive Propagating Torsion

We will consider the Dirac theory with torsion taken to be massive enough to allow the effective approximation.

In effective approximation, the torsion field loses all its propagating properties, with field equations reducing to

$$\begin{array}{cc}& M^2{W}^{\mu }=X{S}^{\mu }\end{array}$$

(41)

so that now torsion can be replaced in terms of the spin.

When this is done in the expression of the energy density tensor (37) remarkable simplifications occur. Taking in particular the purely spinorial contribution, it reads

$$\begin{array}{cc}\hfill & E^{\rho \sigma }={\varphi }^2[2(mcos\beta -{\varphi }^2{X}^2/M^2-s^{\mu }{\nabla }_{\mu }\beta /2-\frac{1}{4}{R}^{\pi \tau \eta }{s}^{\kappa }{\epsilon }_{\kappa \pi \tau \eta })u^{\rho }{u}^{\sigma }+\\ \hfill & +2{\varphi }^2{X}^2/M^2(g^{\rho \sigma }-u^{\rho }{u}^{\sigma })+(s^{\rho }{u}^{\sigma }+s^{\sigma }{u}^{\rho })u^{\mu }{\nabla }_{\mu }\beta /2+s^{\rho }{\nabla }^{\sigma }\beta /2+s^{\sigma }{\nabla }^{\rho }\beta /2+\\ \hfill & +Z_{\mu }{u}_{\pi }{s}_{\tau }({\epsilon }^{\mu \pi \tau \sigma }{u}^{\rho }+{\epsilon }^{\mu \pi \tau \rho }{u}^{\sigma })-\frac{1}{4}{R}_{\pi \tau \eta }{s}_{\kappa }({\epsilon }^{\rho \pi \tau \kappa }{g}^{\sigma \eta }+{\epsilon }^{\sigma \pi \tau \kappa }{g}^{\rho \eta }+\\ & +{\epsilon }^{\pi \tau \eta \mu }{u}_{\mu }{u}^{\sigma }{g}^{\rho \kappa }+{\epsilon }^{\pi \tau \eta \mu }{u}_{\mu }{u}^{\rho }{g}^{\sigma \kappa }\left)\right]\end{array}$$

(42)

which can be worked out in detail in the following way.

Defining the quantities given by

$$\begin{array}{cc}& \mu =E_{\rho \sigma }{u}^{\rho }{u}^{\sigma }\end{array}$$

(43)

$$\begin{array}{cc}& p=-\frac{1}{3}{E}_{\rho \sigma }(g^{\rho \sigma }-u^{\rho }{u}^{\sigma })\end{array}$$

(44)

it is easy to prove that

$$\begin{array}{cc}& E^{\rho \sigma }=\mu u^{\rho }{u}^{\sigma }-p(g^{\rho \sigma }-u^{\rho }{u}^{\sigma })\end{array}$$

(45)

as a general identity. Re-writing the energy density with this structure helps identifying the quantities $\mu$ and p as the internal energy density and the pressure of the field.

As a consequence, in our case we have that

$$\begin{array}{cc}& \mu =2{\varphi }^2(mcos\beta -{\varphi }^2{X}^2/M^2)-\left[2{\varphi }^2(s^{\mu }{\nabla }_{\mu }\beta /2-\frac{1}{2}{\epsilon }^{\kappa \alpha \mu \nu }{s}_{\kappa }{u}_{\alpha }{\nabla }_{\mu }{u}_{\nu })\right]\end{array}$$

(46)

$$\begin{array}{cc}& p=-2{\varphi }^4{X}^2/M^2-\frac{1}{3}\left[2{\varphi }^2(s^{\mu }{\nabla }_{\mu }\beta /2-\frac{1}{2}{\epsilon }^{\kappa \alpha \mu \nu }{s}_{\kappa }{u}_{\alpha }{\nabla }_{\mu }{u}_{\nu })\right]\end{array}$$

(47)

are the internal energy density and pressure of spinors.

Introducing the volume $2{\varphi }^2=1/V$ as well as the internal energy $U=\mu V$ we can write them respectively as

$$\begin{array}{c}U=mcos\beta +3RT-\frac{X^2}{2M^2}\frac{1}{V}\end{array}$$

(48)

$$\begin{array}{c}\left(p+\frac{X^2}{2M^2}\frac{1}{V^2}\right)V=RT\end{array}$$

(49)

in which

$$\begin{array}{cc}& 3RT=-s^{\mu }{\nabla }_{\mu }\beta /2+\frac{1}{2}{\epsilon }^{\kappa \alpha \mu \nu }{s}_{\kappa }{u}_{\alpha }{\nabla }_{\mu }{u}_{\nu }\end{array}$$

(50)

has also been defined. Notice that (49) is exactly the van der Waals equation of state in the case in which $b=0$ and $2a=X^2/M^2$ showing that the torsional effective force is indeed attractive. Also notice that (48) can be recognized as the van der Waals gas internal energy if $C_V=3R$ and $m=\omega$ for small values of the chiral angle. It is essential to remark that in order for (49-48) to be structurally similar to those of a van der Waals gas condition (50) must hold.

The valitidy of (50) can be interpreted as the definition of temperature for the Dirac field, and it can be read as the fact that the internal dynamics of the Dirac field gets contributions from its chiral angle and its vorticity. It is not surprising that the chiral angle, the phase difference between the chiral parts, be tied to the internal dynamics, and thus thermodynamically associated to the concept of temperature. Such an association is also clear in the fact that $mcos\beta$ is another contribution of the chiral angle to the internal energy, which is just the relativistic mass.

We recall to the reader that the association of the chiral angle to temperature, while justified by an interpretation employing the conpcet of internal dynamics, is only the axiomatic type of connection in the sense explained here above. The definition of temperaure assigned by means of the internal energy and its equation (38) is formal and not functional. We have defined T according to (50) with the aim of rendering (38) satisfied but T does not represent a chaotic motion of particles in the kinetic theory of gases.

The definition of temperature as given by (50) seems to us the only way to define something conceptually close to the idea of temperature even for systems that are not constituted by randomly distributed particles.

4. Zero Chiral Angle and Weyssenhoff Fluid

4.1. Non-Relativistic Regime

In [3] and references therein, we have discussed the idea of non-relativistic limit as the regime for which

$$\begin{array}{c}\overrightarrow{u}\to 0\end{array}$$

(51)

$$\begin{array}{c}\beta \to 0\end{array}$$

(52)

characterizing the difference between the two conditions in the fact that, while the first represents the lost motion, the second represents the loss of the dynamical properties that would remain even in rest-frame, thus the intrinsic, internal dynamics. This fits well in the discussion above, where it is even more reasonably justified the fact that, in non-relativistic regime, the temperature (50) would lose all contributions coming from the material distribution.

Therefore, while the pair of conditions (51-52) are the non-relativistic limit, the single condition (52) is thought as internal triviality. Or in other words, when the chiral angle vanishes we lose the internal dynamics. This is also reasonable if we consider that $\beta =0$ means no difference between the two chiral parts. Or that the zitterbewegung effect vanishes as it was discusses in references [5,6].

The condition of internal triviality has also the advantage of being covariant, so it does make sense to see what is going to happen when it is consistently assumed.

4.2. Hydrodynamics with Spin

Assuming $\beta =0$ from the start implies that the bi-linear pseudo-scalar $\Theta =0$ identically. Hence

$$\begin{array}{cc}& M_{ik}{u}^i=0{\Sigma }_{ik}{u}^i=2{\varphi }^2{s}_k\end{array}$$

(53)

$$\begin{array}{cc}& M_{ik}{s}^i=0{\Sigma }_{ik}{s}^i=2{\varphi }^2{u}_k\end{array}$$

(54)

alongside to

$$\begin{array}{cc}& M_{ab}=2{\varphi }^2{u}^j{s}^k{\epsilon }_{jkab}\end{array}$$

(55)

$$\begin{array}{cc}& {\Sigma }_{ab}=2{\varphi }^2{u}_{[a}{s}_{b]}\end{array}$$

(56)

and

$$\begin{array}{cc}& M_{ab}{M}^{ab}=-{\Sigma }_{ab}{\Sigma }^{ab}=8{\varphi }^4\end{array}$$

(57)

$$\begin{array}{cc}& M_{ab}{\Sigma }^{ab}=0\end{array}$$

(58)

as Fierz identities. By employing (4) into (54) one has

$$\begin{array}{cc}& M_{ik}{u}^i=0\frac{1}{2}{\epsilon }_{kiab}{M}^{ab}{u}^i=2{\varphi }^2{s}_k\end{array}$$

(59)

$$\begin{array}{cc}& M_{ik}{s}^i=0\frac{1}{2}{\epsilon }_{kiab}{M}^{ab}{s}^i=2{\varphi }^2{u}_k\end{array}$$

(60)

so that focusing in particular on the first, we can re-write the two expressions according to

$$\begin{array}{cc}& M^{ki}{u}_i=0\end{array}$$

(61)

$$\begin{array}{cc}& M^{[ab}{u}^{c]}={\epsilon }^{abck}{S}_k\end{array}$$

(62)

telling that the momentum is orthogonal to the velocity and that the completely antisymmetric part of $M^{ij}{u}^k$ is the Hodge dual of the spin density axial-vector. It follows that the momentum $M^{ki}$ is the fundamental spin tensor of Weyssenhoff fluids [7,8], where (61) is the consitutive condition and (62) the condition saying that it is only the completely antisymmetric part of the spin density tensor that is excited. This is expected as the Dirac spinor has a completely antisymmetric spin density tensor.

5. Spinlessness and Newton Mechanics

5.1. Classical Limit

At last, we discuss the case of spinlessness. Such a case is obtained in the approximation

$$\begin{array}{c}{s}^a\to 0\end{array}$$

(63)

and it means that we are losing quantum effects. Indeed, if we were not to choose natural units, the spin would be seen to be proportional to ℏ and the limit $\hslash \to 0$ is what would give rise to the classical approximation condition.

Notice also that the validity of the Dirac equation gives

$$\begin{array}{c}{\nabla }_i{S}^i=2m\Theta \end{array}$$

(64)

showing that $\beta \to 0$ is implied by $S^i\to 0$ and stating that there can be no chirality if there is no helicity.

The present limit is therefore compatible with the limit that we discussed in the previous section.

5.2. Point Particle

Let us then re-consider the momentum (36) as well as the energy density tensor (37) in effective approximation and in this limit. We have

$$\begin{array}{cc}& P^{\eta }=(m-2{\varphi }^2{X}^2/M^2)u^{\eta }\end{array}$$

(65)

and

$$\begin{array}{cc}& T^{\rho \sigma }=\frac{1}{4}{F}^2{g}^{\rho \sigma }-F^{\rho \alpha }{F}_{\phantom{\sigma }\alpha }^{\sigma }+2{\varphi }^2(m-2{\varphi }^2{X}^2/M^2)u^{\rho }{u}^{\sigma }+2{\varphi }^4{X}^2/M^2{g}^{\rho \sigma }\end{array}$$

(66)

which next we discuss in view of their conservation laws.

To this purpose set $2{\varphi }^2=\rho$ being $\rho$ the density distribution of the material field. The above become

$$\begin{array}{cc}& P^{\eta }=(m-\rho X^2/M^2)u^{\eta }\end{array}$$

(67)

and

$$\begin{array}{cc}& T^{\rho \sigma }=\frac{1}{4}{F}^2{g}^{\rho \sigma }-F^{\rho \alpha }{F}_{\phantom{\sigma }\alpha }^{\sigma }+\rho (m-\rho X^2/M^2)u^{\rho }{u}^{\sigma }+\frac{1}{2}{\rho }^2{X}^2/M^2{g}^{\rho \sigma }\end{array}$$

(68)

and for them we know that

$$\begin{array}{cc}& {\nabla }_{\alpha }{T}^{\alpha \nu }=0\end{array}$$

(69)

and

$$\begin{array}{cc}& {\nabla }_{\alpha }\left(\rho u^{\alpha }\right)=0\end{array}$$

(70)

must be valid as a consequence of the Dirac spinorial field equations. Taking (67) into (68) and the result into (69) and then employing (70) we arrive at

$$\begin{array}{cc}& \frac{1}{2}{F}_{\alpha \pi }{\nabla }^{\sigma }{F}^{\alpha \pi }+F_{\alpha \pi }{\nabla }^{\pi }{F}^{\sigma \alpha }-{\nabla }^{\eta }{F}_{\eta \alpha }{F}^{\sigma \alpha }+\rho u^{\nu }{\nabla }_{\nu }{P}^{\sigma }-{\nabla }^{\sigma }p=0\end{array}$$

(71)

where the pressure $p=-\frac{1}{2}{\rho }^2{X}^2/M^2$ was used.

By employing now the Maxwell equations (17) we get

$$\begin{array}{cc}& \rho u^{\nu }{\nabla }_{\nu }{P}^{\sigma }={\nabla }^{\sigma }p+q\rho F^{\sigma \alpha }{u}_{\alpha }\end{array}$$

(72)

which is the Newton equation of hydrodynamic motion.

In total absence of torsion, no pressure remains so that it becomes possible to simplify the density on both sides and we reduce to the final

$$\begin{array}{cc}& u^{\nu }{\nabla }_{\nu }{P}^{\sigma }=q{F}^{\sigma \alpha }{u}_{\alpha }\end{array}$$

(73)

as the Newton equation for the motion of material points.

It is important to remark that the Newton law has been obtained without any assumption on localization for the matter distribution. With this we do not mean to imply that matter distributions cannot be localized, but rather that there is no need for this assumption at this stage.

6. Conclusions

In this work, we have considered the Dirac spinor field theory re-formulated in terms of the polar variables given by the $\varphi$ and $\beta$ scalars with the $u_a$ and $s_a$ vectors. After conversion, the full relativistic quantum mechanics turns into a type of hydrodynamics in which $2{\varphi }^2$ is the density distribution and $\beta$ the chiral angle while $u_a$ is the velocity and $s_a$ is the spin. This hydrodynamics is, therefore, an extension of the usual one since not only the density and velocity, but also the chiral angle and spin are present.

However, the general construction can be restricted to the standard hydrodynamics by removing these two extra variables. The general theory, with torsion in its effective approximation, has the same thermodynamic features of a van der Waals gas, with van der Waals pressure due to torsion, always negative since torsion is always attractive, and with temperature and internal energy being tied to the chiral angle. In the limit $\beta \to 0$ (corresponding to the requirement of losing the phase difference between chiral parts) the general theory reduces to that of a Weyssenhoff fluid with completely antisymmetric spin. And for $s_a\to 0$ (corresponding to the condition of non-quantum limit) it reduces to a Newton fluid in presence of pressure due to torsion. By vanishing torsion the usual Newton equation for the motion of material points is eventually recovered.

Aside from allowing us to see that torsion is a form of pressure and that the chiral angle can be interpreted like a type of temperature, the polar re-formulation of spinors allows the relativistic quantum mechanics to convert into a specific hydrodynamics, whose variables may perfectly be visualized, and because of this, better understood.

The challanges of relativistic quantum mechanics have no resolution in a reformulation of the theory alone, and many questions remain still open. Nonetheless, questions that can be answered more easily when made clearer will receive a boost by a Dirac spinor field theory formulated in terms of variables that are visualizable.

In its polar form, the Dirac theory is precisely this.