<<

Article An Alternative to Dark and Dark : Scale-Dependent in Superfluid Theory

Konstantin G. Zloshchastiev

Institute of Systems , Durban University of Technology, P.O. Box 1334, Durban 4000, South Africa; [email protected]

 Received: 29 August 2020; Accepted: 10 October 2020; Published: 15 October 2020 

Abstract: We derive an effective gravitational potential, induced by the quantum wavefunction of a physical vacuum of a self-gravitating configuration, while the vacuum itself is viewed as the superfluid described by the logarithmic quantum equation. We determine that gravity has a multiple-scale pattern, to such an extent that one can distinguish sub-Newtonian, Newtonian, galactic, extragalactic and cosmological terms. The last of these dominates at the largest length scale of the model, where superfluid vacuum induces an asymptotically Friedmann–Lemaître–Robertson–Walker-type , which provides an explanation for the accelerating expansion of the Universe. The model describes different types of expansion mechanisms, which could explain the discrepancy between measurements of the Hubble constant using different methods. On a galactic scale, our model explains the non-Keplerian behaviour of galactic rotation curves, and also why their profiles can vary depending on the . It also makes a number of predictions about the behaviour of gravity at larger galactic and extragalactic scales. We demonstrate how the behaviour of rotation curves varies with distance from a gravitating center, growing from an inner galactic scale towards a metagalactic scale: A squared orbital ’s profile crosses over from Keplerian to flat, and then to non-flat. The asymptotic non-flat regime is thus expected to be seen in the outer regions of large spiral .

Keywords: ; ; superfluid vacuum; emergent spacetime; ; galactic rotation curve; quantum Bose ; logarithmic fluid; logarithmic wave equation

1. Introduction Astronomical observations over many length scales support the existence of a number of novel phenomena, which are usually attributed to dark matter (DM) and (DE). Dark matter was introduced to explain a range of observed phenomena at a galactic scale, such as flat rotation curves, while dark energy is expected to account for cosmological-scale dynamics, such as the accelerating expansion of the Universe. For instance, the ΛCDM model, which is currently the most popular approach used in cosmology and galaxy-scale , makes use of both DE and cold DM concepts [1]. In spite of being a generally successful framework purporting to explain the large-scale structure of the Universe, it currently faces certain challenges [2,3]. There is also growing consensus that a convincing theory of DM- and DE-attributed phenomena cannot be a stand-alone model; but should, instead, be a part of a fundamental theory involving all known interactions. In turn, we contend that formulating this fundamental theory will be impossible without a clear understanding of the dynamical structure of the physical vacuum, which underlies all interactions that we know of. Moreover, this theory must operate at a quantum level, which necessitates us rethinking of the concept of gravity using basic notions of quantum .

Universe 2020, 6, 180; doi:10.3390/universe6100180 www.mdpi.com/journal/universe Universe 2020, 6, 180 2 of 25

One of the promising candidates for a theory of physical vacuum is superfluid vacuum theory (SVT), a post-relativistic approach to high-energy and gravity. Historically, it evolved from Dirac’s idea of viewing the physical vacuum as a nontrivial quantum object, whose phase and derived velocity are non-observable in a quantum-mechanical sense [4]. The term ‘post-relativistic’, in this context, means that SVT can generally be a non-relativistic theory; which nevertheless contains relativity as a special case, or limit, with respect to some dynamical value such as (akin to being a superset of the Newton’s theory of gravity). Therefore, underlying three-dimensional space would not be physically observable until an observer goes beyond the above-mentioned limit, as will be discussed in more detail later in this article. The dynamics and structure of superfluid vacuum are being studied, using various approaches which agree upon the main paradigm (physical vacuum being a background quantum liquid of a certain kind, and elementary being excitations thereof), but differ in their physical details, such as an underlying model of the liquid [5–7]. It is important to with a precise definition of superfluid, to ensure that we avoid the most common misconceptions which otherwise might arise when one attempts to apply superfluid models to astrophysics and cosmology, some details can be found in AppendixA. In fact, some superfluid-like models of dark matter based on classical perfect fluids, scalar field theories or scalar-tensor , turned out to be vulnerable to experimental verification [8]. Moreover, superfluids are often confused not only with perfect fluids, but also with the concomitant phenomenon of Bose-Einstein condensates (BEC), which is another kind of quantum matter occurring in low-temperature condensed matter [9]. However, even though BEC’s do share certain features with superfluids, this does not imply that they are superfluidic in general. In particular, quantum excitations in laboratory superfluids that we know of have dispersion relations of a distinctive shape called the Landau “roton” spectrum. Such a shape of the spectral curve is crucial, as it ensures the suppression of dissipative fluctuations at a quantum level [10,11], which results in inviscid flow [12,13]. If plotted as an excitation energy versus momentum, the curve starts from the origin, climbs up to a local maximum (called the maxon peak), then slightly descends to a local nontrivial minimum (called the “roton” energy gap); then grows again, this time all the way up, to the boundary of the theory’s applicability range. In fact it is not the roton energy gap alone, but the energy barrier formed by the maxon peak and roton minimum in momentum space, which ensures the above-mentioned suppression of quantum fluctuations in quantum liquid and, ultimately, causes its flow to become inviscid. In other words, it is the global characteristics of the dispersion curve, not just the existence of a nontrivial local minimum and related energy gap, which is important for superfluidity to occur. Obviously, these are non-trivial properties, which cannot possibly occur in all quantum and condensates. Further details and aspects are discussed in AppendixA. This paper is organized as follows. Theory of physical vacuum based on the logarithmic superfluid model is outlined in Section2, where we also demonstrate how four-dimensional spacetime can emerge from the three-dimensional dynamics of quantum liquid. In Section3, we derive the gravitational potential, induced by the logarithmic superfluid vacuum in a given state, using certain simplifying assumptions. Thereafter, in Section4, we give a brief physical interpretation of different parts of the derived gravitational potential and estimate their characteristic length scales. In Section5, profiles of induced matter density are derived and discussed for the case of spherical symmetry. Galactic scale phenomena are discussed in Section6, where the phenomenon of galactic rotation curves is explained without introducing any ad hoc. In Section7, we discuss the various mechanisms of the accelerating expansion of the Universe, as well as the cosmological singularity, “vacuum catastrophe” and cosmological coincidence problems. Conclusions are drawn in Section8. Universe 2020, 6, 180 3 of 25

2. Logarithmic Superfluid Vacuum Superfluid vacuum theory assumes that the physical vacuum is described, when disregarding quantum fluctuations, by the fluid condensate wavefunction Ψ(r, t), which is a three-dimensional Euclidean scalar. The state itself is described by a ray in the corresponding Hilbert space, therefore this wavefunction obeys a normalization condition Z hΨ|Ψi = ρ dV = M, (1) V where M and V are the total and volume of the fluid, respectively, and ρ = |Ψ|2 is the fluid mass density. The wavefunction’s dynamics is governed by an equation of a U(1)-symmetric Schrödinger form: " # h¯ 2 −ih¯ ∂ − ∇2 + V (r, t) + F(|Ψ|2) Ψ = 0, (2) t 2m ext where m is the constituent particles’ mass, Vext(r, t) is an external or trapping potential and F(ρ) is a duly chosen function, which effectively takes into account many-body effects inside the fluid. This wave equation can be formally derived as a minimizing condition of an action functional with the following Lagrangian:

ih¯ h¯ 2 L = (Ψ∂ Ψ∗ − Ψ∗∂ Ψ) + |∇Ψ|2 + V (r, t) |Ψ|2 + V(|Ψ|2), (3) 2 t t 2m ext where V(ρ) equals to a primitive of F(ρ) up to an additive constant: F(ρ) = V0(ρ); throughout the paper the prime denotes a derivative with respect to the function’s argument. In this picture, massless excitations, such as , are analogous to acoustic propagating p 0 with velocity cs ∝ |p (ρ)|, where fluid pressure p = p(ρ) is determined via the . For the system (2), both the equation of state and speed of sound can be derived using the fluid-Schrödinger analogy, which was established for a special case in [14], and generalized for an arbitrary F(ρ) in works [7,15]. In a leading-order approximation with respect to the constant, we obtain 1 Z 1 p = − ρF0(ρ) dρ, c2 = ρ|F0(ρ)|, (4) m s m while higher-order corrections would induce Korteweg-type effects, thus significantly complicating the subject matter [15]. Furthermore, it is natural to require that superfluid vacuum theory must recover Einstein’s at a certain limit. One can show that at a limit of low momenta of quantum excitations, often called a “phononic” limit by analogy with laboratory quantum liquids, Lorentz symmetry does emerge. This can be easily shown by virtue of the fluid/gravity analogy [16], which was subsequently used to formulate the BEC-spacetime correspondence [7]; it can also be demonstrated by using dispersion relations [11,17], which are generally become deformed in theories with non-exact Lorentz symmetry [18–21]. This correspondence states that Lorentz symmetry is approximate, while four-dimensional spacetime is an induced phenomenon, determined by the dynamics of quantum Bose liquid moving in Euclidean three-dimensional space. The latter is only observable by a certain kind of observer, a F(ull)-observer. Other observers, R(elativistic)-observers, perceive this superfluid as a non-removable background, which can be modeled as a four-dimensional pseudo-Riemannian manifold. What is the difference between these types of observers? F-observers can perform measurements using objects of arbitrary momenta and “see” the fundamental superfluid wavefunction’s in three-dimensional Euclidean space according to Equation (2) or an analogue thereof. On the other hand, R-observers are restricted to measuring only small-momentum small-amplitude excitations of the background superfluid. This is somewhat Universe 2020, 6, 180 4 of 25 analogous to listening to acoustic waves () in the conventional Bose-Einstein condensates, but being unaware of higher-energy particles such as photons or . According to BEC-spacetime correspondence, a R-observer “sees” himself located inside four-dimensional curved spacetime with a pseudo-Riemannian metric. The latter can be written in Cartesian coordinates as [7]:

 .   2 2 2 . − cs − η (∇S) . −η∇ S ρ   gµν ∝  ··········  , (5) cs  .  −η∇S . I where η = h¯ /m, S = S(r, t) = −i ln (Ψ(r, t)/|Ψ(r, t)|) is a phase of the condensate wavefunction √ written in the Madelung representation, Ψ = ρ exp (iS), and I is a three-dimensional unit matrix. To maintain the correct metric signature in Equation (5), condition |cs| > η |∇S| must be imposed, which indicates that cs is the maximum attainable velocity of test particles (i.e., small-amplitude excitations of the condensate), moving along geodesics on this induced spacetime. Therefore, cs is the velocity of those excitations of vacuum, which describe massless particles in the low-momentum limit, whereas massive test particles move along geodesics of a pseudo-Riemannian manifold with metric Equation (5). According to a R-observer, they are freely falling, independently of their properties including their rest mass. In this approach, we interpret Einstein field equations not as differential equations for an unknown metric; but as a definition for an induced stress-energy tensor, describing some effective matter to which test particles couple. Therefore, this would be the gravitating matter observed by a R-observer. We thus obtain   −1 1 Teµν ≡ κ Rµν(g) − gµνR(g) , (6) 2 = 2 where κ 8πG/c(0) is the Einstein’s . An example of usage of this procedure will be considered in Section 7.1. While Equation (6) is in fact an assumption, it should hold not only under the validity of conventional general relativity, but also in other Lorentz-symmetric theories of gravity which are linear with respect to the Riemann tensor, because the form of Einstein equations is quite universal (up to a conformal transform). For other Lorentz-symmetric theories, whose field equations cannot be transformed into this form, definition (6) can be adjusted accordingly. Furthermore, one can see from Equation (4), that cs contains an unknown function F(ρ). To determine its form, let us recall that one of the relativistic postulates implies that velocity cs should not depend on density, at least in the classical limit. More specifically, at low momenta, 10 −1 this velocity should tend to the value c(0) ≈ c, where c = 2.9979 × 10 cm s is called the in vacuum, for historical reasons. Recalling Equation (4), this requirement can be written as a differential equation [7]: 0 2 ρ|F (ρ)| = mcs ≈ const(ρ), (7) where const(ρ) denotes a function which does not depend on density. The solution of this differential equation is a logarithmic function: F(ρ) = −b ln (ρ/ρ¯), (8) where b and ρ¯ are generally real-valued functions of coordinates. The wave Equation (2) thus narrows down to " # h¯ 2 ih¯ ∂ Ψ = − ∇2 + V (r, t) − b ln (|Ψ|2/ρ¯) Ψ, (9) t 2m ext where b is the nonlinear coupling; b = b(r, t) in general. Correspondingly, Equation (4) yields q p = −(b/m)ρ, cs = |b|/m, (10) Universe 2020, 6, 180 5 of 25 thus indicating that logarithmic Bose liquid behaves like barotropic perfect fluid; but only when one neglects quantum corrections, and assumes classical averaging. This reaffirms the statement made in the previous Section about the place of perfect-fluid models when it comes to gravitational phenomena. The way gravity emerges in the superfluid vacuum picture is entirely different from those models, as will be demonstrated shortly, after we have specified our working model. Some special cases of Equation (9), for example when b → b0 = const, were extensively studied in the past, although not for reasons related to quantum liquids [22,23]. There were also extensive mathematical studies of these equations, to mention just some very recent literature [24–37]. Interestingly, wave equations with logarithmic nonlinearity can be also introduced into fundamental physics independently of relativistic arguments [7,38,39]. This nonlinearity readily occurs in the theory of open quantum systems, quantum and information [40,41]; as well as in the theory of general condensate-like materials, for which characteristic kinetic are significantly smaller than interparticle potentials [42]. One example of such a material is helium II, the superfluid phase of helium-4. For the latter, the logarithmic superfluid model is known to have been well verified by experimental data [10,43]. Among other things, the logarithmic superfluid model does reproduce the sought-after Landau-type spectrum of excitations, discussed in the previous Section; detailed derivations can be found in [10]. One of underlying reasons for such phenomenological success is that the ground-state wavefunction of free (trapless) logarithmic liquid is not a de Broglie plane wave, but a spatial Gaussian modulated by a de Broglie plane wave. This explains the liquid’s inhomogenization followed by the formation of fluid elements or parcels; which indicates that such models do describe fluids, rather than gaseous matter [44–49]. To summarize, a large number of arguments to date, both theoretical and experimental, demonstrate the robustness of logarithmic models in the general theory of superfluidity. In the next Section we shall demonstrate the logarithmic superfluid model’s capabilities when assuming superfluidity of the physical vacuum itself. In what follows, we shall make use of a minimal inhomogeneous model for the logarithmic superfluid which was proposed in [42], based on statistical and thermodynamics arguments. In the F-observer’s picture, its wave equation can be written as " # h¯ 2  q   |Ψ|2  ih¯ ∂ Ψ = − ∇2 + V (r, t) − b − ln Ψ, (11) t 2m ext 0 r2 ρ¯ √ where r = |r| = r · r is a radius-vector’s absolute value, and b0 and q are real-valued constants. For definiteness, let us assume that b0 > 0, because one can always change the overall signs of the nonlinear term F(ρ) and the corresponding field-theoretical potential V(ρ). As always, this wave equation must be supplemented with a normalization condition (1), boundary and initial conditions of a quantum-mechanical type; which ensure the fluid interpretation of Ψ [50]. 2 One can show that nonlinear coupling b = b(r) = b0 − q/r is a linear function of the quantum temperature TΨ, which is defined as a thermodynamic conjugate of entropy sometimes dubbed as the Everett-Hirschman information entropy. The latter can be written as 2 R 2 2 SΨ = −hΨ| ln (|Ψ| /ρ¯)|Ψi = − V |Ψ| ln (|Ψ| /ρ¯) dV, where a factor 1/ρ¯ is introduced for the sake of correct dimensionality, and can be absorbed into an additive constant due to the normalization condition (1). Therefore, one can expect that the thermodynamical parameters

b0 = b0(TΨ), q = q(TΨ) (12) are constant at a fixed temperature TΨ. Thus, for a trapless version of the model (11) we have four parameters, but only two of them, m and ρ¯, are a priori fixed, whereas the other two, b0 and q, can vary depending on the environment. Universe 2020, 6, 180 6 of 25

3. Induced Gravitational Potential Invoking model (11), while neglecting quantum fluctuations, let us assume that physical vacuum is a collective described by wavefunction Ψ = Ψvac(r, t), which forms a self-gravitating configuration with a center at r = 0. Therefore, for this state, the solution of Equation (11) is equivalent to the solution of the linear Schrödinger equation, " # h¯ 2 ih¯ ∂ Ψ = − ∇2 + V (r, t) Ψ, (13) t 2m eff for a of mass m driven by an effective potential

 |Ψ (r, t)|2  V (r, t) = V (r, t) − b ln vac eff ext ρ¯ (14)  q   |Ψ (r, t)|2  = V (r, t) − b − ln vac , ext 0 r2 ρ¯ when written in Cartesian coordinates [42]. If working in curvilinear coordinates, the last formula must be supplemented with terms which arise after separating out the angular variables in the wave equation. In the absence of quantum excitations and other interactions, it is natural to associate this effective quantum-mechanical potential with the only non-removable that we know of: gravity. This interpretation will be further justified in Section4. Therefore, in Cartesian coordinates one can write the induced gravitational potential as

1 1  q   |Ψ (r, t)|2  Φ(r, t) = − V (r, t) = b − ln vac , (15) m eff m 0 r2 ρ¯ where we assume that the background superfluid is trapless, i.e., we set Vext = 0. It should also be remembered that in curvilinear coordinates, this formula must be modified according to the remark after Equation (14); but for now we shall disregard any and rotation. It should be noticed that if one regards this potential as a multiplication operator then its quantum-mechanical average would be related to the Everett-Hirschman information entropy 2 discussed in the previous Section: hΦi ∼ TΨ hΨ| ln (|Ψ| ) |Ψi ∼ TΨSΨ. This not only makes theories of (which are essentially based on the ideas of Bekenstein, Hawking, Jacobson and others) a subset of the logarithmic superfluid vacuum approach, but also endows them with an underlying physical meaning and origin of the entropy implied. We can see that the induced potential maintains its form as long as the physical vacuum stays in the state |Ψvaci. If the vacuum were to transition into a different state, then it would change its wavefunction; hence the induced gravitational potential would also change. We expect that our vacuum is currently in a stable state, which is close to a ground state or at least to a metastable state, with a sufficiently large lifetime. It is thus natural to assume that the state |Ψvaci is stationary and rotationally invariant. As we established earlier, the wavefunction describing such a state should be the solution of a quantum wave equation containing logarithmic nonlinearity. In the case of trivial spatial topology and infinite extent, the amplitude of such a solution is known to be the product of a Gaussian function, which was mentioned in the previous Section, and a conventional quantum-mechanical part, which is a product of an exponential function, power function and a polynomial. Thus we can write the amplitude’s general form as:

χ /2 p  r  0  a2 2 a1 a0  |Ψvac| = ρ¯ P(r) exp − r + r + , (16) `¯ 2`¯2 2`¯ 2 Universe 2020, 6, 180 7 of 25

1/3 where P(r) is a polynomial function, χ0 and a’s are constants, and `¯ = (m/ρ¯) is a classical `¯ characteristic length scale of the logarithmic√ nonlinearity (alternatively, one can choose being equal to the quantum characteristic length, h¯ / mb0, which might be more useful for h¯ -expansion techniques). If quantum liquid occupies an infinite spatial domain then the normalization condition (1) requires

a2 > 0, (17) which is also confirmed by analytical and numerical studies of differential equations with logarithmic nonlinearity of various types [23–25,28,31,35,42]. Both the form of a function P(r) and the values of χ0 and a’s must be determined by a solution of an eigenvalue problem for the wave equation under normalization and boundary conditions. At this stage, those conditions are not yet precisely known; even if they were, we do not yet know which quantum state our vacuum is currently in. Therefore, these constants’ values remain theoretically unknown at this stage, yet can be determined empirically. Furthermore, for the sake of simplicity, let us approximate the power-polynomial term (r/`¯)χ0/2P(r), by the single power function (r/`¯)χ/2, where the constant χ is the best fitting parameter. Therefore, we can approximately rewrite Equation (16) as h a a  r  i |Ψ |2 ≈ ρ¯ exp − 2 r2 + 1 r + χ ln + a , (18) vac `¯2 `¯ `¯ 0 which is more convenient for further analytical studies than the original expression (16). From the empirical point of view, the function (18) can be considered as a trial function, whose parameters can be fixed using experimental data following the procedure we describe below. For the trial solution (18), the normalization condition (1) immediately imposes a constraint for one of its parameters:

(χ+3)/2     −1 Ma2 3 1 a1 3 exp(a0) ≈ Z˘ , + √ Z˘ 2, , (19) 2πm 2 2 a2 2

˘ 2  where we introduced an auxiliary function Z (a, b) = Γ (a + χ/2) 1F1 a + χ/2, b; a1/4a2 , where Γ(a) and 1F1(a, b; z) are the gamma function and Kummer confluent hypergeometric function, respectively. If values of a’s and χ are determined, e.g., empirically, then this formula can be used to estimate the ratio M/m. Furthermore, by substituting the trial solution (18) into the definition (15), we derive the induced gravitational potential as a sum of seven terms:

Φ(r) = Φsmi(r) + ΦRN(r) + ΦN(r) (20) +Φgal(r) + Φmgl(r) + ΦdS(r) + Φ0, Universe 2020, 6, 180 8 of 25 where

χ q ln (r/`¯) L2 ln (r/`¯) Φ (r) = − = −ζ c2 smi , (21) smi m r2 χq b r2 a q 1 L2 Φ (r) = − 0 = −ζ c2 RN , (22) RN m r2 a0q b r2 a q 1 GM Φ (r) = − 1 = − , (23) N m`¯ r r χ b Φ (r) = 0 ln (r/`¯) = c2 χ ln (r/`¯), (24) gal m b a1b0 2 r Φmgl(r) = r = ζa1 cb , (25) m`¯ Lmgl a b r2 Φ (r) = − 2 0 r2 = −c2 , (26) dS `¯2 b 2 m LdS and ! 1  a q  1 q Φ = a b + 2 = a b + (27) 0 m 0 0 `¯2 m 0 0 2 LdS is the additive constant. Here, and throughout the paper, we denote the sign functions by ζ’s: ζα = sign (α), and use the following notations:

r s b0 a1q |χ q| cb = , GM = , Lsmi = , m m`¯ b0 s |a0q| `¯ |q| LRN = , Lmgl = = , (28) b0 |a1| mGM `¯ χ`¯ χq LdS = √ , Lχ = = , a2 a1 mGM where G is the Newton’s gravitational constant as per usual. Furthermore, Lorentz symmetry emerges in the “phononic” low-momentum limit of the theory, as discussed in the previous Section. Therefore, a R-observer would perceive the gravity induced by potential (20) as curved four-dimensional spacetime, which is a local perturbation (not necessarily small) of the background flow metric, such as the one derived in Section 5.3 of [7], see Section 7.1 below. In a rotationally invariant case, the line element of this spacetime can be written in the Newtonian ( ) 2  gauge; if Φ r /c(0) 1, then it can be approximately rewritten in the form " # 2Φ(r) dr2 ds2 ≈ −c2 1 + dt2 + + R2(r)dσ2, (29) (0) 2 + ( ) 2 c(0) 1 2Φ r /c(0) h  i ( ) = + O ( ) 2 ≈ 2 = 2 + 2 2 where R r r 1 Φ r /c(0) r, d σ dθ sin θ dϕ is the line element of a unit two-sphere, and a leading-order approximation with respect to the is implied, as usual. The mapping (29) is valid for regions where the induced metric maintains a signature ‘− + ++’, and its matrix is non-singular. In other regions, such as close vicinities of spacetime singularities or horizons, the relativistic approximation is likely to fall outside its applicability range, thus it should be replaced with the F-observer’s description of reality. The main simplifying assumptions and approximations underlying the derivation of our gravitational potential are summarized and enumerated in the AppendixB. Universe 2020, 6, 180 9 of 25

4. Physical Interpretation It should be noticed that if we did not have a logarithm in the original model (11), then in Equations (14) and (15), then we would not have arrived at the polynomial functions in Equations (20)–(26), which are easily recognizable. This reaffirms our expectations that the underlying model can be successfully confirmed by experiment; but first those functions must be endowed with precise physical meaning. In this Section, we shall assign a physical interpretation to each term of the derived gravitational potential. For the sake of brevity, we shall be omitting an additive constant Φ0, assuming that it is 2 small compared to c(0).

4.1. Potential ΦN and Gravitational We begin with term (23), which has the most obvious meaning. In a non-relativistic picture, it represents Newton’s model of gravity. According to the BEC-spacetime correspondence manifested through the mapping (29), a R-observer can observe an effect of the potential ΦN by measuring probe particles moving along geodesics in the Schwarzschild spacetime:

2 2 2  rH  2 dr 2 2 ds(N) ≈ −c(0) 1 − dt + + r dσ , (30) r 1 − rH/r

= 2 where rH 2GM/c(0) is the Schwarzschild radius. Therefore, in absence of asymptotically non-vanishing terms, M can be interpreted as the gravitational mass of the configuration. This mass can be expressed in terms of superfluid parameters as ( 2a q 1 gravity, r = 1 , sign (a q) = (31) H 2 `¯ 1 −1 anti-gravity, mc(0) thus assigning physical meaning to a combination of parameters a1q/m`¯. In particular, one can see that a sign of the product a1q determines whether the ΦN interaction is attractive (gravity) or repulsive (anti-gravity). For most systems that we know of, anti-gravitational effects have not yet been observed, therefore one can assume that M > 0 or a1q > 0 (32) from now on. Nevertheless, it should be emphasized that the anti-gravity case is not a priori forbidden in superfluid vacuum theory. Indeed, the spacetime singularity occurs at r = 0 in a relativistic picture only, which poses certain issues for a R-observer, especially in the case of anti-gravity when a singularity is not covered by an (the existence of naked singularities is often doubted, on grounds of the cosmic censorship hypothesis). However, a F-observer would see no singular behavior in either case, because the wavefunction Ψvac remains regular and normalizable at each point of space and at any given time—as it should be in a quantum-mechanical theory. This reaffirms the fact that spacetime singularities are an artifact of incomplete information accessible to observers operating with relativistic particles [7]. Thus, the mapping from Equations (23)–(30) can be used to reformulate phenomena in the language of continuum mechanics and the theory of superfluidity; which can resolve certain long-standing problems occurring in the relativistic theory of gravity. For instance, neglecting asymptotically non-flat terms for simplicity, one can view Equations (23), (30) and (31) as the gravitational mass generation mechanism: such mass is not a fundamental notion, but a composite quantum phenomenon induced by the background superfluid’s dynamics (through the elementary inertial mass m and critical density ρ¯), its quantum temperature (through q), and an exponential part of the condensate’s wavefunction (through a1). Such a mechanism can be thus considered as the Universe 2020, 6, 180 10 of 25

quantum-mechanical version of the Mach principle [7]. For example, if either a1 or q vanish, then the system would not possess any gravitational mass, but it still would be gravitating in a non-Newtonian way, if other potentials from Equation (20) are non-zero.

4.2. Potential ΦRN and Abelian Charges Equation (22) represents another potential which can be easily recognized. According to the mapping (29), potential ΦRN is observed by a R-observer as Reissner–Nordström spacetime, when taken together with the ΦN potential:

2 ! 2 r rQ dr ds2 ≈ −c2 1 − H + dt2 + + r2dσ2, (33) (N+RN) (0) 2 2 2 r r 1 − rH/r + rQ/r where rQ is a characteristic length scale: v √ c u 2|a q| r = 2L b = u 0 , (34) Q RN c t 2 (0) mc(0) provided a0q < 0. The Reissner–Nordström metric is known to be associated with the gravitational field caused by a charge related to an abelian group. An example would be an electric charge Q, which is related to the abelian group U(1) of . This charge can be thus revealed through the formula

2 4 2 rQc(0) 2|a q|c Q2 = ≈ 0 , (35) keG keGm where ke is the Coulomb constant. In other words, potential ΦRN describes, together with ΦN, the gravitational field created by an object of a charge Q and gravitational mass M. Thus, from a F-observer’s viewpoint, an electrical charge is not an elementary notion, but a composite quantum phenomenon, induced by the background superfluid’s dynamics (through the elementary inertial mass m), its quantum temperature (through q), and overall constant coefficient of the condensate’s wavefunction (through a0). It should be noted also that ΦRN is a short-range potential, therefore, it becomes substantial only 2 ¯ at those microscopical length scales, of an order rQH = rQ/rH = ` |a0/a1| or below. Since rQ < rH for most objects we know of, we have 0 6 rQH < rH. Thus, those scales would be causally inaccessible to a R-observer, but a F-observer would have no problem accessing them, per usual.

4.3. Potential Φsmi and Strong Gravity As the distance from a gravitating center decreases, it is term (21) which eventually predominates. According to the mapping (29), this potential Φsmi, when taken together with the ΦN and ΦRN potentials, induces spacetime with the line element:   r2 r2   2 ≈ − 2 − rH + Q − W r 2 ds(N+RN+smi) c(0) 1 r r2 ζχq r2 ln `¯ dt dr2 (36) + + r2dσ2, 2 2 2 ¯ 2 1 − rH/r + rQ/r − ζχqrW ln r/` /r where rW is a characteristic length scale: v √ c u 2|χ q| r = 2L b = u . (37) W smi c t 2 (0) mc(0) Universe 2020, 6, 180 11 of 25

The potential Φsmi has a distinctive property: unlike other sub-Newtonian potentials in Equation (20), it can switch between repulsive and attractive regimes, depending on whether the distance is larger or smaller than `¯. The magnitudes of Φsmi and ΦRN become comparable at two values of r, shown by the formula

(±) ¯ 2 2 ¯ rWQ = ` exp (±rQ/rW ) = ` exp (±|a0/χ|), (38)

(−) (+) which indicates that |Φsmi| overtakes |ΦRN| either at r < rWQ or at r > rWQ. If |a0/χ| is large then (+) (−) rWQ is exponentially large and rWQ is exponentially small. Magnitudes of Φsmi and ΦN become comparable at a certain value of r, which is shown by the formula

¯  rWM = −`W W −`/`W 6 `W, (39)

2 where `W ≡ rW /|rH| = |Lχ|, and by W(a) we denoted the Lambert function. The function |Φsmi(r)| becomes larger than |ΦN(r)| at r < rWM 6 |Lχ|. Furthermore, considering Φsmi as a perturbation of ΦN, one can deduce that the effective gravitational coupling,  ln (r/`¯)  G ≈ G 1 + ζ L , (40) eff χq χ r becomes larger as r gets smaller. Here an approximation symbol reminds us that we are working with wavefunction (18), and assume a leading-order approximation with respect to the Planck constant. One can see that gravity naturally becomes stronger at shorter scales, without introducing any additional effects or matter, which suggests the way towards resolving the . Indeed, the latter states the large discrepancy between magnitudes of forces of interactions and classical gravity, but in this case, the gravitational force grows stronger than inverse square as r decreases, as ln r/r3. Moreover, it is possible to have even stronger small-length behaviour in our approach: if one goes beyond a minimal model (11) and generalizes nonlinear coupling to a series: 2 3 4 5 b(r) = b0 − q/r + q1/r + ..., then this will induce terms ln r/r , ln r/r , et cetera. ¯ On the other hand, as r grows, the function Geff converges to G, especially if |Lχ| or 1/` are sufficiently small; which makes Geff approximately constant for a large range of values of r.

4.4. Potential ΦdS and de Sitter Spacetime Let us turn our attention to terms which do not tend to zero at r → ∞. The physical meaning of one of these terms, given by Equation (26), becomes clear upon using the mapping (29): ! r2 dr2 ds2 ≈ −c2 1 − dt2 + + r2dσ2, (41) (dS) (0) 2 2 2 RdS 1 − r /RdS where s c mc2 1 (0) ¯ (0) RdS = √ LdS = ` (42) 2 cb 2a2b0 is a radius of de Sitter horizon. This metric represents de Sitter spacetime (written in static coordinates), which belongs to a class of Friedmann–Lemaître–Robertson–Walker (FLRW) . Indeed, by applying a coordinate 2 2  −1 2 2 −1/2 transformation τ − t = τdS ln 1 − r /RdS , $ = α0 r exp(−t/2τdS) 1 − r /RdS , one can rewrite Equation (41) in isotropic coordinates   2 2 2 2 τ  2 2 2 ds(dS) → −c(0)dτ + α0 exp d$ + $ dσ , (43) τdS Universe 2020, 6, 180 12 of 25

where τdS = RdS/2c(0) and α0 is an integration constant. The physical implications of the term ΦdS will be further examined in Section7.

4.5. Potentials Φmgl and Φgal and Gravity on Astronomical Scales The remaining asymptotically non-vanishing potentials are given by Equations (24) and (25). With respect to the dependence upon radial distance from the gravitating center, they occupy an intermediate place between de Sitter term (26) and Newtonian potential (23). It is thus natural to expect that these terms are responsible for large scale dynamics—from galaxies (a kiloparsec scale) to metagalactic objects, such as voids and (a megaparsec scale). According to the mapping (29), the terms Φgal and Φmgl modify de Sitter metric (41): causing the resulting spacetime to be asymptotically de Sitter only. The line element, which corresponds to the terms ΦdS, Φmgl and Φgal taken together, can be written in static coordinates as     2 ds2 ≈ −c2 1 + β2 ln r + r − r dt2 (gal+mgl+dS) (0) χ `¯ RM 2 RdS dr2 (44) + + r2dσ2, 2 ¯ 2 2 1 + βχ ln r/` + r/RM − r /RdS q √ = 2 = where βχ 2χ b0/mc(0) 2χcb/c(0), and

2 ¯ 2 mc( )` ζ c( ) = 0 = a1 0 RM 2 Lmgl (45) 2a1b0 2 cb is a characteristic length scale constant, which can be positive or negative depending on the sign of a1. Note that in Equation (44), we included the ΦdS-induced (de Sitter) term to remind us that at large r we still have a spacetime of a FLRW type. Interestingly, linear terms in metrics occur also in an alternative theory of gravity, Weyl gravity [51,52]. This coincidence can be explained by conformal symmetry, which often emerges in logarithmic models at the relativistic limit, see for example Section 7.1. However, Weyl gravity does not produce a logarithmic term in metric, while in our theory it is induced by Φgal. This term is responsible for the flat rotation curves phenomenon in galaxies, which will be discussed, along with other astronomical implications of the terms Φmgl and Φgal, in further detail in Section6. It is also useful to know that the metric induced by the term (25) alone,

  2 2 2 r 2 dr 2 2 ds(mgl) ≈ −c(0) 1 + dt + + r dσ , (46) RM 1 + r/RM transforms to the conformally FLRW-type metric

2 " 2 # (1 + $/4RM) a (τ)   ds2 → −c2 dτ2 + M d$2 + $2dσ2 , (47) (mgl) 2 2 (0) 2 2 2 aM(τ)(1 − $/4RM) (1 − $ /16RM)

2 R upon applying the coordinate transformation r = $/(1 − $/4RM) , t = dτ/aM(τ). For R-observers, this represents a surrounding homogeneous and isotropic spacetime with aM(τ) and 2 negative-definite spatial scalar , −1/(2RM) . Finally, one could mention that potentials of type (24) were studied, albeit in the context of a linear Schrödinger equation, in [53,54]. Universe 2020, 6, 180 13 of 25

5. Density of Effective Gravitating Matter In this Section, let us derive spherically-symmetric density profiles of effective gravitating matter ρΦ, which formally corresponds to the superfluid-vacuum induced potential (20). Contrary to the directly observable value of orbital velocity derived in Section6, such density depends more substantially on the choice of an observer. It can be defined either via the Lorentz-covariant definition (6), or via the Poisson equation, which is a non-relativistic version of Equation (6). Correspondingly, one would obtain different results, which will be discussed below. In accordance with the last paragraph of AppendixB, this Section’s computations cannot take into account any secondary induced matter, such as the equilibrium configurations of mass-energy emerging as a result of interaction between scalar and tensor modes of superfluid vacuum’s excitations. One can show that such equilibria do exist, manifesting themselves in a form of general relativistic nonsingular horizon-free stellar-like objects or particle-like Q-balls, therefore, in reality such objects would definitely contribute to the density associated with dark matter. In other words, here we are restricting ourselves to background values of density. In this Section only, we shall be temporarily assuming that M = M(r), otherwise the corresponding contribution to the density profile would be identically zero, and thus non-indicative.

5.1. Galilean Symmetry In this case, one defines an effective matter density by virtue of the Poisson equation. Taking the whole potential (20) and assuming spherical symmetry, we obtain

1  2  ρ (r) ≡ Φ00(r) + Φ0(r) Φ 4πG r

= ρsmi(r) + ρRN(r) + ρN(r) (48)

+ρgal(r) + ρmgl(r) + ρdS(r), where

3χ q  2  r  ρ (r) = 1 − ln , (49) smi 4πGmr4 3 `¯ a q k Q2 ρ (r) = − 0 = e , (50) RN 2πGmr4 4πr4 1 ρ (r) = − M00(r), (51) N 4πr 2 χ b vgal ρ (r) = 0 = , (52) gal 4πGmr2 4πGr2 a1 b0 χ b0 ρmgl(r) = = , (53) 2πGm`¯r 2πGmLχr 3a b 3b ρ (r) = − 2 0 = − 0 = const, (54) dS `¯2 2 2πGm 2πGmLdS whereby the sum of the last three densities can be regarded as a density corresponding to the astronomical-scale “dark matter” and “dark energy”, which will be further justified in Sections6 and7 . Furthermore, these formulae were derived on the assumption that the gravitational G is the same for all length scales; which is valid when any influence from the term (21) can be disregarded. However, this term might cause an additional effect, discussed in Section 4.3: it makes the gravitational coupling constant vary with distance. If this does happen, then in Equations (48)–(54) one should replace G with a running constant given by Equation (40). Universe 2020, 6, 180 14 of 25

5.2. Lorentz Symmetry In this case, one defines effective matter density by virtue of Einstein field equations with the induced stress-energy tensor defined by Equation (6). Using it together with Equation (20) and metric (29), we obtain

1  1  ρ˜ (r) ≡ − Φ0(r) + Φ(r) Φ 4πGr r

= ρ˜smi(r) + ρ˜RN(r) + ρ˜N(r) (55)

+ρ˜gal(r) + ρ˜mgl(r) + ρ˜dS(r) + ρ˜Φ0 (r), where χ q h  r i ρ˜ (r) = 1 − ln , (56) smi 4πGmr4 `¯ a q k Q2 ρ˜ (r) = − 0 = e , (57) RN 4πGmr4 8πr4 1 ρ˜ (r) = M0(r), (58) N 4πr2 χ b h  r i ρ˜ (r) = − 0 1 + ln , (59) gal 4πGmr2 `¯ a1 b0 χ b0 ρ˜mgl(r) = − = − , (60) 2πGm`¯r 2πGmLχr 3a b 3b ρ˜ (r) = 2 0 = 0 = const, (61) dS `¯2 2 4πGm 4πGmLdS ! Φ 1 q ( ) = − 0 = − + ρ˜Φ0 r 2 2 a0b0 2 , (62) 4πGr 4πGmr LdS whereby the sum of last four densities can be regarded corresponding to the astronomical-scale “dark matter” and “dark energy”, which will be further justified in Sections6 and7. Density (62) is a somewhat surprising contribution, because it corresponds to the constant term Φ0 in Equation (20), which is not supposed to affect trajectories; at least, in classical mechanics. Unless its presence can be confirmed by observations, it must be regarded as a gauge term, or as an artifact of approximations underlying Equation (29). Furthermore, a running gravitational coupling constant can not be implemented in the relativistic case as simply as in Section 5.1. To preserve Lorentz invariance, one has to associate this coupling with a four-dimensional scalar; which automatically upgrades general relativity to a scalar-tensor gravity with a non-minimally coupled scalar field. This theory cannot be written by hand, but it must be derived in a way which is similar to that used in Section 7.1. Comparing results of Sections 5.1 and 5.2, one can conclude that the definition of effective matter density is somewhat ambiguous: in particular, it drastically depends on symmetry assumptions. Therefore, further experimental studies should help to empirically establish which symmetry is more appropriate to use when dealing with “dark” phenomena.

6. Galactic Rotation Curves In this Section, we demonstrate how induced gravitational potential can explain various phenomena, which are usually attributed to dark matter. Let us focus on the terms (24) and (25), which were partially discussed in Section 4.5. Because they become significant at a galactic scale and above (i.e., a kiloparsec to megaparsec scale), it is natural to conform them to astronomical observations; such as those of rotation curves in galaxies. In a spherically symmetric case, velocity curves of orbiting with non-relativistic on a plane in a central gravitational potential Φ(r) can be estimated using a simple formula Universe 2020, 6, 180 15 of 25

2 0 v = Rac = R Φ (R), where v is the orbital velocity, ac is the centripetal acceleration, and R is the orbit’s radius. The cylindrically symmetric case can be considered by analogy, by assuming various disk models [55,56]. Considering the terms (24) and (25) in conjunction with the Newtonian term (23), we thus obtain

1/2 n d h io v(R) = R dR ΦN(R) + Φgal(R) + Φmgl(R) (63) q 2 2 = vN + vgal + Φmgl(R), where GM a q 1 v2 = = 1 = −Φ (R), (64) N R m`¯ R N χb v2 = χc2 = 0 = const, (65) gal b m while the contribution from the term (26) is disregarded for now, due the assumed smallness of the 2 “local” 1/RdS; and the contribution from the term (21) is disregarded due the assumed smallness of the corresponding characteristic length, according to discussion in Section 4.3. In the case of a galaxy, the contribution from the Newtonian term ΦN rapidly decreases as R grows. Correspondingly, the main contribution would then come from the second term in a row, Φgal, and then from the third term, Φmgl. From Equation (65) one can see that the contribution from Φgal is constant, which explains the average flatness of galactic rotation curves. Notice that the value of velocity vgal depends on one of the wavefunction parameters χ and one of quantum temperature parameters b0. Both are not a priori fixed parameters of the model, cf. Equations (12) and (18), but vary depending on the environment and conditions: background superfluid gets affected by the gravitational potential it induces, because this potential acts upon the surrounding conventional matter, thus creating density inhomogeneity. Therefore, χ and b0 should generally be different for each galaxy; and hence should vgal be. Similar to the case of vgal, the parameter a1 hence a value Φmgl(R) at a fixed R will also be dependent on the gravitating object they refer to. This potential should usually be negligible on the inner scale length of a galaxy, but as R grows towards the extragalactic length scale, R & 10 kpc, rotation curves should start to deviate from flat: s R v(R) ≈ vgal 1 + , (66) Lχ which can be used for estimating the combination of superfluid vacuum parameters χ`¯/a1 empirically. In cases where the contribution from other terms of the induced potential cannot be neglected, Equation (63) must be generalized to include those too. Possible galactic-scale regions, where this non-flat asymptotics should become visible, depending on a value Lχ, are the outer regions of large spiral galaxies, such as M31 or M33 [57–61], where Φmgl(r) can not only overtake ΦN(r) but also become comparable with Φgal(r).

7. Accelerating Expansion of the Universe This phenomenon is usually explained by introducing exotic forms of relativistic matter, such as dark energy; usually modeled by various long-range scalar fields, which are assumed not to affect the numerous experiments on the Earth. The superfluid vacuum approach offers a simple framework, which can explain the Universe’s expansion as an observer-dependent effect, without involving any matter other than the background superfluid itself. Universe 2020, 6, 180 16 of 25

7.1. Conformally Flat Spacetime and Field Following work [7], let us consider the most simple possible special case: laminar flow of a logarithmic background condensate in a state |Ψ(0)i, described by Equation (9) at b = const, with a constant velocity u(0), if viewed as, from the F-observer’s perspective, an embedding into underlying Euclidean space. On the other hand, what does a R-observer see? Due to a well-known separability property of the logarithmic Schrödinger equation [23,25,28,31,42], the phase of its simplest ground-state solutions is a linear function of a radius-vector:   (0) i ln Ψ(0)(r, t)/|Ψ(0)(r, t)| ∝ u · r + f (t), (67) where u(0) is a constant 3-vector, and f (t) is an arbitrary function of time. In this case, the fluid-Schrödinger analogy confirms that the background condensate does flow with a constant   (0) velocity u = −iη∇ ln Ψ(0)/|Ψ(0)| ∝ u . Recalling Equation (5), it means that the background geometry induced by such solutions is conformally flat:

1   2 ds2 ∝ |Ψ (r, t)|2 −c2dt2 + dr − u(0)dt , (68) (0) ρ¯ (0) s where cs is given by Equation (10), in a leading-order approximation with respect to the Planck constant. Spacetime of a type (68) lies within a large class of manifolds with the vanishing Weyl tensor—a type O in the Petrov classification. This is the class that all FLRW spacetimes belong to, including those which describe the Universe with accelerating expansion—simply written in conformally-flat coordinates instead of comoving ones. Using definition (6), we obtain an induced stress-energy tensor for our system: h  i ˜ λ 1 ˜ λ κTeµν = D ∇µ∇νφ − ∇µφ∇νφ − gµν ∇λ∇ φ + 2 (D − 1)∇λφ∇ φ , (69) where D˜ = D − 2 = 2, ∇ is a covariant derivative with respect to metric g, and by φ we denote the induced scalar field:  2  φ = ln |Ψ(0)(r, t)| /ρ¯ , (70) up to an additive constant. This stress-energy tensor strongly resembles the one occurring in the theory of gravity with a scalar field. One can verify that it can be indeed derived from the following scalar-tensor gravity action functional Z p ˜ h i S˜[g, 6 φ] ∝ dDx −g eDφ R + D˜ (D˜ + 1)(∇φ)2 , (71) where the notation “6 φ” reminds us that the field φ is fixed by the solution of the original quantum wave equation, cf. Equation (70), while the variation of action must be taken with respect to the metric only. In other words, both metric g and dilaton φ are induced by the superfluid vacuum being in a state described by Ψ(0)(r, t). Thus, we have found yet another example of the differences between the F-observer’s and R-observer’s pictures of reality. While the former sees a background quantum fluid flowing with a constant velocity in three-dimensional Euclidean space, the latter observes itself as being inside four-dimensional spacetime governed by a Lorentz-covariant scalar-tensor gravity. An action functional (71) therefore explains why covariant models involving scalars provide a robust description of the large-scale evolution of the Universe, agreeing with current observational data; yet no quanta of relativistic dilaton have thus far been detected. This correspondence also reveals the limitations of the relativistic description itself: if the superfluid vacuum goes into a different quantum state, then one gets a different expression for the induced metric, scalar, stress-energy tensor and covariant action. In fact, for more complicated Universe 2020, 6, 180 17 of 25 superfluid flows, even the condition (67), leading to a conformal flatness, can become relaxed to an asymptotic one. Therefore, depending on the physical configuration (determined by external potential, if any, and boundary conditions), nonlinear coupling behaviour and the quantum state the vacuum is in, small fluctuations and probe particles would obey different covariant actions. Consequently, a R-observer would have to tweak its field-theoretical models by hand; with the unified picture being observable only at the level of a F-observer.

7.2. Cosmological Constant and Local Expansion Mechanism In Section 7.1 we considered an example of the global superfluid flow which would be seen by a relativistic observer as the accelerating expansion of the . What about locally (20), can it cause similar effects? Let us consider once again the induced potentials of Sections 4.4 and 4.5. The de Sitter term from Section 4.4 predominates if we consider the physical vacuum alone, without any generated matter therein. This is perhaps only valid for the early Universe, such as the one which existed during the inflationary epoch. In the current epoch, our spacetime can only be de Sitter asymptotically or even approximately; therefore the potential ΦdS must be negligible, unless one considers very large length scales. For example, if a length scale RdS, which generally depends on a massive object defining the of the function (16), is comparable to a size of the observable Universe (∼10 Gpc), then RdS can be related to the cosmological constant Λ as √ (cos) −56 −2 RdS = 3/Λ, Λ ∼ 10 cm , (72) therefore, the term (26) becomes substantial only at a gigaparsec scale. This relation yields an empirical constraint for a combination of characteristic parameters of superfluid vacuum, average quantum temperature and Ψvac: (cos) (cos) a b 1 2 0 ≈ Λc2 ∼ 10−36 s−2, (73) m`¯2 6 see also the discussion around Equation (12). All this essentially means that Λ is not a fundamental constant of , but a combination of various parameters of superfluid vacuum, including quantum temperature and Gaussian width of the condensate function Ψvac. Therefore, its smallness, sometimes referred as “vacuum catastrophe” [1], can easily be explained by the fact that: either average quantum temperature across the Universe, or wavefunction’s width, or both, are sufficiently small compared to mc2 and `¯, respectively; thus resulting in the overall smallness of the ratio on the left-hand side of Equation (73). The line element, which results from taking ΦdS, Φmgl, Φgal, Φ0 and ΦN together, can be written in static coordinates as     2 ds2 ≈ −c2 1 + δ − rH + β2 ln r + r − r dt2 (cos+N) (0) 0 r χ `¯ RM 2 RdS (74) dr2 + + r2dσ2, 2 ¯ 2 2 1 + δ0 − rH/r + βχ ln r/` + r/RM − r /RdS

= 2 = where δ0 2Φ0/c(0). In this metric, the Schwarzschild term ensures that spacetime singularity at r 0 is “dressed” by the black hole horizon, while at large r we still have spacetime of a FLRW type. While from the viewpoint of a F-observer, no spacetime singularity would pose a problem, because quantum wavefunction remains regular and normalizable at any non-negative value of r. This reaffirms our earlier statement that cosmological singularity is an artifact of the low-momentum approximation of superfluid vacuum [7,62]. According to Equation (47), the linear potential term (25), when taken alone, induces the universal 2 acceleration c(0)/2RM, occurring due to the spatial curvature, when seen by a R-observer in its own local static coordinate system. This contributes to the Hubble expansion induced by the quadratic Universe 2020, 6, 180 18 of 25 term (26). Thus, in the relativistic picture, the non-small fluctuation of superfluid vacuum produces an effect at the center of a gravitating configuration; and therefore contributes to the explicit rotational of the stars inside this configuration, which can be seen as a consequence of curved spacetime.

7.3. Expansion Mechanisms and Cosmological Coincidences Comparing Sections 7.1 and 7.2, one can see that they describe different expansion mechanisms. The mechanism of Section 7.1 occurs due to the global laminar flow of background superfluid, assumed to be laminar, which is “seen” by a R-observer as a FLRW-type spacetime; the resolution of various cosmological problems related to this mechanism was discussed in Section 5 of Ref. [7]. On the other hand, in Section 7.2, expansion is explained as a cumulative effect from terms in metric (74), which do not vanish at spatial infinity, induced by a “local” wavefunction associated with a gravitating configuration or body. This wavefunction can be regarded as a fluctuation (not necessarily small) of the global wavefunction from Section 7.1. The interplay between these mechanisms depends on the length scales of the quadratic and linear terms, RdS and RM. Unless a cumulative expansion effect from asymptotically non-vanishing terms taken together is, by some extraordinary coincidence, exactly equal to the expansion due to the global flow mechanism, its rate must be different from that of the global flow-induced expansion. The occurrence of an additional expansion mechanism, at the scale of a , such as our Virgo or Laniakea, could explain the remarkable discrepancy between measurements of the Hubble constant using different methods from those based on the whole Universe expansion, such as cosmic microwave background radiation (CMB). Among non-CMB methods one could mention Cepheid calibration, time-delay cosmography, and geometric distance measurements to megamaser-hosting galaxies [63–65]. From a theoretical point of view, a scenario with different expansion rates seems slightly more plausible, because it does not require an explanation why accelerations from different mechanisms should be exactly equal to each other (this coincidence should not be confused with the conventional cosmological coincidence which we will discuss next). Furthermore, our approach offers a simple explanation of the cosmological coincidence problem itself [1], in both the simplified and quantitative versions thereof. The simplified formulation of the cosmological coincidence states that if dark matter and dark energy were different kinds of matter, then during the Universe’s evolution they should have evolved independently of each other; therefore their distributions would be uncorrelated by now—which does not seem to be the case. Superfluid vacuum theory trivially resolves this paradox: because “dark matter” and “dark energy” are actually induced phenomena and manifestations of the same object, superfluid vacuum, they cannot be independent from each other. The quantitative formulation of the coincidence problem is an explanation requirement for why the ratio of DM- and DE-associated densities is of order one, despite the reasons given in the simplified formulation. In our approach, we can regard all terms of the induced potential, which do not vanish at spatial infinity, as being associated with “dark” effects; but inside this group we cannot unambiguously separate DM-attributed effects from DE ones. For example, potential (25) is intermediate between de Sitter and logarithmic, thus it affects both galactic rotation curves and Hubble expansion, cf. Sections6 and 7.2. If, for simplicity, we consider only the local expansion mechanism of Section 7.2 and omit the contribution from Φ0, then the cosmological coincidence can be reformulated as the following condition: Ωe gal + Ωe mgl = O(1), (75) Ωe mgl + Ωe dS where Ωe ’s are average values corresponding to densities from Section 5.2. The numerator of this ratio represents average density of effective “dark matter”, while the denominator represents effective “dark energy” density; or at least the predominating proportions thereof. Universe 2020, 6, 180 19 of 25

In general, this condition simply imposes yet another constraint for the parameters of the theory. It is trivially satisfied if the involved parts of the gravitational potential, hence the associated densities are of the same order of magnitude, if averaged on a large scale. Moreover, even if Ωe gal or Ωe dS are much smaller than the remaining involved densities, but the value Ωe mgl is substantial, then the relation (75) still holds, due to the presence of Ωe mgl in both parts of the ratio.

8. Conclusions Working within the framework of the post-relativistic theory of physical vacuum, based on the logarithmic superfluid model, we derived induced gravitational potential, corresponding to a generic quantum wavefunction of the vacuum. This mechanism is radically different from the one used in models of relativistic classical fluids and fields, which are based on modifying the stress-energy tensor in Einstein field equations. The form of such a wavefunction is motivated by ground-state solutions of quantum wave equations of a logarithmic type. Such equations find fruitful applications in the theory of strongly-interacting quantum fluids, and have been successfully applied to laboratory superfluids [10,43,49]. We note that, in principle, one is not precluded from adding other types of nonlinearity, such as polynomial ones, into the condensate wave equations, but the role of logarithmic nonlinearity is crucial. Thus, we used a logarithmic superfluid model with variable nonlinear coupling, because it accounts for an effect of the environment in a more realistic way than the logarithmic model with a constant coupling. As a result, for the trapless version of our model, we have four parameters, but only two of them are a priori fixed, whereas the other two can vary, depending on the quantum thermodynamic properties of the environment under consideration. Additionally, a number of parameters come from the wavefunction solution itself. Those are not independent parameters of the theory, but functions thereof. Because we do not yet know the exact form of the superfluid wavefunction, see remarks at the end of Section3, we leave those parameters to be empirically estimated, or bound, at the stage of current knowledge. It turns out that gravitational interaction has a multiple-scale structure in our theory: induced potential is dominated by different terms at each length scale; such that one can distinguish sub-Newtonian, Newtonian (inverse-law), galactic (logarithmic-law), metagalactic (linear-law), and cosmological (square-law) parts. A relativistic observer, who operates with low-momentum small-amplitude fluctuations of superfluid vacuum, observes this induced potential by measuring the trajectories of probe particles moving along geodesics in induced four-dimensional pseudo-Riemannian spacetime. The metric of the latter is determined by virtue of the BEC-spacetime correspondence and fluid-Schrödinger analogy, applied jointly. The sub-Newtonian part of the induced gravitational potential is defined as one which grows faster than the inverse law, as distance tends towards zero. It can be naturally divided into the following two parts. One part has an inverse square law behaviour, and thus can be associated with the gravitational field caused by a U(1) gauge charge, such as an electric charge. On a relativistic level, it is described by Reissner–Nordström spacetime. The other part has ‘inverse square times logarithm’ law behaviour, which might become substantial at both ultra-short and macroscopic distances, depending on the values of the corresponding parameters. If it “survives” at macroscopic distances, then it upgrades Newton’s gravitational constant to a function of length, such that gravity has both strong and weak regimes. With the potential or spacetime metric in hand, one can, in principle, assign effective fictitious matter density to our potential, which corresponds to “dark matter” and “dark energy”. This can be done in two ways: either by Einstein field equations in a relativistic case, or the Poisson equation in a non-relativistic one. It should be noted that the resulting density in each of the cases can be modified, depending on whether the gravitational constant is considered to be running or not. This will require more verification from future experimental and theoretical studies. Universe 2020, 6, 180 20 of 25

Furthermore, on a galactic scale and above, the potential is dominated by non-Newtonian terms, which do not vanish at spatial infinity. This explains the non-Keplerian behaviour of rotation curves in galaxies, which is often attributed to dark matter. Our model, not only explains the average flatness of galactic rotation curves, but also makes a number of new predictions. One of them is the approximately linear law behavior of gravitational potential on a metagalactic scale, which is an intermediate scale between galactic distances and the size of the observable universe. This should partially affect galactic rotation curves too: as the distance from the gravitating center grows further towards the metagalactic length scale, a squared velocity’s profile asymptotically changes from being flat towards linear, cf. Equation (66). On the other hand, at the largest length scale, the induced potential displays square law behaviour. If the quadratic term is negative-definite, then the corresponding metric describes (asymptotically) de Sitter space, merely written in static coordinates. Taken together with the contribution from the linear potential term, this explains the accelerating expansion of the corresponding spacetime region, which is usually associated with dark energy. Such expansion could supplement the “global” one, caused by laminar flow of background logarithmic superfluid absent any other matter, which induces a FLRW-type spacetime. The occurrence of more than one type of expansion mechanism, could be responsible for the discrepancy between measurements of the Hubble constant using different methods. The relevant problems, such as smallness of cosmological constants and cosmological coincidence, were also discussed. To conclude, we used the BEC-spacetime correspondence and fluid-Schrödinger analogy to argue that the description of reality and fundamental symmetry crucially depend on the choice of an observer. We demonstrated that both dark matter and dark energy are related phenomena, and different manifestations of the same object, superfluid vacuum, which acts by inducing both gravitational potential and spacetime.

Funding: This research was funded by Department of Higher Education and Training of South Africa and in part by National Research Foundation of South Africa. Acknowledgments: Proofreading of the manuscript by P. Stannard is greatly appreciated. Conflicts of Interest: The author declares no conflict of interest.

Abbreviations The following abbreviations are used in this manuscript:

BEC Bose-Einstein condensate CDM CMB Cosmic microwave background DE Dark energy DM Dark matter dS de Sitter FLRW Friedmann–Lemaître–Robertson–Walker F-observer Full observer RN Reissner–Nordström R-observer Relativistic observer SVT Superfluid vacuum theory

Appendix A. Superfluid or Not When it comes to astrophysics and cosmology, including dark matter and dark energy applications, the notion of superfluid is probably the most misunderstood and misused term of all. In addition to having a catchy, but not exactly informative, prefix, the first thing which causes misunderstanding of the term is inviscid flow. This, the most striking feature of superfluids, is often regarded as the Universe 2020, 6, 180 21 of 25 definition thereof. As a result, superfluids are often confused with perfect fluids, which have no viscosity by definition, and can be easily implemented into Einstein field equations by virtue of a stress-energy tensor borrowed from classical fluid mechanics. However, perfect fluids, being classical hydrodynamical objects by construction, can not properly reflect the essentially quantum nature of the superfluidic matter that we know of. Therefore, perfect or non-quantum inviscid fluids can be used as a crude approximation, at best, of superfluids. Additionally, superfluids are often also confused with the concomitant phenomenon of Bose-Einstein (BEC), which is another kind of quantum matter occurring in low-temperature condensed matter, such as cold cesium in a trap. However, in the superfluid helium phase, for example, the BEC’s content comprises only about ten per cent. Therefore, this condensate alone cannot fully account for -free flow and the other distinctive features of superfluids. Even though Bose-Einstein condensates do share certain features with superfluids, this does not imply that they are superfluidic in general. In particular, quantum excitations in laboratory superfluids are known to have dispersion relations of a distinctive shape called the Landau “roton” spectrum. This shape of the spectral curve is crucial, as it ensures the suppression of dissipative fluctuations at a quantum level [10]. If plotted as an excitation energy versus momentum, the curve starts from the origin, climbs up to a local maximum (called the maxon peak), then slightly descends to a local nontrivial minimum (called the “roton” energy gap); and then grows again, this time all the way up, to the boundary of the theory’s applicability range, cf. a curve in the Figure 1a from [11]. In fact it is not the roton energy gap alone, but the energy barrier formed by the maxon peak and roton minimum in momentum space, which ensures the above-mentioned suppression of quantum fluctuations in quantum liquid and, ultimately, causes its flow to become inviscid. In other words, it is the global characteristics of the curve, not just the existence of a nontrivial local minimum and related energy gap, which are required for superfluidity to occur. Obviously, Landau’s shape is a non-trivial property which cannot possibly occur in all quantum liquids. The final reason for the misuse of the term ‘superfluid’ is the extensive, but not always careful, utilization of relativistic scalar field models in astrophysics and cosmology. Historically, the four-dimensional scalar field came about as a bold extrapolation of a non-relativistic wavefunction into the realm of Lorentz-symmetric theories. However, when it comes to quantum liquids and condensates, there is a fundamental difference between relativistic scalar field and condensate wavefunction, which makes a correspondence between them far from isomorphic. The fluid condensate wavefunction obeys both a normalization condition (1) and a wave Equation (2), therefore, it is a three-dimensional Euclidean scalar related to a ray in the associated Hilbert space. One can see that Equations (1)–(3) are essentially non-relativistic and three-dimensional. While one can still make the last two relativistic, by replacing derivative parts with the Lorentz-covariant analogues thereof, Equation (1) strongly violates the Lorentz invariance. This condition requires the foliation of a spacetime manifold into (3 + 1)-dimensional spacelike hypersurfaces, and makes mass-energy a three-dimensional scalar, not a time component of a four-dimensional vector. As a result, relativistic scalar fields models offer a useful, but approximate description of superfluidic phenomena (which is valid for small wavefunction amplitudes and low momenta of excitations, running ahead), while the rigorous extrapolation of BEC and superfluid notions into the relativistic domain requires special treatment. Therefore, it is not surprising that some “superfluid” models of dark matter; which are based on classical perfect fluid models, scalar field theories or scalar-tensor theories of gravity, turn out to be vulnerable to experimental verification [8]. Universe 2020, 6, 180 22 of 25

Appendix B. Assumptions and Approximations Let us summarize and enumerate the main simplifying assumptions and approximations underlying the derivation of our gravitational potential. We used them to keep our calculations as analytical and non-perturbative as possible, which is crucial for an essentially nonlinear theory, such as ours. First, it should be emphasized that superfluid vacuum theory, even when narrowed down to its logarithmic version, is a framework which potentially contains a set of models. Our chosen underlying quantum superfluid model, defined by Equation (11), is a minimal one. It can easily be expanded by adding polynomial terms |Ψ|k to the logarithmic nonlinearity in wave equations, to make it describe the phenomena in a more precise way. The reason for this is that the logarithmic nonlinearity is a leading-order approximation for the condensate-like matter, as discussed in [42]; but non-logarithmic terms can also come into play: one example is to be found in [43]. The nonlinear coupling function b = b(r, t) can also be made more detailed, to account for the thermodynamic environment in a more realistic way. Second, for a chosen model, the signs of coupling parameters can be changed, which is often equivalent to changing the topological sector a model belongs to. Because we do not precisely know the topological structure of the physical vacuum we live in, it essentially doubles the number of candidate models to be tested empirically. For example, changing the sign of the field-theoretical potential; or, alternatively, the overall sign of the nonlinear coupling b; switches between the topologically trivial, which is considered in this paper, and the topologically non-trivial sectors. In a logarithmic liquid model with constant nonlinear coupling, switching to the topologically non-trivial sector of the theory means that wavefunction changes from a droplet-like non-topological soliton type to a bubble-like topological soliton type, as discussed in [47,49]. Consequently, the forms of a trial wavefunction and induced potential in Section3 would also change. Third, we have avoided the problem of ambiguity in choosing trapping potential and boundary conditions by making the former identically zero, and the latter of a conventional quantum-mechanical type on an infinite spatial domain of trivial topology R3. These are simplifying assumptions which are yet to be proven to work, otherwise they must be replaced with something more sophisticated. Fourth, even if we have chosen our model correctly, there still remains the ambiguity of how to determine the state of the vacuum—is it in a ground state, an excited but metastable state, a pure or mixed state, a superposition of states, or even in a quantum transition between states? Any new empirical information about this could drastically change the assumptions underlying the computations in Section3. Fifth, the technical approximation which led us from Equations (16)–(18) might oversimplify the picture. It is suitable for the purposes of this study, but should be modified in more precise considerations. Finally, throughout the paper, we consider the physical vacuum alone: assuming that its small excitations, which would be observed by a R-observer as relativistic matter with deformed dispersion relations [11,17–21], do not back-react, for example, via interaction between scalar and tensor modes. This is obviously an over-simplified picture of reality, which is sufficient for our current purpose, but unlikely to be valid in general.

References

1. Carroll, S.M. The cosmological constant. Living Rev. Rel. 2001, 4, 1. [CrossRef] 2. Bullock, J.S.; Boylan-Kolchin, M. Small-scale challenges to the ΛCDM paradigm. Ann. Rev. Astron. Astrophys. 2017, 55, 343–387. [CrossRef] 3. Tulin, S.; Yu, H.-B. Dark matter self-interactions and small scale structure. Phys. Rept. 2018, 730, 1–57. [CrossRef] 4. Dirac, P.A.M. Is there an aether? Nature 1951, 168, 906–907. [CrossRef] 5. Volovik, G.E. The Universe in a Helium Droplet; Oxford University Press: Oxford, UK, 2009. Universe 2020, 6, 180 23 of 25

6. Huang, K. A Superfluid Universe; World Scientific: Hackensack, NJ, USA, 2016. 7. Zloshchastiev, K.G. Spontaneous symmetry breaking and mass generation as built-in phenomena in logarithmic nonlinear quantum theory. Acta Phys. Polon. 2011, 42, 261–292. [CrossRef] 8. Lisanti, M.; Moschella, M.; Outmezguine, N.J.; Slone, O. The inconsistency of superfluid dark matter with dynamics. arXiv 2019, arXiv:1911.12365. 9. Pethick, C.J.; Smith, H. Bose-Einstein Condensation in Dilute ; Cambridge University Press: Cambridge, UK, 2008. 10. Zloshchastiev, K.G. Volume element structure and roton-maxon- excitations in superfluid helium beyond the Gross-Pitaevskii approximation. Eur. Phys. J. B 2012, 85, 273. [CrossRef] 11. Zloshchastiev, K.G. Superfluid vacuum theory and deformed dispersion relations. Int. J. Mod. Phys. A 2020, 35, 2040032. [CrossRef] 12. Kapitsa, P.L. Viscosity of liquid helium below the λ-point. Nature 1938, 141, 74. [CrossRef] 13. Allen, J.F.; Misener, A.D. Flow of liquid helium II. Nature 1938, 141, 75. [CrossRef] 14. De Martino, S.; Falanga, M.; Godano, C.; Lauro, G. Logarithmic Schrödinger-like equation as a model for magma transport. Europhys. Lett. 2003, 63, 472–475. [CrossRef] 15. Zloshchastiev, K.G. Matrix logarithmic wave equation and multi-channel systems in fluid mechanics. J. Theor. Appl. Mech. 2019, 57, 843–852. [CrossRef] 16. Barceló, C.; Liberati, S.; Visser, M. Analogue gravity. Living Rev. Rel. 2005, 8, 12. [CrossRef][PubMed] 17. Zloshchastiev, K.G. Vacuum Cherenkov effect in logarithmic nonlinear quantum theory. Phys. Lett. A 2011, 375, 2305–2308. [CrossRef] 18. Schreck, M. Vacuum for Lorentz-violating . Phys. Rev. D 2017, 96, 095026. [CrossRef] 19. Nassif, C.; dos Santos, R.F.; Amaro, A.C., Jr. Lorentz violation with an invariant minimum speed as foundation of the tachyonic inflation within a Machian scenario. Phys. Dark Universe 2018, 22, 116–126. 20. Lecian, O.M. Alternative uses for quantum systems and devices. Symmetry 2019, 11, 462. [CrossRef] 21. Pu, J.; Li, G.P.; Jiang, Q.Q.; Zu, X.T. Deformed dispersion relation constraint with 1S − 2S transition. Chin. Phys. C 2020, 44, 014001. [CrossRef] 22. Rosen, G. Particlelike solutions to nonlinear complex scalar field theories with positive-definite energy densities. J. Math. Phys. 1968, 9, 996–998. [CrossRef] 23. Bialynicki-Birula, I.; Mycielski, J. Nonlinear wave mechanics. Ann. Phys. (N. Y.) 1976, 100, 62–93. [CrossRef] 24. Scott, T.C.; Shertzer, J. Solution of the logarithmic Schrödinger equation with a Coulomb potential. J. Phys. Commun. 2018, 2, 075014 . [CrossRef] 25. Alves, C.O.; de Morais Filho, D.C. Existence and concentration of positive solutions for a Schrödinger logarithmic equation. Z. Angew. Math. Phys. 2018, 69, 144. [CrossRef] 26. Alves, C.O.; de Morais Filho, D.C.; Figueiredo, G.M. On concentration of solution to a Schrödinger logarithmic equation with deepening potential well. Math. Meth. Appl. Sci. 2019, 42, 4862–4875. [CrossRef] 27. Wang, Z.-Q.; Zhang, C. Convergence from power-law to logarithm-law in nonlinear scalar field equations. Arch. Rational Mech. Anal. 2019, 231, 45–61. [CrossRef] 28. Bao, W.; Carles, R.; Su, C.; Tang, Q. Regularized numerical methods for the logarithmic Schrödinger equation. Numer. Math. 2019, 143, 461–487. [CrossRef] 29. Chen, S.; Tang, X. Ground state sign-changing solutions for elliptic equations with logarithmic nonlinearity. Acta Math. Hungar. 2019, 157, 27–38. [CrossRef] 30. Wen, L.; Tang, X.; Chen, S. Ground state sign-changing solutions for Kirchhoff equations with logarithmic nonlinearity. . J. Qual. Theory Differ. Equ. 2019, 47, 1–13. [CrossRef] 31. Li, H.; Zhao, X.; Hu, Y. Numerical solution of the regularized logarithmic Schrödinger equation on unbounded domains. Appl. Numer. Math. 2019, 140, 91–103. [CrossRef] 32. Zhang, L.; Hou, W. Standing waves of nonlinear fractional p-Laplacian Schrödinger equation involving logarithmic nonlinearity. Appl. Math. Lett. 2020, 102, 106149. [CrossRef] 33. Mohammadi, M. An energetically stable Q-ball solution in 3 + 1 dimensions. Phys. Scr. 2020, 95, 045302. [CrossRef] 34. Alves, C.O.; Ji, C. Multiple positive solutions for a Schrödinger logarithmic equation. Discret. Contin. Dyn. Syst. 2020, 40, 2671–2685. [CrossRef] Universe 2020, 6, 180 24 of 25

35. Shertzer, J.; Scott, T.C. Solution of the 3D logarithmic Schrödinger equation with a central potential. J. Phys. Commun. 2020, 4, 065004. [CrossRef] 36. Boudjeriou, T. On the diffusion p(x)-Laplacian with logarithmic nonlinearity. J. Elliptic Parabol. Equ. 2020, 1–12. [CrossRef] 37. Boudjeriou, T. Global existence and blow-up for the fractional p-Laplacian with logarithmic nonlinearity. Mediterr. J. Math. 2020, 17, 162. [CrossRef] 38. Zloshchastiev, K.G. Logarithmic nonlinearity in theories of quantum gravity: Origin of time and observational consequences. Grav. Cosmol. 2010, 16, 288–297. [CrossRef] 39. Scott, T.C.; Zhang, X.; Mann, R.B.; Fee, G.J. Canonical reduction for dilatonic gravity in 3 + 1 dimensions. Phys. Rev. D 2016, 93, 084017. [CrossRef] 40. Yasue, K. of nonconservative systems. Annals Phys. 1978, 114, 479–496. [CrossRef] 41. Brasher, J.D. Nonlinear wave mechanics, information theory, and thermodynamics. Int. J. Theor. Phys. 1991, 30, 979–984. [CrossRef] 42. Zloshchastiev, K.G. On the dynamical nature of nonlinear coupling of logarithmic quantum wave equation, Everett-Hirschman entropy and temperature. Z. Naturforsch. A 2018, 73, 619–628. [CrossRef] 43. Scott, T.C.; Zloshchastiev, K.G. Resolving the puzzle of sound propagation in liquid helium at low temperatures. Low Temp. Phys. 2019, 45, 1231–1236. [CrossRef] 44. Avdeenkov, A.V.; Zloshchastiev, K.G. Quantum Bose liquids with logarithmic nonlinearity: Self-sustainability and of spatial extent. J. Phys. B At. Mol. Opt. Phys. 2011, 44, 195303. [CrossRef] 45. Bouharia, B. Stability of logarithmic Bose-Einstein condensate in harmonic trap. Mod. Phys. Lett. B 2015, 29, 1450260. [CrossRef] 46. Zloshchastiev, K.G. Stability and metastability of trapless Bose-Einstein condensates and quantum liquids. Z. Naturforsch. A 2017, 72, 677–687. [CrossRef] 47. Zloshchastiev, K.G. Nonlinear wave-mechanical effects in Korteweg fluid magma transport. Europhys. Lett. (EPL) 2018, 122, 39001. [CrossRef] 48. Kraiev, M.; Domina, K.; Kraieva, V.; Zloshchastiev, K.G. Soliton solutions of logarithmic wave equation and their application for polycrystalline metals. J. Phys. Conf. Ser. 2019, 1416, 012020. [CrossRef] 49. Zloshchastiev, K.G. Temperature-driven dynamics of quantum liquids: Logarithmic nonlinearity, phase structure and rising force. Int. J. Mod. Phys. B 2019, 33, 1950184. [CrossRef] 50. Rylov, Y.A. Spin and as attributes of ideal fluid. J. Math. Phys. 1999, 40, 256–278. [CrossRef] 51. Mannheim, P.D.; Kazanas, D. Exact vacuum solution to conformal Weyl gravity and galactic rotation curves. Astrophys. J. 1989, 342, 635–638. [CrossRef] 52. Mannheim, P.D. Linear potentials and galactic rotation curves. Astrophys. J. 1993, 419, 150–154. [CrossRef] 53. Znojil, M.; Semorádová, I. Quantum square well with logarithmic central spike. Mod. Phys. Lett. A 2018, 33, 1850009. [CrossRef] 54. Znojil, M.; Semorádová, I. Log-anharmonic oscillator and its large-N solution. Mod. Phys. Lett. A 2018, 33, 1850223. [CrossRef] 55. Toomre, A. On the distribution of matter within highly flattened galaxies. Astrophys. J. 1963, 138, 385–392. [CrossRef] 56. Casertano, S. Rotation curve of the edge-on NGC 5907: Disc and halo . Mon. Not. R. Astron. Soc. 1983, 203, 735–747. [CrossRef] 57. Rubin, V.C.; Ford, W.K., Jr.; Thonnard, N. Rotational properties of 21 SC galaxies with a large range of and radii, from NGC 4605 (R=4kpc) to UGC 2885 (R=122kpc). Astrophys. J. 1980, 238, 471–487. [CrossRef] 58. Corbelli, E.; Salucci, P. The extended rotation curve and the of M33. Mon. Not. R. Astron. Soc. 2000, 311, 441–447. [CrossRef] 59. Chemin, L.; Carignan, C.; Foster, T. HI kinematics and dynamics of Messier 31. Astrophys. J. 2009, 795, 1395–1415. [CrossRef] 60. Corbelli, E.; Lorenzoni, S.; Walterbos, R.; Braun, R.; Thilker, D. A wide-field HI mosaic of Messier 31. II. The disk warp, rotation and the dark matter halo. Astron. Astrophys. 2010, 511, A89. [CrossRef] 61. Kam, S.Z.; Carignan, C.; Chemin, L.; Foster, T.; Elson, E.; Jarrett, T.H. HI kinematics and mass distribution of Messier 33. Astron. J. 2017, 154, 41. [CrossRef] Universe 2020, 6, 180 25 of 25

62. Zloshchastiev, K.G. Resolving cosmological singularity problem in logarithmic superfluid theory of physical vacuum. J. Phys. Conf. Ser. 2020, 1557, 012038. [CrossRef] 63. Freedman, W.L.; Madore, B.F.; Gibson, B.K.; Ferrarese, L.; Kelson, D.D.; Sakai, S.; Mould, J.R.; Kennicutt, R.C., Jr.; Ford, H.C.; Graham, J.A.; et al. Final results from the key project to measure the Hubble constant. Astrophys. J. 2001, 553, 47–72. [CrossRef] 64. Chen, G.C.-F.; Fassnacht, C.D.; Suyu, S.H.; Rusu, C.E.; Chan, J.H.H.; Wong, K.C.; Auger, M.W.; Hilbert, S.;

Bonvin, V.; Birrer, S.; et al. A SHARP view of H0LiCOW: H0 from three time-delay systems with adaptive optics imaging. Mon. Not. R. Astron. Soc. 2019, 490, 1743–1773. [CrossRef] 65. Pesce, D.W.; Braatz, J.A..; Reid, M.J.; Riess, A.G.; Scolnic, D.; Condon, J.J.; Gao, F.; Henkel, C.; Impellizzeri, C.M.V.; Cuo, C.Y.; et al. The Megamaser Cosmology Project. XIII. Combined Hubble constant constraints. Astrophys. J. Lett. 2020, 891, L1. [CrossRef]

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

c 2020 by the author. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).