The Halo Effect and Quantum Vortices. Not So Dark with Alena Tensor

1. Introduction

Research on the dark sector has been ongoing for many years [1] and there is still no theoretical consensus [2] or convincing experimental evidence regarding its nature [3]. We have not found dark matter signals for the WIMP models [4], axions/ALPs [5], SIDM [6], despite new experimental approaches [7], experiments in underground detectors [3], LZ experiments [4], XENON-nT [8] or SuperCDMS [9], and the latest observational data (e.g., "Hubble tension", ”Sigma-8 tension”) make the issue of the dark sector even more puzzling [10]. The dominant theoretical model is still $\mathsf{\Lambda }CDM$ with dark matter haloes [11], although alternative theories such as MOND/relativistic generalizations [12,13], dark photons [14], TeVeS/related constructions [12], f(R) [15], black holes [16], or the recently popular ”emergent/entropic gravity” [17] have achieved some success at selected scales. However, all these approaches experience ups and downs depending on the subsequent observational data [4,8,10,18]. There is also considered the possibility of hybrid models (e.g. a superfluid DM combining features of MOND and DM) [19], or even incorrect/weakened estimates of baryon masses and systematic errors in the M/L estimate [20], and the interaction of dark matter and energy [21] but so far this entire massive effort by the scientific community has not yielded a definitive conclusion. We have more certainty about dark energy because the universe is definitely expanding, and there are ample, multiple, independent confirmations, eg. from SN Ia, BAO, and CMB [11,22]. It is still uncertain whether the cosmological constant is indeed constant, although recent DESI (BAO) analyses indicate the possibility of w(z) dynamics [23], and the Euclid mission is rapidly delivering high-quality weak lensing maps and surveys that will be crucial for measuring the dark energy equation of state [24].

The Alena Tensor is a relatively young field of research and has not had much relevance to the dark sector until now. Previous work has focused on developing a dual description for physical systems with matter and fields in which the metric tensor is not a feature of spacetime but only a method of describing it. Its aim was to provide a smooth transition between curvilinear description consistent with GR [25] and a flat (classical and quantum [26]) description for simple cases with dust [27], which was analyzed mainly for physical systems with electromagnetic fields [28]. This paper will demonstrate that extending the Alena Tensor to the general case yields results, that provide a representation of the spherical dark matter halo and, naturally, provide an explanation for the cosmological constant. The obtained picture of GR equations leads to a potential explanation of the dark matter phenomenon based on the energy associated with the rotation itself and for about 100 diverse galaxies from the SPARC database [29] gives a fairly good approximation of the observed rotation curves.

In the first part of the paper, an introduction to the Alena Tensor model for dust will be presented, then this approach will be extended to the general form of the matter energy-momentum tensor, leading to the "dark matter halo effect". The quantum description resulting from the proposed approach will also be analyzed and tested to generate elementary particles. The results will be analyzed and discussed and shown to lead to the conclusions described in the abstract.

2. The Extended Alena Tensor approach

In this chapter, the conclusions reached so far regarding the Alena Tensor will be recalled, the notation will be introduced, and the reasoning from the previous articles will be generalized to all gauge fields and all forces acting in the physical system. The obtained results will be then used to formulate conclusions regarding the dark matter and compared with observational data. The author uses the Einstein summation convention, metric signature $(+,-,-,-)$ and commonly used notations.

2.1. Transforming Curved Path into Geodesic for Dust

As a first step, one may generalize the solution proposed in [27] in such a way that the forces resulting from all gauge fields are related to the metric tensor of curved spacetime.

One may begin the reasoning by introducing tensor $F^{\mu \nu \alpha \beta }$ defined in terms of the gauge field tensors ${\mathbb{F}}_A^{\mu \nu }$ for each gauge group A, and the stress-energy tensor ${{\rm Y}}^{\mu \nu }$ for such generalized field

$$F^{\mu \nu \alpha \beta }\equiv \sum _A{\mathbb{F}}_A^{\mu \nu }\otimes {\mathbb{F}}_A^{\alpha \beta };{{\rm Y}}^{\mu \nu }\equiv F^{\mu \alpha \nu \beta }{g}_{\alpha \beta }-g^{\mu \nu }{\displaystyle \frac{1}{4}}{F}^{\alpha \gamma \beta \delta }{g}_{\alpha \beta }{g}_{\gamma \delta }$$

(1)

and denote the invariant of this field as $p_{\mathsf{\Lambda }}$. Following reasoning from [27] let this field invariant be defined dually as follows

$$p_{\mathsf{\Lambda }}\equiv {\displaystyle \frac{1}{4}}{F}^{\alpha \gamma \beta \delta }{g}_{\alpha \beta }{g}_{\gamma \delta }\equiv p_o{\mathbb{k}}^2;\mathbb{k}\equiv {\mathbb{k}}_{\mu \nu }{g}^{\mu \nu }$$

(2)

where $p_o$ is certain constant (or simply invariant, independent of the metric) and where ${\mathbb{k}}^{\mu \nu }$, as raised in [26], in this approach is a metric tensor describing a curved spacetime in which all motion occurs along geodesics. By making variation on $-p_{\mathsf{\Lambda }}$ with respect to metric $g_{\mu \nu }$ (Hilbert’s method) one obtains the energy-momentum tensor of the field from (1) expressed dually as

$${{\rm Y}}^{\mu \nu }=p_{\mathsf{\Lambda }}\left({\displaystyle \frac{4}{\mathbb{k}}}{\mathbb{k}}^{\mu \nu }-g^{\mu \nu }\right)$$

(3)

Such approach, exactly as shown in [27], establishes a relationship between the field and the metric tensor ${\mathbb{k}}^{\mu \nu }$ and in the spacetime considered as described by the metric tensor $g^{\mu \nu }\to {\mathbb{k}}^{\mu \nu }$, one obtains $\mathbb{k}=4$ what yields that in curvilinear description of the system, energy-momentum tensor of the field ${{\rm Y}}^{\mu \nu }$ vanishes, maintaining continuity of function. As shown in previous articles, this causes the presence of a field in curved spacetime to manifest itself solely through curvature, which replaces the four-force densities in flat spacetime.

In flat spacetime, one may assume that the equations of motion for the gauge fields A are satisfied, thus one obtains gauge four-currents $J_A^{\nu }\equiv D_{\mu }{\mathbb{F}}_A^{\mu \nu }$. Therefore, the total density of the Yang-Mills four-forces [30,31]$f_{YM}^{\nu }$ is

$$f_{YM}^{\nu }={\partial }_{\mu }{{\rm Y}}^{\mu \nu }=\sum _A{J}_A^{\alpha }{\mathbb{F}}_{A\alpha }^{\nu }$$

(4)

where the self-interactions of gauge fields in a non-Abelian theory are reducing. Then, following the reasoning presented in [27] one may define coefficient ${\chi }_m\equiv {\displaystyle \frac{\varrho c^2}{p_{\mathsf{\Lambda }}}}$ where $\varrho$ represents matter density and is associated with the translational current ${\chi }_m{U}^{\mu }$. In the Alena Tensor approach, the existence of matter is thus a manifestation of the existence of fields what ensures, that without fields matter does not exist.

This allows to define the Lagrangian ${\mathcal{L}}_{dust}$ and obtain from it the strest-energy tensor (Alena Tensor) for the system with dust $T_{dust}^{\mu \nu }$ by variation on the metric

$${\mathcal{L}}_{dust}\equiv p_{\mathsf{\Lambda }}(1-{\chi }_m)\to T_{dust}^{\mu \nu }\equiv \varrho U^{\mu }{U}^{\nu }-(1-{\chi }_m){{\rm Y}}^{\mu \nu }$$

(5)

In accordance with [28], assuming ${\varrho }_o$ as rest mass density, four-momentum density is defined as $\varrho U^{\alpha }\equiv {\varrho }_o\gamma U^{\alpha }$ what takes into account motion and Lorentz contraction of the volume. Total translational four-force density acting on matter is therefore defined as

$$f^{\nu }\equiv {\partial }_{\mu }\varrho U^{\mu }{U}^{\nu }=\varrho U^{\mu }{\partial }_{\mu }{U}^{\nu }=\varrho {\displaystyle \frac{d{U}^{\nu }}{d\tau }};{\partial }_{\mu }\varrho U^{\mu }=0$$

(6)

As shown in [26], the above amendment introduces a natural property concerning curved spacetime, assuming that for dust, geodetic motion is expected

$$\begin{array}{c}\hfill {\partial }_{\alpha }\varrho U^{\alpha }=0\to U_{,\alpha }^{\alpha }=-{\displaystyle \frac{d\gamma }{dt}}\to U_{;\alpha }^{\alpha }=0\end{array}$$

(7)

$$\begin{array}{c}\hfill U^{\alpha }{U}_{;\alpha }^{\beta }=0\to {\displaystyle \frac{D{U}^{\beta }}{D\tau }}=0;{\left(\varrho U^{\alpha }{U}^{\beta }\right)}_{;\alpha }=0\end{array}$$

(8)

Now it can be noticed, that in flat spacetime the four-divergence of the above tensor $T_{dust}^{\mu \nu }$ can be interpreted as the density of the four-forces acting on matter $f^{\nu }$ reduced by the density of the field-related four-forces

$$f_{field}^{\nu }\equiv {\partial }_{\mu }\left[(1-{\chi }_m){{\rm Y}}^{\mu \nu }\right]=(1-{\chi }_m)f_{YM}^{\nu }+f_{gr}^{\nu };f_{gr}^{\nu }\equiv {{\rm Y}}^{\mu \nu }{\partial }_{\mu }(1-{\chi }_m)=-{{\rm Y}}^{\mu \nu }{\partial }_{\mu }{\chi }_m$$

(9)

As shown in [26], $f_{gr}^{\nu }$ can be associated with the existence of gravity in the system, while $f_{rr}^{\nu }\equiv {\chi }_m{f}_{YM}$ behaves as a radiation-reaction force, reducing the value of forces due to the field and upholding the conservation of energy, ensuring that the increasing energy density associated with matter $\varrho c^2$ does not exceed the total energy density $p_{\mathsf{\Lambda }}$ available in the system

$${\partial }_{\mu }{T}_{dust}^{\mu \nu }=f^{\nu }-f_{field}^{\nu }=f^{\nu }-f_{YM}^{\nu }-f_{gr}^{\nu }+f_{rr}^{\nu }$$

(10)

where $\underset{{\chi }_m\to 1}{lim}{f}_{field}^{\nu }=0$. Presented approach also indicates the anstaz for the Kerr-Schild type metrics for curved spacetime

$${\mathbb{k}}^{\mu \nu }\equiv \sum _i{c}_i{l}_i^{\mu }{l}_i^{\nu }+{\displaystyle \frac{\mathbb{k}}{4}}{\eta }^{\mu \nu };0={\eta }_{\mu \nu }{l}_i^{\mu }{l}_i^{\nu }$$

(11)

where $l_i^{\mu }$ are null vectors and $c_i$ are related coefficients.

In the next section, the above model will be expanded to include rotation-related components, which will prove crucial for describing dark sector phenomena and allow to obtain a description that agrees with the observational results.

2.2. Rotational Energy

It can be noticed that the radiation reaction force should take into account the total energy associated with the body, so in addition to the energy associated with the translational motion, it seems necessary to take into account the rotational energy.

One may thus introduce a projector ${\Delta }^{\mu \nu }$, flow vorticity tensor ${\omega }^{\mu \nu }$, positive coefficient ${\chi }_{\omega }$ equal to the rotational energy ${\mathcal{E}}_{rot}$ up to $p_{\mathsf{\Lambda }}$ and some metric independent auxiliary $\alpha$ with the dimension of the square of time

$$\begin{array}{c}\hfill {\Delta }^{\mu }{}_{\nu }\equiv g^{\mu }{}_{\nu }-{\displaystyle \frac{1}{c^2}}{U}^{\mu }{U}_{\nu };{\omega }_{\mu \nu }\equiv {\Delta }_{\mu }{}^{\alpha }\Delta _{\nu }{}^{\beta }\nabla _{[\alpha }{U}_{\beta ]};{\chi }_{\omega }\equiv {\displaystyle \frac{\alpha }{2}}{\omega }^{\mu \nu }{\omega }_{\mu \nu };{\mathcal{E}}_{rot}\equiv p_{\mathsf{\Lambda }}{\chi }_{\omega }\end{array}$$

(12)

Defining Lagrangian density ${\mathcal{L}}_T$ for the whole system and tensor $S^{\lambda \mu \nu }$ to describe emerging boundary terms

$$\begin{array}{c}\hfill {\mathcal{L}}_T=p_{\mathsf{\Lambda }}\left(1-{\chi }_{\omega }-{\chi }_m\right);S^{\lambda \mu \nu }\equiv p_{\mathsf{\Lambda }}\alpha \left(U^{\lambda }{\omega }^{\mu \nu }-U^{(\mu }{\omega }^{\nu )\lambda }\right)\end{array}$$

(13)

and introducing ${\mathsf{\Xi }}^{\mu \nu }$ tensor related to rotational properties of the system

$$\begin{array}{c}\hfill -{\mathsf{\Xi }}^{\mu \nu }\equiv \alpha p_{\mathsf{\Lambda }}\left({\omega }^{\mu \gamma }{\omega }^{\nu }{}_{\gamma }-{\displaystyle \frac{1}{2}}{\Delta }^{\mu \nu }{\omega }^{\alpha \beta }{\omega }_{\alpha \beta }\right)+{\nabla }_{\lambda }{S}^{\lambda \left(\mu \nu \right)}\end{array}$$

(14)

one obtains Alena Tensor $T^{\mu \nu }$ for the system derived with help of variational method on ${\mathcal{L}}_T$ in the form

$$\begin{array}{c}\hfill T^{\mu \nu }=\varrho U^{\mu }{U}^{\nu }-{\mathsf{\Xi }}^{\mu \nu }-\left(1-{\chi }_{\omega }-{\chi }_m\right){{\rm Y}}^{\mu \nu };T^{\mu \nu }{g}_{\mu \nu }=\varrho c^2-{\mathcal{E}}_{rot}\end{array}$$

(15)

Considering description in curved spacetime, described by the metric tensor ${\mathbb{k}}^{\mu \nu }$, the field tensor ${{\rm Y}}^{\mu \nu }$ vanishes, the system tensor reduces to the form $T_{matt}^{\mu \nu }\equiv \varrho U^{\mu }{U}^{\nu }-{\mathsf{\Xi }}^{\mu \nu }$ and its vanishing four-divergence means that any deviations from the geodesic motion with $a^{\mu }\equiv U^{\nu }{\nabla }_{\nu }{U}^{\mu }$ are compensated by rotation related forces. Using the standard kinematic decomposition one may calculate

$$\begin{array}{c}\hfill \epsilon \equiv {\displaystyle \frac{1}{c^2}}{T}_{matt}^{\mu \nu }{U}_{\mu }{U}_{\nu }=\varrho c^2+2{\mathcal{E}}_{rot};T_{matt}^{\mu \nu }{U}_{\nu }=\left(\varrho c^2+2{\mathcal{E}}_{rot}\right)U^{\mu }-\alpha p_{\mathsf{\Lambda }}\left(a_{\nu }{\omega }^{\nu \mu }+{\displaystyle \frac{c^2}{2}}{\Delta }^{\mu }{}_{\nu }\nabla _{\lambda }{\omega }^{\lambda \nu }\right)\end{array}$$

(16)

where the element in brackets in last equation represents in fact the purely spatial vorticity divergence. Assuming the classical definition of the energy flux $q^{\alpha }$ one also gets

$$\begin{array}{c}\hfill q^{\alpha }\equiv {\displaystyle \frac{1}{c^2}}{\Delta }^{\alpha }{}_{\mu }T_{matt}^{\mu \nu }{U}_{\nu }=-\alpha p_{\mathsf{\Lambda }}\left(a_{\nu }{\omega }^{\nu \mu }+{\displaystyle \frac{c^2}{2}}{\Delta }^{\mu }{}_{\nu }\nabla _{\lambda }{\omega }^{\lambda \nu }\right)\end{array}$$

(17)

Introducing classical shear tensor ${\sigma }^{\mu \nu }$ and effective vortex stress tensor ${\tau }^{\mu \nu }$ as

$$\begin{array}{c}\hfill {\sigma }^{\mu \nu }\equiv {\Delta }^{\mu \alpha }{\Delta }^{\nu \beta }{\nabla }_{(\alpha }{U}_{\beta )}-{\displaystyle \frac{1}{3}}{\nabla }_{\alpha }{U}^{\alpha }{\Delta }^{\mu \nu };{\tau }^{\mu \nu }\equiv {\displaystyle \frac{p_{\mathsf{\Lambda }}\alpha }{2}}\left({\sigma }^{\mu }{}_{\lambda }\omega ^{\nu \lambda }+{\sigma }^{\nu }{}_{\lambda }\omega ^{\mu \lambda }\right)\end{array}$$

(18)

one may thus rewrite $T_{matt}^{\mu \nu }$ in curved spacetime as

$$\begin{array}{c}\hfill T_{matt}^{\mu \nu }={\displaystyle \frac{\epsilon }{c^2}}{U}^{\mu }{U}^{\nu }+{\displaystyle \frac{1}{c^2}}\left(U^{\mu }{q}^{\nu }+U^{\nu }{q}^{\mu }\right)-{\mathcal{E}}_{rot}{\Delta }^{\mu \nu }-{\tau }^{\mu \nu }\end{array}$$

(19)

One may notice, that the system has a built-in anisotropic stress described by ${\tau }^{\mu \nu }$, but its source is not viscosity, but the coupling between shear and vorticity (between flow deformation and local spin angular momentum).

Considered in flat spacetime such approach introduces additional four-force density $f_{\mathsf{\Xi }}^{\nu }$ acting on matter and also changes $f_{gr}^{\nu }$ and the radiation reaction $f_{rr}^{\nu }$ to the form

$$\begin{array}{c}\hfill f_{\mathsf{\Xi }}^{\nu }\equiv {\partial }_{\mu }{\mathsf{\Xi }}^{\mu \nu };f_{rr}^{\nu }=\left({\chi }_m+{\chi }_{\omega }\right)f_{YM}^{\nu };f_{gr}^{\nu }={{\rm Y}}^{\mu \nu }{\partial }_{\mu }\left(1-{\chi }_{\omega }-{\chi }_m\right)\end{array}$$

(20)

One may now consider the impact of the above expansion of the Alena Tensor on quantum equations and GR equations.

2.2.1. Noether Tensor and Quantum Interpretation

Introducing for simplicity $\xi \equiv 1-{\chi }_{\omega }-{\chi }_m$ one may in flat spacetime define Belinfante superpotential $Q^{\lambda \mu \nu }$ referring to (13) and improvement $I^{\mu \nu }$ resulting from the variable $\xi$ as

$$\begin{array}{c}\hfill Q^{\lambda \mu \nu }\equiv {\displaystyle \frac{1}{2}}\left(S^{\mu \nu \lambda }+S^{\nu \mu \lambda }-S^{\lambda \mu \nu }\right);I^{\mu \nu }\equiv {\partial }_{\lambda }{Q}^{\lambda \mu \nu }+{\eta }^{\mu \nu }\left({\partial }_{\lambda }ln\xi \right){\mathbb{C}}^{\lambda }-({\partial }^{(\mu }ln\xi ){\mathbb{C}}^{\nu )}\end{array}$$

(21)

where based on [32] approach, $I^{\mu \nu }$ takes into account virial currents, where

$$\begin{array}{c}\hfill {\mathbb{C}}^{\mu }\equiv {\mathsf{\Pi }}^{\mu \alpha }{U}_{\alpha }-\sum _A{\mathbb{F}}_A^{\mu \lambda }{\mathbb{A}}_{\lambda }{}_A;{\mathsf{\Pi }}^{\mu \alpha }\equiv {\displaystyle \frac{\partial {\mathcal{E}}_{rot}}{\partial \left({\partial }_{\mu }{U}_{\alpha }\right)}}=\alpha p_{\mathsf{\Lambda }}{\Delta }^{\alpha }{}_{\rho }\Delta _{\beta }{}^{\mu }\omega ^{\rho \beta }\end{array}$$

(22)

As will be shown, the combination ${\mathbb{C}}_{\mu }/\xi$ defines the total Noether current associated with the underlying symmetry. It remains gauge invariant and includes both matter and field contributions, ensuring a consistent formulation of the continuity equation.

It may be noticed, that ${\partial }_{\mu }{I}^{\mu \nu }=0$ becomes zero on-shell, using the equations of motion and smooth boundary conditions. Since $Q^{\lambda \mu \nu }$ is antisymmetric in ${\partial }_{\mu }{\partial }_{\lambda }$ and partial derivatives commute in flat spacetime, therefore ${\partial }_{\mu }{\partial }_{\lambda }{Q}^{\lambda \mu \nu }=0$. One may thus simplify ${\partial }_{\mu }{I}^{\mu \nu }$ to the form of

$$\begin{array}{c}\hfill {\partial }_{\mu }{I}^{\mu \nu }=({\partial }^{\nu }ln\xi ){\partial }_{\mu }\left({\displaystyle \frac{{\mathbb{C}}^{\mu }}{\xi }}\right)+({\partial }_{\mu }ln\xi )\left({\partial }^{\mu }\left({\displaystyle \frac{{\mathbb{C}}^{\nu }}{\xi }}\right)-{\partial }^{\nu }\left({\displaystyle \frac{{\mathbb{C}}^{\mu }}{\xi }}\right)\right)\end{array}$$

(23)

The equations of motion for $U_{\mu }$ and gauge fields ${\mathbb{A}}_{\mu A}$ derived from the Lagrangian ${\mathcal{L}}_T=\xi p_{\mathsf{\Lambda }}$ imply the following two relations

$$\begin{array}{c}\hfill {\partial }_{\mu }\left({\displaystyle \frac{{\mathbb{C}}^{\mu }}{\xi }}\right)=0;{\partial }^{\nu }{\mathbb{C}}^{\lambda }-{\partial }^{\lambda }{\mathbb{C}}^{\nu }=\left({\partial }^{\nu }ln\xi \right){\mathbb{C}}^{\lambda }-\left({\partial }^{\lambda }ln\xi \right){\mathbb{C}}^{\nu }\\ \hfill \to {\partial }^{\nu }\left({\displaystyle \frac{{\mathbb{C}}^{\lambda }}{\xi }}\right)-{\partial }^{\lambda }\left({\displaystyle \frac{{\mathbb{C}}^{\nu }}{\xi }}\right)=0\to {\displaystyle \frac{{\mathbb{C}}^{\mu }}{\xi }}\equiv {\partial }^{\mu }{\displaystyle \frac{\mathsf{\Phi }}{\kappa }}\end{array}$$

(24)

The first relation is the Euler equation for the component $U_{\mu }$ (with ${\mathsf{\Pi }}^{\mu \alpha }$ as the momentum conjugate to ${\partial }_{\mu }{U}_{\alpha }$) after taking into account the fact that ${\mathcal{L}}_T$ has variable $\xi$. The second relation follows from the gauge field equations ${\partial }_{\mu }\left(\xi {\mathbb{F}}_A^{\mu \nu }\right)={\mathbb{J}}_A^{\nu }$ and from the definition of the "virial" part ${\sum }_A{\mathbb{F}}_A^{\mu \lambda }{\mathbb{A}}_{\lambda A}$ in ${\mathbb{C}}^{\mu }$, which ensures the proper gauge covariance.

The condition ${\partial }^{[\nu }({\mathbb{C}}^{\lambda ]}/\xi )=0$ guarantees the local existence of a scalar potential $\mathsf{\Phi }$ such that ${\displaystyle \frac{{\mathbb{C}}^{\mu }}{\xi }}={\partial }^{\mu }{\displaystyle \frac{\mathsf{\Phi }}{\kappa }}$, where $\kappa$ was introduced to maintain dimensional consistency. In spaces with non-trivial topology this potential can be multi-valued, which corresponds to quantized vortices and circulations. In this representation ${\mathbb{C}}^{\mu }/\xi$ plays the role of an ordinary Noether current of $\mathsf{\Phi }$, since the above relations follow from the Euler-Lagrange equations for $U_{\mu }$ and ${\mathbb{A}}_{\mu A}$ derived from ${\mathcal{L}}_T$.

One may now introduce Noether Lagrangian density ${\mathcal{L}}_N$ as

$${\mathcal{L}}_N[U,\mathbb{A},\mathsf{\Phi }]\equiv {\mathcal{L}}_T[U,\mathbb{A}]+{\displaystyle \frac{{\mathbb{C}}^{\mu }}{\xi }}[U,\mathbb{A}]{\partial }_{\mu }\mathsf{\Phi }\to {\mathsf{\Pi }}_{\mathsf{\Phi }}^{\mu }={\displaystyle \frac{\partial {\mathcal{L}}_N}{\partial \left({\partial }_{\mu }\mathsf{\Phi }\right)}}={\displaystyle \frac{{\mathbb{C}}^{\mu }}{\xi }}$$

(25)

The field $\mathsf{\Phi }$ is a local Noether phase (not a gauge phase) associated with the conserved total current ${\mathbb{C}}_{\mu }/\xi$. It introduces no new degrees of freedom but parametrizes the underlying $U\left(1\right)$ symmetry, so that the continuity equation can be written in terms of the phase gradient ${\partial }_{\mu }\mathsf{\Phi }$. From the above one obtains the Noether tensor in canonical form as

$$T_N^{\mu \nu }={\mathsf{\Pi }}^{\mu }{}_{\alpha \left(U\right)}\partial ^{\nu }{U}^{\alpha }+\sum _A{\mathsf{\Pi }}^{\mu A}{}_{\alpha \left(A\right)}\partial ^{\nu }{\mathbb{A}}^{\alpha A}+{\displaystyle \frac{{\mathbb{C}}^{\mu }}{\xi }}{\partial }^{\nu }\mathsf{\Phi }-{\eta }^{\mu \nu }{\mathcal{L}}_N$$

(26)

By the Rosenfeld-Belinfante identity, which equates the Hilbert energy-momentum tensor obtained by metric variation with the Belinfante-improved Noether tensor (on-shell), one obtains Alena Tensor as

$$\begin{array}{c}\hfill T^{\mu \nu }=T_N^{\mu \nu }+I^{\mu \nu }=\varrho U^{\mu }{U}^{\nu }-{\mathsf{\Xi }}^{\mu \nu }-\xi {{\rm Y}}^{\mu \nu }\end{array}$$

(27)

This construction yields a symmetric and conserved energy-momentum tensor, suitable for quantization. In the quantum description $T^{\mu \nu }$ is promoted to an operator whose expectation values reproduce the classical results. One may thus denote $D_{\mu }={\partial }_{\mu }-i{\sum }_A{g}_A{\mathbb{A}}_{\mu A}{T}_A$ with $g_A$ denoting the gauge couplings, ${\mathbb{A}}_{\mu A}$ the corresponding gauge fields, and $T_A$ the generators of the gauge group in the representation of the matter field. In flat spacetime, one therefore obtains the Alena Tensor $T^{\mu \nu }$ as

$$\begin{array}{c}\hfill T^{\mu \nu }={\displaystyle \frac{i}{4}}\overline{\psi }{\gamma }^{(\mu }{\overleftrightarrow{D}}^{\nu )}\psi -{\mathsf{\Xi }}^{\mu \nu }-\xi {{\rm Y}}^{\mu \nu }\end{array}$$

(28)

where $\overleftrightarrow{D}$ denotes the two-way derivative. It may be calculated, that the inclusion of the rotational energy term allows to construct below quantum-effective Lagrangian density, giving the expected structure of the energy-momentum tensor and providing a simple quantum interpretation

$$\begin{array}{ccc}\hfill {\mathcal{L}}_{\mathrm{eff}}& =& \underset{{\mathcal{L}}_{\mathrm{scalar}\left(\mathrm{vortex}\right)}}{\underbrace{{\displaystyle \frac{\hslash c}{2}}\left({\partial }_{\mu }\rho \right)\left({\partial }^{\mu }\rho \right)+{\displaystyle \frac{\hslash c}{2}}{\rho }^2\left({\partial }_{\mu }\varphi \right)\left({\partial }^{\mu }\varphi \right)}}+\underset{{\mathcal{L}}_{\psi }}{\underbrace{\overline{\psi }\left(i\hslash c{\gamma }^{\mu }{D}_{\mu }-m{c}^2\right)\psi }}\hfill \\ & +& \underset{{\mathcal{L}}_{\omega }+{\mathcal{L}}_{\mathrm{spin}-\mathrm{vortex}}}{\underbrace{{\displaystyle \frac{e^{-A\left(\rho \right)}}{2c^2\kappa }}{\omega }_{\mu \nu }{\omega }^{\mu \nu }-g\hslash \overline{\psi }{\Sigma }_{\mu \nu }\psi {\omega }^{\mu \nu }}}\hfill \end{array}$$

(29)

where

• $\varphi \left(x\right)$ - vortex phase field (action phase). Its gradient ${\partial }_{\mu }\varphi$ represents the generalized four-momentum flow associated with the vortex structure.
• $\rho \left(x\right)$ - amplitude of the complex condensate $\phi =\hslash c\rho e^{i\varphi }$. It determines the vortex core profile and sets the symmetry-breaking scale.
• ${\omega }_{\mu \nu }$ - vorticity tensor of the underlying medium. In this Lagrangian it is treated as an independent antisymmetric field capturing local rotational structure.
• ${\Sigma }_{\mu \nu }$ - spin generator in the fermionic representation ${\Sigma }_{\mu \nu }=\frac{i}{4}[{\gamma }_{\mu },{\gamma }_{\nu }]$.
• $e^{A\left(\rho \right)}$ - plays the role of the a dimensionless state-dependent stiffness function, encoding the effective elastic response of the vortex condensate, where it is assumed for calculation simplicity $p_{\mathsf{\Lambda }}\alpha =e^{-A\left(\rho \right)}/\left(c^2\kappa \right)$
• g - dimensionless spin-vorticity coupling constant, determining the strength of the interaction between fermionic spin and the vortex background.

The Madelung decomposition $\phi =\hslash c\rho e^{i\varphi }$ introduces an effective phase field whose gradient reproduces the conserved current $\hslash c{\rho }^2{\partial }^{\mu }\varphi$ so that the hydrodynamic Noether current defined in the classical theory ${\partial }^{\mu }{\displaystyle \frac{\mathsf{\Phi }}{\kappa }}$ is represented in the effective description by the amplitude-phase variables through the identification ${\partial }^{\mu }\mathsf{\Phi }=\kappa \hslash c{\rho }^2{\partial }^{\mu }\varphi$. In this sense, $\varphi$ does not introduce an additional dynamical degree of freedom, but parametrizes the same conserved current encoded in the underlying field texture, ensuring quantized circulation and capturing the spin-vorticity structure of the medium.

2.2.2. General Relativity Interpretation

One may now repeat the reasoning from [27] and define the generalized Ricci and Einstein tensors, where the ˜ sign indicates normalization with the constant $\kappa /2={\displaystyle \frac{4\pi G}{c^4}}$, and additional tensor ${\Theta }^{\mu \nu }$ as

$$\begin{array}{ccc}\hfill {\tilde{R}}^{\mu \nu }& \equiv & 2T_{matt}^{\mu \nu }+2{\chi }_{\omega }{{\rm Y}}^{\mu \nu }+\left(p_{\mathsf{\Lambda }}+{\mathcal{E}}_{rot}-\varrho c^2\right)g^{\mu \nu }\hfill \end{array}$$

(30)

$$\begin{array}{ccc}\hfill \tilde{R}& \equiv & {\tilde{R}}^{\mu \nu }{g}_{\mu \nu }=4p_{\mathsf{\Lambda }}+2{\mathcal{E}}_{rot}-2\varrho c^2\hfill \end{array}$$

(31)

$$\begin{array}{ccc}\hfill {\tilde{G}}^{\mu \nu }& \equiv & {\tilde{R}}^{\mu \nu }-{\displaystyle \frac{\tilde{R}}{2}}{\displaystyle \frac{4}{\mathbb{k}}}{\mathbb{k}}^{\mu \nu }={\tilde{R}}^{\mu \nu }-{\displaystyle \frac{\tilde{R}}{2}}{g}^{\mu \nu }-{\displaystyle \frac{\tilde{R}}{2p_{\mathsf{\Lambda }}}}{{\rm Y}}^{\mu \nu }\hfill \\ & =& 2T_{matt}^{\mu \nu }-\left(2-{\chi }_m-{\chi }_{\omega }\right){{\rm Y}}^{\mu \nu }-p_{\mathsf{\Lambda }}{g}^{\mu \nu }\hfill \end{array}$$

(32)

$$\begin{array}{ccc}\hfill {\Theta }^{\mu \nu }& \equiv & -\left({\chi }_m+{\chi }_{\omega }\right){{\rm Y}}^{\mu \nu }\to {\tilde{G}}^{\mu \nu }+p_{\mathsf{\Lambda }}{g}^{\mu \nu }=2T^{\mu \nu }+{\Theta }^{\mu \nu }\hfill \end{array}$$

(33)

where the last equality holds in any considered spacetime and for $G^{\mu \nu }={\displaystyle \frac{\kappa }{2}}{\tilde{G}}^{\mu \nu }$, $\mathsf{\Lambda }={\displaystyle \frac{\kappa }{2}}{p}_{\mathsf{\Lambda }}$ becomes the classical GR equation in curved spacetime (in curved spacetime ${{\rm Y}}^{\mu \nu }$ and ${\Theta }^{\mu \nu }$ vanish).

The Lagrangian ${\mathcal{L}}_{\Theta }$ for ${\Theta }^{\mu \nu }$ may be obtained the same way as in [27] with the use of the interpolating path method $g^{\mu \nu }\left(\lambda \right)=(1-\lambda ){\mathbb{k}}^{\mu \nu }+\lambda g^{\mu \nu }$. Using this method one obtains

$$\sqrt{-g}{\mathcal{L}}_{\Theta }\equiv {\displaystyle \frac{1}{2}}{\int }_0^{1}d\lambda \sqrt{-g\left(\lambda \right)}{\Theta }_{\mu \nu }\left[g^{\mu \nu }\left(\lambda \right)\right]{\partial }_{\lambda }{g}^{\mu \nu }\left(\lambda \right)$$

(34)

Since the variation of the functional is located on the boundary $\lambda =1$, thus ${\Theta }_{\mu \nu }={\displaystyle \frac{2}{\sqrt{-g}}}{\displaystyle \frac{\delta \left(\sqrt{-g}{\mathcal{L}}_{\Theta }\right)}{\delta g^{\mu \nu }}}$.

In the last equation of (33) considered in curved spacetime $(g^{\mu \nu }\to {\mathbb{k}}^{\mu \nu })$ field invariant $p_{\mathsf{\Lambda }}$ acts as double the vacuum energy density (vacum pressure). This equation may be derived from the Lagrangian density ${\mathcal{L}}_G$ in the form

$${\mathcal{L}}_G\equiv 2{\mathcal{L}}_T+{\mathcal{L}}_{\Theta }={\mathcal{L}}_{matt}+{\displaystyle \frac{1}{2}}\tilde{R}-2{\mathcal{E}}_{rot};{\mathcal{L}}_{matt}\equiv {\mathcal{L}}_{\Theta }-p_{\mathsf{\Lambda }}\left({\chi }_m+{\chi }_{\omega }\right)$$

(35)

where the variation by Hilbert’s method on ${\mathcal{L}}_{matt}$ gives $T_{matt}^{\mu \nu }\equiv \varrho U^{\mu }{U}^{\nu }-{\mathsf{\Xi }}^{\mu \nu }$.

The equation (32) considered in curved spacetime simplifies to ${\tilde{G}}^{\mu \nu }={\tilde{R}}^{\mu \nu }-{\displaystyle \frac{\tilde{R}}{2}}{g}^{\mu \nu }$ where ${{\rm Y}}^{\mu \nu }$ and thus also ${\Theta }^{\mu \nu }$ vanishes. It is also worth noting that in flat spacetime the generalized Einstein tensor is associated with the four-divergence of ${\Theta }^{\mu \nu }$

$${\partial }_{\mu }{\tilde{G}}^{\mu \nu }={\partial }_{\mu }{\Theta }^{\mu \nu }=f_{gr}^{\nu }-f_{rr}^{\nu }$$

(36)

so the curvature it describes in curved spacetime replaces this four-force density, where $f_{gr}^{\nu }$ is related to gravity and $f_{rr}^{\nu }$ is the density of radiation-reaction four-force. The presence of the radiation-reaction force has already been discussed in previous works [25], and it now prevents the matter energy and rotational energy from increasing beyond the maximum energy density $p_{\mathsf{\Lambda }}$ available in the system.

It is worth noting that internal energy density and rotational energy, essentially exhaust the possible forms of energy that can be attributed to material bodies (other forms of energy, e.g., chemical energy, can be treated as their components) which could be present in the radiation reaction force. This means that the model proposed here seems complete (with the possible extension of $\alpha$ to a tensor form for more complex systems) and should allow for reproducing the results obtained from GR, as well as reproducing observational results that are inconsistent (such as the dark sector) with currently used interpretation of GR.

Alena Tensor approach therefore allows to look at Einstein’s equations in a new light and analyze the possibilities of explaining the dark sector in a consistent mathematical framework that allows analysis in both flat and curved spacetime. Importantly, it is also possible to analyze the system using a quantum approach (in the description for flat spacetime) and to use standard tools of continuum mechanics for continuous media in flat and curved spacetime, where the description of the behavior of matter has been separated into effects related to fundamental interactions $f_{YM}^{\nu }$, gravity and radiation reaction $f_{gr}^{\nu }-f_{rr}^{\nu }$, and forces related to the distribution of matter $f_{\mathsf{\Xi }}^{\nu }$.

3. Results

The following section will present the results of applying the Alena Tensor model both to cosmological objects and to describe quantum vortices.

3.1. The Halo Effect

The obtained results de facto means, that Alena Tensor ensures correct operation of the standard continuum mechanics equations and GR equations (Euler equations, EOS, TOV, first integrals for rotating stars, etc.), with assumption that the energy density used in them is $\epsilon =\varrho c^2+2{\mathcal{E}}_{rot}$ and pressure is equal to ${\mathcal{E}}_{rot}$. In the next steps, this approach will be analyzed to show that it leads to consistency with the observational data.

It’s worth starting with a simple approximation. Denoting $u_{rot}$ as rotational velocity and assuming

$$p_{\mathsf{\Lambda }}\alpha ={\displaystyle \frac{1}{c^2\kappa }};{\varrho }_{rot}\equiv {\displaystyle \frac{2{\mathcal{E}}_{rot}}{c^2}}={\displaystyle \frac{{\omega }^{\alpha \beta }{\omega }_{\alpha \beta }}{8\pi G}}\to \underset{r\to \infty }{lim}{\varrho }_{rot}={\displaystyle \frac{u_{rot}^2}{4\pi G{r}^2}}$$

(37)

one may notice, that $p_{\mathsf{\Lambda }}\alpha$ in (14) plays the role of the density of the moment of inertia, while ${\varrho }_{rot}$ increases the body’s effective mass within its own frame. This would allow to consider galaxies as continuous media, where the effective mass $M_{eff}$ and its density ${\varrho }_{eff}$ responsible for gravity ${\varrho }_{eff}\equiv {\displaystyle \frac{\epsilon }{c^2}}=\varrho +{\varrho }_{rot}$ from (16) increases with the galactic disk size and angular velocity, causing the halo effect.

For far regions, denoting $M_b$ as baryonic mass, for spherical symmetry one obtains from Poisson’s equation simple linear ODE in the Newtonian limit. In the far regions it could determine a constant rotation speed and might be used to measure of deviation from the vacuum solution.

Vacuum solution in curved spacetime from (15) and (16) yields

$$\begin{array}{c}\hfill \tilde{R}=0;p_{\mathsf{\Lambda }}=2\epsilon \to \varrho c^2=-3{\mathcal{E}}_{rot};p_{\mathsf{\Lambda }}=-2{\mathcal{E}}_{rot};\mathsf{\Lambda }=-\kappa {\mathcal{E}}_{rot};0={\tau }^{\mu \nu }=U^{(\mu }{q}^{\nu )};{\mathcal{L}}_T=0\end{array}$$

(38)

This means that Kepplerian profiles are still possible for systems that can be approximated by a vacuum solution or does not rotate $({\mathcal{E}}_{rot}\approx 0)$.

Going into a more detailed analysis, one may consider $\left(0={\tau }^{\mu \nu }=U^{(\mu }{q}^{\nu )}\right)$ as a simple ideal fluid system with pressure $p\equiv {\mathcal{E}}_{rot}$ according to (16) and (19). In the GR equations one may consider above with an axisymmetric, spherical metric

$$\begin{array}{ccc}\hfill d{s}^2& =& N^2{c}^{2}d{t}^2-A^2\left(d{r}^2+r^{2}d{\theta }^2\right)-B^2{r}^2{sin}^2\theta {\left(d\varphi -\omega dt\right)}^2\hfill \\ \hfill T_{ideal}^{\mu \nu }& \equiv & \left(\varrho +{\displaystyle \frac{2p}{c^2}}\right)U^{\mu }{U}^{\nu }-p{\Delta }^{\mu \nu };\mathsf{\Omega }\equiv {\displaystyle \frac{d\varphi }{dt}}\hfill \end{array}$$

(39)

Analyzing Euler’s energy and momentum equations in ${\nabla }_{\mu }{T}_{ideal}^{\mu \nu }=0$ one may notice, that $u_{ZAMO}$ is not a geodetic movement, and on the equator it takes the value

$$\begin{array}{c}\hfill {\displaystyle \frac{u_{ZAMO}^2}{c^2}}={\displaystyle \frac{{\left(Br\right)}^2}{N^2}}{\left(\mathsf{\Omega }-\omega \right)}^2=r{\partial }_{r}lnN+{\gamma }_p^2-1;{\gamma }_p^2\equiv 1-r{\displaystyle \frac{p^{\prime }}{\varrho c^2+3p}}\end{array}$$

(40)

where ${\gamma }_p$ coefficient determines the deviation from the geodetic. One may therefore define pressure according to conclusions from previous section as follows, what yields

$$\begin{array}{c}\hfill \kappa p\equiv {\displaystyle \frac{u_{ZAMO}^2}{r^2{c}^2}}\to {\displaystyle \frac{d{M}_{eff}}{dr}}=4\pi r^2\varrho \left(r\right)+{\displaystyle \frac{u_{ZAMO}^2}{G}}\to {\displaystyle \frac{d}{dr}}{u}_{ZAMO}^2=4\pi Gr\varrho \left(r\right)\end{array}$$

(41)

In obtained picture the velocity increase depends solely on the baryon mass distribution, while the flattening of the tail is maintained by the rotational energy. This precisely corresponds to the expected behavior of a dark matter "halo."

The introduction of $q^{\mu }$ and ${\tau }^{\mu \nu }$ into a system can be approximated by defining parameter $\chi \left(r\right)$ changing oryginal p used for isotropic model

$$\begin{array}{c}\hfill p_r\equiv \alpha p;p_{\theta }\equiv \beta p;p_{\varphi }\equiv p\to \chi =\alpha +\beta ;3p_{\chi }\to (1+\chi )p;2p\to \chi p\end{array}$$

(42)

In practice, even a constant $\chi$ should be sufficient for analyzing the fit of galaxy rotation curves. Using a constant $\chi$ also provides a simpler ODE and the ability to quickly perform preliminary fits of $\chi$ to observational data for large amounts of data.

$$\begin{array}{ccc}\hfill {\displaystyle \frac{d{M}_{eff}}{dr}}& =& 4\pi r^2\varrho \left(r\right)+{\displaystyle \frac{\chi }{2}}{\displaystyle \frac{u_{ZAMO}^2}{G}}\to {\displaystyle \frac{d}{dr}}{u}_{ZAMO}^2=4\pi Gr\varrho \left(r\right)+\left({\displaystyle \frac{\chi }{2}}-1\right){\displaystyle \frac{u_{ZAMO}^2}{r}};\hfill \\ \hfill \underset{\varrho \to 0}{lim}{u}_{ZAMO}^2& \propto & r^{{\displaystyle \frac{\chi }{2}}-1}\hfill \end{array}$$

(43)

Below one may find the expected course of the rotation curves depending on the assumed constant $\chi$. The calculations used an averaged Hernquist bulge baryon model ${\varrho }_{bulge}$ and a ”spherical proxy” of the exponential disk ${\varrho }_{disc}$

$$\begin{array}{ccc}\hfill {\varrho }_{bulge}& =& {\displaystyle \frac{M_{bulge}}{2\pi }}{\displaystyle \frac{a}{r{(r+a)}^3}};M_{bulge}={10}^{10}{M}_{\odot };a=1kpc\hfill \end{array}$$

(44)

$$\begin{array}{ccc}\hfill {\varrho }_{disc}& =& {\displaystyle \frac{M_{disc}}{4\pi R_d^2}}{\displaystyle \frac{e^{-r/R_d}}{r}};R_d=3kpc\hfill \end{array}$$

(45)

with total baryon density ${\varrho }_b\left(r\right)={\varrho }_{bulge}+{\varrho }_{disc}$, standard G value, anisotropy and energy stream simulated by constant $\chi$. As can be seen from the graph, the increasing anisotropy towards the outskirts of the galaxy $\chi (r=0)=0;\underset{r\to \infty }{lim}\chi \left(r\right)=const$ would allow the graph to align with the expected curve shapes for spiral galaxies.

Figure 1. Approx . rotation curves in Alena Tensor model

Figure 1. Approx . rotation curves in Alena Tensor model

As it appears, constant $\chi$ actually allows to tune the rotation velocity distribution for some part of galaxies what gives an overview of the method and may help in further analysis and tuning of $\chi \left(r\right)$ function to achieve full agreement with observations. The A presents the results of tuning the $\chi$ constant for about 100 galaxies from the SPARC catalog. The results seem encouraging, and worth further analysis with the $\chi \left(r\right)$ function or with the full $T^{\mu \nu }$ tensor representation. For most galaxies the fit is very good, and for the rest it is obtainable by introducing a simple function $\chi \left(r\right)$ which provides $\chi \left(r\right)\approx 0$ for small r (chaotic motions in the center, no ordered rotation) to a stabilized rotation in the outskirts with a fixed $\chi$.

The $\chi$ parameter was adjusted for each galaxy by iteratively solving the radial motion differential equation resulting from (43), with the condition of normalizing the rotation rate in the outer disk ${\langle V_{\chi }\rangle }_{\mathrm{outer}}={\langle V_{\mathrm{obs}}\rangle }_{\mathrm{outer}}$. The entire procedure, including reading the rotation data from the SPARC catalog, interpolating the baryonic component $V_{\mathrm{bar}}\left(r\right)$, numerically solving the equation for $V_{\chi }(r,\chi )$ and saving the resulting plots in PDF format, was fully automated in the Mathematica script, and run in an environment with the SPARC source data. The script is included in the supplementary files.

3.2. Quantum Vortices and Elementary Particles

The obtained Alena Tensor equations allow for the analysis of vortices in a quantum model, it is therefore worth analyzing their consequences. In the next steps it will be shown that by treating the Lagrangian (29) as an effective field model and choosing an appropriate stiffness function $A\left(\rho \right)$, both Yukawa-type mass generation and Higgs-type potential structure can emerge naturally.

The vorticity tensor enters the Lagrangian (29) only through the elastic term and the spin-vorticity interaction. Its algebraic equation of motion yields

$${\displaystyle \frac{\partial {\mathcal{L}}_{\mathrm{eff}}}{\partial {\omega }^{\mu \nu }}}=0\to {\omega }_{\mu \nu }=c^2\kappa e^{A\left(\rho \right)}g\hslash \overline{\psi }{\Sigma }_{\mu \nu }\psi$$

(46)

so that eliminating ${\omega }_{\mu \nu }$ produces the effective fermion term

$${\mathcal{L}}_f\equiv -{\displaystyle \frac{g^2{\hslash }^2{c}^2\kappa e^{A\left(\rho \right)}}{2}}\left(\overline{\psi }{\Sigma }^{\mu \nu }\psi \right)\left(\overline{\psi }{\Sigma }_{\mu \nu }\psi \right)=-{\mathcal{E}}_{rot}$$

(47)

Evaluated at the homogeneous equilibrium, this yields a dynamically generated fermion mass

$$m_{\mathrm{eff}}={\displaystyle \frac{g^2{\hslash }^2\kappa e^{A\left(\rho \right)}}{2}}\langle \Sigma \rangle ;\langle \Sigma \rangle \equiv {\displaystyle \frac{\langle \left(\overline{\psi }{\Sigma }^{\mu \nu }\psi \right)\left(\overline{\psi }{\Sigma }_{\mu \nu }\psi \right)\rangle }{\langle \overline{\psi }\psi \rangle }}$$

(48)

In this picture, the mass generation results from the equilibrium between the scalar amplitude $\rho$ and the spin-vorticity background, rather than from a fundamental Yukawa coupling where the condensate $\langle \Sigma \rangle$ plays the role of an order parameter. For a stationary vortex configuration $\phi =\rho e^{i\varphi }$, the scalar sector of (29) leads to the standard energy functional

$$E[\rho ,\varphi ]={\displaystyle \frac{\hslash c}{2}}\int d^{3}x\left[{(\nabla \rho )}^2+{\rho }^2{(\nabla \varphi )}^2\right]$$

(49)

Euler-Lagrange equations dedived from (29) are

$$\hslash c{\partial }_{\mu }{\partial }^{\mu }\rho -\left[\hslash c\rho \left({\partial }_{\mu }\varphi \right)\left({\partial }^{\mu }\varphi \right)-{\displaystyle \frac{A^{\prime }\left(\rho \right)e^{-A\left(\rho \right)}}{2c^2\kappa }}{\omega }_{\alpha \beta }{\omega }^{\alpha \beta }\right]=0;{\partial }_{\mu }\left({\rho }^2{\partial }^{\mu }\varphi \right)=0$$

(50)

and show, that the amplitude adjusts to the local phase gradient. In the presented model, $\rho$ has the dimension of the inverse of the length while Standard Model vev $v_{\mathrm{SM}}$ is expressed as energy. Considering the equations in the limit $\rho =v$ one may thus assume that $\hslash cv=v_{\mathrm{SM}}\sqrt{2}$. Since the effective Yukawa parameter can be extracted directly from (48) therefore in order to reproduce the Standard Model it schould satisfy $m_{\mathrm{eff}}{c}^2=y_{\mathrm{f}}{\displaystyle \frac{\hslash c}{2}}v$, thus

$${\displaystyle \frac{\hslash c}{2}}{y}_{\mathrm{f}}\equiv {\left.{\displaystyle \frac{\partial m_{\mathrm{eff}}\left(\rho \right)c^2}{\partial \rho }}\right|}_{\rho =v}=m_{\mathrm{eff}}{c}^2{A}^{\prime }\left(v\right)\to {\displaystyle \frac{\hslash c{y}_{\mathrm{f}}}{2m_{\mathrm{eff}}{c}^2}}={\displaystyle \frac{1}{v}}\to A^{\prime }\left(v\right)={\displaystyle \frac{1}{v}}$$

(51)

Assuming the stationary energy-minimizing configuration (static state)

$$\rho \left(x\right)=v=\mathrm{const};{\partial }_i\rho =0;\varphi \left(x\right)=\mu ct;{\partial }_i\varphi =0$$

(52)

the Euler-Lagrange equation for $\rho =v$ reduces therefore to

$${\mathcal{E}}_{rot}\left(v\right)=m_{\mathrm{eff}}{c}^2\langle \overline{\psi }\psi \rangle ={\displaystyle \frac{e^{-A\left(v\right)}}{2c^2\kappa }}{\omega }_{\alpha \beta }{\omega }^{\alpha \beta }=\hslash c{\displaystyle \frac{v{\mu }^2}{A^{\prime }\left(v\right)}}=\hslash c{v}^2{\mu }^2$$

(53)

Under suitable assumptions on the stiffness function $A\left(\rho \right)$ and on the spin-vorticity condensate encoded in ${\omega }_{\alpha \beta }{\omega }^{\alpha \beta }$, this relation plays the role of a gap equation, equivalent to the condition $\partial V_{\mathrm{eff}}{\left(\rho \right)/\partial \rho |}_{\rho =v}=0$ for an effective Higgs-like potential. The competition between ${\mathcal{E}}_{rot}$ and $\hslash c{v}^2{\mu }^2$ might therefore reproduce the characteristic symmetry-breaking structure of the Higgs potential.

Assuming that the effective potential $V_{\mathrm{eff}}\left(\rho \right)$ is indeed of Higgs type, the mass of the small radial (Higgs-like) excitation around the equilibrium value $\rho =v$, using emergent mass formula (48) and using the condensate identity (46), may be obtained in a compact expression

$$m_H^2\equiv {\displaystyle \frac{\hslash }{c^3}}{\left.{\displaystyle \frac{{\partial }^2{V}_{\mathrm{eff}}\left(\rho \right)}{\partial {\rho }^2}}\right|}_{\rho =v}={\displaystyle \frac{\hslash }{c}}{m}_{\mathrm{eff}}\langle \overline{\psi }\psi \rangle \left[A^{″}\left(v\right)+{\left(A^{\prime }\left(v\right)\right)}^2-{\displaystyle \frac{A^{\prime }\left(v\right)}{v}}\right]={\displaystyle \frac{\hslash }{c}}{m}_{\mathrm{eff}}\langle \overline{\psi }\psi \rangle A^{″}\left(v\right)$$

(54)

Therefore (54) may be further simplified, imposing a specific form on the function $A\left(\rho \right)$ as below, for the preservation of the Yukawa condition and the Higgs-like potential

$$A\left(\rho \right)\equiv ln{\displaystyle \frac{\rho }{v}}+{\displaystyle \frac{\lambda }{2v^2}}{(\rho -v)}^2\to {\mathcal{E}}_{rot}\left(v\right)=m_{\mathrm{eff}}{c}^2\langle \overline{\psi }\psi \rangle ={\displaystyle \frac{m_H^2{c}^3}{\hslash }}{\displaystyle \frac{v^2}{\lambda -1}}$$

(55)

Since in the considered model, the relation between Standard Model vev $v_{\mathrm{SM}}$ and v is $\hslash cv=v_{\mathrm{SM}}\sqrt{2}$, therefore introducing the Standard Model Higgs self-coupling as ${\lambda }_{\mathrm{SM}}$ and substituting $m_{\mathrm{eff}}{c}^2$ from (51) one obtains

$$\hslash cv=v_{\mathrm{SM}}\sqrt{2}={\displaystyle \frac{m_H{c}^2}{\sqrt{{\lambda }_{\mathrm{SM}}}}}\to {\displaystyle \frac{y_{\mathrm{f}}}{2}}·\langle \overline{\psi }\psi \rangle ={\displaystyle \frac{{\lambda }_{\mathrm{SM}}}{\lambda -1}}·v^3$$

(56)

Therefore, making the natural assumption $\langle \overline{\psi }\psi \rangle =v^3$, one may separate condensate from Yukawa parameter as follows

$$m_{\mathrm{eff}}=m_H{\displaystyle \frac{\sqrt{{\lambda }_{\mathrm{SM}}}}{\lambda -1}}$$

(57)

replacing the Yukawa parameter $y_{\mathrm{f}}$ with the parameter $\lambda$ describing the quantum vortex. In this way the above approach allows to model elementary particles as stable quantum vortices, reproducing the successful Yukawa and Higgs mechanism, where all stability requirements follow directly from the field equations and the structure of the effective potential.

In the resulting framework, elementary particles arise as stable, self-sustained vortices of the underlying vacuum field, whose rotational structure is captured by the antisymmetric vorticity tensor ${\omega }_{\mu \nu }$. This interpretation is fully consistent with the general paradigm of topological solitons, where localized and finite-energy excitations of a continuous field behave as particle-like objects, as in the Skyrme model [33], the Faddeev-Niemi hopfion model [34,35], and the superfluid-vacuum approach of Volovik [36]. In presented case, the spin-vorticity coupling and the rotational energy functional play the role of stabilizing terms, ensuring that the vortex configuration propagates as a massive fermionic quasiparticle, closely analogous to knotted vortical solitons in classical and quantum field theories [34,37].

One may notice, that the effective Lagrangian (29) naturally supports four distinct conserved charges, arising from its symmetries and from the structure of the vorticity sector

• Phase (Noether) charge $Q_{\mathsf{\Phi }}=\int d^{3}x{j}^0;j^{\mu }\equiv {\rho }^2{\partial }^{\mu }\varphi$ originating from the global shift symmetry $\mathsf{\Phi }\to \mathsf{\Phi }+\mathrm{const}$. It corresponds to the conserved circulation associated with the phase field.
• Topological vortex number $N_v={\displaystyle \frac{1}{2\pi }}{\oint }_{\mathcal{C}}\nabla \varphi ·d\overrightarrow{\ell }\in \mathbb{Z}$ defined for static configurations with nontrivial winding of the phase $\varphi$ around the vortex core. This integer counts the number of $2\pi$ windings.
• Spin-vorticity charge $Q_{\mathrm{sv}}=\int d^{3}x{\partial }_{\mu }\left(g\hslash {\omega }^{\mu 0}\right)$ where the vorticity tensor ${\omega }_{\mu \nu }$ satisfies the algebraic field equation ${\omega }_{\mu \nu }=c^2\kappa e^{A\left(\rho \right)}g\hslash \overline{\psi }{\Sigma }_{\mu \nu }\psi$. This charge reflects the conserved flow associated with the spin-vorticity coupling term $g\hslash \overline{\psi }{\Sigma }_{\mu \nu }\psi {\omega }^{\mu \nu }$.
• Hopf (linking) charge $Q_H={\displaystyle \frac{1}{32{\pi }^2}}\int d^{3}x{ϵ}^{ijk}{\mathcal{A}}_i{\mathcal{F}}_{jk}$ defined when the dual vorticity vector ${\omega }_i\equiv \frac{1}{2}{ϵ}_{ijk}{\omega }_{jk}$ is normalized to a unit field $\overrightarrow{n}\left(x\right)=\overrightarrow{\omega }/\left|\overrightarrow{\omega }\right|$, with ${\mathcal{F}}_{ij}={\partial }_i{\mathcal{A}}_j-{\partial }_j{\mathcal{A}}_i$ denoting the pullback of the area form on $S^2$. This integer-valued invariant characterizes the knotting and linking of vorticity lines.

The presence of these conserved quantities enables a deeper analysis of the model in the context of reproducing known structural features of elementary particles, including the possibility of identifying particle families with different topological sectors.

It may be also noticed, that assuming non-zero m, variation of (29) with respect to $\overline{\psi }$ gives the modified Dirac equation

$$\left(i\hslash c{\gamma }^{\mu }{D}_{\mu }-m{c}^2-g\hslash {\Sigma }^{\mu \nu }{\omega }_{\mu \nu }\right)\psi =0$$

(58)

where the last term is a local spin-vorticity coupling formally analogous to the Mashhoon effect [38,39]. For the normalization of ${\omega }_{\mu \nu }$ assumed in (12) this correspondence fixes $g={\displaystyle \frac{1}{2}}$. This thus allows the obtained model to be used to describe the atom and for other classical applications.

4. Discussion and Conclusions

It is worth discussing the conclusions of this article by dividing them into issues concerning GR/Cosmology and quantum issues.

4.1. Discussion and Conclusions Regarding GR and Cosmology

As seen in the above article, supplementing the Alena Tensor with the energy associated with the rotation of bodies naturally leads to the creation of halo effects, known from dark matter studies. Preliminary analysis allows for a fairly good match of this effect to observational results, although this obviously requires further development and verification for a larger number of cosmological objects. Importantly, the proposed approach does not require modification of the GR equations, but rather fits naturally into the applied GR equations and continuum mechanics. Since the source is described in this approach by $\epsilon =\varrho c^2+2{\mathcal{E}}_{rot}$, this means that the observed increase in effective mass also affects gravitational lensing to an extent precisely corresponding to the increase in effective mass by the energy associated with (in this case - rotational) ”dark matter”. This is precisely what is obtained in observations [40,41].

The proposed solution fits quite well with the research direction represented by [42,43,44] and also [45] (including baryotropy), who investigated anisotropic fluid in cosmology and its potential connections with the dark sector. However, it complements these studies with the natural halo effect resulting directly from the GR equations for the Alena Tensor. The proposed model also expands and, in a sense, substantiates the hypothesis posed by C. Rourke [46], complementing the research [47,48,49,50] with a justification for linking rotation with the halo effect. The idea that rotation-related effects can mimic dark matter is not new, but Alena Tensor gives it some additional structure, making it a direct consequence of a coherent mathematical model.

Importantly, the Alena Tensor also provides a natural interpretation of dark energy. The value of $p_{\mathsf{\Lambda }}$ is an invariant of the field tensor and becomes constant (or, at least, metric-independent invariant) in curvilinear description. In a sense, a nonzero value of $\mathsf{\Lambda }$ can therefore be interpreted as a scale of deviation from pure wave solutions, without matter (for example, for the electromagnetic field, $p_{\mathsf{\Lambda }}=0$ would mean that the electric and magnetic fields are equal, so the solutions must be pure electromagnetic waves). Since the value of $p_{\mathsf{\Lambda }}$ measured in flat spacetime is $p_{\mathsf{\Lambda }}=p_o{\left({\mathbb{k}}_{\mu \nu }{\eta }^{\mu \nu }\right)}^2$, it is a measure of the "flatness" of spacetime, or more precisely, a measure of how much the metric tensor for the curvilinear description deviates from the Minkowski tensor. This interpretation seems particularly interesting in the context of the works [51,52], because it strengthens and details the conclusions described therein, providing a geometric, anisotropic source that can be interpreted as a specific backreaction mechanism leading to acceleration.

This approach could be applied to many other continuous systems (e.g., stars or black holes) and seems worth to describe the extreme in which $\varrho =0$ and all the energy in the system is rotational energy. Although at first glance it seems absurd, one may notice, that the source in vacuum solution (38) is indeed solely rotational energy. As shown in 3.2, replacing rotation with vorticity and treating it as a consequence of the circulating field creates the possibility to model elementary particles as quasi-stable systems of three-dimensional vortices.

4.2. Discussion and Conclusions Regarding Quantum Issues

In quantum theory, the same structure of the energy-momentum tensor gives rise to topological vortices whose stability determines states with interpretable mass. Obtained interpretation of the Higgs mechanism creates the possibility of further analysis of the properties of elementary particles depending on properties of the associated vortices, because in the obtained picture the particle masses result from the self-consistent spin condensation (gap-equation), not from the external Yukawa condition. Hence, the model admits particle-like, topologically stabilized vortices, naturally emerging from the improved Belinfante tensor structure.

It is thus possible, taking into account the conclusions from these analyses, to further develop the idea of rotation in quantum systems. In quantum field theory systems with significant local vorticity, the inclusion of a rotational-energy contribution in the energy-momentum tensor implies a new spin-vorticity coupling term, which may lead to measurable spin polarization of fermions as predicted e.g. in heavy-ion physics [53]. Such additional rotational energy terms effectively modify the Hamiltonian density and thus alter propagators or dispersion relations of excitations in the medium, analogous to how vorticity-induced spin currents appear in hydrodynamic and condensed-matter contexts as in [54]. Since in considered approach the modified tensor now couples orbital angular momentum and intrinsic spin via the rotational energy density, therefore the renormalization group flows and transport coefficients (like shear viscosity or spin diffusion) may pick up new contributions, offering a potential probe of vorticous quantum media such as the quark-gluon plasma [55].

The extended Alena Tensor naturally suggests a broad research programme connecting microscopic particle physics, topological soliton theory and relativistic hydrodynamics. The modified Dirac equation (58) exhibits a local coupling of the intrinsic spin to the vorticity tensor ${\omega }_{\mu \nu }$. In the nonrelativistic limit this reduces to a Mashhoon-like term $-\mathsf{\Omega }·(\mathbf{L}+\mathbf{S})$, in agreement with earlier analyses of inertial couplings in the Dirac theory [56,57,58]. This provides a direct physical interpretation of the Alena vorticity field as an observable rotational background.

The emergence of spin-vorticity coupling suggests several measurable signatures. These include possible energy-level shifts in atomic systems, modifications of fermion propagators in rotating backgrounds, and a connection to the experimentally observed spin polarization in the quark-gluon plasma (QGP) [59,60]. Since heavy-ion experiments have already demonstrated substantial vorticity and spin alignment of $\mathsf{\Lambda }$ hyperons, the Alena Tensor offers a concrete framework to reinterpret such data and potentially extract constraints on the coupling parameter g.

The effective theory possesses four distinct conserved quantities: the phase charge $Q_{\mathsf{\Phi }}$, the vortex number $N_v$, the spin-vorticity charge $Q_{\mathrm{sv}}$ and the Hopf charge $Q_H$. These invariants classify vortex solutions according to their topology and encode information about stability and possible decay channels. Such a structure closely parallels the theory of topological solitons, including Skyrme-type models [33], Faddeev-Niemi hopfions [34] and modern knot-like soliton constructions [61,62]. In this context, the Alena Tensor model suggests a possible reinterpretation of elementary particles as quantized vortex excitations distinguished by their topological invariants. For instance, fermions may correspond to stable spin-carrying vortices, while different particle families could be captured by distinct Hopf sectors. This motivates the construction of fully three-dimensional vortex solutions with prescribed $(Q_{\mathsf{\Phi }},N_v,Q_{\mathrm{sv}},Q_H)$ and the analysis of their energetics and stability.

Eliminating the antisymmetric field ${\omega }_{\mu \nu }$ leads to an effective four-fermion interaction term and a gap equation for the modulus $\rho$. This mechanism resembles dynamical mass generation in Nambu-Jona-Lasinio models [63], composite Higgs scenarios [64] and superfluid-vacuum analogies [36]. In this picture, particle masses arise not from fundamental Yukawa couplings but from vortex stability and the curvature of the effective potential encoded in $A\left(\rho \right)$. The Higgs particle appears as a radial excitation of the condensate, and its mass follows from differentiating the effective potential. This approach suggests a unified mechanism for fermion and Higgs masses, and establishes a link between soliton stability and the observable mass spectrum. Determining the function $A\left(\rho \right)$ from either phenomenological constraints or first-principles requirements is therefore a key target for future investigations.

The three directions outlined above: spin-vorticity dynamics, topological classification of vortex excitations and dynamical mass generation, form a coherent research strategy in which microscopic and macroscopic aspects of the Alena Tensor are treated on equal footing. Future work may include, e.g.:

• computing atomic and astrophysical signatures of the modified Dirac equation,
• constructing fully nonlinear three-dimensional vortex solitons with conserved topological charges,
• deriving fermion mass spectra from vortex equilibrium and comparing with Standard Model data,
• exploring possible vortex-induced contributions to anisotropic stress in neutron stars [65,66],
• studying hydrodynamic transport and spin polarization in vortical media [67,68].

Equally interesting direction of further analysis could be e.g. the use of the possibilities of quantum description of the dark sector in the Alena Tensor model, for further development of works such as [69]. It also seems that describing matter (e.g. a neutron star, as in [70,71]) using the mechanism proposed here for GR, would be the simplest way to confirm or falsify the Alena Tensor, due to the high symmetry of such a solution. However, all these analyses deserve separate articles.

In conclusion, it remains an open question whether the Alena Tensor is a correct way to describe physical systems, but this paper shows that, beyond the compliances with available knowledge achieved so far, it naturally leads to the existence of halo effects, interpretation of dark energy and modeling quantum vortices. Taken together, these results indicate that the extended Alena Tensor offers a unified geometric and topological framework that connects elementary particles, relativistic fluids and large-scale astrophysical structures. This opens a way for theoretical and phenomenological studies, extending far beyond the cosmological applications emphasized in the present work.

The author hopes that the results obtained in this paper will facilitate further use and development of the discussed approach and, potentially, many similar concepts. It also seems that further analysis of Alena Tensor may provide useful descriptions of the transformation between curved and flat spacetime and bring new insights that will contribute to a better understanding of issues related to the broadly understood unification of physical theories.

5. Statements

All data that support the findings of this study are included within the article (and any supplementary files).

During the preparation of this work the author did not use generative AI or AI-assisted technologies, except for continuous learning.

Author did not receive support from any organization for the submitted work.

Author have no relevant financial or non-financial interests to disclose.