Inhomogeneous spacetimes in Weyl integrable geometry with matter source

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


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, which is why 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 [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 [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 [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, whose value determines the asymptotic spatial inhomogeneities in various physical quantities. It 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 [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 Λ-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 [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 Λ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 interesting in determining 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 µν is a metric tensor with Levi-Civita connection Γ κ µν , then in Weyl Integrable theory the manifold is supported by the set g µν ,Γ κ µν whereΓ κ µν is the Levi-Civita connection for the conformally related metricg µν = φg µν , where φ 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 [34][35][36]. In [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 [38]. In [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 [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 [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 [44], while higher-or lower-dimensional gravitational models were studied in [40,42,43]. In [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 [46][47][48][49].
The plan of this paper is as follows. In Section 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 Section 4. In Sections 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 Section 6 we discuss our results and draw our conclusions.

WEYL INTEGRABLE GRAVITY
Weyl geometry is an extension of Riemannian geometry, specified by a metric tensor g µν and a gauge vector field ω µ . The covariant derivative∇ µ is defined by the (Weyl) affine connectionΓ κ µν , with the property, from which we infer that the gauge vector ω µ plays a significant role in the geometry. Specifically, the Weyl affine connectionΓ κ µν is related to the Christoffel symbols Γ κ µν of the metric tensor g µν as follows: The curvature tensor in the Weyl geometry is defined as where, by using (2), we observe that in generalR κλµν is not symmetric as it is in the case of Riemannian geometry.
In this work, we are interesting in the case whereR κλµν has the same symmetric properties as in Riemannian geometry. This is true when ω µ is a gradient vector, which means that there exists a scalar φ such that ω µ = φ ,µ . 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 µν 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 tensorR µν is related to the Riemannian Ricci tensor R µν by, where the Ricci scalar, in a four-dimensional manifold, is written as [45]

Gravitational Action Integral
We define the simple gravitational action Integral, which includes the Weyl Ricci scalarR and the field φ ;µ , as where ξ 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],G whereG µν is the Weyl Einstein tensor. Moreover, variation with respect to the scalar field φ gives that is, a Klein-Gordon equation of the form, where ∇ µ denotes the Riemannian covariant derivative. The gravitational field equations (8) can be rewritten by using the Riemannian Einstein tensor G µν as follows in which the new constant λ is defined as 2λ ≡ 4ξ − 3. Consequently, the field equations (11) for λ > 0 are those of general relativity 1 with a massless scalar field. A new possibility is introduced when λ < 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.

The presence of dust
When a pressureless fluid dust source is included, the gravitational field equations become where T while the conservation equation for the matter field reads∇ ν T (m)µν = 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 ρ m . The interaction between scalar field and dust fluid has been proposed as a potential mechanism to explain the cosmic coincidence problem. [50][51][52][53]. Various interaction models have been studied before in the literature; for instance, see: [54][55][56][57][58], and references therein.

INHOMOGENEOUS SPACETIMES
The gravitational model that we have considered in Weyl Integrable geometry is equivalent to that of general relativity and with an effective energy-momentum tensor, where the field equations are of the form and T µν is the effective energy-momentum tensor. It consists of a massless scalar field φ, and an additional fluid source interacting with the field φ. In particular, µν +T µν , whereT µν describes the energy-momentum tensor of a pressureless fluid, i.e.T µν = T (m) µν , or of the second scalar field, ψ; that is, where µν is the energy-momentum tensor of the massless scalar field, However, as we discussed before the continuous equation ∇ ν T µν = 0, provides ∇ ν T (φ)µν +T µν = 0, that is, In this work, we assume that the underlying spacetime is described by the inhomogeneous and anisotropic diagonal line element in which A = A (t, r, y, z) and B = B (t, r, y, z). 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 (t, y, z), the exact solutions belong to the Kantowski-Sachs (-like) family, while for B ,r = 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 [59]. 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 [59] when the scalar field is massless, we recover the exact solution of the vacuum Weyl Integrable geometry for λ > 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 µν in (14) to be homogeneous such that the FLRW limit to be provided. Thus, φ = φ (t) while for the matter source we have ρ m = ρ m (t). The latter follows easily, if we rewrite the energy-momentum tensor T µν 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.
( [59]). 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) , so that the line element is written as [59]: and the curvature, K, of the two-dimensional surface of constant curvature {y − z} 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 α (t) , β (t) 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 [59]: ∂C (r, y, z) ∂r 2 dr 2 + e 2C(r,y,z) dy 2 + dz 2 .

Kantowski-Sachs spacetimes
The unknown scale factors of the Kantowski-Sachs spacetime (18) are given by the following system, 2β 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.

FLRW spacetimes
Analogously, the unique scale factor for the FLRW (-like) spacetime (19) is given by the (modified) Friedmann equations and the scalar field φ satisfies the Klein-Gordon equation, while the conservation equation for the dust fluid source iṡ At this point, we remark that for λ = 0, only the vacuum solutions of general relativity are recovered, while the Szekeres spacetimes are recovered when φ = φ 0 and λ → ∞. This is reminiscent of the range of the constant Brans-Dicke parameter, ω, in scalar-tensor theory such, where the limit of general relativity to be recovered as ω → ∞ [60]. This family of spacetimes includes also the inhomogeneous Lemaître-Tolman-Bondi (LTB) spacetimes [61].
In the following we show the analytic solution for the inhomogeneous FLRW (-like) spacetime.

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 =ȧ a is the Hubble function and Φ =φ. We continue by defining the new variables {R, Θ} 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λ = 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.

DYNAMICAL ANALYSIS
The field equations (14) with time-derivatives can be written in a covariant form by using the kinematic variables for the observer: the volume expansion rate θ = 3H, , the shear scalar σ, the electric part of the Weyl tensor E , and the components of the effective fluid energy density ρ and pressure p.
In Weyl Integrable theory with a dust fluid source the effective energy density and pressure are ρ = e − φ 2 ρ m +ρ φ , p = p φ , in which ρ φ = λ 2φ 2 and p φ = λ 2φ 2 . In addition, from equation (36a) we can write the equivalent systeṁ In the following, we rewrite the field equations (36a)-(36e) by 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.

Dimensionless variables
We define the new expansion-normalised dimensionless variables In the new variables, the Szekeres system becomes Kantowski-Sachs Unstable with (first integral) constraint equation where the prime derivative is defined by Ω ′ m ≡ dΩm dτ , where τ = ln a and a(τ ) is the geometric mean expansion scale factor (ȧ/a = H). Moreover, the parameter for the equation of state of an effective fluid source, w tot = p ρ , is expressed in terms of the dimensionless variables as

Stationary points
The set of stationary points, P, have coordinates P = (Ω m (P) , x (P) , β (P) , α (P)), 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 = (0, 0, 0, 0) describes an empty isotropic universe, with spatial curvature Ω R (P 1 ) = −1 and a parameter for the equation of state w tot (P 1 ) = − 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 (P 1 ) = −1, e 2 (P 1 ) = −1, e 3 (P 1 ) = −1 and e 4 (P 1 ) = − 2 3 , hence P 1 is an attractor.
describes a Kantowski-Sachs universe where Ω R (P 9 ) = . The eigenvalues are calculated numerically, from which we infer that point P 9 is a saddle point.
The above results are summarized in Table I. 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 λ.

Past attractors
When analyzing the dynamics of the system (39) towards the past, it is convenient to make a time reversal τ → −τ . In this case, we have the same points as before, but there is an overall change of sign in the eigenvalues. Then, the correspond to the past attractors of the original one. We study the points P 3 with coordinates (Ω m , x, β, α) = 0, 0, − 1 3 , 0 , Ω R = 0, and P 4 with coordinates (Ω m , x, β, α) = 0, 0, 1 3 , 2 9 , Ω R = 0, corresponding to Bianchi I models, and we show they are unstable for the original system using the center manifold theorem.