On the Effects of Non-metricity in an Averaged Universe

In the covariant averaging scheme of macroscopic gravity, the process of averaging breaks the metricity of geometry. We reinterpret the back-reaction within macroscopic gravity in terms of the non-metricity of averaged geometry. This interpretation extends the effect of back-reaction beyond mere dynamics to kinematics of geodesic bundles. With a 1+3 decomposition of the spacetime, we analyse how geometric flows are modified by deriving the Raychaudhuri and Sachs equations. We also present the modified forms of Gauss and Codazzi equations. Finally, we derive an expression for the angular diameter distance in Friedmann Lema\^itre Robertson Walker universe and show that non-metricity modifies it only through the Hubble parameter. Thus, we caution against overestimating the influence of back-reaction on the distances.


Introduction
The universe is traditionally modelled by employing cosmological solutions to the Einstein field equations.In these models, the matter distribution is taken to be averaged over cosmological scales, and hence the Einstein tensor needs to be averaged as well [1,2].However, averaging tensor fields on a differentiable manifold is not straightforward and is a well-known problem in relativistic cosmology [3,4].The difficulty arises because of the lack of a mathematically rigorous integration operation on non-Euclidean geometries.Moreover, the non-linearity of the Einstein field equations exacerbates the problem since the averaged dynamics of a given metric are different from the dynamics of an averaged metric, i.e., ⟨G[g]⟩ ̸ = G[⟨g⟩].This is referred to as the averaging problem in the literature.There is a long history of approaches towards this problem (see [3,5,6] and the references within).In all the approaches, however, the problem is effectively addressed by modifying the dynamics of the averaged spacetime.This modification may take different forms in various approaches, and is generically referred to as cosmological back-reaction [5,6].One such approach is macroscopic gravity [7][8][9][10] which offers a non-perturbative and fully covariant solution to the averaging problem.
Macroscopic gravity calculates space-time averages of arbitrary tensor fields through averaging bivectors and Lie dragging of the averaging regions.When this averaging scheme is applied to general relativity (GR), the Riemannian1 structure of the geometry is broken.The averaged geometry (the geometry at the macroscopic scales) is equipped with a metric-incompatible but symmetric connection that arises from averaging.This leads to a non-vanishing non-metricity but a vanishing torsion.The dynamics of macroscopic gravity is given by taking a space-time average of the Einstein field equations.The back-reaction in these dynamics is also given in terms of the non-metricity, as has been noted in [11,12].
The aim of this work is to analyse the effects of non-metricity in macroscopic gravity on various aspects of relativistic cosmology.We do this by treating macroscopic gravity as a metric-affine theory.In such a theory, the geometry is endowed with a metric and an independent connection which is precisely how macroscopic gravity proposes the averaged geometry to be.This means that the back-reaction (arising due to non-metricity) modifies the kinematics in the averaged universe.This work emphasises this by demonstrating its effects on geometric flows, which in turn affects the distance measures in cosmology.
We perform a 1+3 decomposition of the macroscopic geometry.This enables us to derive the kinematics of geodesic flows.We show that both the Raychaudhuri and Sachs optical equations have contributions from the non-metricity.In other words, back-reaction in macroscopic gravity modifies not only the dynamics but also the kinematics.However, it turns out that the non-metricity has no effect on the cross sectional area of geodesic bundles.We use the Sachs equation to derive an expression for the angular distances in an arbitrary space-time.As a quick application of this geometric setup, we present the case of a spatially flat Friedman-Lemaître-Robertson-Walker (FLRW) universe and discuss how the non-metricity affects the distance measurements.We also show that the spatial hypersurfaces have a non-Riemannian geometry and derive the modified Gauss and the Codazzi equations.
The organisation of the paper is as follows: in section 2, we present the basic formalism of macroscopic gravity and reinterpret it as a metric-affine theory.In section 3, we derive the structure equations for geometric flows.Then, in section 4, we discuss the geometry of spatial hypersurfaces and derive the Gauss-Codazzi equations.In section 5, we derive an expression for the angular distances.We show that the back-reaction affects the distances only through the dynamics (Hubble expansion).Finally, in section 6, we discuss the results in this paper and give a few concluding remarks.
The notation used is as follows: Greek indices run from 0 to 3 while Latin indices from 1 to 3. The partial derivative is denoted by ∂ while ∇ and ∇ denote the covariant derivatives with respect to the non-Riemannian and the Levi-Civita connections, respectively.Averaged quantities are denoted using angular brackets, ⟨. ..⟩.The indices within (. . . ) and [. . .] are are symmetrised/anti-symmetrised; and underlined indices are not included in (anti-)symmetrisation.The sign convention followed is the 'Landau-Lifshitz Space-like Convention (LLSC)' [13].

Basics of Macroscopic Gravity Formalism
In this section we briefly present the basic formalism of macroscopic gravity.We use the same terminology as used in [7][8][9]11].We will reinterpret certain quantities within the macroscopic gravity formalism in terms of more familiar geometric objects in a metric-affine gravity.Macroscopic gravity employs a covariant averaging procedure to average out the (pseudo-)Riemannian geometry and the field equations of general relativity [7-11, 14, 15].Given a geometric object, p α β (x), defined on an n-dimensional differentiable metric manifold (M, g), the space-time averaged value of this object is defined as, where, the average has been taken over a compact region ) is called the bilocal extension of the object p α β (x); and the objects A α µ ′ (x, x ′ ) and A ν ′ β (x ′ , x) are called the bilocal averaging operators.The average is taken over finite volume regions in the manifold (i.e., V Σ is finite).Hence, well defined microscopic objects average out without the integral diverging in equation (2.1).In the limit of vanishing averaging volume, the average of the object is equal to the object itself [7][8][9].Further, there is freedom to choose these averaging operators according to the problem at hand.However, all such operators must follow the following properties: coincidence limit, idempotency, existence of an inverse, and factorisation (see [7,8] for details).Moreover, the averaging procedure is constructed in such a way that the average of a quantity is uniquely defined and does not depend on the choice of the averaging operators [7][8][9].
The average of the microscopic Riemann curvature tensor is written as R α βρσ and behaves like a curvature tensor [7][8][9]11].Further, bilocal objects defined as, F α βρ (x, ρ , behave as connection coefficients at x and hence, can be considered to be the bilocal extension of the microscopic connection coefficients, Γ α βρ .They have the following coincidence limit, lim The average of the connection coefficients, x ′ , are the macroscopic affine connection coefficients [7,8].Then, one cam construct a curvature tensor, M α βρσ , with these macroscopic connection coefficients.This will not be the same as the average of the microscopic curvature tensor, but rather is related to it through, (2. 3) The additional terms on the right hand side are arising from the fact that the curvature is a non-linear function of the connection coefficients.
One can calculate another connection, Π α βρ , associated with the curvature tensor R α βρσ .Therefore, in the macroscopic space-time, there are two connections associated with two curvature tensors which are related to each other by equation (2.3).An affine deformation tensor, A α βρ = Γ α βρ − Π α βρ , can be defined that measures the difference in the two connections.Such a deformation tensor exists due to the metric incompatibility of the connection, Π α βρ .That is, for macroscopic metric, g αβ , we have, where | denotes the covariant derivative with respect to the connection Π α βρ .Further, looking at the construction of both the curvature tensors in terms of their respective connections, equation (2.3) takes the form, where ∥ denotes the covariant derivative with respect to the metric-compatible (Levi-Civita) connection, ⟨Γ⟩, of the macroscopic manifold.The tensor, P, is dubbed as the gravitational polarisation tensor [7], and is related to the the bilocal objects F as, (2.6) The tensor, P, although not a curvature tensor, follows the algebraic and differential properties of one.
From equations (2.3) and (2.6), we can see the that it measures the difference between the two curvature tensors -a curvature deformation tensor, Both the curvature tensors, M and ⟨R⟩, follow the same algebraic identities and also the differential Bianchi identities with respect to their corresponding connections.The former ensures that the connection, Π, even though metric incompatible, is symmetric.Moreover, averaging of the differential Bianchi identities for the microscopic curvature tensor gives the following relations, In order to evaluate the above equations, a splitting rule is required for the products in the second and third terms.This is constructed using an object called the connection correlation, which is defined as [7,8], (2.9) Using equations (2.6) and (2.9), we see that, (2.10) In summary, when the averaging procedure of macroscopic gravity (2.1) is applied to the 4-dimensional (pseudo-)Riemannian manifold of GR, one finds that at macroscopic scales, the manifold is characterised by a metric, two symmetric affine connections 2 ravity theories with multiple connections both in the Riemannian and non-Riemannian framework are already present in the literature [16,17]., two curvature tensors associated with these connections and a correlation tensor that needs to be constructed in order to solve differential Bianchi identities for these curvature tensors.

The Field Equations of Macroscopic Gravity
The average of Einstein field equations with a cosmological constant, yields the field equations of macroscopic gravity, which are given by [7,8], where, νγ] , P = P µ µ , P µν = P ϵ µϵν , and ⟨T⟩ is the averaged energy- momentum tensor (usually taken to be that of a perfect fluid in cosmology).Therefore, upon averaging, the field equations get modified and now include additional terms given (partially) in terms of nonmetricity of the space-time.

Non-metricity of the Averaged Geometry
As discussed above, applying the averaging procedure of macroscopic gravity to GR leads to the interesting result that at macroscopic scales, the manifold is non-Riemannian.That is, the manifold has a Riemannian structure only at small (microscopic) scales.At such scales, GR is a valid theory of gravity.However, when one is concerned with large scales (like in cosmology), an averaging has to be performed on the geometry, and in the process, the Riemannian structure of the manifold is broken.In other words, for cosmological analyses, one must work with a theory of gravity that employs a non-Riemannian geometrical framework.The theory of macroscopic gravity is precisely such a theory.
According to the macroscopic gravity formalism, the non-Riemannian (macroscopic) manifold is characterised by a symmetric but metric-incompatible connection, given by Π.The covariant derivative of the metric (of the macroscopic manifold) with respect to this connection is not zero.The non-metricity tensor is then defined as, Using this, the general connection, Π, can be decomposed as, where, Γ is the Riemannian part (Levi-Civita connection) 3 , and D is the so-called disformation tensor that encapsulates the deviation from the Riemannian nature of geometry.The disformation is given in terms of the non-metricity as, Further, the Ricci identity is given by,  2 G 3 The Levi-Civita connection can be written in terms of metric as, where, R ρ σαβ , is the Riemann curvature tensor defined in terms of the general connection of the non-Riemannian geometry, Using equation (2.14), it can be shown that, where, R ρ σαβ is the Riemannian part of the curvature tensor which is the same as M of the previous section, and ∇ is the covariant derivative both defined with the Levi-Civita connection.One can define a new quantity that captures the deviation of the curvature tensor from its Riemannian counterpart, It is not difficult to see that the so-called affine deformation tensor, A, and the polarisation tensor, P, of macroscopic gravity are, respectively, the disformation tensor, D, and the object, N, of the metricaffine gravity literature.Henceforth, we will use these symbols to denote the affine deformation tensor and polarisation tensor since we are treating the theory of macroscopic gravity as a metric-affine theory.

Structure Equations for Geometric Flows
Since the geometry at the macroscopic scales is no longer Riemannian, determining the geodesic curves (curves with extremum length) is not trivial.In Riemannian geometry, the autoparallel curves with respect to the Levi-Civita connection and also the geodesic curves.Then, the geodesic equation is given by, However, in general, autoparallel curves (of a non-Riemannian geometry) are not necessarily geodesics.In fact, since in non-Riemannian geometry, there exist additional geometric fields (like non-metricity) the paths traced by freely moving physical particles (the geodesics) will depend on whether these particles interact with these 'new' geometric fields [18].In cosmology, one takes the content of the universe to be a perfect fluid or dust.Since such matter does not interact with non-metricity, despite the averaged geometry being non-Riemannian, the path traced is taken to be Riemannian geodesics given by equation (3.1).
We will now derive equations governing the evolution of the congruences of geodesics on the averaged (non-Riemannian) geometry.These are sometimes called kinematic equations.Given a congruence, one can characterise its flow in terms of the evolution of the irreducible components of its velocity gradient.Physically, these equations characterise the flow of the fluid content of the Universe (for the timelike curves) and bundle of photons (for the null curves).We will derive (in detail) the equations for the trace part of the velocity gradient (also called expansion) for both the timelike and null curves.This will give us the macroscopic gravity analogue of the Raychaudhuri and Sachs equations [19,20].

The Raychaudhuri Equation
To derive the Raychaudhuri equation, we begin with a congruence of timelike geodesics.Let u α be the vector field tangent to the congruence (the velocity vector).In cosmology, this is simply the comoving velocity of the (fluid) content of the Universe, u α = (1, 0, 0, 0).The norm of such a velocity is normalised as, u α u α = −1 [21].We are interested in the transverse evolution of the congruence and hence we define a transverse metric in the following manner [21], The 3-metric, h αβ , is completely transverse to the congruence, i.e., we have, u α h αβ = 0 = h αβ u β .It is easy to see that, h α α = 3.This splits the space-time in temporal and spatial parts, referred to as 1 + 3 decomposition.All the quantities transverse to the congruence reside on a spatial hypersurface orthogonal to it.Here, we will investigate the transverse properties of the congruence while in the next section we present the properties of spatial hypersurfaces.
Using the 3-metric, one can separate out the transverse and longitudinal parts of the velocity gradient as follows, where, ∇ represents the covariant derivative with respect to the non-Riemannian connection, Π.It is easy to check that, u α (T) ∇ α u β = 0 = (T) ∇ α u β u β .
The velocity gradient can be decomposed into its irreducible components as follows, where, θ is the trace part known as expansion, ω is the anti-symmetric part known as rotation or vorticity, and σ is the symmetric traceless part known as shear.
Using the expression for the expansion scalar, we can write, where, we have used, u α ∇ α u β = D β αρ u α u ρ which follows from equation (3.1).One can decompose the first term above into a Riemannian and a non-Riemannian part.The Riemannian part would be the usual expansion scalar, θ = ∇ α u α .Then, we have, Taking a derivative on both sides gives us the Raychaudhuri equation in the averaged geometry, where, the terms with a bar over them are Riemannian quantities.The last term on the right hand side is the contribution from averaging (the back-reaction) given solely on terms of the non-metricity of the averaged space-time.

The Sachs Optical Equation
The Sachs optical equation is the null counterpart of the Raychaudhuri equation.To derive it, we begin with considering a congruence of null geodesics with tangent vector, k α .The norm of the tangent is, k α k α = 0.In this case, the transverse metric is given by, where, n α is an auxiliary null vector satisfying the conditions, n α n α = 0 and n α k α = −1.Once again, one can check, k α H αβ = 0 = H αβ k β and h α α = 2.The steps involved in deriving the Sachs equation are analogous to previous section.We define the transverse part of the velocity gradient as, This can be can be decomposed into irreducible components as follows, ) Using the definition of the expansion scalar, we get, where, we have used, The first term can again be decomposed into Riemannian and non-Riemannian parts.The Riemannian part is the usual null expansion scalar, Θ = ∇ α k α .This gives us, Taking a derivative on both sides gives the Sachs optical equation in the averaged geometry, Therefore, the expansion of both the timelike and null geodesics is affected by the non-metricity that arises due to averaging.The non-metricity in macroscopic gravity partially constitutes cosmological back-reaction.Hence, we have shown here that the back-reaction not only modifies the dynamics but also the kinematics in the averaged Universe.This is expected since in addition to the field equations, the process of averaging in macroscopic gravity also modifies the underlying geometric structure.

Expansion Scalar and Cross Section of Null Bundles
It is useful (especially in cosmology) to show how the expansion scalar is related to the cross sectional area of the congruence.A formal way of doing this can be found in [22].We will follow an analogous treatment here.Consider a cross section around a point P on a curve in the congruence.This can be constructed by taking a set of points, P ′ , on the neighbouring curves that correspond to the same parameter value as that of P. Let this set of points be called, S P .To see the change in cross section, we need to compare S P to a similar set of neighbouring points at some other point, Q, along the same curve.Let us call this S Q .Since the points on S P (and S Q ) are also points on neighbouring curves, they can be used to label the curves, thus giving us a consistent coordinate system, y A (A = 2, 3), on the cross section.Then, one can construct a 2-metric on the cross section as, Using this, it can be shown that (see [22]) the fractional rate of change of the cross sectional area is, Calculating the right hand side of the above equation, we get, Using equation (3.17), this becomes, 1 2 Therefore, the fractional rate of change of the cross sectional area of the null geodesic congruence in the averaged geometry is equal to the Riemannian part of the expansion scalar, An important point to note here is that this is turning out to be the case since we are dealing with geodesic curves (k α ∇ α k β = 0).For non-geodesic curves one should expect contributions from nonmetricity 4 .The fractional rate of change of cross sectional area is closely related to the angular diameter distances in cosmology.The consequence of equation (3.22) is that the non-metricity will contribute to the distances only through the dynamics.

The Geometry of Hypersurfaces
In this section, we will continue with the 1 + 3 decomposition of the space-time presented in the last section.However, here we will present the equations governing the geometric properties of the spatial hypersurfaces.This approach is commonly employed in GR to split the curvature and the field equations in temporal and spatial parts.This leads to the well known Gauss-Codazzi equations and the definitions of the shift and lapse functions.Such a splitting is particularly useful in probing various physical situations in numerical relativity.
Let us consider a spacelike hypersurface normal to a timelike geodesic.The 3-metric, h αβ , in equation (3.2) serves as the metric on this hypersurface.This is sometimes referred to as the induced metric.In the signature adapted in this paper, this is a positive definite non-Riemannian metric.To analyse the geometric properties of the hypersurface, one needs to take the spatial projection of the space-time objects.This is done using the induced metric in a similar manner as equation (3.3).

The Non-Metricity of the Hypersurface
The spatial projection of the covariant derivative is defined as, This serves as the covariant derivative on the hypersurface.Using this definition, we can write, Therefore, we see that the spatial covariant derivative of the 3-metric is not zero.In other words, the spatial connection is not metric compatible 5 .This means that the hypersurface retains the non-Riemannian nature.We define the the non-metricity of the hypersurface as, (3) Similarly, one can evaluate that, 1 2 Note, however, that, (3) ∇ ρ h α β = 0.The extrinsic curvature of the hypersurface is defined as, where, u is a timelike vector normal to the hypersurface i.e., u α K αβ = 0 = K αβ u β .Due to the nonmetricity of the hypersurface, we have, where, K α β = g ρβ K αρ .Since the connection is torsion-less (and hence, so is the connection on the hypersurface), we have, K αβ = K βα .

The Gauss Equation
The Gauss equations arises from evaluating the Ricci identity on the hypersurface.It gives an expression for the Riemannian curvature of the full space-time in terms of its components projected onto the hypersurface.To derive the Gauss equation, we begin with defining the Riemann tensor on the hypersurface in the following manner, where, n is some vector tangential to the hypersurface, i.e., n α u α = 0. Evaluating the left side of the above equation, we get, Using equations (2.17), (4.4), (4.5), and recognising, n ϕ = h ϕ σ n σ , we finally get the Gauss equation, Contracting the appropriate indices, we can write similar equations for the Ricci tensor and Ricci scalar, The above identities are useful in performing a 1 + 3 splitting of the field equations.

The Codazzi Equation
The Gauss equation relates the Riemann curvature tensor with all components projected onto the hypersurface to the extrinsic curvature.Similar to this, one can project only three components of the Riemann tensor onto the hypersurface and one component along a timelike vector, u.This quantity also be written in terms of the extrinsic curvature.This is done through the Codazzi equation.
To derive the Codazzi equation, we start with evaluating, Anti-symmetrising ρ and α, the first term in the brackets on the right hand side vanishes.The second term turns into the Ricci identity (2.17).Taking care of raising and lowering of the indices introduces a few non-metricity terms, finally giving, This is the Codazzi equation.Contracting the indices ρ and β in the above equation, we get, Using the 1 + 3 covariant formalism presented here, it will be straightforward to decompose the four dimensional averaged geometry into foliations and fibrations as is frequently done.The definition of the lapse function and the shift vector remain the same and the 1 + 3 decomposition of the metric retains the same form [24].

FLRW Cosmology in Macroscopic Gravity
The field equations of macroscopic gravity, obtained by averaging the Einstein field equations, contain different traces of the connection correlation and the non-Riemannian part of the curvature tensor.These additional terms in equation (2.12) are the result of a change in the geometric structure of the space-time due to averaging.These terms constitute precisely what is referred to as the cosmological back-reaction in the literature.Considering a particular averaged geometry (i.e., taking an ansatz on the metric of the averaged space-time), one can find the back-reaction in terms of the metric coefficients.This has been previously done for the FLRW cosmological model [25][26][27].Here, we will build upon the analyses in [25][26][27] and discuss how the non-metricity affects the distance measurements in cosmology.

Dynamics and Expansion History
Let us begin by considering the metric at the macroscopic scales is the spatially flat FLRW metric, ds 2 = −dt 2 + a 2 (t) dr 2 + r 2 dθ 2 + r 2 sin 2 θdϕ 2  (5.1) Considering the averaged matter distribution to be a perfect fluid, we have, where, ρ ≡ ρ(t) and p ≡ p(t) are the energy density and pressure for the averaged matter-energy content, respectively; and u α = (1, 0, 0, 0) is the average 4-velocity of the fluid.Then, writing equation (2.12) explicitly gives us the modified Friedmann equations [25,26], where, Q is a constant coming from the back-reaction term and 'dot' represents the derivative with respect to the cosmological time, t.One can immediately notice that the back-reaction term evolves as a curvature.The first equation above can be written in the form, where, H(t) = ȧ a is the Hubble parameter, we have defined, ρ Λ = Λ 8π and ρ Q = − Q 8πa 2 , and we have added the subscript m to differentiate the matter density from other energy densities.Further, defining the associated density parameters as, Ω X = 8πρ X 3H 2 , we can write, Multiplying both sides by a factor of H 2 H 2 0 , we get, where, is the 'non-metricity density' today, and we have used, a = 1 1+z .The above expression gives the expansion history of the universe taking into account the effect due to non-metricity at the macroscopic scales.One can see that depending on the sign of the constant, Q, a spatially flat FLRW Universe may evolve identical to a spatially curved one.

Angular Diameter Distance
In observational cosmology, one deals with large amounts photons, and therefore, it is the averaged photon paths that should be used to calculate distances.In the context of macroscopic gravity, this means the paths of photon bundles (null geodesics) on the macroscopic geometry should be used to determine the distance measures in cosmology.The angular diameter distance, D A , is related to the cross sectional area, A, of null congruence as [28], Using this, one can write, 1 A Using equation (3.23), this becomes, Taking a derivative on both sides, one gets, (5.11) where, in the second equality we have used the Sachs equation (3.18).
The above equation holds for arbitrary space-times.For a spatially flat FLRW space-time, we have, Ω = 0 = Σ, and the 4-velocity of null geodesics is given by [29], Therefore, we see that the non-metricity affects the distance measurements only through its presence in the Hubble parameter.That is, even though the cross sectional area of null congruences only depends on the Riemannian part of the expansion scalar, since photon trajectories are also affected by the Hubble parameter (H), the distance measurements do get modified by the non-metricity.In figure 1, we plot the expansion history and the angular distance to show the modifications caused by small values of the non-metricity parameter.
The effects of the non-metricity (back-reaction) become more pronounced for large redshifts.For small negative values of the non-metricity parameter, the Hubble parameter is becoming larger than its ΛCDM value.This could prove useful in the context of Hubble tension [30] in which the cosmic microwave data suggests a smaller value of the Hubble constant than the local observations.However, it is worth noting that the back-reaction parameter has been previously constrained to have small positive values [31,32].

Summary and Discussion
In macroscopic gravity, the averaging procedure breaks the metricity of geometry, and the cosmological back-reaction can be given in terms of the non-metricity of the space-time.Since the non-metricity results from a change in the geometric structure of the space-time, the effects of back-reaction are not limited to dynamics.In this paper, we analysed the effects of non-metricity on the geometric properties of the averaged universe.The primary goal of this paper was to demonstrate how non-metricity modifies kinematics by deriving the Raychaudhuri and Sachs equations.However, the non-Riemannian nature of space-time raises questions about the trajectories of physical particles.It has been argued that only matter that possesses a 'microstructure' (properties like hypercharge and dilation) couples with the non-Riemannian geometric fields and thus will be affected by them [18].In the metric-affine gravity literature, such matter has been termed as a 'hyperfluid' [33,34].Ordinary matter (e.g.perfect fluid or dust used in cosmology) does not possess such microstructure.Hence, despite the averaged geometry being non-Riemannian, we took the paths traced by physical particles to be simply Riemannian geodesics.
Having clarified this issue, we presented a 1 + 3 decomposition of the averaged space-time geometry within macroscopic gravity.We started with deriving the Raychaudhuri and Sachs equations which govern the expansion of the geodesic flows.Physically, these equations characterise the flow of the cosmic fluid and photon bundles, respectively.Both of these equations were modified due to nonmetricity.However, we found that only the Riemannian part of the expansion scalar was related to the fractional rate of change of cross sectional area of the geodesic bundles.Then, we presented geometric properties of the spatial hypersurfaces orthogonal to these geodesics.Using the definition of the spatial covariant derivative, we showed that the spatial hypersurfaces are also non-Riemannian in nature.This became even more apparent in the expressions for the extrinsic curvature.Using the Ricci identity for hypersurfaces, we derived the Gauss-Codazzi equations.As expected, they showed us that nonmetricity modifies the Riemann curvature projected onto the hypersurfaces.
Finally, we applied the results derived in this paper to observational cosmology.Using the Sachs optical equation, we derived the angular diameter distance formula.We showed that, despite the exis-tence of a non-metricity, the angular diameter distances do not explicitly get modified by it.This was a consequence of the cross sectional area depending only on the Riemannian part of the expansion scalar.Therefore, one needs to be careful when using distances for observational cosmology within macroscopic gravity.Including the averaged Ricci tensor directly in the Sachs equation might result in an incorrect evaluation of the effects of back-reaction on distances in cosmology.However, it is important to note that the non-metricity does modify the distances implicitly through the dynamics (Hubble parameter).We plotted the expansion history and the angular distance formula to demonstrate the deviation from their ΛCDM counterparts for small values of the back-reaction.We note that a small negative value of the back-reaction parameter could help alleviate the Hubble tension.This warrants further investigations into the averaging problem.