Quaternion Algebra on 4D Spacetime: The Lorentz Gauge Is a Door to the Invisible World of Dark Matter and Dark Energy

1. Introduction

The cosmological principle states that the spatial distribution of matter in the universe is homogeneous and isotropic when viewed on a sufficiently large scale. In this view, let us first evaluate the mass density of the visible universe. With its diameter d equal to about 28.5 gigaparsecs [1], or $8.8\times {10}^{26}\mathrm{m}$, its volume is as follows $\Delta V=(\pi /6)d^3\approx 3.57\times {10}^{80}{\mathrm{m}}^3$. Hereafter we will adhere to SI units. The total mass of visible matter in the universe according to different sources [2,3,4] is about ${10}^{53}$ to $2\times {10}^{53}$ kg. We take the average $M=1.5\times {10}^{53}$ kg. From here we may evaluate the mass density of the visible matter filling the visible universe

$${\rho }_{\mathrm{visible}}=\frac{M}{\Delta V}\approx 5.6\times {10}^{-28}\mathrm{kg}·{\mathrm{m}}^{-3}.$$

(1)

Planck Collaboration [5], Planck Collaboration [6], Planck Collaboration [7] has shown from the measurements of the CMB anisotropies a good consistency with the standard spatially flat 6-parameter $\Lambda$CDM cosmology. It follows from here that one can take Euclidean flat geometry, where the average density is close to the critical density [8,9]:

$${\rho }_{\mathrm{crit}}=\frac{3}{8\pi G}{H}_0^2\approx {10}^{-26}\mathrm{kg}·{\mathrm{m}}^{-3}.$$

(2)

Here $G\approx 6.7\times {10}^{-11}{\mathrm{m}}^3{\mathrm{s}}^{-2}{\mathrm{kg}}^{-1}$ is the gravitational constant. The Hubble constant $H_0$ expressed in units adopted in astrophysics [10] is about $73\mathrm{km}·{\mathrm{Mpc}}^{-1}·{\mathrm{s}}^{-1}$. Or it is $2.36\times {10}^{-18}{\mathrm{s}}^{-1}$ being represented in SI unit [11]. Now, we can estimate the presence of visible matter and the matter composing the critical density (dark matter and energy) in the observed universe as a percentage

$$\frac{{\rho }_{\mathrm{crit}}}{{\rho }_{\mathrm{visible}}+{\rho }_{\mathrm{crit}}}·100\%\approx 96\%$$

(3)

So, the visible (baryon) matter is only about $4\%$ in the visible universe. The remaining $96\%$ are the dark mass and dark energy. It is exactly what our instruments do not detect, except that the dark matter has gravity.

There is a reason to believe that the dark matter and dark energy are basically Bose-Einstein condensates [8,12,13,14,15,16,17] representing the superfluid quantum medium [18,19,20,21]. We call it as the superfluid quantum space-time [22,23], since this medium fills all the space everywhere densely, being the physical vacuum. Albareti et al. [24,25] consider that the large-scale structure of physical vacuum may manifest very much in the same way as cold dark matter does. One can notice that the above medium is an updated version of the ether, which was unfairly rejected but then accepted again, thanks in particular to the authorities of Einstein [26] and Dirac [27].

Qualitatively, the distribution of matter in the universe can be represented in pie charts, as shown in Figure 1.

The first pie chart (a) in Figure 1 shows the presence of dark energy, dark matter, and visible matter that make up the universe. It can be seen that, against the background of dark energy and dark matter, visible matter is only a small part. Figuratively speaking, if dark matter and energy are the bulk of the boundless ocean, then visible matter appears to be a fleeting foam on it.

The following pie chart (b) shown in Figure 1 demonstrates the presence of various atomic elements within this fleeting baryon foam. It can be seen that these elements are mainly hydrogen and helium. Only hydrogen, the carrier of protons, and helium, the carrier of both protons and neutrons, account for up to 98% of all visible matter, which washes the entire universe in a thin layer. The remaining two percent of visible matter is accounted for by all other atoms heavier than helium.

As for the third diagram represented in Figure 1, it shows what chemical and biochemical components the living brain consists of. This group includes such elements as sodium, potassium, nitrogen, carbon, calcium, oxygen, iron, zinc, selenium, magnesium, and other higher elements that are the basis of living matter. As seen, the organization of the living brain needs many chemicals, which are a tiny percentage of visible matter. It means that living, intelligent matter is an exceptionally rare phenomenon in the universe. But it is this intelligent matter that realizes the greatness of the universe.

The article is organized as follows: In Sec. 2 we apply the quaternion generating operators to the energy-momentum tensor written on the same quaternion basis. Note that the energy-momentum tensor is defined up to two scalar fields. One scalar field, the massless field, relates to the EM field. The other field, the mass-loaded field, deals with the mechanical part of the tensor. To arrive at the gravitomagnetic equations, we should accept the Lorentz gauge as cutting off these scalar fields. Our attention here is on the invisible world that arises due to applying the Lorentz gauge. We declare in sec. 3 that these scalar fields deal with dark energy and dark matter. The manifestation of dark matter will be confirmed by its contribution to the stabilization of the rotational velocities of spiral galaxy arms. With that aim in mind, we first derive, in sec. 4, the gravitmagnetic equations. Further, we get the Navier-Stoces equation in sec. 5, and next we obtain the vorticity equation and its solutions for the orbital speed, subsets. 5.1 and 5.2. The result is a flat profile of orbital speed stabilized by applied dark matter chunks. Concluding remarks are given in sec. 6.

2. Quaternion Framework of the Space-Time Continuum

In mathematics, quaternions $\mathbf{1},\mathbf{i},\mathbf{j},\mathbf{k}$ are a commutative number system that extend complex numbers. Quaternions were first described in print by Linda Rodriguez about in 1840 and independently discovered by the Irish mathematician Sir William Rowan Hamilton [28] in 1843 and applied to mechanics in three-dimensional space. Here is what James Clerk Maxwell [29], the discoverer of electromagnetic theory, wrote about the discovery of the quaternion calculus: “The invention of the calculus of quaternions is a step towards the knowledge of quantities related to space which can only be compared, for its importance, with the invention of triple coordinates by Descartes. The ideas of this calculus, as distinguished from its operations and symbols, are fitted to be of the greatest use in all parts of science.” There is, however, an amazing compatibility between the coordinate framework of 4D spacetime and the quaternions, $\mathbf{1},\mathbf{i},\mathbf{j},\mathbf{k}$ (${\mathbf{i}}^{\mathbf{2}}={\mathbf{j}}^{\mathbf{2}}={\mathbf{k}}^{\mathbf{2}}=\mathbf{ijk}=-\mathbf{1}$) described by William Rowan Hamilton. Nevertheless, instead of these numbers, further we introduce the quaternion basis represented by four $4\times 4$ matrices with real-valued components [23,30,31,32,33,34,35,36,37,38]:

$$\begin{array}{ccc}& & {\eta }_0=\left(\begin{array}{cccc}1& 0& 0& 0\\ 0& 1& 0& 0\\ 0& 0& 1& 0\\ 0& 0& 0& 1\end{array}\right),{\eta }_z=\left(\begin{array}{cccc}0& 0& 0& -1\\ 0& 0& 1& 0\\ 0& -1& 0& 0\\ 1& 0& 0& 0\end{array}\right),\hfill \\ & & {\eta }_x=\left(\begin{array}{cccc}0& -1& 0& 0\\ 1& 0& 0& 0\\ 0& 0& 0& 1\\ 0& 0& -1& 0\end{array}\right),{\eta }_y=\left(\begin{array}{cccc}0& 0& -1& 0\\ 0& 0& 0& -1\\ 1& 0& 0& 0\\ 0& 1& 0& 0\end{array}\right).\hfill \end{array}$$

(4)

These four $4\times 4$-matrices, ${\eta }_0$, ${\eta }_x$, ${\eta }_y$, ${\eta }_z$, are isomorphic to four $2\times 2$ matrices ${\sigma }_0$, ${\sigma }_x$, ${\sigma }_y$, ${\sigma }_z$ [39,40,41,42,43,44]. The latter matrices belonging to the group SU(2) are the main elements used by Roger Penrose at formulations of his twistor theory [45,46,47,48].

The matrices ${\eta }_0$, ${\eta }_x$, ${\eta }_y$, ${\eta }_z$ obey the following multiplication rules:

$${\eta }_x^2={\eta }_y^2={\eta }_z^2={\eta }_z{\eta }_y{\eta }_x=-{\eta }_0.$$

(5)

All mathematical relations of the quaternion group can be obtained from these algebraic equations. Note that the orders of multiplication $\mathbf{ijk}=-\mathbf{1}$ and ${\eta }_z{\eta }_y{\eta }_x=-{\eta }_0$ are opposite to each other. This chirality violation was introduced when studying the passage of polarized neutrons through magnetic fields [30], since the correspondence of the direction of the magnetic field and the spin of the neutron was observed. The recovery of symmetry is achievable by changing the sign at any matrix ${\eta }_x$, ${\eta }_y$, or ${\eta }_z$.

Let us now define the generating differential operators [23,35,36,37,49], as follows:

$$\begin{array}{ccc}\hfill \mathcal{D}& =& \mathbf{i}{c}^{-1}{\partial }_t{\eta }_0+{\partial }_x{\eta }_x+{\partial }_y{\eta }_y+{\partial }_z{\eta }_z,\hfill \end{array}$$

(6)

$$\begin{array}{ccc}\hfill {\mathcal{D}}^{\mathbf{T}}& =& \mathbf{i}{c}^{-1}{\partial }_t{\eta }_0^{\mathbf{T}}+{\partial }_x{\eta }_x^{\mathbf{T}}+{\partial }_y{\eta }_y^{\mathbf{T}}+{\partial }_z{\eta }_z^{\mathbf{T}}\hfill \\ \hfill & =& \mathbf{i}{c}^{-1}{\partial }_t{\eta }_0-{\partial }_x{\eta }_x-{\partial }_y{\eta }_y-{\partial }_z{\eta }_z,\hfill \end{array}$$

(7)

which are needed for reproducing the space-time dynamics of a medium described by the energy-momentum density. In particular, the D’Alembert wave operator for the case of the negative metric signature $\{-1,+1,+1,+1\}$ reads

$${\mathcal{D}}^{\mathbf{T}}\mathcal{D}=\mathcal{D}{\mathcal{D}}^{\mathbf{T}}=(-c^{-2}{\partial }_t^2+{\partial }_x^2+{\partial }_y^2+{\partial }_z^2){\eta }_0=\square$$

(8)

Now let us turn to the consideration of the medium evolving in this 4D space-time. We begin with terms $ϵ$ and $\overrightarrow{p}$ representing the kinetic energy ${ϵ}_0$ and momentum ${\overrightarrow{p}}_0$ of matter enveloped by the electromagnetic potentials $\Phi$ and $\overrightarrow{A}$:

$$\left\{\begin{array}{c}ϵ={ϵ}_0-{\rho }_q\Phi =\rho ·(\gamma \left(\mathit{v}\right)m{c}^2-e\Phi ).\\ \\ \overrightarrow{p}={\overrightarrow{p}}_0-{\rho }_{q}c\overrightarrow{A}=c\rho ·(\gamma \left(\mathit{v}\right)m\overrightarrow{\mathit{v}}-e\overrightarrow{A}),\end{array}\right.$$

(9)

Here $\gamma \left(\mathit{v}\right)={(1-{\mathit{v}}^2/c^2)}^{-1/2}$ is the Lorentz factor. At the non-relativistic limit, $\mathit{v}\to 0$, it tends to the unit. The density $\rho$ in these equations multiplied by m or by e (mass and charge of the unit carrier – particle) returns the mass density ${\rho }_M$ and the charge density ${\rho }_q$, respectively. Here $M=mN$ and $q=eN$ are the rough mass and charge of the material body containing N particles within the elementary volume $\Delta V$ under consideration. Note that the particles within this volume make chaotic walks, obeying the uncertainty principle at collisions [50,51]. In the first approximation, these walks permit description by the Wiener process [52,53]. As a result, a blurred cloud of the particles arising permits us to describe their evolution by a wave function.

Now, we introduce two gauge invariant scalar fields [36]. One scalar field, ${\psi }_0$, relates to the massless electromagnetic medium loaded by the electromagnetic potentials $\Phi$ and $\overrightarrow{A}$. Carriers here are EM quanta with zero rest mass. The other scalar field, ${\varphi }_1$, deals with the mass loaded medium which has the kinetic energy ${ϵ}_0$ and momentum ${\overrightarrow{p}}_0$ describing a state of this medium. The other scalar field, ${\varphi }_1$, deals with the mass-loaded medium. It has the kinetic energy ${ϵ}_0$ and momentum ${\overrightarrow{p}}_0$ describing the state of this medium. Let us rewrite Eq. (9) in light of the above said

$$\begin{array}{ccc}\hfill \widetilde{ϵ}& \leftarrow & \stackrel{ϵ}{\overbrace{({ϵ}_0-{\rho }_e\Phi )}}-c^{-1}{\partial }_t({\varphi }_1-{\varphi }_0),\hfill \\ \hfill \\ \hfill \tilde{\overrightarrow{p}}& \leftarrow & \underset{\overrightarrow{p}}{\underbrace{({\overrightarrow{p}}_0-{\rho }_{e}c\overrightarrow{A})}}+\nabla ({\varphi }_1-{\varphi }_0).\hfill \end{array}$$

(10)

The energy-momentum tensor written in the quaternion basis looks as follows

$$T=\mathbf{i}ϵ{\eta }_0+p_x{\eta }_x+p_y{\eta }_y+p_z{\eta }_z-{\mathcal{D}}^{\mathbf{T}}\varphi .$$

(11)

Here $ϵ$ and $\overrightarrow{p}=(p_x,p_y,p_z)$ are the generalized energy and momentum densities (9) and $\varphi ={\varphi }_0-{\varphi }_1$ is an arbitrary scalar field, having dimensionality of $\mathrm{Joule}/{\mathrm{m}}^2=\mathrm{kg}/{\mathrm{s}}^2$. Observe that, $\varphi$ has the same dimension, applied by Fedi [54] for the explanation of the ether being the Newtonian fluid manifesting this effect at the relativistic speeds.

Let us now apply the generating operator $\mathcal{D}$ to the energy-momentum tensor (11):

$$\begin{array}{ccc}\hfill \mathcal{D}T& =& \left(-\frac{1}{c}\frac{\partial ϵ}{\partial t}-\frac{\partial p_x}{\partial x}-\frac{\partial p_y}{\partial y}-\frac{\partial p_z}{\partial z}-\square \varphi \right){\eta }_0\hfill \\ \hfill & +& \left(\mathbf{i}\left(\frac{\partial ϵ}{\partial x}+\frac{\partial p_x}{c\partial t}\right)-\left(\frac{\partial p_z}{\partial y}-\frac{\partial p_y}{\partial z}\right)\right){\eta }_x\hfill \\ \hfill & +& \left(\mathbf{i}\left(\frac{\partial ϵ}{\partial y}+\frac{\partial p_y}{c\partial t}\right)-\left(\frac{\partial p_x}{\partial z}-\frac{\partial p_z}{\partial x}\right)\right){\eta }_y\hfill \\ & +& \left(\mathbf{i}\left(\frac{\partial ϵ}{\partial z}+\frac{\partial p_z}{c\partial t}\right)-\left(\frac{\partial p_y}{\partial x}-\frac{\partial p_x}{\partial y}\right)\right){\eta }_z\hfill \end{array}$$

(12)

We proclaim that the multiplier of the unit quaternion ${\eta }_0$ returns the Lorentz gauge, namely:

$$\frac{1}{4}\mathrm{trace}\mathcal{D}T=-\frac{1}{c}\frac{\partial ϵ}{\partial t}-\frac{\partial p_x}{\partial x}-\frac{\partial p_y}{\partial y}-\frac{\partial p_z}{\partial z}-\square \varphi =0.$$

(13)

So that only the expressions of the quaternions ${\eta }_x$, ${\eta }_y$, ${\eta }_z$ remain, which do not contain terms with $\varphi$. As a result, we may rewrite equation (12) in the following manner

$${\mathbb{F}}_{\mathcal{E}\mathcal{B}}=\mathcal{D}·T=\left(\begin{array}{ccccc}& 0& {\mathcal{B}}_x-\mathbf{i}{\mathcal{E}}_x& {\mathcal{B}}_y-\mathbf{i}{\mathcal{E}}_y& {\mathcal{B}}_z-\mathbf{i}{\mathcal{E}}_z\\ & -{\mathcal{B}}_x+\mathbf{i}{\mathcal{E}}_x& 0& -{\mathcal{B}}_z+\mathbf{i}{\mathcal{E}}_z& {\mathcal{B}}_y-\mathbf{i}{\mathcal{E}}_y\\ & -{\mathcal{B}}_y+\mathbf{i}{\mathcal{E}}_y& {\mathcal{B}}_z-\mathbf{i}{\mathcal{E}}_z& 0& -{\mathcal{B}}_x+\mathbf{i}{\mathcal{E}}_x\\ & -{\mathcal{B}}_z+\mathbf{i}{\mathcal{E}}_z& -{\mathcal{B}}_y+\mathbf{i}{\mathcal{E}}_y& {\mathcal{B}}_x-\mathbf{i}{\mathcal{E}}_x& 0\end{array}\right).$$

(14)

This form is the gravitomagnetic tensor represented in the quaternion basis.

In the quaternion basis the electromagnetic tensor has the same form [42]. It is instructive to compare this electromagnetic tensor represented in the quaternion basis with a tensor written in the usual form in the theory of electromagnetism

$$F_{\mu \nu }={\partial }_{\mu }{A}_{\nu }-{\partial }_{\nu }{A}_{\mu }=\left(\begin{array}{ccccc}& 0& -E_x& -E_y& -E_z\\ & E_x& 0& -B_z& B_y\\ & E_y& B_z& 0& -B_x\\ & E_z& -B_y& B_x& 0\end{array}\right).$$

(15)

The complex gravitomagnetic field $\overrightarrow{\mathcal{B}}-\mathbf{i}\overrightarrow{\mathcal{E}}=(\overrightarrow{\Omega }-{\rho }_{q}c\overrightarrow{B})-\mathbf{i}(\overrightarrow{\Xi }+{\rho }_q\overrightarrow{E}).$ in Eq. (14), as follows from computations in Eq. (12), reads:

$$\left\{\begin{array}{ccc}& \overrightarrow{\mathcal{B}}=[\nabla \times \overrightarrow{p}],& \overrightarrow{\mathcal{E}}={\displaystyle \frac{\partial \overrightarrow{p}}{c\partial t}}+\nabla ϵ,\\ & \overrightarrow{B}=[\nabla \times \overrightarrow{A}],& \overrightarrow{E}={\displaystyle -\frac{\partial \overrightarrow{A}}{c\partial t}}-\nabla \Phi .\end{array}\right.$$

(16)

Here $\overrightarrow{E}$ and $\overrightarrow{B}$ are electric and magnetic fields. By analogy, $\overrightarrow{\Xi }$ and $\overrightarrow{\Omega }$ are gravitational and torsion fields. But more definitely, they determine irrotational and solenoidal fields. We will raise this topic further when we discuss Helmholtz’s theorems relating to such flows.

3. The Arbitrary Scalar Field Represents a Dark Energy Ocean

Let us study the arbitrary scalar field $\varphi ={\varphi }_0+{\varphi }_1$ in more detail. We rewrite Eq. (13):

$${\square \varphi }_1+{\square \varphi }_0+{\partial }_{\mu }{\epsilon }^{\mu }=0.$$

(17)

Here $\mu =0,1,2,3$ and ${\epsilon }^{\mu }=(ϵ,p_x.p_y,p_z)$. The scalar product of the 4-gradient operator ${\partial }_{\mu }$ by the 4-vector ${\epsilon }^{\mu }$ defines the 4-divergence ${\partial }_{\mu }{\epsilon }^{\mu }=(c^{-1}{\partial }_tϵ+{\partial }_x{p}_x+{\partial }_y{p}_y+{\partial }_z{p}_z)$. Thus, it is a 4D divergence of the energy-momentum of our world acting on the invisible world. Let us evaluate the scalar field $\varphi ={\varphi }_0+{\varphi }_1$ represented in Eq. (17). Note first that since $\epsilon$ has the dimension J·m⁻³ then the term ${\partial }_{\mu }{\epsilon }^{\mu }$ has the dimension J·m⁻⁴. From here it follows that $\varphi$ has the dimension J·m⁻² = kg·s⁻². The scalar field $\varphi$ as this estimation indirectly indicates, relates to the physical vacuum that behaves as a non-Newtonian (dilatant) fluid as shear stress increases [54,55]. This means that the fluid viscosity grows with the increasing material body velocity moving through the fluid (it is equivalent to its inert mass growing). The fluid "hardens" with increasing shear stress.

The scalar field ${\varphi }_0$ relates to the massless electromagnetic field. While ${\varphi }_1$ deals with the gravitomagnetic field loaded by mass. If so, then let us rewrite the first term and in Eq. (17) as follows:

$${\square \varphi }_1=\underset{\left(\mathrm{a}\right)}{\underbrace{\left(\square +{\left(\frac{mc}{\hslash }\right)}^2\right){\varphi }_1}}-\underset{\left(\mathrm{b}\right)}{\underbrace{{\left(\frac{mc}{\hslash }\right)}^2{\varphi }_1}}$$

(18)

We declare that the part covered by brace (a) vanishes:

$$\begin{array}{ccc}\hfill & & \stackrel{\left({\mathrm{a}}_0\right)}{\overbrace{\left({\hslash }^2\square +m^2{c}^2\right){\varphi }_1}}=\left({\hslash }^2{\gamma }_{\mu }{\gamma }^{\mu }{\partial }^{\mu }{\partial }_{\mu }+{\left(mc\right)}^2\right){\varphi }_1\hfill \\ & =& \stackrel{\left({\mathrm{a}}_1\right)}{\overbrace{\underset{\mathrm{Dirac}\mathrm{operator}+}{\underbrace{(\mathbf{i}\hslash {\gamma }_{\mu }{\partial }^{\mu }-mc)}}\underset{\mathrm{Dirac}\mathrm{operator}-}{\underbrace{(-\mathbf{i}\hslash {\gamma }^{\mu }{\partial }_{\mu }-mc)}}{\psi }_1}}=0.\hfill \end{array}$$

(19)

The first equation covered by brace (a₀) is the Klein-Gordon equation that obeys the "on-shell condition" [38]${(\hslash \overrightarrow{k})}^2+m^2{c}^2=p^2+m^2{c}^2=0$. Further, we rewrite this equation by extracting two Dirac operators multiplied by each other, as shown in the bracket (a₁). Here, the multiplication of these operators yields the sum $-i\hslash {\gamma }_{\mu }{\partial }^{\mu }mc+i\hslash {\gamma }^{\mu }{\partial }_{\mu }mc$. Note that ${\gamma }^{\mu }{\eta }_{\mu \nu }{\eta }^{\nu \mu }{\partial }_{\mu }={\gamma }_{\nu }{\partial }^{\nu }$, where ${\eta }_{\mu \nu }$ is the Minkowski metric with signature $(+,-,-,-)$ and ${\eta }_{\mu \nu }{\eta }^{\nu \mu }=1$. Thus, $-i\hslash {\gamma }_{\mu }{\partial }^{\mu }+i\hslash {\gamma }_{\nu }{\partial }^{\nu }$ vanishes.

This means that the multiplication of two Dirac operators under brace (a₁) describes the action of a coupled pair of particle and antiparticle, each with a half-integer spin. It is a boson with an integer spin. Its behavior is described by the Klein-Gordon equation. This equation states that (a) both fermions rotate about the common mass center; (b) the bonded fermions have zero charge; (c) the bonded fermions have integer spin equal to either 0 or ±1; (d) and they have equal masses, Figure 2. There are possible spontaneous transitions from bound fermions to massive scalar waves described by the following equation:

$$-\frac{1}{c^2}\frac{{\partial }^2}{\partial t^2}\Psi +{\nabla }^2\Psi =-\frac{m^2{c}^2}{{\hslash }^2}\Psi ,$$

(20)

and vice versa. The mass shell is a hyperboloid in the energy-momentum space, which is described by solving the following equation:

$$\mathcal{E}=\pm c\sqrt{m^2{c}^2+{\left|\overrightarrow{p}\right|}^2}.$$

(21)

This formula states that the energy and momentum of the particle and antiparticle are interconnected through the Einstein relation, and they lie on the mass surface.

Let’s return to Eq. (18). A hidden term covered by bracket (b) can be accepted as the dark matter manifestation. From the computations made above, these unknown particles are bosons composed of coupled pairs of particles and antiparticles. Following Occam’s Razor principle, we will not invent exotic particles but focus on baryon particles, like widespread atomic hydrogen [56,57]. Such a Bose-Einstein condensate of baryon-antibaryon matter exists at a temperature of T = 2.75 K. It is the cold dark matter [14,21,24,25]. Because of the Meissner effect, it is not visible by EM radiation.

Evaluations say that the critical density, ${\rho }_{\mathrm{crit}}$, can be achieved if 18 proton masses are contained within the volume of $\Delta V=3$ m³, namely, $18·(1.67\times {10}^{-27}\mathrm{kg})/\left(3{\mathrm{m}}^3\right)\approx {10}^{-26}\mathrm{kg}/{\mathrm{m}}^3$, compare with Eq. (2). For the visible matter, only one proton mass is in the volume, about $\Delta V=3$ m³. We find $1.67\times {10}^{-27}\mathrm{kg}/\left(3{\mathrm{m}}^3\right)\approx 5.57\times {10}^{-28}\mathrm{kg}/{\mathrm{m}}^3$, compare with Eq. (1).

Eighteen invisible proton masses versus one visible proton mass are seen to define a general picture of the gravitational influence of dark matter on visible baryon matter. The problem is to define the mechanism leading to the mechanism of vortex capture by invisible matter, visible as they rotate around a common center. Surprisingly, such a capture mechanism was experimentally studied by [58,59] in terrestrial laboratory conditions. Further, we will consider such a capture mechanism. But before that, we need to get the vorticity equation that describes a vortex exchanging with the quantum ether through viscosity fluctuations.

4. Equations of the Real Gravitomagnetic Field

Let us now return to the end of Sect. 2, to Eq, (14). Our aim here is to obtain the vorticity equation and to show how dark matter influences the rotation of vortex tails.

First, we note that the gravitomagnetic fields stem from the force density tensor, ${\mathbb{F}}_{\Xi \Omega }$, after applying the transposed generating differential operator ${\mathcal{D}}^{\mathbf{T}}$ [36,38]. For this reason, first we define a 4D current of the force density that takes into account gravitation influence, electromagnetic, and other accompanying forces acting on the medium within the volume $\Delta V$ under consideration. We define the density distribution term of all forces acting on the medium enclosed within this unit volume as follows:

$$\wp =\frac{1}{4\pi }(\nabla ·\overrightarrow{F}).$$

(22)

The divisor $4\pi$ is introduced to emphasize the commonality with Maxwell’s EM equations. Here $\overrightarrow{F}$ considers external forces acting on the medium inside the volume $\Delta V$ and its internal reaction on these forces with pressure gradients arising in this medium [42]. In addition, we must consider the shear viscosity due to friction of the medium layers induced by the pressure gradients. We could omit this effect from our consideration in the first approximation. Note, however, that because of this neglect, we can face unwanted singularities in subsequent calculations. The following step defines the 3D current density as $\overrightarrow{\Im }=\overrightarrow{\mathit{v}}\wp$. Here $\overrightarrow{\mathit{v}}$ is the velocity of the unit volume $\Delta V$ of the medium in question. Finally, we define the 4D current density as follows:

$$\mathbb{J}=\mathbf{i}c\wp {\eta }_0+{\Im }_x{\eta }_x+{\Im }_y{\eta }_y+{\Im }_z{\eta }_z.$$

(23)

The continuity equation for this current is

$$\frac{1}{4}\mathrm{trace}\mathcal{D}\mathbb{J}=\frac{\partial \wp }{\partial t}+(\nabla ·\overrightarrow{\Im })=0.$$

(24)

This is a manifestation of Newton’s third law.

By applying the transposed generating differential operator ${\mathcal{D}}^{\mathbf{T}}$ to the force density tensor (14) and equating it to the 4D current density $\mathbb{J}$ we come a set of gravitomagnetic equations:

$${\mathcal{D}}^{\mathbf{T}}·{\mathbb{F}}_{\mathcal{E}\mathcal{B}}={\mathcal{D}}^{\mathbf{T}}\mathcal{D}·T=\square T=\frac{4\pi }{c}\mathbb{J}.$$

(25)

One can see first that here we get the wave equation for the energy-momentum density tensor subjected to the external action of 4D density current $\mathbb{J}$. Amazing complementarity takes place between Eqs. (13) and (25) if we use the term $\mathbb{J}$ accepted in Eq. (23). Note that Eq. (25) describes wave oscillations of the visible matter, while Eq. (13) does the same with the invisible matter and energy.

Now revealing in detail the product ${\mathcal{D}}^{\mathbf{T}}·{\mathbb{F}}_{\mathcal{E}\mathcal{B}}$ we obtain the following expanded equation, where its terms have been sorted as coefficients under the quaternions ${\eta }_0$, ${\eta }_x$, ${\eta }_y$, ${\eta }_z$:

$$\begin{array}{ccc}\hfill {\mathcal{D}}^{\mathbf{T}}·{\mathbb{F}}_{\mathcal{E}\mathcal{B}}& =& \{-\underset{=0}{\underbrace{({\partial }_x{\mathcal{B}}_x+{\partial }_y{\mathcal{B}}_y+{\partial }_z{\mathcal{B}}_z)}}+\mathbf{i}\underset{=4\pi \wp }{\underbrace{({\partial }_x{\mathcal{E}}_x+{\partial }_y{\mathcal{E}}_y+{\partial }_z{\mathcal{E}}_z)}}\}{\eta }_0\hfill \\ \hfill & +& \left\{\left(\underset{4\pi {\Im }_x/c}{\underbrace{-\frac{1}{c}{\partial }_t{\mathcal{E}}_x-({\partial }_y{\mathcal{B}}_z-{\partial }_z{\mathcal{B}}_y)}}\right)+\mathbf{i}\left(\underset{=0}{\underbrace{-\frac{1}{c}{\partial }_t{\mathcal{B}}_x+({\partial }_y{\mathcal{E}}_z-{\partial }_z{\mathcal{E}}_y)}}\right)\right\}{\eta }_x\hfill \\ \hfill & +& \left\{\left(\underset{4\pi {\Im }_y/c}{\underbrace{-\frac{1}{c}{\partial }_t{\mathcal{E}}_y-({\partial }_z{\mathcal{B}}_x-{\partial }_x{\mathcal{B}}_z)}}\right)+\mathbf{i}\left(\underset{=0}{\underbrace{-\frac{1}{c}{\partial }_t{\mathcal{B}}_y+({\partial }_z{\mathcal{E}}_x-{\partial }_x{\mathcal{E}}_z)}}\right)\right\}{\eta }_y\hfill \\ \hfill & +& \left\{\left(\underset{4\pi {\Im }_z/c}{\underbrace{-\frac{1}{c}{\partial }_t{\mathcal{E}}_z-({\partial }_x{\mathcal{B}}_y-{\partial }_y{\mathcal{B}}_x)}}\right)+\mathbf{i}\left(\underset{=0}{\underbrace{-\frac{1}{c}{\partial }_t{\mathcal{B}}_z+({\partial }_x{\mathcal{E}}_y-{\partial }_y{\mathcal{E}}_x)}}\right)\right\}{\eta }_z\hfill \\ & =& \mathbf{i}4\pi \wp {\eta }_0+\frac{4\pi }{c}({\Im }_x{\eta }_x+{\Im }_y{\eta }_y+{\Im }_z{\eta }_z).\hfill \end{array}$$

(26)

On the left side of this equation, braces enveloping terms equal those being represented on the right side when their real and imaginary values permit it. Otherwise, they are zero. By gathering real and imaginary parts at the quaternions ${\eta }_0$, ${\eta }_x$, ${\eta }_y$, ${\eta }_z$, we obtain the gravitomagnetic equations

$$\begin{array}{ccc}& & (\nabla ·\overrightarrow{\mathcal{B}})=0,\hfill \end{array}$$

(27)

$$\begin{array}{ccc}& & (\nabla ·\overrightarrow{\mathcal{E}})=4\pi \wp ,\hfill \end{array}$$

(28)

$$\begin{array}{ccc}& & {\displaystyle [\nabla \times \overrightarrow{\mathcal{E}}]-\frac{1}{c}\frac{\partial }{\partial t}\overrightarrow{\mathcal{B}}}=0,\hfill \end{array}$$

(29)

$$\begin{array}{ccc}\hfill & & {\displaystyle [\nabla \times \overrightarrow{\mathcal{B}}]+\frac{1}{c}\frac{\partial }{\partial t}\overrightarrow{\mathcal{E}}=-\frac{4\pi }{c}\overrightarrow{\Im }.}\hfill \end{array}$$

(30)

As noted above, the gravitomagnetic fields $(\overrightarrow{\mathcal{E}},\overrightarrow{\mathcal{B}})$ are represented by the superposition of the gravitotorsion fields $(\overrightarrow{\Xi },\overrightarrow{\Omega })$ and the electromagnetic fields $(\overrightarrow{E},\overrightarrow{B})$. It can be noted here that these equations are ancestral equations leading to sets of wave equations, vorticity equations, and Hamilton-Jacobi equations describing the motion of material bodies in ideal and viscous media, as well as the quantum mechanical equations [23,34,36,37,42,60,61,62,63]. Here, we omit this detailed analysis. Instead, we will immediately proceed to the last part - the analysis of the Navier-Stokes equation, opening a path to the vorticity equations.

The problem here is to determine parameters that can lead to the appearance of a flat orbital velocity profile. Looking ahead, it can be noted that the critical parameter in this matter may be the viscosity of the medium. Since the medium is a superfluid with zero viscosity, the viscosity dispersion around zero can be such a critical parameter.

5. Navier-Stokes Equation: Forward to the Vorticity Equation

Let us look on Eq. (28). By rewriting the expression $(\nabla ·\overrightarrow{\mathcal{E}})=4\pi \wp$ in details we have

$$\nabla \left(\frac{\partial \overrightarrow{p}}{c\partial t}+\nabla ϵ-\overrightarrow{F}\right)=0\to \frac{\partial \overrightarrow{p}}{c\partial t}+\nabla ϵ=\overrightarrow{F}+{\overrightarrow{C}}_{*}.$$

(31)

The term ${\overrightarrow{C}}_{*}$ is an arbitrary solenoid force density. It disappears when the gradient operator is applied to this equation. Here, the force density $\overrightarrow{F}$ comes from Eq. (22). It is the sum of all force densities acting on the medium under consideration within a unit volume. They are (a) external forces; (b) internal forces represented by the pressure gradients arising in the medium under the action of the external forces; and (c) forces leading to the sliding of the layers of the medium on each other (they can lead to heat losses).

Considering the above, we get the equation for incompressible fluids, $\nabla ·\overrightarrow{\mathit{v}}=0$ [64], divided preliminarily by the density distribution $\rho$:

$$m\left(\frac{\partial \overrightarrow{\mathit{v}}}{\partial t}+(\overrightarrow{\mathit{v}}·\nabla \overrightarrow{\mathit{v}})\right)=-\nabla (U+Q)+\nu \left(t\right){\nabla }^2\overrightarrow{\mathit{v}}+{\overrightarrow{\Omega }}_{*}.$$

(32)

This is the modified Navier-Stokes equation. Here ${\overrightarrow{\Omega }}_{*}={\overrightarrow{C}}_{*}/\rho$ and we have (a) the external force is conservative and can be represented as a negative gradient of the external potential energy U, and (b) the pressure gradient can be rewritten by adding an extra term, i.e., the pressure multiplied by the gradient of logarithm from the density distribution. As a result, it leads to a gradient of the quantum potential Q. In this equation, we represent the external potential and the quantum potential as the sum under the gradient. And (c) the viscosity coefficient $\nu \left(t\right)$ is due to the layers sliding because of the pressure gradient arising within the medium.

It would seem that if the medium is superfluid, then the viscous term can be zeroed out in the first approximation. However, as practice has shown, sooner or later, when analyzing such an equation, you will have to face singularities. To get around these obstacles, a scientist usually shifts the poles of a Green function to the negative real half-plane. In fact, he manually enters the attenuation. We will not reject viscosity. However, we believe that the average viscosity vanishes over time, but its variance is not zero.

The term $\nu =\mu /{\rho }_M$ is the kinematic viscosity coefficient. Its dimension is [m²$·$s⁻¹] such as that of the diffusion coefficient. However, the medium, which fills all spaces densely, is superfluid [18,20,65,66]. Therefore, such a medium is considered to have a zero viscosity coefficient, $\nu =0$. However, as was said above we assume that the viscosity, adhesion, of the medium is zero on average over time, but its dispersion is not zero [36,42,60,67]:

$$\langle \nu \left(t\right)\rangle =0_{+},\langle \nu \left(t\right)\nu \left(0\right)\rangle >0.$$

(33)

These expressions describe an exchange of energy with the superfluid ether medium, the invisible world, as described in Sec. 3. The first relation states that, on average, there is no exchange with the invisible world in time. From the second expression, it follows that the exchange occurs regularly, both with the return of energy and with its acquisition randomly. It should be noted that since the viscosity coefficient depends on time but does not contain spatial variables, the invisible medium is non-local and consists of the Bose-Einstein condensates.

Before considering the vorticity equation, we apply the Helmholtz decomposition theorem of velocity into two components, irrotational velocity and solenoidal velocity:

$$\overrightarrow{\mathit{v}}={\overrightarrow{\mathit{v}}}_I+{\overrightarrow{\mathit{v}}}_{\mathrm{o}},$$

(34)

where subscripts I and O indicate the existence of irrotational and orbital (solenoidal, rotational, tangential) fluid velocities. We write the subscript O (orbital) instead of S (solenoidal) because the letter S relates to the notation of the scalar field and the action function S.

The component irrotational, ${\overrightarrow{\mathit{v}}}_I$, and orbital, ${\overrightarrow{\mathit{v}}}_{\mathrm{o}}$, velocities satisfy the following conditions:

$$\left\{\begin{array}{cc}(\nabla ·{\overrightarrow{\mathit{v}}}_I)\ne 0,& [\nabla \times {\overrightarrow{\mathit{v}}}_I]=0,\\ \\ (\nabla ·{\overrightarrow{\mathit{v}}}_{\mathrm{o}})=0,& [\nabla \times {\overrightarrow{\mathit{v}}}_{\mathrm{o}}]=\overrightarrow{\omega }.\end{array}\right.$$

(35)

Now, let us apply the curl operator to Eq. (32). We obtain [61]

$$\frac{\partial \overrightarrow{\omega }}{\partial t}+(\overrightarrow{\omega }·\nabla )\overrightarrow{\mathit{v}}=\underset{\left(\mathrm{a}\right)}{\underbrace{\nu \left(t\right){\nabla }^2\overrightarrow{\omega }}}+{\omega }_{*}.$$

(36)

Here ${\overrightarrow{\omega }}_{*}=[\nabla {\overrightarrow{\Omega }}_{*}]$ is a vorticity resulting from applying the curl operator to the arbitrary external solenoidal force ${\overrightarrow{\Omega }}_{*}$. For simplicity, we consider it equal to zero. The term covered by brace (a) describes the dissipation of energy stored in the vortex.

Figure 3. The vortex shown in the cylindrical coordinate system extends from minus infinity to plus infinity [36]. The vortex wall is a boundary where the orbital speed, ${\mathit{v}}_o$, reaches maximal values; a is the radius of the vortex tube, $\overrightarrow{\omega }$ is the vorticity. Profiles of the orbital velocity ${\overrightarrow{\mathit{v}}}_o$ at different N [36]: (a) $N=1$, normal profile of the Lamb-Oseen vortex orbital velocity ${\mathit{v}}_{\mathrm{o}}$ [68]); (b) $N\gg 1$, profile of the non-decreasing solenoidal velocity ${\mathit{v}}_{\mathrm{s}}$ with distance as for rotating spiral galaxies, for example [69].

Figure 3. The vortex shown in the cylindrical coordinate system extends from minus infinity to plus infinity [36]. The vortex wall is a boundary where the orbital speed, ${\mathit{v}}_o$, reaches maximal values; a is the radius of the vortex tube, $\overrightarrow{\omega }$ is the vorticity. Profiles of the orbital velocity ${\overrightarrow{\mathit{v}}}_o$ at different N [36]: (a) $N=1$, normal profile of the Lamb-Oseen vortex orbital velocity ${\mathit{v}}_{\mathrm{o}}$ [68]); (b) $N\gg 1$, profile of the non-decreasing solenoidal velocity ${\mathit{v}}_{\mathrm{s}}$ with distance as for rotating spiral galaxies, for example [69].

By omitting this term (i.e., assuming $\nu =0$), we proceed to the action of the Helmholtz theorem [70]: (i) if fluid particles form, in any moment of time, a vortex line, then the same particles support the vortex line both in the past and in the future; (ii) ensemble of the vortex lines traced through a closed contour forms a vortex tube; (iii) intensity of the vortex tube is constant along its length and does not change in time. The vortex tube drawn in Figure 3 either (a) goes to infinity by both endings; or (b) these endings lean on boundary walls containing the fluid; or (c) these endings are locked to each other, forming a vortex ring.

Now we will not remove the viscosity. Instead, we hypothesize that the viscosity fluctuates around zero. Its expectation is zero. These fluctuations represent the exchange of stored energy in the vortex with zero-point vacuum oscillations. Since $\langle \nu \left(t\right)\rangle =0_{+}$, the vortex does not completely disappear. It is a long-lived object.

Here, we consider a redial component of the vortex dynamics. Let us look at the vortex tube in its cross-section, which is oriented along the z-axis and its center is placed in the origin of the plane $(x,y)$, Figure 3. Eq. (36), written in this cross-section, appears as follows:

$$\frac{\partial \omega }{\partial t}=\nu \left(t\right)\left(\frac{{\partial }^2\omega }{\partial r^2}+\frac{1}{r}\frac{\partial \omega }{\partial r}\right).$$

(37)

A general solution to this equation has the following view of vorticity [36]:

$$\omega (r,t)=\frac{\Gamma }{4\Sigma (\nu ,t,\sigma )}exp\left\{{\displaystyle -\frac{r^2}{4\Sigma (\nu ,t,\sigma )}}\right\},$$

(38)

and orbital speed

$${\mathit{v}}_{\mathrm{o}}(r,t)=\frac{1}{r}\underset{0}{\overset{r}{\int }}\omega (r^{\prime },t)r^{\prime }d{r}^{\prime }=\frac{\Gamma }{2r}\left(1-exp\left\{{\displaystyle -\frac{r^2}{4\Sigma (\nu ,t,\sigma )}}\right\}\right).$$

(39)

Here $\Gamma$ is the integration constant having dimension [length²/time]. And the denominator $\Sigma (\nu ,t,\sigma )$ looks as

$$\Sigma (\nu ,t,\sigma )=\underset{0}{\overset{t}{\int }}\nu \left(\tau \right)d\tau +{\sigma }^2.$$

(40)

Here $\sigma$ is an integration constant. It is such that the denominator must be positive regardless of how the viscosity coefficient fluctuates, whether the fluctuations look like white noise, or whether they contain periodic components.

5.1. Lamb-Oseen Vortex, Pancake Vortex, and Rogue Wave

Note that for viscous hydrodynamic flows ($\nu >0$), the above formulas with $\sigma \ge 0$ simulate a decaying line vortex because of nonzero viscosity

$$\omega (r,t)=\frac{\Gamma }{4(\nu t+{\sigma }^2)}exp\left\{{\displaystyle -\frac{r^2}{4(\nu t+{\sigma }^2)}}\right\},$$

(41)

$${\mathit{v}}_{\mathrm{o}}(r,t)=\frac{\Gamma }{2r}\left(1-exp\left\{{\displaystyle -\frac{r^2}{4(\nu t+{\sigma }^2)}}\right\}\right).$$

(42)

At $\sigma =0$ this solution is that of the Lamb-Oseen vortex [68]. A normal profile of the orbital velocity ${\mathit{v}}_{\mathrm{o}}$ of the Lamb-Oseen vortex rotating about the central axis of the vortex tube is shown in Figure 4a. One can see that the velocity grows from the central axis up to the vortex wall. Next, it sharply decreases as one moves away from the origin. At $\sigma \ne 0$ the orbital velocity also grows from the central axis up to the vortex wall. But next, it begins monotonically to decrease the slower, the larger $\sigma$, as the distance terns to infinity, Figure 4b. This indicates the path of the search for the flattening of the orbital velocity profile in the region far from the origin.

When we reject $\nu$ in Eqs. (41) and (42) these solutions degenerate into the so-called long-lived Gaussian pancake vortices [71,72,73]. Note that such long-lived Gaussian pancake formations arise because of the inserted $\sigma$ parameter as an arbitrary integration constant, which was not originally provided in the equation (36). We had to introduce this parameter to prevent singularities that may occur due to fluctuations in the viscosity coefficient.

There are phenomena observed in ocean hydrodynamics that show pancake vortex formations. These are killer waves (wandering waves, monster waves, white waves, rogue waves, freak waves, crazy wave) – giant single waves [74,75] arising in the ocean, 20–30 meters high (and sometimes larger), having uncharacteristic behavior for sea waves. An important circumstance is the suddenness of their appearance and disappearance. We can hypothesize that the appearance and disappearance of killer waves are caused by fluctuations of clumps of dark matter, which induce giant wandering waves on the expanses of movable ocean water. Note that the mass of dark matter is five times greater than that of visible matter. Ancient civilizations may have learned about this phenomenon by erecting monumental structures such as pyramids.

5.2. Rotation of Spiral Galactic Arms

Let us return to Eqs. (38)-(40). Because the viscosity coefficient $\nu \left(t\right)$ fluctuates about zero according to our assumption, Eq. (33), the integral in Eq. (40) returns zero. For that reason the solution written in Eq. (39) takes an exceptionally simple form

$${\mathit{v}}_{\mathrm{o}}(r,t,\sigma )=\frac{\Gamma }{2r}\left(1-exp\left\{{\displaystyle -\frac{r^2}{4{\sigma }^2}}\right\}\right).$$

(43)

The coherent Gaussian vortices, being time-independent objects, allow superposition:

$${\mathit{v}}_s(r,t)=\sum _{n=1}^N{\mathit{v}}_{\mathrm{o}}(r,t,{\sigma }_n),$$

(44)

where ${\sigma }_n$ grows with increasing n. For simplicity, we consider the following simple dependence:

$${\sigma }_n=\sigma ·n=0.1n,n=1,2,\cdots ,N.$$

(45)

Figure 3 shows the orbital speed as a function of the distance r for two different cases: $N=1$, and $N\gg 1$, respectively. The profile of the speed grows monotonically as r increases from the center of rotation. In the case of $N=1$, the speed after reaching the vortex wall begins to drop monotonously, Figure 3, curve (a). For a sufficiently large N, the speed behind the vortex wall comes to a flat level, Figure 3, curve (b). This case agrees well with the astronomical observations of the orbital speeds of rotating spiral galaxies [76]. Figure 5 shows the orbital speeds for different N as functions of time.

The fact that time is not explicitly represented in the formulas (43)-(44) is because fluctuations of the viscosity coefficient, $\nu \left(t\right)$, are defined as the white noise (33). In principle, we can define the existence of some harmonic oscillations [77] against white noise, which simulate the breathing of a galaxy. Figure 6 shows such an imitation of breathing. It is visible as ripples of the flat plateau on the curve.

This breathing is caused by the exchange of energy with the physical vacuum at ultra-low frequencies owing to the term ${\Sigma }_n\left(t\right)$. At these ultra-low frequencies, this exchange occurs via virtual particle-antiparticle pairs, i.e., proton-antiproton pairs, as stated above. This manifestation stems from dark matter [69]. The ultra-low frequencies are perhaps due to the slow shifting of the dark matter clumps as the galaxy evolves further. Such a shifting of the dark matter one can guess is a consequence of the galaxy evolution processes.

As follows from numerous observations, the flat profile of orbital velocities during the rotation of spiral galaxies around their central core can be provided by additional invisible matter, called dark matter [69,78]. De Blok et al. [76] have presented a rich set of rotation curves, Figure 7, supporting the tendency of the orbital velocities of rotating galaxy arms to possess a flat profile. Thirty diagrams showing the orbital velocity ${\mathit{v}}_{\mathrm{o}}$ (km·s⁻¹) versus the distance r (in kiloparsec) demonstrate different profiles, where do not always follow up well with the monotonous outlet of curves on to a flat profile.

Sbitnev and Fedi [11] have approximated these curves using the formulas given in Eqs. (43)-(45). Only instead of ${\sigma }_n=\sigma ·n$ we have assigned ${\sigma }_n=\sigma$ being fitted for each $n=0,1,2,\cdots$ individually. The task was to draw the curve passing within confidence intervals so that it came to experimentally measured points as closely as possible. Fitted curves are drawn in red in Figure 7 against the black-white data [76] representing the orbital speed of galaxies measured by [72]. These curves fit the confidence intervals quite well. However, an important remark according to this fitting is that the fitting coefficients $\sigma$ for low $n=1,2,3,\cdots$ have a large scatter. This indicates that the blocks of dark matter are not ordered but rather chaotically organized. It appears that the observed galaxies are in a transitional stage of evolution. Most likely, there is a dynamic evolution occurring here - a change in the motion of the components of the galaxy (in this case, clumps of dark matter, perhaps, together with visible matter as well).

6. Conclusion

Quaternions written out as four matrices of rank 4x4 are natural variables describing the algebra of the 4D spacetime having 10 degrees of motion freedom. Within the framework of the quaternion basis superimposed on 4D space-time, the energy-momentum density tensor is recorded, which describes the current state of visible matter in this 4D space-time. Note that the density tensor consists of the superposition of two parts. The first part describes the mechanic-rotational state of the visible matter, loaded by brute mass. The second part deals with the electromagnetic state of this matter, loaded by distributed charges. By applying differential generating operators to this energy-momentum density tensor, we disclose a whole set of equations describing the behavior of our visible world.

One point that the author has tried to avoid until now concerns an unclear description of the energy-momentum density tensor. In principle, it can be determined up to an arbitrary scalar field, which subsequently does not affect the results of calculations in any way. Cutting off this scalar field is achieved by imposing a Lorentz gauge condition. After performing that operation, you can forget about the existence of such a scalar field.

Observe, however, that the Lorentz gauge represents the Klein-Gordon equation for this arbitrary scalar field, where the right part is represented by partial derivatives of the energy-momentum density tensor. While the left part contains two scalar fields subjected to the action of the D’Alembert wave operator. These two scalar fields relate to the mechanical-rotational and electromagnetic energy-momentum density tensors. The first scalar field is a mass-loaded one. We can add and subtract a certain mass-carrying term and thereby derive the Klein-Gordon equation loaded with mass in explicit form. This equation describes the coexistence of particles and antiparticles on a mass shell. The second scalar field relating to the electromagnetic field is mass-free. This gives rise to the following assumption: that the Lorentz gauge offers a description of the invisible world, represented by dark matter and dark energy. In fact, this invisible world is the superfluid quantum ether that many scientists write about.

The interaction of this ether with our visible world is possible due to the fluctuation of viscosity, which is zero on average over time. Note, first, that solutions of the vorticity equation contain both the viscosity coefficient multiplied by time and a vague integration parameter having the dimension [meter2]. Usually, this parameter is discarded as unnecessary. The remaining solution is known as the Lamb-Oseen vortex, which rapidly relaxes to zero as time and distance go on. However, the integration parameter turns out to be decisive in describing the so-called long-lived Gaussian pancake formations, which represent dangerous rogue waves in the ocean. If under terrestrial conditions the lifetime of such killer waves is limited to the coastlines of continents, then in the gaseous atmosphere of Jupiter their lifetime can be unlimited. An important point for us is not only the extension of the lifetime of such Gaussian pancake formations, but also the flattening of their profiles with distance.

We argue that the inclusion of this vague integration parameter in the vortex solution is necessary to prevent the appearance of singular jumps that could arise due to fluctuations in the viscosity of the medium. The emergence of this vague integration parameter is due to the dark matter chunks exchanging energy fluctuations with the visible matter through the viscosity fluctuations

Because of the fluctuation in viscosity, chunks of dark matter pull visible matter behind them. Can it be observed in terrestrial experimental conditions? Perhaps a similar pulling up of a non-rotated driven disk by the drive disk has been observed by Prof. Samokhvalov [79]. The author, having studied the works of Samokhvalov, proposed to investigate the effect of this unknown force, pulling an unstoppable, driven disk, with the involvement of neutron interferometry [36].

The meaning of neutron interferometry is to repeat the Aharonov-Bohm experiment [80] to register the effect of the electromagnetic vector potential "immured" in a cylinder shielding magnetic fields. The diagram of the neutron interference experiment is shown in Figure 8. It is expected that if the capture of a stationary slave disk is caused by the reorganization of clumps of dark matter, the measurement result will show the appearance of asymmetry in interference fringes.