Careful measurements have proven that there must be attractive forces in the galaxies that cannot be explained by the action of mass of the stellar objects. It has been suggested that these forces emanate from matter which does not emit light and therefore one has called the source of these forces as dark matter. Despite 20 years of searching, one has not found matter, and other explanations have not been satisfactory.
We can imagine that the forces do not come from localizable objects, but stem from the geometric structure of space. This concept is in perfect harmony with the principles of general relativity. Gravity, tidal force, and also rotational forces are usually described by the properties of space. In addition, pressure and mass density are derived from geometrical quantities.
Our new approach is based on mathematical methods already hidden in Schwarzschild’s paper on the interior solution (1916). We have worked out these methods in  and   and have also successfully applied them to the gravitational collapse  of the interior Schwarzschild solution. Now we want to apply the procedure to a cosmological model. We call our model P-model. This model and the Schwarzschild model both are based on the static de Sitter metric as seed metric. In the first case, this metric is a metric on a pseudo-hyper sphere, in the second case it is a metric on a cap of a hypersphere. From both metrics repulsive forces are derivable, which can be interpreted in the dS case as the cause of an expansion. We describe a process that switches these forces from repulsive to attractive. Although the elaboration of the P-model requires a considerable amount of mathematics, the underlying physical principles of the Schwarzschild world are well-known, but little respected in the literature.
We note that carrying out this procedure is only possible within the framework of the tetrad calculus. Coordinate systems play a minor role, they are only used for basic mathematical operations. In Sec. 2 we briefly outline the static de Sitter model, we explain the use of 4-bein systems, the Ricci-rotation coefficients, and the graded derivative. We discuss the representation of the Einstein field equations in tetrad form and point out that this is useful for many cosmological models and advisable for their clarity. In Sec. 3, we first discuss a static P-model in order to clarify the basic structure of such models. In Sec. 4 and Sec. 5 we develop the expanding P-model and discuss it in Sec. 6.
2. Preliminary Remarks
The dS model in its static version is based on a pseudo-hyper sphere with the radius , the polar coordinates , and the coordinate time t. The metric of the dS model in the canonical form is
We identify the quantities and as the recession velocity of the galaxies and the assigned Lorentz factor. From the metric (2.1) we read the 4-bein system
We use it to calculate the Ricci-rotation coefficients
which we decompose into the radial and the two lateral parts
By means of the unit vectors
we continue to decompose into vector quantities
Finally, we get from (2.3) the geometric quantities
At the same time, we have executed the operation
with m as a tetrad index and i as a coordinate index ( ). The lateral field quantities B and C describe the curvatures of the greater circles and the latitudinal circles of the hyper sphere. The quantity points away from any point on the pseudo-hyper sphere.
are the subequations of Einstein’s field equations. The first relation in (2.8) is the Friedman equation. is a submatrix of the tetrad metric . The use of the graded derivatives 
proves to be very advantageous for the representation of field equations. For the Ricci one has
The structures (2.8) and (2.10) can be used for spherically symmetric systems, static, expanding, or collapsing systems. If we reinterpret the cosmological constant in the field equations of the dS model, the stress-energy-momentum tensor has the form
by using (2.6)-(2.9). p is the pressure and the matter density of the cosmic fluid. Consequently, is the equation of state of the dS cosmos
In one of our papers  we have critically examined the dS cosmos concerning its interpretation. Our P-model is based on the dS model. The above structures are the source for further considerations.
3. The Static P-Model
As a first step towards a generalized model which admits attractive forces, we examine a static cosmological model by introducing the projector technology borrowed from the interior Schwarzschild solution. We do not attribute physical significance to the model. It aims to explain the use of the projectors which we introduce in order to transform repulsive forces into attractive forces. It is the starting point for a more sophisticated model.
We start from the previously described dS model. But now we put
with reference to (2.7). is a space-dependent function which we investigate as follows. We call it as a projector.
First, we recalculate the subequations of Einstein’s field equations. For the lateral subequations it is sufficient to calculate their timelike components. With and we get
and a similar equation for the quantity C. Finally one has
The radial field quantities (3.1) are treated according to
Therein the first line of (2.8) was used. To get a typical Friedman equation, we put zero the last brackets in the above expression and we get the important relation
and finally a relation whose structure is familiar to us from (2.8). We realize that we get the equations (2.8) of the dS metric for .
With (3.2) and (3.4) the Ricci tensor
can be calculated and also the Ricci scalar
and finally the Einstein tensor
Therefore the Einstein field equations result in
For one again has the pressure and the mass density (2.11) of the dS model. For later use we note
and we finally get the equation of state
which reduces for to the equation of state (2.12) of the dS model.
The conservation laws
for this static model are
and it has to be examined whether the pressure gradient leads to a contradiction.
If we differentiate (3.5)
we see immediately with (3.7) that (3.8) is satisfied. Thus, the introduction of the function has proved to be useful and opens the way to cosmological models in which the pressure is position-dependent.
4. The Expanding P-Model
In this section we want to show that the static P-model can be extended to an expanding one. We address the problem without reference to conventional cosmology. We avoid using an expanding metric and the scale factor. Nor do we refer to a comoving coordinate system, but we use orthogonal local reference systems. We start from the static de Sitter model. Its metric is the seed metric for our expanding model. In this case, the radius of curvature of the universe will be time-dependent. Then at any moment the expanding cosmos is a snapshot of the dS-cosmos. The coordinate system of the dS-cosmos is completely sufficient to carry out simple mathematical operations. We develop the main part of the theory in the tetrad calculus. This provides reference systems in which the quantities of the model have physically interpretable components. We use comoving and non-comoving reference systems. Between them a Lorentz transformation mediates, which includes as a parameter the rate of expansion of the cosmos.
In an earlier paper  we have extended the dS model by introducing a second velocity. We again apply this method, but we abandon the constraint which applies to the dS cosmos in order to obtain a genuine expanding model. The double-velocity ansatz is
Therein and is connected to the polar angle of a pseudo-hyper sphere. is the proper time of an observer who participates in the expansion of the cosmos. The pseudo-hyper sphere with the time-dependent radius forms the skeleton of the cosmological model. r, as in previous models, is the radial coordinate of the non-comoving observer. is, like in the ‘expanding’ dS cosmos, the speed of a fictitious freely falling observer. and are to be understood in analogy to the definition of , where is another time-dependent parameter. Thus, we can assume that in a fictive cosmos, which is preliminary to the physical cosmos, has a similar geometric meaning as has in the dS cosmos. The validity range of the quantities is
With the ansatz (4.1) we have established the method of double velocities. The Lorentz transformation
mediates between the observer system m’ comoving with the expansion of the cosmos and the non-comoving observer system m. v is the recession velocity and the associated Lorentz factor.
The velocities in (4.1) are supposed to satisfy Einstein’s addition law
The composition of the speeds can be taken from the drawing Figure 1.
Figure 1. Composition of the velocities.
The system m’’ is a fictitious reference system in a cosmos that can be represented by a pseudo-hyper sphere with the radius . At any given moment, a snapshot of the cosmos corresponds to the static dS cosmos. This enables us to use the coordinate system of the dS cosmos for basic mathematical operations. The physical system m’ is the reference system which is in the empirical space and follows an expansion that is less than that of the free-fall and m the non-comoving reference system. In order to describe all states of motion, we complement (4.2) with
Furthermore, the Lorentz relations
With the definitions in (4.1) one obtains for the recession velocity
From (4.1) it can be seen that accepts the value 1 at the equator of the pseudo-hyper sphere. That is the speed of light in the natural measuring system. In  we have shown that the velocity can only reach asymptotically the speed of light, i.e. it can come infinitely close to the value 1 only after an infinitely long time. Furthermore, we can see from Einstein’s addition law of velocities (4.3) that the recession velocity is lower than the one of free fall
Thus, we have envisaged a cosmological model whose rate of expansion is slower than the ones of the FRW models and whose rate of expansion can be manipulated with the parameter . Since the recession velocity does not reach the speed of light and definitely cannot exceed it, there is no violation of the special theory of relativity and no violation of causality nor is there any galactic island formation.
For the further development of the theory some basic relations are necessary. We are already familiar with
referring to the geometry of the hyper sphere and the Lorentz transformations (4.2) and (4.4).
In calculating the changes of the two velocities and , and finally of the recession velocity v, we must proceed with care. Both quantities contain, apart from the radial coordinate, the variable quantities and . Therefore we define the two quantities
With this and with (4.5) and (4.9) we get
From the above relations, it is obvious that the change of the two speeds has two causes. First, the speeds depend on the position. The recession velocity and thus the redshift of the light of a galaxy is higher, the farther the galaxy is dislocated from the observer. Secondly, the recession velocity depends on the expansion of the cosmos, i.e., on the quantities and . So we can expect accelerations of expansion, which is suggested by recent observations. We will take a closer look at this problem later.
Now the necessity arises to develop the properties of the quantities (4.10) from the geometric and the kinematic structure of the model. A look at Figure 1 shows us that for an observer with the velocity the system m’ is the rest frame. Therefore we demand that in analogy to the non-comoving observers of the dS cosmos 
applies. Thus, we expect that from the view of the physical observer m’ no change in time of the velocity can be experienced. With this restrictive condition one can calculate . This restriction simplifies the structure of the model considerably. From (4.12) and (4.11) and using the Lorenz relations we obtain
For the velocity , similar constraints would bring us back to the dS cosmos. In order to determine the quantity from the structure of the model, further considerations are necessary. First of all, it is evident that the comoving observers cannot ascertain any position dependency of . The model is homogeneous and its curvature is everywhere the same for . So we can put .
We do not derive our model by solving Einstein’s field equations under specific conditions, but we make assumptions about a particular geometry, we derive from them all necessary quantities. Then we check, whether they satisfy Einstein’s field equations. Since a homogeneous world model is based on a spherical space, some of the required quantities are already familiar to us. Thus, we know the expressions for the two lateral field quantities B and C in the system at rest
We do not derive the field quantities of the comoving system from an expanding metric, but with a Lorentz transformation from the metric of the static seed system. With
we obtain the quantities of the comoving system
From these relations we realize that the spatial components of these quantities do not appear to be flat. This was the case for the FWR models and the simple subluminal model, i.e. for all models that expand in free fall. For the special case we get , and thus
The spatial parts and of the above equation now appear to be flat and correspond to the quantities of a free-falling system, as one can easily convince oneself. As we have repeatedly stated, this does not mean that the space is flat, but according to Einstein’s elevator principle  , no forces act on the free-falling observers.
Let us take another look at the relations (2.4)-(2.7) taking into account (3.1). To complete the model we need to know the quantity . It is decisive for the development of the universe and enters into the Friedman equation.
The quantity , which is closely related to the quantity , can however be deduced without further calculation, if we demand that the expansion in all three spatial directions is equal. Thus, after looking at (4.14) one has
For further considerations, we turn to the conservation law. The stress-energy-momentum tensor of the model in the comoving system should have the form of the dS model, if we assume that the matter density has the usual form . For the pressure, we expect a more complicated expression, as the preliminary study of Sec. 3 shows.
From the conservation law of the moving system with and
using (4.15). On the other hand
follows from the relations . By comparison with (4.16) we get
and with (3.6)
Therein is the projector already used in the last section, which will be discussed in more detail below. Knowing according to (4.15), we have derived the general form of and at the same time we have allowed the model to expand according to the definition (4.10). Thus,
For we get back the dS model. Then one has and . For the two systems m’ and m one finally has
Now we are also able to specify completely the velocity changes. If we substitute the above relations into (4.11) we get
The 4th components of the velocity changes do not vanish. Thus, the P-model allows an acceleration of the cosmos.
So far we have developed some quantities necessary for the construction of the theory without much mathematical effort. We now want to put the model on a safe mathematical basis. To this end, we deal in detail with the transformation to the comoving and non-comoving observer systems and with the inhomogeneous transformation law of the Ricci-rotation coefficients. The field quantities are constituents of these Ricci-rotation coefficients, which transform from the non-comoving system m to the comoving system m’ according to
The last term in the above relation is called the Lorentz term. This formula which is fundamental in the tetrad calculus must be generalized to the P-cosmology so that the Ricci-rotation coefficients of the auxiliary model can be used to generate the Ricci-rotation coefficients A of the P-model. First of all, we again use the comoving system and now we define most generally the projector . This projector contains only the diagonal components
For the P-cosmology, we modify the inhomogeneous transformation law (4.21) to
With this we have first accessed the seed metric and we have transformed it into a comoving system. At the same time we have with the help of the projector created a new geometry.
However, it can be seen that this modified transformation law has no effect on the lateral field quantities B and C. They transform homogeneously and are not changed by the projector. They are given in the comoving system by (4.14). Thus, we can concentrate on the quantity U.
If the Ricci-rotation coefficients of the seed model are designated by in the comoving system, one first obtains
is contained. Finally one has
Since the Lorentz transformation is to be interpreted as a pseudo rotation, which takes place in the [1,4]-slice of the surface, the inhomogeneous term can be simplified with to
From (4.21) only Einstein’s elevator law in cosmology
With a look at (4.2) we calculate with (4.26)
and evaluate this with the relation
With (4.20), (4.12) we get
We define three quantities
where G describes the change of , l the change of , and f the change of . Finally, one has
and the Lorentz term for the inverse transformation
If we take (4.19) into account we can write
In (4.15) we have deduced that the quantity must have the same value as the quantities and in order to correctly represent the expansion scalar. Exactly this value would have to be calculated with a Lorentz transformation from the static auxiliary system. We therefore rely on the relation (4.27) and we calculate the desired quantity using together with
and indeed we get the expression predicted by (4.15).
Now the 1st component of the quantity has to be calculated. It is
Thus, we have completely calculated the quantity
Rather new to FRW models is the radial component in ‘U. Similar to the static dS cosmos a radial force occurs which can be attractive or repulsive depending on the sign of . We realize that with the projector one can manipulate the structure of the geometry. We will come back to this later.
After we have succeeded in deriving a complete set of field quantities from the geometry, we want to check whether these quantities satisfy the field equations and the conservation laws and which equation of state results for the cosmos.
It can be seen that we will come across a term in the progress of the calculations. The analogy to (3.3)
will prove. After some calculation we find the Friedman equation
quite analogous to (3.4). Now the lateral field equations have to be calculated. With the definition of the graded derivatives
we get the remaining subequations of Einstein’s field equations
which, according to the scheme (2.10), we compose to the Ricci tensor and to the Ricci scalar
The result is the stress-energy-momentum tensor
From the conservation law results for the pressure
From (4.40) we get
thus showing that the ansatz (4.34) is correct. From the conservation law also follows
With the value (4.40) for we have previously inferred the relation (4.17). The relation (3.6) is valid for the expanding model, as well as the equation of state (3.7) of the cosmic fluid
So far, we have worked very successfully with the function without specifying it. From (4.43) one recognizes that is closely connected with the equation of state of the cosmos. By restructuring one gets
In the literature the relation
is common. With this one can write
It can be seen that the function allows a class of different cosms, depending on the choice of the equation of state. We will discuss some specific cases later on.
In addition, we want to note that the P-geometry can also be derived from the freely falling system of the seed metric. One has for the inhomogeneous transformation law
With this and the Lorentz relations one can easily obtain the already known quantities of the comoving system m’. However, it must be noted that a return to the dS auxiliary model is possible only if the conditions for this model are satisfied, viz , .
5. The P-Model in the Non-Comoving System
The cosmologists endeavor to find field quantities and field equations also for the non-comoving system. Although we believe that a representation in the non-comoving system does not provide any new insights, we nevertheless show that such a representation is possible and that the P-model proves to be extremely robust to transformations.
First, the meaning of a non-comoving observer must be clarified. In the literature, one finds the demand , i.e. the fixing of the radial coordinate r of the higher-dimensional embedding space. However, it is obvious that such an observer does not comove with the expansion but still moves with respect to an adjacent observer, e.g. to an observer located on the arbitrarily chosen pole of the pseudo-hyper sphere.
For the distance of the non-comoving observer from the pole one has specifically
In an expanding cosmos its radius depends on time, and its polar angle also changes if one holds the radial coordinate r tight at a fixed position. Thus, the condition is quite artificial for an expanding cosmos.
In order to complement the P-model in the sense above, field quantities and field equations still have to be revealed in the non-comoving system. The Lorentz transformation is given by (4.2). The lateral field quantities transform homogeneously and can be taken from (2.7). However, more attention is given to the radial field quantities which transform inhomogeneously according to
we have already calculated in (4.30). One now has to perform
and with this after some calculation the Friedman equation in the non-comoving system
Comparison with (4.35) shows that the U-equation is form invariant. With (5.3) we recognize that the radial force U of the non-comoving system does not reduce to that of the static auxiliary model. That would only be the case if , i.e., if .
The lateral field equations again are formulated with the graded derivative. According to (5.3) these equations now contain terms with
From this one gains
In the non-comoving system, the stress-energy-momentum tensor has the form
Now all that one has to do is to show how the quite different expressions in (5.5) and (5.6) interrelate. First, we take from (4.45)
Furthermore, we get with (4.19) and (5.3)
Thus, it can be seen that the stress-energy-momentum tensor (5.6) takes the expected values in the non-comoving system.
6. Discussion of the Model
In the previous section we have set up a cosmological model with the help of a projector, which is formulated very generally. However, there is still the possibility of giving the model a specific appearance by specializing the projector. Now we assign some specific values to the projector. We recall that is a function and a specification can only be made after the model has been worked out.
According to (4.19) one obtains , , i.e. a model which is geometrically represented by a pseudo-hyper sphere of constant radius. In addition, if is one obtains the dS cosmos, i.e. the static model from which our considerations started.
For a comoving observer in the dS cosmos the polar angle changes on the pseudo-hyper sphere. There is a change in position of these observers by a geometrically induced drifting.
However, if one has , one obtains the extended dS cosmos, which we have discussed in detail in  . This model is also static, i.e. has constant radius of the pseudo-hypersphere. However, observer crowds drift apart at a speed less than that of free fall.
The motion of the galaxies in this model is accelerated. The Hubble law deviates from its linear form. This can affect the calculated redshift values. If we take the acceleration into account it may be possible to better adapt the calculated values to the measured ones.
Both models, with or without , have the same equation of state. As a result, no matter flow or energy transport can be detected in the comoving system. Both models have the speed of light as the highest recession velocity, which, however, can only be reached asymptotically. After an infinitely long time, the recession velocity approaches infinitely close the speed of light. The location that can only be reached by the galaxies asymptotically, i.e. after an infinite time, is the cosmic horizon. This is the equatorial surface of the hyper sphere, related to the pole of which the observer of the galaxies is located. Any observer can define himself as a pole on the hyper sphere and has his own individual cosmic horizon.
Both models have repulsive forces responsible for the galaxy drift. These forces are directly derived from the dS ansatz.
With one obtains the expanding subluminal cosmos  . Repulsive forces act in the system in rest of this model, the comoving system is free of forces. The subluminal cosmos expands in free fall; the Einstein elevator principle applies  . In the non-comoving system the metric has a form from which can be read the curvature parameter . This means that the space is positively curved and finite. After a transformation into the comoving system, the metric has the form . This does not mean that the space is flat now. It rather means that the free-falling observers do not experience forces, they recognize the space to be locally flat, in the same way as the free-falling observers do in the Schwarzschild field.
In this model, the polar angle on the pseudo-hyper sphere does not change in time for the comoving observer. With one finds
According to  is in this case and one finally has for the recession velocity
The subluminal model is closely related to the model by Melia. Melia’s calculated redshift values fit the observed values much better than those of FRW models. Therefore, our subluminal model is very close to Nature.
With results an extended subluminal model which does not prove to be particularly useful.
The pressure in this model is positive. By some authors the equation of state is attributed to an ultra-relativistic gas. What is interesting about this model is that in the comoving system, the system in which we live, an attractive force
occurs. It corresponds to the repulsive force of the extended dS model. The corresponding 4th components of the field quantities are identical in both models.
This new force depends on the quantity . If one puts this quantity zero, this force vanishes. Then the cosmos expands in free fall. The new attractive force is based on two geometric features of the model: the implementation of a second speed that reduces the recession velocity, and the operation of the projector. The latter changes the sign in (4.33). The attractive force stems from geometry and may explain the dark matter without our having to assume new types of matter as the source of this force. Thus, dark matter would be a property of space.
In the general P-model the polar angle changes on the pseudo-hyper sphere for the comoving observer. If one recalls the relation and one now obtains with the relations
The change of the polar angle has two causes: the geometric structure of the model with the forces that occur, and the expansion of the cosmos. With (4.18) the above expression is reduced to
Depending on the choice of the cases discussed above are included in this relation. For the change of has as the cause the geometric structure on the pseudo-hyper sphere and for no temporal change of takes place.
Let us return to the problem of the recession velocity. Combining (4.11) with (4.19) one has
From this it is clear that only for , i.e. for the subluminal model, results. Since, by definition, is zero, the recession velocity is also temporally constant. For other values for acceleration occurs. It turns out that both the generalized dS model and the P-model are exact solutions to Einstein’s field equations that allow accelerations of the recession velocity. Furthermore, one gets for and from the relation and thus also in accordance with the equation of state .
By recent measurements of the redshift of galaxies one has found that the expansion of the universe is likely to be accelerated. We therefore want to work out how the acceleration for individual models in different frames of reference is obtainable. In the textbooks it is mentioned that the acceleration is not covariant, i.e. it does not transform like vectors between differently moving frames of reference. The problem is rarely discussed in more detail. But we can say that we have already outlined the problem, although it has not been explicitly stated. Accelerations are identical to forces on the unit mass. We have already discussed in detail how the forces behave in the various frames of reference. The forces are the field quantities and are components of the Ricci-rotation coefficients. The accelerations thus transform according to the inhomogeneous transformation law of the Ricci-rotation coefficients.
we go into the equations of motion for the comoving observer , represented in the frame of reference m
and thus we get the inhomogeneous transformation law of the Ricci-rotation coefficients. In a more conventional way the above relation is written as
The total acceleration of a galaxy relative to the non-comoving system is composed of a kinematic component resulting from the motion of the reference systems relative to each other
and a geometric component originating from the basic structure of the model. The latter can be positive, negative, or zero, depending on the value assigned to the projector. ‘L is the Lorentz term of the inhomogeneous transformation law of the Ricci-rotation coefficients, which we have got to know in (4.29).
In addition, if one multiplies the last but one equation with the rest mass of the galactic matter, one finally has
The kinematic expression contained therein corresponds to the definition of the Lorentz force in the theory of electrodynamics. Finally, one has
in referring to classical mechanics. The geometric term and the kinematic term have, according to Sec. 4, the form
The total force is repulsive for , i.e. for the extended dS model, but the kinematic term has the opposite sign and decreases the repulsive effect. For the simple dS model with , the two terms cancel each other. The comoving observers are in free fall. For , for the subluminal model, is . Thus, no forces act on the comoving observers. For both forces have an attractive effect and add up to
and suggest an effect that can be attributed to dark matter.
From the 4th component of the equation of motion one obtains with the intermediate result
and after multiplication with the proper mass
one has calculated the work that must be done in transporting the galaxies in the time interval .
A strict geometric treatment of cosmology allows for the ability to modify the underlying geometry to reveal certain physical properties of the cosmos. Thus, the calculated recession velocity of the galaxies can be reduced, enforcing an acceleration of these and implementing attractive forces that can be attributed to dark matter. Such forces can also be responsible for ensuring that the internal structures of the galaxies are stable and not subject to the expansion of the universe.
Although the model proposed here is not very close to Nature, it contains mathematical methods that prove to be useful in the construction of cosmological models.
Because of the homogeneity of the universe, the new attractive forces are active throughout the universe and are not limited to galaxies, although this would be desirable. Here one could pursue an idea that was initiated by Einstein and Strauss  and was taken up by some other authors: The insertion of matter condensates into the homogeneous cosmic fluid.
In additions, the approach with seems to be too simple. A look at the interior Schwarzschild solution, from which we have borrowed the projector technology, shows that only a more complicated function for leads to a functioning model. Nevertheless, it is encouraging that three known cosmological models can be obtained from the P-model, i.e. the dS model, the extended dS model, and the subluminal model. The P-model leaves open the possibility of designing the projector so that the model takes up new structures that may possibly explain the properties of our world.