Below is an interesting article from December 9, 2020, on dark matter (which was discussed above) as Gravitational Solitons in the weak field limit.

Interestingly, Nathan Seiberg and Edward Witten ( both mentioned in Section 4 and both are at The Institute for Advanced Study at Princeton, NJ) posited the full set of corrections for SU(2) Yang–Mills theory was calculated in "

Electric – magnetic duality, monopole condensation, and confinement in N=2 supersymmetric Yang–Mills theory," in the process creating a subject that is today known as

Seiberg–Witten theory. They extended their calculation to SU(2) gauge theories with fundamental matter in

Monopoles, duality and chiral symmetry breaking in N=2 supersymmetric QCD. These results were later extended for various gauge groups and matter contents, and the direct gauge theory derivation was also obtained in most cases. For gauge theories with gauge group U(N) the Seiberg-Witten geometry has been derived from gauge theory using

Nekrasov partition functions in 2003 by

Nikita Nekrasov and

Andrei Okounkov and independently by

Hiraku Nakajima and

Kota Yoshioka.

M-theory is a theory in physics that unifies all consistent versions of superstring theory. Edward Witten first conjectured the existence of such a theory at a string-theory conference at the University of Southern California in the spring of 1995.

Hartmann

12/11/2020

**Dark Matter as Gravitational Solitons in the Weak Field Limit **
by

Torsten Asselmeyer-Maluga
1,*,†

and

Jerzy Król
2,†

1

German Aerospace Center (DLR), Rosa-Luxemburg-Str. 2, D-10178 Berlin, Germany

2

Chair of Cognitive Science and Mathematical Modeling, University of Information Technology and Management, ul. Sucharskiego 2, 35-225 Rzeszw, Poland

*

Author to whom correspondence should be addressed.

†

These authors contributed equally to this work.

Universe 2020, 6(12), 234;

https://doi.org/10.3390/universe6120234
Received: 16 October 2020 / Revised: 22 November 2020 / Accepted: 7 December 2020 / Published: 9 December 2020

(This article belongs to the Special Issue

Geometric and Topological Models of Dark Matter and Dark Energy)

Download PDF
**Abstract**
In this paper, we will describe the idea that dark matter partly consists of gravitational solitons (gravisolitons). The corresponding solution is valid for weak gravitational fields (weak field limit) with respect to a background metric. The stability of this soliton is connected with the existence of a special foliation and amazingly with the smoothness properties of spacetime. Gravisolitons have many properties of dark matter, such as no interaction with light but act on matter via gravitation. In this paper, we showed that the gravitational lensing effect of gravisolitons agreed with the lensing effect of usual matter. Furthermore, we obtained the same equation of state w=0 as matter.

Keywords:

dark matter;

gravitational solitons;

exotic smoothness of 4-manifolds
**1. Introduction**
In the 1930s, astronomers, among them the Swiss astrophysicist Fritz Zwicky [

1,

2], made observations which seem to imply that the orbits of stars in various galaxies did not have the expected relation between speed and distance from the galactic center. However, it took more than 30 years, until 1975, for careful studies of several different galaxies to confirm the results of Zwicky. Using new techniques, the astronomer Vera Rubin [

3] showed that the speeds of stars in orbit in every galaxy tended to be independent of the distance from the center of the galaxy. However, this observation is in sharp contrast to the theory: If the visible mass of galaxies is responsible for star acceleration, then the speeds should vary as 1/r−−−√, with respect to the distance. With this contradiction in mind, astronomers and physicists quickly realized that almost all (90%) of the mass of galaxies resides in an invisible halo of unknown matter, called dark matter, sticking out into space for a distance around 10 times the visible galactic radius. Dark matter is matter that emits or reflects minimal to no light, but does have a gravitational influence. Evidence for dark matter is present in

- the motion of stars in galaxies,
- the orbits of galaxies in galaxy clusters,
- the temperature of intracluster gas in galaxy clusters and
- the gravitational lensing of distant galaxies.

The appearance of large-scale structure (e.g., the distribution of galaxies) is very hard to understand without Dark matter. This large scale structure is particularly encoded into the microwave background as measured by the satellites like Cosmic Background Explorer (COBE), Wilkinson Microwave Anisotropy Probe (WMAP) or PLANCK [

4,

5]. One way to accommodate this is to go to a dark matter model in which you have cold dark matter (CDM) to act as a seed for galaxy formation. But there are also results which contrast with the CDM hypothesis like the distribution of satellite galaxies around the Milky Way [

6], the absence of dark matter in our sun system [

7] and the existence of a vast polar structure of satellite galaxies around the Milky Way [

8]. Dark matter can be divided into cold, warm, and hot categories, according to the free streaming length (FSL) indicating how far corresponding objects moved due to random motions in the early universe. Warm dark matter comprises particles with an FSL comparable to the size of a protogalaxy. Predictions based on warm dark matter are similar to those for cold dark matter on large scales, but with less small- scale density perturbations. A conjectured candidate is the sterile neutrino. The warm dark matter model is considered to be a better fit to observations. However, in both cases, there is currently no experimental detection of dark matter particles. Hot dark matter consists of particles, like the neutrino, whose FSL is much larger than the size of a protogalaxy. However, current limits on the neutrino masses by using cosmic microwave background data induces the result that the neutrinos cannot be dark matter. Here we will follow a different path. To explain all these properties, we need a kind of matter which cannot be compressed, couples to gravity only and is stable. In this paper we present such a substance: a gravitational soliton (in the weak field limit—it will be explained later). It is a special very stable structure of the spacetime appearing in two types, the gravisoliton S and the anti-gravisoliton A. Furthermore, there is no coupling to light and therefore only small acoustic oscillations (caused by possible collisions of SS and AA) are allowed.

Here we will consider gravisolitons as another component of dark matter. Loosely speaking, Gravisolitons are topological defects of the space like a density wave of space. They behave like matter but with different scattering properties. At first we will study the properties of the spacetime which is needed to generate gravitational solitons. Then we will connect the appearance of gravisolitons with a special codimension-1 foliation (also known as helical wobble, see [

9]). In

Section 4 we will discuss how a gravisoliton can be generated by matter. At the same time, this mechanism uncovers the topological reason for the appearance of gravisolitons: the underlying spacetime must admit an exotic smoothness structure. Dark matter was found by the gravitational lensing effect. Therefore, we will study the deflection of light rays by the gravisoliton in

Section 5. Interestingly, we will find no difference to ordinary (baryonic) matter. Then we will discuss a coherent model of spacetime with all necessary properties. Here, we will see that all relevant spacetime with exotic smoothness structure are fulfilling these properties. Finally, in

Section 7 we will also discuss the state equation for gravisolitons to be w=0, i.e., no difference to usual matter. The discussion of strong field limit is reserved to a separate paper.

**2. Spacetime Properties for the Generation of Gravitational Solitons**
Dark matter has many properties known for matter itself, but it has also some other specific properties, like the very weak interaction to usual matter, which leads to the idea that dark matter consists only partly of (particle-like) matter (WIMPS etc.). Here we will follow this idea, i.e., we will consider dark matter as a special property of the space. It is just a trial which, nevertheless, is aiming toward deciding important question namely whether and how the idea can be realized by a spacetime itself. What properties of spacetime are required to this end? In the following sections we will describe these properties. First let us consider the following ansatz for the metric:

ds2=f(t,z)(dz2−dt2)+gab(t,z)dxadxb

(1)

for a,b=1,2 and with the non-diagonal 2-metric gab. The non-diagonality of gab is very important, it is responsible for the difference between non-linear solutions (for non-diagonal gab) and linear solutions (for diagonal gab). The linear case was first considered by Einstein and Rosen in 1937 (gravitational cylindrical waves). It is convenient to introduce the coordinates

θ=12(z−t)ζ=12(z+t)

transforming dz2−dt2 into dθdζ and α=det(gab)−−−−−−−−√ for the determinant. Then the Einstein equation Rab=0is represented by the matrix equation

∂∂θ(α∂g∂ζg−1)+∂∂ζ(α∂g∂θg−1)=0

(2)

and we obtained the additional conditions from the equations R00=0, R33=0 and R03=0. For the following analysis, we will refer to the first two chapters of the book [

10] where the solution is presented using the Inverse Scattering Method (ISM). The method itself is complicated and cannot be presented here completely. We will try to give some background here to make the paper as self-contained as possible. At first we must explain the wording weak field limit. A solution of the non-linear equation (

2) can be obtained by the ISM as perturbation with respect to a background metric. Therefore, this solution is only valid for weakly gravitational fields. Then the solution is a small perturbation, i.e., we will only obtain a gravisoliton for weak gravitational fields. Or we will say that the gravisoliton is obtained in the weak field limit. The strong field limit will be considered in a separate paper.

Main idea relies on a proper introduction of a background metric. Let a particular solution g0 serves as a background metric so that metric g(1) derived from the ISM method is a soliton solution with respect to g0. g0 can be chosen as a diagonal metric (for a reason discussed in

Section 4)

g0=diag(αeu0,αe−u0)

where α=detg0−−−−−−√and

∂2α∂θ∂ζ=0∂∂θ(α∂u0∂ζ)+∂∂ζ(α∂u0∂θ)=0

follows from (

2). Now we will get a metric g from the background metric g0by using ISM. Main idea of the ISM is the introduction of two differential operators, the L-A or Lax pair, so that the eigenfunctions of these operators determine the new solution out from the old solution. In this case, we started with the constant diagonal (background) metric. Then the 1-soliton solution is given by

g(1)=1μ1⋅cosh(ρ1)((μ21eρ1+α2e−ρ1)eu0α2−μ21α2−μ21(μ21eρ1+α2e−ρ1)e−u0)

where μ1,ρ1are integration constants. It is interesting to note that the off-diagonal elements are the topological charge of the soliton. Unfortunately, from the physics point of view this solution is difficult to interpret. Therefore, we will study invariant reparameterizations to study the invariant properties of the solution. With the notation

∂∂ζ(ln(g11α−1))=R1cos(γ+ω2)∂∂θ(ln(g11α−1))=R2cos(γ−ω2)g11α−1∂∂ζ(g12g−111)=R1sin(γ+ω2)g11α−1∂∂θ(g12g−111)=R2sin(γ−ω2)

(3)

we have a new parameterization so that R1,R2 and ω are an invariant parameterization (not depending on arbitrary linear transformations of the coordinates x1,x2) for the metric g (there is also a system of two equations to determine γ which depends only on R1,R2 and ω). Now we will make our main assumption:

α=const.⟶R1=R2=const.

which will be motivated later. Then we obtain the following equation for the variable ω

∂2ω∂ζ∂θ=R1R2sinω

which is the famous sine-Gordon equation. This equation admits soliton solutions of the form (now in coordinates (z,t))

ω(z,t)=4⋅arctan(emC(z−vt+zo))

(4)

with velocity v=−2αμ1(α2+μ21)−1, C=(1−v2)−1 and mass m=R1=R2 which are stable over the time. The constant μ1 is an integration constant from the ISM. The solution depends on the sign of α2−μ21 which is the topological charge of the soliton. The topological charge is defined as the difference ω(+∞,t)−ω(−∞,t) in units of 2π. Usually, the solution (

4) of the sine-Gordon equation describes a soliton where the sign of α2−μ21(topological charge) determines whether it is a gravisoliton (positive sign) or anti-gravisoliton (negative sign). Or, the gravisoliton S has topological charge +1 whereas the anti-gravisoliton A has −1. The sine-Gordon soliton is a non-linear wave which can be described by a rotating vector. Then the sign of the topological charge gives the direction of the rotation (with respect to the propagation direction). The difference between S and A defines the behavior with respect to the scattering. The scattering of solitons like SS or AA shows an amazing behavior: there is a repulsive force so that these solitons bounce back by a collision. In contrast the soliton SA can be represented by a bound state (a so-called breather solution), so there is an attractive force in this system. By considering bound states of solitons (like SA but technically generated by a Bäcklund transformation) one obtains 2-gravisolitons which form only stable stationary solution (including an axisymmetry). These 2-solitons in a Minkowski background are equivalent to a Kerr-NUT or Schwarzschild solution [

10]. Then one obtains an attraction of a test particle by a 2-gravisoliton (or 2-anti-gravisoliton). Summarizing this approach, we need to construct a spacetime

- with metric (1): ds2=f(t,z)(dz2−dt2)+gab(t,z)dxadxb where the 2-metric gab is parameterized by the t,z—coordinates only,
- the 2-metric gab is non-diagonal, i.e., gab≠0 for a≠b and
- det(gab)=const., i.e., the area of the surface defined by gab is fixed.

How to fulfill these conditions? At first one needs a spacetime which splits into a (non-trivial, curved) surface and a (conformally flat) (1+1)—dimensional subspace (first condition). The surface is locally generated by two independent vector fields (the dual to dx1,dx2) which are never orthogonal to each other (second condition). Furthermore this surface has a finite volume (third condition). The first and third conditions can be simply realized for instance by spaces like R1,1×S for every compact closed surface S. The problem is only with the second condition. In principle, we need a complicated codimension-1 foliation of a 3-dimensional manifold (formed by x1,x2,z) so that the three generating vector fields are never orthogonal to each other. In the next section we will construct this foliation for the 3-sphere but also for all other compact 3-manifolds.

**3. A Codimension-1 Foliation Realizing the Conditions for Gravitational Solitons**
In the previous section we discussed gravitational solitons, called gravisolitons, which are connected with the metric (

1). Main characteristics of this metric is the appearance of a non-diagonal 2-metric gab. Then the non-diagonal term of the 2-metric can be interpreted as non-orthogonal vector fields in this space. With other words, we must construct a space or spatial subspace of the spacetime which is generated by non-orthogonal vector fields. A direct way would be the foliation of the space by surfaces described by the 2-metric gab. Does such foliation exist? And if it exists, what are conditions guaranteeing this which must be implied on the foliation and the spacetime?

First, we will give a flavor of the main idea without the technical details. Our task is now to construct a foliation of a 3-manifold into surfaces so that the three generating vector fields are never orthogonal to each other. This problem was solved by Thurston [

9] who constructed a codimension-1 foliation of the 3- sphere which can be easily extended to all 3-manifolds. The details of the complicated construction are not important for the following discussion and can be found in [

11]. Main idea is the usage of the group PSL(2,R). Then the vector fields of the foliation are the left-invariant vector fields of the Lie algebra of this group. Furthermore, we will construct an invariant of the foliation (Godbillon–Vey invariant) which is directly connected to the non-orthogonal vector fields. Differently expressed: the vector fields are non- orthogonal if this invariant is non-zero. The reader not interested in the technical details of this construction can switch to the next section keeping in mind that there is such a foliation with the suitable properties.

A foliation (M,F) of a manifold M is an integrable Sub-bundle F⊂TM of the tangent bundle TM. The leaves L of the foliation (M,F) are the maximal connected submanifolds L⊂M with TxL=Fx for all x∈L. The number dimM−dimL is the codimension of the foliation. A codimension-1 foliation on a 3-manifold M can be constructed by a smooth 1-form ω fulfilling the integrability condition dω∧ω=0. This construction can be found in [

11]. In [

12] the authors analyzed the local structure of the foliation: the three vector fields (forming the frame of the foliation) are not orthogonal to each other. This foliation is constructed by using the group PSL(2,R), the isometry group of the hyperbolic plane H2. Let T,N,Z be the frame formed by three vector field dual to the one-forms ω,η,ξ as a basis for the tangent bundle TM (remember that the tangent bundle of every 3-manifold is trivial, see [

13]). Now one considers the vector fields as left-invariant vector fields of the group PSL(2,R) (forming the basis of the Lie algebra) fulfilling the commutator relations [Z,N]=Z, [N,T]=T and [Z,T]=N. As shown in [

12], the vector fields are not orthogonal to each other leading to a (local) non-zero torsion τ≠0. For codimension-1 foliations of 3-manifolds, there exists the Godbillon–Vey invariant which we are going to describe now. Recall that codimension-1 foliation on M is defined by a PSL(2,R)—invariant integrable one-form ω (dω∧ω=0) and we define another one-form η by dω=−η∧ω. Then the Godbillon–Vey class gv=η∧dηis a closed form (by using the integrability) which is not exact. Furthermore, the integral

GV(M)=∫Mη∧dη

(5)

is a topological invariant of the foliation, known as Godbillon–Vey number GV(M). From the physics point of view, it is the abelian Chern–Simons functional. The mathematical construction of the foliation can be found in [

9,

14] including a global calculation of the Godbillon–Vey invariant. However, here we will concentrate on the local expression for the Godbillon–Vey invariant. Let κ,τ be the curvature and torsion of a normal curve to the foliation, respectively. Furthermore, let T,N,Z be the frame formed by three vector field dual to the one-forms ω,η,ξ and let lTbe the second fundamental form of the leaf. Then the Godbillon–Vey class is locally given by

η∧dη=κ2(τ+lT(N,Z))ω∧η∧ξ

(6)

where τ≠0 for PSL(2,R) invariant foliations i.e., [Z,N]=Z, [N,T]=T and [Z,T]=N. Therefore at least two vector fields are not orthogonal to each other and the corresponding metric gab is not diagonal. Finally, let Mbe a compact 3-manifold with codimension-1 foliation of non-zero Godbillon–Vey number GV(M)≠0. Then the spacetime M×R (or M×[0,1] for finite time) has all properties to admit gravitational solitons (in the weak field limit).

**4. How to Generate Gravitational Solitons**
In the previous section we connected the gravitational soliton with a foliation of the 3-space. The result can be expressed by the simple statement: if the space admits a codimension-1 foliation with non-zero Godbillon–Vey invariant (

5), locally given by (

6), then, there is a gravitational soliton (gravisoliton) in space described by the metric (

1). Therefore, if we can get a non-zero Godbillon invariant then a gravisoliton is generated. Starting point is a general but arbitrary foliation of the space into surfaces. In [

15], we described a strong relation (using the work of [

16]) between a surface and a spinor ψ. We will give a different meaning to the Godbillon–Vey invariant. In principle, this invariant is generated by a 3-form η∧dη with the real 1-form η. Now we interpret η as a U(1) gauge field A which is an imaginary valued 1-form i.e., A=iη. Then we couple this gauge field η to the spinor ψ. The Godbillon–Vey invariant is also known as abelian Chern–Simons form. Then the action of a spinor coupled to the gauge field given by the Chern–Simons form is a dynamical theory to generate a foliation of the space (see the Dirac–Chern–Simons functional below). The foliation is defined by the gradient flow. However, the gradient flow is equivalent to the Seiberg–Witten equations, a theory to detect different smoothness structures on 4-manifolds. A non-trivial solution of this gradient system is given by ψ≠0,η≠0, which is equivalent to the existence of a non-trivial Godbillon–Vey invariant which means that a gravisoliton is being generated. However, this solution gives also a non-trivial solution of the Seiberg–Witten equation for the 4-dimensional spacetime. This solution is a necessary condition for an exotic smoothness structure of the spacetime. A smoothness structure is the maximal smooth atlas of a manifold, i.e., the collection of smooth charts covering the manifold. In principle, the smoothness structure is the main topological information about the smooth manifold. Loosely speaking, all properties of differentiable quantities are determined by the smoothness structure, even one cannot define a differential equation on a manifold without fixing the smoothness structure. In dimension 4 there are manifolds which are homeomorphic to each other but not diffeomorphic, i.e., these manifolds differ by the smoothness structure. Main results of this section are the connection between the existence of gravisolitons in space and the smoothness structure of the spacetime.

Now we go into the details of this construction. In [

15] we presented a geometrical/topological model for matter. In general, for a manifold M with boundary ∂M=Σ one has the action (see [

17])

SEH(M)=∫MRg√d4x+∫ΣHh−−√d3x

where H is the mean curvature of the boundary with metric h. For the following discussion, we consider the boundary term

SEH(Σ)=∫ΣHh−−√d3x

(7)

along the boundary Σ (a 3-manifold). Following [

15], the action (

7) over a 3-manifold Σ is equivalent to the Dirac action of a spinor over Σ. Interestingly, this relation is induced by a stronger relation between surfaces and spinors (see [

18] for a complete discussion). Main result of [

15] is the following relation between the corresponding Dirac operators

DMΦ=DΣψ−Hψ

(8)

where DΣ or DM denote the Dirac operators on the 3-manifold Σ or 4-manifold M, respectively. Using the work of [

16], the spinor ϕ directly defines the embedding (via an integral representation) of the 3-manifold. Then the restricted spinor Φ|Σ=ϕ is parallel transported along the normal vector and Φ is constant along the normal direction. However, then the spinor Φmust fulfill

DMΦ=0

(9)

i.e., Φis a parallel spinor. Finally, we get

DΣψ=Hψ

(10)

and

∫ΣHh−−√d3x=∫Σψ¯¯DΣψh−−√d3x

(11)

Therefore, we get a direct relation between spinors and geometry. Now we will use this relation to get a connection to the foliation.

For every codimension-1 foliation, it is known: If the 1-form ω defines the leaf via the equation ω=const. then the dual of the 1-form η defines a vector normal to the leaf. However, then the spinor Φ (representing the leaf) is constant along the normal direction fulfilling the relation (

9) i.e., Φ is a parallel spinor. Furthermore, in [

15] we showed that a fermion is given by a knot complement admitting a hyperbolic structure. Then calculations in [

19] implied that the particular knot is only important for the dynamical state (like the energy or momentum) but not for charges, flavors etc.

Now we will use this relation between spinors, surfaces and the foliation to reinterpret the Dirac equation DMΦ=0along the normal direction as covariant constant spinor

DΣηψ=0

where the one-form η is interpreted as abelian gauge field (with structure group U(1). By definition, the covariant constant 1-form ω

Dηω=dω+η∧ω=0

defines the foliation, because the integrability condition ω∧dω=0 is automatically fulfilled. Then we reinterpret the Godbillon–Vey invariant gv=η∧dη of the foliation as abelian Chern–Simons form for the abelian gauge field η. Now we use the relation between the spinor and the surface (as given by ω=const.) to define a similar relation as Dηω=0: the spinor ψ is covariant constant with respect to ηor

DΣηψ=0.

To take the special properties of the foliation into account, we must consider a non-zero Godbillon–Vey invariant at the same time. That means that for a fixed foliation, the coupling between the abelian gauge field η and the spinor ψto the Dirac–Chern–Simons action functional

SDCS=∫Σ(ψ¯¯DΣηψh−−√d3x+η∧dη)

on the 3-manifold is constant with the critical points at the solution

DΣηψ=0dη=τ(ψ,ψ)

where τ(ψ,ψ) is the unique quadratic form for the spinors locally given by ψ¯¯γμψ. If one has a spacetime Σ×I, then we consider the solution which is translationally invariant. Alternatively, it is a spacetime with foliation induced by the foliation of Σextended by a translation. Mathematically we must consider the gradient flow

ddtηddtψ==dη−τ(ψ,ψ)DΣηψ

From every solution of this gradient system, one can construct the corresponding foliation with Godbillon–Vey invariant. It is interesting to note that a solution with vanishing Godbillon–Vey invariant leads also to vanishing spinor. In [

20,

21] this system was shown to be equivalent to the Seiberg–Witten equation for Σ×I by using an appropriate choice of the SpinC structure. Then we have the result that a non-trivial foliation together with the existence of fermions induce a non-trivial solution of the gradient system which results in a non-trivial solution of the Seiberg–Witten equations. However, this non-trivial solution (i.e., ψ≠0,η≠0) is a necessary condition for the existence of an exotic smoothness structure (but not sufficient).

Here we get an unexpected relation between the existence of gravisolitons and topological (or better differential-topological) properties of the spacetime. This property of the spacetime is called exotic smoothness in mathematics (see the book [

22] for an overview).

Again, exotic smoothness denotes a different (non-diffeomorphic) smooth atlas of a 4-manifold. Two 4-manifolds with different smoothness structures are topologically equivalent (homeomorphic) but smoothly different. From the physics point of view, two non-diffeomorphic spacetimes are different if they represent (physically) different systems.

Finally, we are concluding that the foliation above is connected with the smoothness properties of the spacetime. Furthermore, if one starts choosing an exotic smoothness structure instead of the standard one for the spacetime then, there gravisolitons must exist automatically. By using the relation between spinors and surfaces, there is a natural coupling between a spinor and the gravisoliton, or the gravisoliton is generated by a spinor (electrically charged). The description here is only the first sign of how the appearance of gravisolitons could be understood in terms of physical spinor fields. More work is needed to understand it completely.

**5. Gravitational Lensing for Gravitational Solitons**
In this section, we will investigate the gravitational lensing property of the gravitational solitons. To study the effect of these solitons on the propagation of light, we will introduce the optical metric and compare it with the metric of the gravitational soliton. Unexpectedly, we will obtain the same effect as for usual matter. Therefore, one cannot see a difference between the gravitational soliton and matter by using gravitational lensing. This section is important to understand that gravitational solitons cannot be distinguished from matter by gravitational lensing.

Dark matter was indirectly recognized by gravitational lensing. Therefore we are going to study the deflection of light induced by gravitational solitons. The idea is to introduce a metric, the optical metric, which describes this deflection of light.

The optical metric with the line element ds2opt in Schwarzschild coordinates (t,r,ϑ,ϕ)is usually given by

ds2opt=−exp(2A(r))dt2+exp(2B(r))dr2+r2(dϑ2+sin2ϑdϕ2)

with some functions A(r),B(r). To determine these two functions, one must put it in Einstein’s field equations. Here one uses an energy-momentum tensor Tμν=diag(ρ,p,p,p) in the local flat metric, where ρ=ρ(r) denotes the density and p=p(r)the pressure of the lens model. Then one obtains

exp(−2B(r))=1−2μ(r)r

with

μ(r)=4πG∫0rρ(r′)r'2dr′

where ρ(r) denotes the energy density. However, the solution (

4) implied that the dr2part of the corresponding gravitational soliton solution is in a good approximation constant, so that we get

exp(−2B(r))=1−2μ(r)r=const.

However, then we will obtain

μ(r)∼r

and the energy density of a gravitational soliton is given by

ρ(r)∼1r2

i.e., we obtain the same gravitational lensing effect as the singular isothermal sphere serving as a lensing model for galaxies. In [

23], the deflection angle was determined from the above relation. Let us explain briefly how the deflection angle in the solitonic model is obtained. Let σ2be the amplitude of the soliton, then we set

ρ(r)=σ22πGr2

implying μ(r)=2σ2r. This is known as the solution of the Tolman– Oppenheimer–Volkoff equation in the nonrelativistic limit (i.e., first order in σ2). The Gaussian curvature vanishes for r>0 in this limit, so that the equatorial plane in the optical metric is a cone with a singular vertex at r=0. The deficit angle of the conical optical metric is the deflection angle δ. Thus, for our model the deflection angle is given by

δ=4πσ2

in agreement with the gravitational lensing of galaxies. It is interesting to note that this lensing is similar to the lensing of a cosmic string. Finally, gravitational solitons have the same effect as usual matter and gravitational lensing cannot distinguish between the gravitational soliton and matter.

**6. A Coherent Model for a Spacetime with These Properties**
In the

Section 4, we found an unexpected relation between gravisolitons and (differential-)topological properties of the spacetime. In particular, the existence of gravisolitons is connected with a class of complicated foliations (with non-trivial Godbillon–Vey invariant) leading to a coupling of the gravisoliton to a spinor. Then the last relation (spinor-gravisoliton coupling) can be fulfilled for a class of spacetimes, known as spacetime with exotic smoothness structure. In this section, we will construct a suitable spacetime and study its properties. Finally, we will see that for this class of spacetimes we obtain the remaining conditions which are needed to construct the gravisoliton (i.e., α=det(g)=const. and the non-diagonality of the 2-metric).

In a series of papers we have described the class of spacetimes with exotic smoothness structure. The influence of the smoothness structure on the topological properties of the spacetime is tremendous. Let us consider the spacetime given by the smooth manifold S3×R. In the standard smoothness structure on S3×R one has a global foliation of the spacetime into slides S3×{t} for every t∈R (or by copies of S3). In contrast to that, the 4-spacetime S3×Rwith exotic smoothness structure is given by a sequence …→Σ1→Σ2→… of topological transitions where all Σi have the same homology like S2×S1. In the limit of infinite transitions, one gets a manifold which is topologically equivalent to the 3-sphere. We remark that this is a general rule for all exotic versions of spacetimes N×R (with N a compact 3-manifold). To understand the general structure of the Σn, we must understand the general decomposition of 3-manifolds into pieces.

Every 3-manifold can be decomposed into prime 3-manifolds according to

K1#⋯#Kn1#n2S1×S2#n3S3/Γ

where Γ is a finite subgroup of SO(4). Every (irreducible) 3-manifold Ki can be also split into a hyperbolic and a graph manifold. As shown in [

18] and further extended in [

15,

19], the 3-manifolds Ki are the matter part. It consists of hyperbolic 3-manifolds as fermions and graph manifolds (torus bundles) for the interaction (gauge bosons). The further details of this model are not important here, but we want to call the reader’s attention to the sum operation #. By definition we have for the connecting sum K1#K2the expression

K1#K2=(K1∖D3)∪S2(S2×[0,1])∪S2(K2∖D3)

with the tube S2×[0,1] to identify K1∖D3 and K2∖D3 along its common boundary S2. Therefore we obtained an additional manifold S2×[0,1], called sphere bundle, between the matter Ki. The corresponding spacetime version of these sphere bundles admit the topology S2×[0,1]2. Obviously, we have the following properties:

- S2×[0,1]2 admits a Lorentz metric with coordinates (x1,x2) for the sphere S2 and (z,t) for [0,1]2 and
- the line element is given by

ds2=f(t,z)(dz2−dt2)+gab(t,z)dxadxb

in agreement with (1).

The next two properties of the spacetime are directly related to the exotic smoothness structure on it. However, from the point of view of geometry, exotic smoothness is connected with hyperbolic geometry (see for instance [

24,

25] and the result stating that 4-manifolds with positive scalar curvature in the standard smoothness structure fail to be paired with hyperbolic geometry). In [

26] we analyzed this situation and explained the reason: exotic smoothness structures enforce the appearance of saddle points (in the Morse-theoretic sense, i.e., non-canceling 2-/3-handle pairs) carrying a hyperbolic geometry. The corresponding analysis was extended in [

27,

28]. This fact has two important consequences: hyperbolic 3- and 4-manifolds (with finite volume, i.e., the neighborhood of the saddle points) have the property called Mostow rigidity. As shown by Mostow [

29] and extended by Prasad [

30], every hyperbolic n—manifold n>2 (with finite volume) has the following property: Every diffeomorphism (especially every conformal transformation) of a hyperbolic n—manifold (with finite volume) for n>2 is induced by an isometry. Therefore, one cannot scale a hyperbolic 3- and 4-manifold. Then the volume vol() and the curvature (both combined into the Chern–Simons invariant) are topological invariants. Secondly, the foliation of this hyperbolic 4-manifold part admits a non-trivial Chern–Simons invariant and also a non-trivial Godbillon–Vey invariant (foliations of type III) [

31]. Both consequences imply the following properties:

- The background for the sphere bundle has a constant scale-invariant volume (by Mostow rigidity) from the matter, i.e., det(g)=const.
- The metric g is non-diagonal (otherwise the Godbillon–Vey invariant of the foliation vanishes, see [32] and Section 3and Section 4).

All four properties together are the prerequisites to obtain the model of the gravisoliton as dark matter.

**7. The Equation of State for Gravisolitons**
In the astrophysical context, the scaling behavior of the energy density ρ(a) depends on the scale factor a(t). Now we consider the fluid equation

ρ˙+3a˙a(ρ+p)=0

and the equation of state p=w⋅ρ. Usually, one assumes a constant w so that one gets

ρ∼a−3(1+w)

For gravisolitons, the energy E is constant. In the previous section we discussed that the volume of the sphere S2in the sphere bundle S2×[0,1] (representing the gravisoliton) is fixed. Furthermore, we know from exotic smoothness structure that this sphere bundle is the equator region of a 3-sphere which carries a hyperbolic structure. By Mostow rigidity (see

Section 6), the volume of S2×[0,1]is also fixed, i.e., gravisolitons do not scale with the expansion of the universe. Then we obtain

ρsoliton=∑iNλEiV

where N is the number of gravisolitons in our cosmic horizon, λ is the linear density of gravisolitons, and Ei is the energy of each gravisoliton. Then we can get the dependence of ρon the scale a

ρsoliton(a)∼1a3=a−3

and then

3=3(1+wsoliton)⟹wsoliton=0

Then we have the following consequences:

(i)

The state equation of the gravisoliton (as dark matter) is p=0 or w=0.

(ii)

Gravisolitons interact only via gravity with matter by an attractive force.

(iii)

Gravisolitons have two possible scattering behaviors: two gravisolitons scatter so that it looks repulsively. In contrast, gravi- and anti-gravisolitons scatter so that it looks attractively.

In [

33], the state equation for dark matter was investigated by using the gravitational lensing. For two clusters, the Coma and the CL0024 cluster (see Figure 2 in [

33]), the variation of the equation of state parameter w over the radius of the cluster was studied. In contrast to the expected value at w=0, the both curves are located at the value w=−1/3 but the error bars are large enough to get w=0 also. Newer measurements [

34] got the value w=0. In [

35] the equation of state for dark matter was studied. In particular, the inverse cosmic volume law for Dark Matter was tested by allowing its equation of state to vary independently for different times values in the cosmic history (i.e., eight redshift bins in the range z=105 and z=0). Again, the equation of state w=0 was also confirmed during the cosmic evolution. According to consequence (i) above, our results agree with these measurements. Therefore, the detection of gravisolitons is rather complicate.

However, what is different for gravisolitons? From the mathematics point of view, the corresponding foliation leading to gravisolitons (see

Section 3) admits torsion. Here we expect that the polarization properties of the light after the gravitational lensing must be different. We conjecture that the light must be circularly polarized. However, what about the strong field limit? One of the main results is the close relation between the smoothness properties and the existence of gravisolitons. As discussed in our previous work, the space admits a fractal structure. Gravisolitons will follow this structure and the distribution of dark matter looks like a fractal which will be shown in a forthcoming paper.

**Author Contributions**
Conceptualization, T.A.-M. and J.K.; methodology, T.A.-M.; validation, T.A.-M. and J.K.; formal analysis, T.A.-M.; investigation, T.A.-M. and J.K.; writing—original draft preparation, T.A.-M.; writing—review and editing, T.A.-M. and J.K. All authors have read and agreed to the published version of the manuscript.

**Funding**
This research received no external funding.

**Acknowledgments**
We thank the anonymous referees for helpful remarks. In the course of the revision we found an error which we correct now. Furthermore, we acknowledge the critical remarks and questions which increased the readability of this paper.

**Conflicts of Interest**
The authors declare no conflict of interest.

**References**
- Zwicky, F. Die Rotverschiebung Von Extragalaktischen Nebeln. Helv. Phys. Acta 1933, 6, 110–127. [Google Scholar]
- Zwicky, F. On the Masses of Nebulae and of Clusters of Nebulae. Astrophys. J. 1937, 86, 217–246. [Google Scholar] [CrossRef]
- Rubin, V.; Thonnard, W.K.J.; Ford, N. Rotational Properties of 21 Sc Galaxies with a Large Range of Luminosities and Radii from NGC 4605 (R = 4 kpc) to UGC 2885 (R = 122 kpc). Astrophys. J. 1980, 238, 471–487. [Google Scholar] [CrossRef]
- Komatsu, E.; Smith, K.M.; Dunkley, J.; Bennett, C.L.; Gold, B.; Hinshaw, G.; Jarosik, N.; Larson, D.; Nolta, M.R.; Page, L.; et al. Seven-year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Interpretation. Astrophys. J. Suppl. 2011, 192, 18. [Google Scholar] [CrossRef]
- Ade, P.A.R.; Aghanim, N.; Armitage-Caplan, C.; Arnaud, M.; Ashdown, M.; Atrio-Barandela, F.; Aumont, J.; Baccigalupi, C.; Banday, A.J.; Barreiro, R.B.; et al. Planck 2013 Results. XVI. Cosmological Parameters. Astron. Astrophys. 2014, 571, A16. [Google Scholar]
- Metz, M.; Kroupa, P.; Theis, C.; Hensler, G.; Jerjen, H. Did the Milkay Way Dwarf Satellites Enter the Halo as a Group? Astrophys. J. 2009, 697, 269–274. [Google Scholar] [CrossRef]
- Moni Bidin, C.; AndCarraro, G.; Méndez, R.; Smith, R. Kinematical and Chemical Vertical Structure of the Galactic Thick Disk II. A Lack of Dark Matter in the Solar Neighborhood. Astrophys. J. 2012, 751, 30. [Google Scholar] [CrossRef]
- Pawlowski, M.; Pflamm-Altenburg, J.; Kroupa, P. The VPOS: A Vast Polar Structure of Satellite Galaxies, Globular Clusters and Streams Around the Milky Way. Mon. Not. R. Astron. Soc. 2012, 423, 1109–1126. [Google Scholar] [CrossRef]
- Thurston, W. Noncobordant Foliations of S3. Bull. Am. Math. Soc. 1972, 78, 511–514. [Google Scholar] [CrossRef]
- Belinski, V.; Verdaguer, E. Gravitational Solitons; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
- Asselmeyer-Maluga, T. Smooth Quantum Gravity: Exotic Smoothness and Quantum Gravity. In At the Frontiers of Spacetime: Scalar-Tensor Theory, Bell’s Inequality, Mach’s Principle, Exotic Smoothness; Asselmeyer-Maluga, T., Ed.; Springer: Stockholm, Switzerland, 2016. [Google Scholar]
- Reinhart, B.; Wood, J. A Metric Formula for the Godbillon-Vey Invariant for Foliations. Proc. AMS 1973, 38, 427–430. [Google Scholar] [CrossRef]
- Milnor, J.; Stasheff, J. Characteristic Classes; Ann. Math. Studies, 76; Princeton Univ. Press: Princeton, NJ, USA, 1974. [Google Scholar]
- Tamura, I. Topology of Foliations: An Introduction; Translations of Math. Monographs Vol. 97; AMS: Providence, RI, USA, 1992. [Google Scholar]
- Asselmeyer-Maluga, T.; Brans, C. How to Include Fermions Into General Relativity by Exotic Smoothness. Gen. Relativ. Grav. 2015, 47, 30. [Google Scholar] [CrossRef]
- Friedrich, T. On the Spinor Representation of Surfaces in Euclidean 3-Space. J. Geom. Phys. 1998, 28, 143–157. [Google Scholar] [CrossRef]
- Gibbons, G.; Hawking, S. Action Integrals and Partition Functions in Quantum Gravity. Phys. Rev. D 1977, 15, 2752–2756. [Google Scholar] [CrossRef]
- Asselmeyer-Maluga, T.; Rosé, H. On the Geometrization of Matter by Exotic Smoothness. Gen. Relativ. Grav.2012, 44, 2825–2856. [Google Scholar] [CrossRef]
- Asselmeyer-Maluga, T. Braids, 3-Manifolds, Elementary Particles: Number Theory and Symmetry in Particle Physics. Symmetry 2019, 11, 1298. [Google Scholar] [CrossRef]
- Morgan, J.; Szabo, Z.; Taubes, C. A Product formula for the Seiberg-Witten invariants and the Generalized Thom Conjecture. J. Diff. Geom. 1996, 44, 706–788. [Google Scholar] [CrossRef]
- Morgan, J.; Szabo, Z.; Taubes, C. Product formulas along T3 for Seiberg-Witten invariants. Math. Res. Lett. 1997, 4, 915–929. [Google Scholar] [CrossRef]
- Asselmeyer-Maluga, T.; Brans, C. Exotic Smoothness and Physics; WorldScientific Publ.: Singapore, 2007. [Google Scholar]
- Gibbons, G.; Werner, M. Applications of the Gauss-Bonnet theorem to gravitational lensing. Class. Quant. Grav.2008, 25, 235009. [Google Scholar] [CrossRef]
- LeBrun, C. Four-Manifolds Without Einstein Metrics. Math. Res. Lett. 1996, 3, 133–147. [Google Scholar] [CrossRef]
- LeBrun, C. Weyl Curvature, Einstein Metrics, and Seiberg-Witten Theory. Math. Res. Lett. 1998, 5, 423–438. [Google Scholar] [CrossRef]
- Asselmeyer-Maluga, T.; Król, J. Inflation and Topological Phase Transition Driven by Exotic Smoothness. Adv. High Energy Phys. 2014, 867460. [Google Scholar] [CrossRef]
- Asselmeyer-Maluga, T.; Krol, J. A topological model for inflation. arXiv 2018, arXiv:1812.08158. [Google Scholar]
- Asselmeyer-Maluga, T. Hyperbolic Groups, 4-Manifolds and Quantum Gravity. J. Phys. Conf. Ser. 2019, 1194, 012009. [Google Scholar] [CrossRef]
- Mostow, G. Quasi-conformal mappings in n-space and the rigidity of hyperbolic space forms. Publ. Math. IHES1968, 34, 53–104. [Google Scholar] [CrossRef]
- Prasad, G. Strong Rigidity of Q-Rank 1 Lattices. Inv. Math. 1973, 21, 255–286. [Google Scholar] [CrossRef]
- Asselmeyer-Maluga, T.; Król, J. Abelian Gerbes, Generalized Geometries and Foliations of Small Exotic R4. arXiv2014, arXiv:0904.1276. [Google Scholar]
- Hurder, S.; Katok, A. Secondary Classes and Trasnverse Measure Theory of a Foliation. Bull. Am. Math. Soc.1984, 11, 347–349. [Google Scholar] [CrossRef]
- Serra, A.L.; Dominguez Romero, M.J.L. Measuring the Dark Matter Equation of State. Mon. Not. R. Astron. Soc.2011, 415, L74–L77. [Google Scholar] [CrossRef]
- Sartoris, B.E. CLASH-VLT: Constraints on the Dark Matter Equation of State from Accurate Measurements of Galaxy Cluster Mass Profiles. Astrophys. J. Lett. 2014, 783, L11. [Google Scholar] [CrossRef]
- Kopp, M.; Skordis, C.; Thomas, D.; Ilić, S. Dark Matter Equation of State through Cosmic History. Phys. Rev. Lett.2018, 120, 221102. [Google Scholar] [CrossRef]

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

© 2020 by the authors. 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/).