Abstract
We investigate the existence of inhomogeneous exact solutions in Weyl Integrable theory with a matter source. In particular, we consider the existence of a dust fluid source while for the underlying geometry we assume a line element which belongs to the family of silent universes. We solve explicitly the field equations and we find the Szekeres spacetimes in Weyl Integrable theory. We show that only the isotropic family can describe inhomogeneous solutions where the LTB spacetimes are included. A detailed analysis of the dynamics of the field equations is given where the past and future attractors are determined. It is interesting that the Kasner spacetimes can be seen as past attractors for the gravitation models, while the unique future attractor describes the Milne universe similar with the behaviour of the gravitational model in the case of General Relativity.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
Analytical and exact solutions play a significant role in the study of gravitational physics. The existence of exact spacetimes is essential in order to understand the physical properties and the nature of the physical space. Inhomogeneous and anisotropic exact spacetimes that have zero magnetic Weyl tensor are very useful in gravitation and cosmology. They include an important family of spacetimes known as the Szekeres universes. The Szekeres spacetimes are the most general cosmological exact solutions of general relativity with a pressureless fluid source [1, 2]. They possess no symmetries but the spatial three-slices have a special geometrical structure. In the Szekeres spacetimes, information does not propagate via gravitational or sound waves, so, they are also known as ’silent’ universes [5].
Szekeres spacetimes are inhomogeneous universes that do not admit any vector field isometry. Moreover, the rotation and acceleration of the fluid source must be identically zero, and the pressure constant. In practice, this means the only inhomogeneous matter sources allowed are dust, with or without a cosmological constant. While, in general, the spacetimes are anisotropic – which means that the shear is non-zero and the expansion rate is non-zero. The inhomogeneous Szekeres spacetimes are classified into two families: the inhomogeneous Kantowski–Sachs (-like) spacetimes and the inhomogeneous FLRW (-like) spacetimes.
There are applications of Szekeres spacetimes in gravitational physics and cosmology [1, 6]. A complete description of the scalar polynomial curvature singularities in both classes of Szekeres solution have been established, and they are velocity-dominated. They have Newtonian counterparts and contain no gravitational waves [4, 10]. In addition, the asymptotic behavior in the distant future has been analyzed [3, 7].
A more general gravitational collapse, known as quasi-spherical by using the Szekeres spacetimes was studied in Ref. [8], where it was found that a strong radial increase in the density, the fluid heralds the onset of a naked singularity. The matter distribution in Szekeres spaces has a dipolar character [9], while there is no gravitational radiation emission from the inhomogeneous moving dust [10], for other applications of Szekeres spacetimes in gravitational physics we refer the reader to [11,12,13,14] and references therein. Tilted Szekeres models were studied in Ref. [15] where it was found that vorticity follows the congruence of the fluid world lines. Recently, the frame rotation of the Szekeres spacetimes which relates the cosmological solutions with the quasi-spherical exact solutions was studied in Ref. [16].
The quasi-spherical Szekeres dust solutions are a generalization of the spherically symmetric Lemaître–Tolman–Bondi dust models where the spherical shells of constant mass are not concentric. A coordinate-independent analysis of the dynamics of the spherically symmetric Lemaître–Tolman-Bondi cosmologies, emphasizing their relation to the Friedmann Lemaître cosmologies was given in [17]. In general, it was shown that ever-expanding Lemaître–Tolman–Bondi cosmologies isotropize at late times, approaching the de Sitter universe, or the Milne universe, depending on whether or not a cosmological constant is present. For the analysis, a dimensionless scalar is introduced to represent the ratio of the Weyl and Ricci curvatures. In all cases, there is a finite limit at late times, its value determines the asymptotic spatial inhomogeneities in various physical quantities. The Lemaître–Tolman–Bondi cosmologies for which the initial singularity is isotropic were also identified. The collapsing quasi-spherical Szekeres dust solution, where an apparent horizon covers all shell-crossings that will occur, can be considered as a model for the formation of a black hole. The apparent horizon can be detected by a Cartan invariant [18]. In the former reference, solutions of this sort are reviewed together with their spin coefficients and curvature scalars in the Newman–Penrose formalism. The Cartan–Karlhede algorithm is used to generate the minimal set of extended Cartan invariants. Cartan scalars are compared with the kinematic scalars [19] and q-scalars [20], which are two well-known sets of scalars used to characterize Szekeres solutions.
Inhomogeneous spacetimes can been seen as limits of FLRW spacetimes with inhomogeneous perturbations, such a comparison between non-spherical Szekeres spaces and the dynamics of cosmological perturbation theory was performed in Ref. [21]. Specifically, it was proved that the linearised Szekeres evolution equations and their solutions fully coincide with the corresponding equations of the linear cosmological perturbation theory and their solutions in the isochronous comoving gauge. Moreover, the conservation of the curvature perturbation holds for the appropriate linear approximation of the exact Szekeres fluctuations in \(\Lambda \)-cosmology, while the different collapse morphologies of Szekeres models yield different growth factors to those that follow from the analysis of redshift space distortions.
There are various generalizations of the Szekeres exact solutions where additional matter sources contribute to the gravitating matter [2]. Indeed, the first generalization presented by Szafron in Ref. [22], where the dust fluid source was replaced by a perfect fluid with non-zero pressure, leading to the Szekeres–Szafron spacetimes. The cosmological constant term was introduced by Barrow et al. [23] where the inhomogeneous analogue of the \(\Lambda \)CDM model was derived. Other kinds of matter source have been introduced, such as heat flow, electromagnetic field, viscosity, and an aether field in the context of Einstein-aether theory [24,25,26,27,28,29,30,31].
In this work, we are interested in determining the exact inhomogeneous spacetimes in Weyl Integrable theory [32]. A Weyl manifold is a conformal manifold equipped with a connection which preserves the conformal structure and is torsion-free. In Weyl Integrable theory the connection structure is related to the Levi-Civita connection, to which it differs by a scalar field of the conformal metric. Specifically, if \(g_{\mu \nu }\) is a metric tensor with Levi-Civita connection \(\Gamma _{\mu \nu }^{\kappa }\), then in Weyl Integrable theory the manifold is supported by the set \(\left\{ g_{\mu \nu },\tilde{\Gamma }_{\mu \nu }^{\kappa }\right\} \) where \({\tilde{\Gamma }}_{\mu \nu }^{\kappa }\) is the Levi-Civita connection for the conformally related metric \(\tilde{g}_{\mu \nu }=\phi g_{\mu \nu }\), where \(\phi \) is a scalar field. An important characteristic of the Weyl Integrable theory is that it is in agreement with current astronomical and other observations [33].
Physical consequences of Weyl invariant theories are discussed, e.g., in Refs. [34,35,36]. In Ref. [37], it is discussed whether or not a general Weyl structure is a suitable mathematical model of spacetime. In this regard, it was found that a Weyl integrable spacetime is the most general structure suitable to model spacetime. The well-posedness of the Cauchy problem for particular kinds of geometric scalar-tensor theories of gravity, which are based on a Weyl integrable spacetime, is given in Ref. [38]. In Ref. [39], a formulation of general relativity on a Weyl-integrable geometry which contains cosmological solutions, exhibiting acceleration in the present cosmic expansion, is studied. The conditions for accelerated expansion of the universe are derived there. A particular solution for the Weyl scalar field describing a cosmological model for the present time is obtained in concordance with the data-combination Planck + WP + BAO + SN. In Ref. [40], the evolution of 4-, 5- and 6-dimensional cosmological models based on the integrable Weyl geometry are considered numerically both for empty spacetime and for scalar field with non-minimal coupling with gravity. In Ref. [41], the motion of massless particles on the background of a toroidal topological black hole is analyzed in the context of conformal Weyl gravity. Null geodesics are found analytically in terms of the Jacobi elliptic functions.
There are various exact solutions of the field equations in Weyl Integrable theory. Vacuum cosmological models were studied in Ref. [44], while higher- or lower-dimensional gravitational models were studied in Refs. [40, 42, 43]. In Ref. [45], the authors studied gravitational models in Weyl Integrable theory with matter source, an electromagnetic field, and additional scalar field. In these models an interaction between the scalar field of the Weyl theory and the matter sources is introduced, when the field equations are written in the covariant form using the tensor quantities of general relativity. Inhomogeneous models in Weyl Integrable models were also studied in Refs. [46,47,48,49], while spherically symmetric solutions can be found in Ref. [50]. In addition in Ref. [50] the authors discuss the similarities and the differences of the Weyl Integrable theory with the Brans–Dicke theory.
The plan of this paper is as follows. In Sect. 2 we present the basic properties and definitions of the gravitational field equations in Weyl Integrable theory. We rewrite the field equations in a way that is equivalent in form to general relativity and show that a scalar field is introduced in the field equations, and we discuss the case where an ideal gas contributes in the gravitational model. For the underlying geometry, we consider the line element which belong to the silent universe class, and describes the Szekeres spacetimes in general relativity. Exact solutions of the field equations are presented in Sect. 4. In Sects. 4 we perform a detailed analysis of the field equations in order to understand the past and future evolution of the cosmological solutions. Finally, in Sect. 5 we discuss our results and draw our conclusions.
2 Weyl integrable gravity
Weyl geometry is an extension of Riemannian geometry, specified by a metric tensor \(g_{\mu \nu }\) and a gauge vector field \(\omega _{\mu }\). The covariant derivative \({\tilde{\nabla }}_{\mu }\) is defined by the (Weyl) affine connection \({\tilde{\Gamma }}_{\mu \nu }^{\kappa },\) with the property,
from which we infer that the gauge vector \(\omega _{\mu }~\) plays a significant role in the geometry. Specifically, the Weyl affine connection \(\tilde{\Gamma }_{\mu \nu }^{\kappa }\) is related to the Christoffel symbols \(\Gamma _{\mu \nu }^{\kappa }\) of the metric tensor \(g_{\mu \nu }\) as follows:
The curvature tensor in the Weyl geometry is defined as
where, by using (2), we observe that in general \({\tilde{R}} _{\kappa \lambda \mu \nu }\) is not symmetric as it is in the case of Riemannian geometry.
In this work, we are interested in the case where \({\tilde{R}}_{\kappa \lambda \mu \nu }\) has the same symmetric properties as in Riemannian geometry. This is true when \(\omega _{\mu }\) is a gradient vector, which means that there exists a scalar \(\phi \) such that \(\omega _{\mu }=\phi _{,\mu }\). In addition, in this case, length variations are integrable along a closed path. This specific theory is known as Weyl Integrable geometry. Moreover, there exists a conformal map which relates the metric tensor \(g_{\mu \nu }\) of a Riemannian space into that of Weyl integrable space, which means that a Weyl integrable space is also conformally a Riemann space.
In Weyl integrable geometry, the Ricci tensor \({\tilde{R}}_{\mu \nu }\) is related to the Riemannian Ricci tensor \(R_{\mu \nu }\) by,
where the Ricci scalar, in a four-dimensional manifold, is written as [45]
2.1 Gravitational action integral
We define the simple gravitational action Integral, which includes the Weyl Ricci scalar \({\tilde{R}}\) and the field \(\phi _{;\mu }\), as
where \(\xi \) is an arbitrary coupling constant. At this point, we remark that
Variation with respect to the metric tensor of the action integral, \(S_{W}\), provides the gravitational field equations [45],
where \({\tilde{G}}_{\mu \nu }\) is the Weyl Einstein tensor. Moreover, variation with respect to the scalar field \(\phi \) gives
that is, a Klein–Gordon equation of the form,
where \(\nabla _{\mu }\) denotes the Riemannian covariant derivative.
The gravitational field equations (8) can be rewritten by using the Riemannian Einstein tensor \(G_{\mu \nu }\) as follows
in which the new constant \(\lambda \) is defined as \(2\lambda \equiv 4\xi -3\).
Consequently, the field equations (11) for \(\lambda >0\) are those of general relativityFootnote 1 with a massless scalar field. A new possibility is introduced when \(\lambda <0,\) which correspond to the addition of a massless phantom scalar field.
Until now, we have considered the case of vacuum. Now, we present the field equations in the presence of a matter source. Specifically, we consider the cases where a pressureless dust fluid source contributes to the gravitational field equations.
2.2 The presence of dust
When a pressureless fluid dust source is included, the gravitational field equations become
where \(T_{\mu \nu }^{\left( m\right) }=e^{-\frac{\phi }{2}}\rho _{m}u_{\mu }u_{\nu }\), while the Klein–Gordon equation (10) becomes
while the conservation equation for the matter field reads \(\tilde{\nabla }_{\nu }T^{\left( m\right) \mu \nu }=0\).
We can see that there is a coupling between the scalar field and the dust fluid; hence, the coupling and effective pressure term can depend on the energy density \(\rho _{m}\). The interaction between scalar field and dust fluid has been proposed as a potential mechanism to explain the cosmic coincidence problem [51,52,53,54]. Various interaction models have been studied before in the literature; for instance, see [55,56,57,58,59], and references therein.
3 Inhomogeneous spacetimes
The gravitational model that we have considered in Weyl Integrable geometry is equivalent to that of general relativity with an effective energy-momentum tensor, where the field equations are of the form
and \(T_{\mu \nu }~\)is the effective energy-momentum tensor. It consists of a massless scalar field \(\phi \), and an additional fluid source interacting with the field \(\phi \). In particular, \(T_{\mu \nu }=T_{\mu \nu }^{\left( \phi \right) }+{\hat{T}}_{\mu \nu }\), where \({\hat{T}}_{\mu \nu }\) describes the energy-momentum tensor of a pressureless fluid, i.e. \({\hat{T}}_{\mu \nu }=T_{\mu \nu }^{\left( m\right) },\) or of the second scalar field, \(\psi \); that is,
where \(T_{\mu \nu }^{\left( \phi \right) }\) is the energy-momentum tensor of the massless scalar field,
However, as we discussed before the continuity equation \(\nabla _{\nu }T^{\mu \nu }=0\), provides \(\nabla _{\nu }\left( T^{\left( \phi \right) \mu \nu }+\hat{T}^{\mu \nu }\right) =0\), that is, \(\nabla _{\nu }\left( T^{\left( \phi \right) \mu \nu }\right) =Q,~\nabla _{\nu }\left( {\hat{T}}^{\mu \nu }\right) =-Q\), where \(Q=Q\left( x^{\mu }\right) \), is the interacting term.
In this work, we assume that the underlying spacetime is described by the inhomogeneous and anisotropic diagonal line element
in which \(A=A\left( t,r,y,z\right) \) and \(B=B\left( t,r,y,z\right) \). The functional forms of the two scale factors A, B are determined by the solution of the field equations (14).
In the case of general relativity, for \(B_{,r}=0\), i.e. \(B=B\left( t,y,z\right) \), the exact solutions belong to the Kantowski–Sachs (-like) family, while for \(B_{,r}\ne 0\) the resulting spacetimes are inhomogeneous and isotropic. Generalization of the Szekeres spacetimes with a purely time-dependent scalar field have been studied before in Ref. [60]. In particular, a quintessence scalar field was considered and the scalar field should be homogeneous. Although the density can be inhomogeneous in metric (17), the pressure must be homogeneous [2].
Hence, from the results for the vacuum solution of [60] when the scalar field is massless, we recover the exact solution of the vacuum Weyl Integrable geometry for \(\lambda >0\). However, in the presence of an additional matter source, as we have here because of the existence of the interaction term, the analytic solutions will be different.
In the following section, we proceed with the presentation of the analytic solutions for the field equations (14), where the underlying spacetime is described by the line element (17).
We require the pressure term of the effective energy-momentum tensor \(T_{\mu \nu }\) in (14) to be homogeneous such that the FLRW limit to be provided. Thus, \(\phi =\phi \left( t\right) \) while for the matter source we have \(\rho _{m}=\rho _{m}\left( t\right) \). The latter follows easily, if we rewrite the energy-momentum tensor \(T_{\mu \nu }\) such that to define a new pressure component in order to eliminate the interaction term. The steps that we follow to solve the field equations are similar to those taken in Ref. [60]. Thus, we omit the presentation and go directly to the main results.
Similarly, in the case of the homogeneous scalar field, we find that the Szekeres-like solutions in the Weyl Integrable theory are classified into two classes of solutions, (A) the inhomogeneous Kantowski–Sachs family of solutions and the (B) inhomogeneous FLRW (-like) solutions.
For the Kantowski–Sachs family of solutions the unknown functions in the line element (17) are \(A\left( t,r,y,z\right) =\alpha \left( t\right) \) and \(B\left( t,r,y,z\right) =\beta \left( t\right) \left( c_{1} uv+c_{2}u+c_{3}v+c_{4}\right) ,\) where the \(y=u+v,~z=i\left( u-v\right) \), so that the line element is written as [60]:
and the curvature, K, of the two-dimensional surface of constant curvature \(\left\{ y-z\right\} \) to be related with the constants \(c_{1},\) \(c_{2},~c_{3}\) and \(c_{4}\) as follows, \(K=c_{1}c_{4}-c_{2}c_{3}\). The unknown time-dependent functions \(\alpha \left( t\right) ,\beta \left( t\right) \) are determined by a set of differential equations that will be presented in the following sections.
The second family of solutions which correspond to the inhomogeneous FLRW-like spacetimes are described by the line element [60]:
The function \(C\left( r,y,z\right) \) is now given by \(C\left( y{,}z\right) {=}{-}2\ln ( \gamma _{1}\left( r\right) \) \(uv+\gamma _{2}\left( r\right) \) \(u+\gamma _{3}\left( r\right) v+\gamma _{4}\left( r\right) )\) , where the functions \(\gamma _{1}\left( r\right) ,~\gamma _{2}\left( r\right) ,~\gamma _{3}\left( r\right) ,\) and \(\gamma _{4}\left( r\right) \) are constrained by \(k=\gamma _{1}\left( r\right) \gamma _{4}\left( r\right) -\gamma _{2}\left( r\right) \gamma _{4}\left( r\right) \), where k is the spatial curvature of the FLRW-like spacetime. The scale factor \(\alpha \left( t\right) \) is given by the generalized Friedmann equations in Weyl Integrable geometry given below.
3.1 Kantowski–Sachs spacetimes
The unknown scale factors of the Kantowski-Sachs spacetime (18) are given by the following system,
while the equation of motion for the scalar field and the matter source are given by,
where overdot means total derivative with respect to the variable t.
3.2 FLRW spacetimes
Analogously, the unique scale factor for the FLRW (-like) spacetime (19) is given by the (modified) Friedmann equations
and the scalar field \(\phi \) satisfies the Klein–Gordon equation,
while the conservation equation for the dust fluid source is
At this point, we remark that for \(\lambda =0\), only the vacuum solutions of general relativity are recovered, while the Szekeres spacetimes are recovered when \(\phi =\phi _{0}\) and \(\lambda \rightarrow \infty \). This is reminiscent of the range of the constant Brans–Dicke parameter, \(\omega \), in scalar-tensor theory such, where the limit of general relativity to be recovered as \(\omega \rightarrow \infty \) [61]. This family of spacetimes includes also the inhomogeneous Lemaître–Tolman–Bondi (LTB) spacetimes [62].
In the following we show the analytic solution for the inhomogeneous FLRW (-like) spacetime.
3.3 Inhomogeneous analytic solution
Now let us consider the case when the spatial curvature is zero, i.e. \(k=0\). The gravitational field equations can be rewritten in an equivalent form,
where \(H=\frac{{\dot{a}}}{a}\) is the Hubble function and \(\Phi ={\dot{\phi }}\).
We continue by defining the new variables \(\left\{ R,\Theta \right\} \) which are given by the point transformation
Therefore, the field equations (29), (30) in the new coordinates are,
from which it follows that the general algebraic solution expressed in parametric form is
where \(I_{0}\) is constant. In the special case where \(6\lambda =1\), the generic algebraic solution follows
We continue our analysis by studying the dynamics of the field equations, specifically, the ones of the (Weyl) Szekeres system.
The phase space portrait of the field equations (32), (33) is presented in Fig. 1.
4 Dynamical analysis
The field equations (14) with time-derivatives can be written in a covariant form using the kinematic variables for the observer: the volume expansion rate \(\theta =3H\), the shear scalar \(\sigma \), the electric part of the Weyl tensor E, and the components of the effective fluid energy density \(\rho \) and pressure p.
In particular, the field equations are then expressed as follows [19, 63]
with the constraint equation,
where \(^{\left( 3\right) }R\) is the spatial curvature of the three-dimensional hypersurfaces. The latter system is known as the Szekeres–Szafron system and has been widely studied in the literature [5, 64,65,66].
In Weyl Integrable theory with a dust fluid source the effective energy density and pressure are \(\rho =e^{-\frac{\phi }{2}}\rho _{m}+\rho _{\phi }~,~p=p_{\phi },~\)in which \(\rho _{\phi }=\frac{\lambda }{2}{\dot{\phi }}^{2} ~\)and\(~p_{\phi }=\frac{\lambda }{2}{\dot{\phi }}^{2}\). In addition, from Eq. (36a) we can write the equivalent system
In the following, we rewrite the field equations (36a)–(36e) using expansion-normalized variables to determine the stationary points of the dynamical system. We remark that every stationary point corresponds to an exact solution of the field equations, which can describe a specific epoch provided by the dynamics of the system. The stability of the stationary points is also determined, which is needed to determine the past and future evolution of the solutions provided by the stationary points.
4.1 Dimensionless variables
We define the new expansion-normalised dimensionless variables
In the new variables, the Szekeres system becomes
with (first integral) constraint equation
where the prime derivative is defined by \(\Omega _{m}^{\prime }\equiv \frac{d\Omega _{m}}{d\tau }\), where \(\tau =\ln a\) and \(a(\tau )\) is the geometric mean expansion scale factor (\({\dot{a}}/a=H)\). Moreover, the parameter for the equation of state of an effective fluid source, \(w_{tot}=\frac{p}{\rho }\), is expressed in terms of the dimensionless variables as
4.2 Stationary points
The set of stationary points, \(\mathbf {P,}\) have coordinates\(~{\mathbf {P}} =\left( \Omega _{m}\left( {\mathbf {P}}\right) ,x\left( {\mathbf {P}}\right) ,\beta \left( {\mathbf {P}}\right) ,\alpha \left( {\mathbf {P}}\right) \right) \), and the physical properties of the exact solutions at these points for the four-dimensional dynamical system (39a)–(39b) are presented below.
Point \(P_{1}=\left( 0,0,0,0\right) \) describes an empty isotropic universe, with spatial curvature \(\Omega _{R}\left( P_{1}\right) =-1\) and a parameter for the equation of state \(w_{tot}\left( P_{1}\right) =-\frac{1}{3}\). From the latter, we infer that the exact solution at the point \(P_{1}\) is the Milne universe. In order to infer the stability of the exact solution at point \(P_{1}\), we determine the eigenvalues of the linearized system around \(P_{1}\). They are \(e_{1}\left( P_{1}\right) =-1\),\(~e_{2}\left( P_{1}\right) =-1\), \(e_{3}\left( P_{1}\right) =-1\) and \(e_{4}\left( P_{1}\right) =-\frac{2}{3}\), hence \(P_{1}\) is an attractor.
Point \(P_{2}=\left( 0,0,\frac{1}{6},0\right) \) has physical quantities \(\Omega _{R}\left( P_{2}\right) =-\frac{3}{4}\). This point describes an anisotropic Kantowski–Sachs universe. The eigenvalues of the linearized system are \(e_{1}\left( P_{2}\right) =0\), \(e_{2}\left( P_{2}\right) =-\frac{3}{4}\), \(e_{3}\left( P_{2}\right) =-\frac{1}{2}\) and \(e_{4}\left( P_{2}\right) =\frac{3}{2}\), from which we infer that \(P_{2}\) is a saddle point, that is, the exact solution at this point is unstable.
Point \(P_{3}=\left( 0,0,-\frac{1}{3},0\right) \) describes a vacuum Bianchi I universe, \(\Omega _{R}\left( P_{3}\right) =0,\) and more specifically, the Kasner universe. The eigenvalues of the linearized system are \(e_{1}\left( P_{3}\right) =0\),\(~e_{2}\left( P_{3}\right) =3\), \(e_{3}\left( P_{3}\right) =3\) and \(e_{4}\left( P_{3}\right) =6\), hence \(P_{3}\) is a source.
Point \(P_{4}=\left( 0,0,\frac{1}{3},\frac{2}{9}\right) \) describes a vacuum Kasner universe, \(\Omega _{R}\left( P_{4}\right) =0\), and the exact solution is unstable. The eigenvalues are \(e_{1}\left( P_{4}\right) =0\) ,\(~e_{2}\left( P_{4}\right) =2\), \(e_{3}\left( P_{4}\right) =3\) and \(e_{4}\left( P_{4}\right) =5\).
Point \(P_{5}=\left( 0,0,-\frac{1}{12},\frac{1}{32}\right) \) gives \(\Omega _{R}\left( P_{5}\right) =-\frac{15}{6}\) which means that the exact solution at the point describes a Kantowski–Sachs universe. The eigenvalues of the linearized system are \(e_{1}\left( P_{5}\right) =-\frac{15}{8}\) ,\(~e_{2}\left( P_{5}\right) =-\frac{3}{4}\), \(e_{3}\left( P_{5}\right) =-\frac{5}{8}\) and \(e_{4}\left( P_{5}\right) =\frac{3}{4}\), which means that \(P_{5}\) is a saddle point.
Points \(P_{6}^{\pm }{=}\left( 0,x,\pm \sqrt{1{-}\lambda x^{2}},\frac{1}{9}\left( 1-\lambda x^{2}{\pm }\sqrt{1{-}\lambda x^{2}}\right) \right) \) are surfaces in the phase space where \(\Omega _{R}\left( P_{6}\right) =0\), which means that the points describe Bianchi I spacetimes. The points are real when \(1-x^{2}\lambda \ge 0\). In the limit where \(x^{2}=\frac{1}{\lambda }\) the solution reduces to that of isotropic FLRW spacetime with a stiff fluid source. The eigenvalues of the linearized system are \(e_{1}\left( P_{6}\right) =0,~e_{2}\left( P_{6}\right) =\frac{6+\sqrt{6}x}{2},~\)
Point \(P_{7}=\left( 1-\frac{1}{6\lambda },-\frac{1}{\sqrt{6}\lambda },0,0\right) \), describes a FLRW spacetime, \(\Omega _{R}\left( P_{7}\right) =0\), where the equation of state parameter for the effective fluid is \(w_{tot}\left( P_{7}\right) =\frac{1}{6\lambda }\). This point is physically acceptable when \(\lambda >\frac{1}{6}\), which means that \(0<w_{tot}\left( P_{7}\right) <1\). The eigenvalues of the linearized system are \(e_{1}\left( P_{7}\right) =1+\frac{1}{2\lambda }\),\(~e_{2}\left( P_{7}\right) =-\frac{3}{2}+\frac{1}{4\lambda }\), \(e_{3}\left( P_{7}\right) =\frac{1}{4}+\frac{5}{72\lambda }+\frac{\sqrt{2916\lambda ^{2}+1044\lambda -191}}{72\lambda } \), \(e_{4}\left( P_{7}\right) =\frac{1}{4}+\frac{5}{72\lambda }-\frac{\sqrt{2916\lambda ^{2}+1044\lambda -191}}{72\lambda }\). Whence it follows that the exact solution at point \(P_{7}\) is unstable.
Point \(P_{8}=\left( -\frac{8}{3}\lambda ,\sqrt{\frac{2}{3}},0,0\right) \) is physical acceptable for \(-\frac{3}{8}<\lambda <0\ \). It describes a FLRW spacetime with spatial curvature \(\Omega _{R}\left( P_{8}\right) =-1-2\lambda ,\) which is always negative for the accepted values of \(\lambda \). The eigenvalues of the linearized system are derived to be, \(e_{1}\left( P_{8}\right) =-1+i\sqrt{2\lambda }\),\(~e_{2}\left( P_{8}\right) =-1-i\sqrt{2\lambda }\), \(e_{3}\left( P_{8}\right) =-\frac{1}{3}-\frac{8}{9}\lambda +\frac{1}{9}\left( \sqrt{64\lambda ^{2}+156\lambda +63}\right) \), \(e_{4}\left( P_{8}\right) =-\frac{1}{3}-\frac{8}{9}\lambda -\frac{1}{9}\left( \sqrt{64\lambda ^{2}+156\lambda +63}\right) \) from which we conclude that the exact solution at \(P_{8}\) is always unstable.
Point \(P_{9}{=}\left( -\frac{3\lambda \left( 2\lambda {+}5\right) }{2\left( \lambda {+}2\right) ^{2}},\frac{1}{2{+}\lambda }\sqrt{\frac{3}{2}},-\frac{1}{3}+\frac{1}{4+2\lambda },\frac{3+4\lambda \left( \lambda +2\right) }{24\left( \lambda +2\right) ^{2}}\right) \) describes a Kantowski–Sachs universe where \(\Omega _{R}\left( P_{9}\right) =-\frac{3\left( 2\lambda +1\right) \left( 2\lambda +5\right) }{4\left( \lambda +2\right) ^{2}}\). The point is physical acceptable for \({-}\frac{23{-}\sqrt{273}}{16}{\le }\lambda <0\) and \(-\frac{5}{2} \le \lambda <-\frac{23+\sqrt{273}}{16}\). The eigenvalues are calculated numerically, from which we infer that point \(P_{9}\) is a saddle point.
The above results are summarized in Table 1. In Fig. 2, the qualitative behaviour of the equation of state parameter \(w_{tot}\) is presented. Moreover, two-dimensional phase portraits for the dynamical system (39a)–(39b) are presented in Figs. 3 and 4 where \(P_{1}\) is the unique attractor. The plots are for positive and negative values of the coupling parameter \(\lambda \).
4.3 Past attractors
When analyzing the dynamics of the system (39a)–(39d) towards the past, it is convenient to make a time reversal \(\tau \mapsto -\tau \). In this case, we have the same points as before, but there is an overall change of sign in the eigenvalues. Then, the possible late-time attractors of the new system, given by
correspond to the past attractors of the original one. We study the points \(P_{3}\) with coordinates \((\Omega _{m},x,\beta ,\alpha )=\left( 0,0,-\frac{1}{3},0\right) \), \(\Omega _{R}=0\), and \(P_{4}\) with coordinates \((\Omega _{m},x,\beta ,\alpha )=\left( 0,0,\frac{1}{3},\frac{2}{9}\right) \), \(\Omega _{R}=0\), corresponding to Bianchi I models, and we show they are unstable for the original system using the center manifold theorem (CMT). The detailed analysis of the CMT for these two points is presented in Appendix A.
5 Conclusions
In this work we found exact inhomogeneous spacetimes which generalize the Szekeres universes into the Weyl integrable theory. Specifically, we assume that the scalar field which defines the Weyl affine connection to be homogeneous such that the limit of FLRW exists. In such scenario, the only inhomogeneous spacetimes are those which belong to the FLRW (-like) solutions included in the family of LTB spacetimes. On the other hand, the Kantowski–Sachs family of solutions is homogeneous and anisotropic. For the inhomogeneous spacetimes, we were able to write in terms of quadratics the generic solution of the field equations.
In order to understand the dynamics and the evolution of the gravitational model we performed a detailed study of the past and future attractors. In particular, we defined Hubble-normalized dimensionless variables. The field equations admit three stationary points which describe a spatially flat, an open, and a closed FLRW space where only the closed FLRW spacetime can be a future attractor, which gives the Milne universe. The other two isotropic solutions correspond to saddle points.
In addition, three homogeneous Kantowski–Sachs spacetimes are supported by the field equations which correspond to saddle points. There are three points which describe Bianchi I spacetimes; the exact solution at one of these points describes a Bianchi I spacetime with a stiff fluid, while the other two points describe vacuum Kasner solutions. The points which describe the Kasner solutions are sources while the third point is a saddle point.
For the sources we performed a detailed study on the past-system in order to investigate if the points are attractors for the past-system. Indeed with the application of the center manifold theorem we were able to prove that the Kasner solutions are past attractors of the field equations for \(\lambda >0\).
Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors’ comment:This is a Theoretical Research Project.]
Notes
We consider the signature of the metric to be \(\left( -,+,+,+\right) \).
References
P. Szekeres, Commun. Math. Phys. 41, 55 (1975)
A. Krasiński, Inhomogeneous Cosmological Models (Cambridge U.P, Cambridge, 1997)
W.B. Bonnor, Mon. Not. R. Astron. Soc. 167, 55 (1974)
A.N. Golubiantnikov, L.M. Truskinovskii, Prikl. Matem. Mekhan. 45, 956 (1981)
M. Bruni, S. Matarrese, O. Pantano, Astrophys. J. 445, 958 (1995)
J.D. Barrow, J. Silk, Astrophys. J. 250, 432 (1981)
S. Goode, J. Wainwright, Phys. Rev. D 26, 3315 (1982)
P. Szekeres, Phys. Rev. D 12, 2941 (1975)
D.J. Shaw, J.D. Barrow, Phys. Rev. D 73, 123506 (2006)
W.B. Bonnor, Commun. Math. Phys. 51, 191 (1976)
B.K. Berger, D.M. Eardley, D.W. Olson, Phys. Rev. D 16, 3086 (1977)
G.M. Covarrubias, J. Phys. A Math. Gen. 13, 3023 (1980)
P.S. Joshi, A. Królak, Class. Quantum Gravity 13, 3069 (1996)
W.B. Bonnor, N. Tomimura, Mon. Not. R. Astron. Soc. 175, 85 (1976)
L. Herrera, A. Di Prisco, J. Ibanez, J. Carot, Phys. Rev. D 76, 044003 (2012)
C. Hellaby, Class. Quantum Gravity 34, 145006 (2017)
J. Wainwright, S. Andrews, Class. Quantum Gravity 26, 085017 (2009)
A. Coley, N. Layden, D. McNutt, Gen. Relativ. Gravit. 51, 164 (2019)
J. Wainwright, G.F.R. Ellis, Dynamical Systems in Cosmology (Cambridge U.P, Cambridge, 2005)
R.A. Sussman, K. Bolejko, Class. Quantum Gravity 29, 065018 (2012)
R.A. Sussman, J.C. Hidalgo, I.D. Gaspar, G. German, Phys. Rev. D 95, 064033 (2017)
D.A. Szafron, J. Math. Phys. 18, 1673 (1977)
J.D. Barrow, J. Stein-Schabes, Phys. Lett. A 103, 315 (1984)
S.W. Goode, J. Wainwright, Gen. Relativ. Gravit. 18, 315 (1986)
N. Tomimura, Nuovo Cimento B 42, 1 (1977)
J.A.S. Lima, M.A.S. Nobre, Class. Quantum Gravity 7, 399 (1990)
J.A.S. Lima, M.R. Garica Maia, Phys. Lett. A 110, 366 (1985)
D.C. Motta, N. Tomimura, Astrophys. Space Sci. 165, 237 (1990)
N. Tomimura, D.C. Motta, Astrophys. Space Sci. 165, 243 (1990)
J.D. Barrow, A. Paliathanasis, EPJC 79, 379 (2019)
A. Paliathanasis, Class. Quantum Gravity 37, 105008 (2020)
C. Romero, J.B. Fonseca-Neto, M.L. Pucheu, Class. Quantum Gravity 29, 155015 (2012)
V. Perlick, Class. Quantum Gravity 8, 1369 (1991)
I. Quiros, arXiv:gr-qc/0011056
I. Quiros, arXiv:1401.2643
I. Quiros, R. De Arcia, arXiv:1811.02458
R. Avalos, F. Dahia, C. Romero, Found. Phys. 48, 253 (2018)
R. Avalos, I. Lobo, T. Sanomiya, C. Romero, J. Math. Phys. 59, 062502 (2018)
R. Aguila, J.E. Madriz Aguilar, C. Moreno, M. Bellini, Eur. Phys. J. C 74, 3158 (2014)
M. Konstantinov, V. Melnikov, Int. J. Mod. Phys. D 4, 339 (1995)
J. Villanueva, F. Tapia, M. Molina, M. Olivares, Eur. Phys. J. C 78, 853 (2018)
J.E. Madriz Aguilar, C. Romero, Int. J. Mod. Phys. A 24, 1505 (2009)
J.E. Madriz Aguilar, C. Romero, J.B. Fonseca Neto, T.S. Almeida, J.B. Formiga, Class. Quantum Gravity 32, 215003 (2015)
K.A. Bonnikov, M. Yu, Konstantinov and V.N. Melnikov, Gravit. Cosmol. 1, 60 (1995)
J.M. Salim, S.L. Sautu, Class. Quantum Gravity 13, 353 (1996)
J.M. Salim, S. Sautu, Class. Quantum Gravity 15, 203 (1998)
J.M. Salim, S. Sautu, Class. Quantum Gravity 16, 3281 (1999)
M. Israelit, Found. Phys. 29, 1303 (1999)
S.E. Perez-Bergliaffa, K.E. Hibrerd, Int. J. Mod. Phys. D 8, 705 (1999)
T.S. Almeida, M.L. Puckeu, C. Romero, J.B. Formiga, Phys. Rev D 89, 064047 (2014)
L. Amendola, Phys. Rev. D 62, 043511 (2000)
L. Amendola, C. Quercellini, Phys. Rev. D 68, 023514 (2003)
D. Pavón, W. Zimdahl, Phys. Lett. B 628, 206 (2005)
S. del Campo, R. Herrera, D. Pavón, Phys. Rev. D 78, 021302 (2008)
L. Amendola, G. Camargo Campos, R. Rosenfeld, Phys. Rev. D 75, 083506 (2007)
D. Pavón, B. Wang, Gen. Relativ. Gravit. 41, 1 (2009)
L.P. Chimento, Phys. Rev. D 81, 043525 (2010)
F. Arevalo, A.P.R. Bacalhau, W. Zimdahl, Class. Quantum Gravity 29, 235001 (2012)
W. Yang, S. Pan, J.D. Barrow, Phys. Rev. D 97, 043529 (2018)
J.D. Barrow, A. Paliathanasis, EPJC 78, 767 (2018)
V. Faraoni, Phys. Rev. D 59, 084021 (1999)
K. Bolejko, A. Krasiński, C. Hellaby, M.-N. Célérier, Structures in the Universe by Exact Methods (Cambridge U. Press, Cambridge, 2010)
H. van Elst, G.F.R. Ellis, Class. Quantum Gravity 13, 1159 (1996)
J. Libre, C. Valls, Phys. Lett. A 383, 301 (2019)
A. Paliathanasis, P.G.L. Leach, Phys. Lett. A 381, 1277 (2017)
A. Gierzkiewicz, Z.A. Golda, J. Nonlinear Math. Phys. 24, 494 (2016)
Acknowledgements
AP & GL were funded by Agencia Nacional de Investigación y Desarrollo - ANID through the program FONDECYT Iniciación Grant no. 11180126. Additionally, GL is supported by Vicerrectoría de Investigación y Desarrollo Tecnológico at Universidad Catolica del Norte. Ellen de los Milagros Fernández Flores is acknowledged for proofreading. JDB is supported by the STFC of the United Kingdom.
Author information
Authors and Affiliations
Corresponding author
Appendix A: CMT for points \(P_{3}\) and \(P_{4}\)
Appendix A: CMT for points \(P_{3}\) and \(P_{4}\)
1.1 Analysis of \(P_{3}\)
Introducing the coordinate transformations,
the equilibrium point \(P_{3}\) is translated to the origin and the linearization matrix is transformed to its canonical real Jordan form.
Therefore, we obtain the equivalent dynamical system to (41), defined by
The eigen system of the origin is
That is, the center manifold of the origin is tangent to the u-axis, and it is given locally by a graph
which satisfies the differential equations
where
Using Taylor expansions, we propose as Ansätze,
Substituting in (A3)–(A5), and equating the coefficients of equal powers of u, we obtain
Therefore,
and we have the parametrization,
where we choose \(\lambda v_{2}\ge 0\).
The dynamics on the center manifold of the origin are dictated by a gradient-like equation \(u^{\prime }=-\nabla U(u)\). For \(\lambda >0\), \(\omega =u\sqrt{\lambda }\), the equation transforms to
for which the origin is a degenerated minimum.
For \(\lambda <0\), \(\omega =u\sqrt{-\lambda }\), the equation transforms to
for which the origin is a degenerated maximum.
Therefore, for \(\lambda >0\) (respectively, \(\lambda <0\)) the center manifold, and hence, the origin of the dynamical system is a local attractor (respectively, a saddle). the origin of the dynamical system is a local attractor of system is \(P_{3}\), and for \(\lambda <0\) it is a saddle point.
Now, we take the time reversal back and work in terms of t. We deduce:
The solutions can be expressed as:
where
Then,
Finally,
and
where \(c_{3}\) and \(c_{4}\) are integration constants, and we set \(\lambda v_{2}=-K_{0}\varepsilon ^{14}\), for a positive constant \(K_{0}\). For \(\lambda >0\), \(\theta (t)\rightarrow \frac{1}{t-t_{0}}\) as \(t\rightarrow 0\) (\(\tau \rightarrow -\infty \)). Hence, \(P_{3}\) it is associated with (an anisotropic) initial singularity.
1.2 Analysis of \(P_{4}\)
Introducing the coordinate transformation
the equilibrium point \(P_{4}\) is translated to the origin and the linearization matrix is transformed to its canonical real Jordan form.
Therefore, we obtain the equivalent dynamical system to (41), defined by
The eigen system of the origin is
That is, the center manifold of the origin is tangent to the u-axis, and it is given locally by a graph
which satisfies the differential equations
where
Using Taylor expansion we propose as Ansätze
Hence it follows
Therefore,
and we have the parametrization
where we choose \(\lambda v_{2}\ge 0\).
The dynamics on the center manifold of the origin are dictated by a gradient-like equation \(u^{\prime }=-\nabla U(u)\). For \(\lambda >0\), \(\omega =u\sqrt{\lambda }\), the equation transforms to
for which the origin is a degenerated minimum. For \(\lambda <0\), \(\omega =u\sqrt{-\lambda }\), the equation transforms to
for which the origin is a degenerated maximum. Therefore, for \(\lambda >0\) (respectively, \(\lambda <0\)) the center manifold, and hence, the origin of the system, is a local attractor (respectively, a saddle). In the original variables mean that for \(\lambda >0\) the past attractor of the dynamical system is \(P_{4}\), and for \(\lambda <0\) is a saddle point. That is exactly the same dynamics as for \(P_{3}\). However, as we will see shortly, the physical solution, although it is Bianchi I, has a different asymptotic expansion.
Now, we take the time reversal back and work in terms of t. Hence
As before,
however for this point
where \(c_{3}\) and \(c_{4}\) are integration constants, and we set \(\lambda v_{2}=-K_{0}\varepsilon ^{14}\), for a positive constant \(K_{0}\). As for \(P_{3} \), for \(\lambda >0\), \(\theta (t)\rightarrow \frac{1}{t-t_{0}}\) as \(t\rightarrow 0\) (\(\tau \rightarrow -\infty \)). Hence, \(P_{4}\) it is associated with (an anisotropic) initial singularity. However, as per the physical solution referred to, this is a different solution with different asymptotic expansions for \(\sigma ,{E}\).
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Funded by SCOAP3
About this article
Cite this article
Paliathanasis, A., Leon, G. & Barrow, J.D. Inhomogeneous spacetimes in Weyl integrable geometry with matter source. Eur. Phys. J. C 80, 731 (2020). https://doi.org/10.1140/epjc/s10052-020-8277-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epjc/s10052-020-8277-z