Geodesic deviation, Raychaudhuri equation, Newtonian limit, and tidal forces in Weyl-type $f(Q,T)$ gravity

We consider the geodesic deviation equation, describing the relative accelerations of nearby particles, and the Raychaudhuri equation, giving the evolution of the kinematical quantities associated with deformations (expansion, shear and rotation) in the Weyl type $f(Q,T)$ gravity, in which the nonmetricity $Q$ is represented in the standard Weyl form, fully determined by the Weyl vector, while $T$ represents the trace of the matter energy-momentum tensor. The effects of the Weyl geometry and of the extra force induced by the nonmetricity-matter coupling are explicitly taken into account. The Newtonian limit of the theory is investigated, and the generalized Poisson equation, containing correction terms coming from the Weyl geometry, and from the geometry matter coupling, is derived. As a physical application of the geodesic deviation equation the modifications of the tidal forces, due to the nonmetricity-matter coupling, are obtained in the weak field approximation. The tidal motion of test particles is directly influenced by the gradients of the extra force, and of the Weyl vector. As a concrete astrophysical example we obtain the expression of the Roche limit (the orbital distance at which a satellite begins to be tidally torn apart by the body it orbits) in the Weyl type $f(Q,T)$ gravity.


I. INTRODUCTION
The twentieth century has seen the birth of General Relativity (GR) and of Quantum Mechanics (QM), which are considered as the most two successful theories describing the nature and properties of the physical world, on scale ranging from the microscopic to the cosmological one. General Relativity, a geometrical theory of gravity, one of the fundamental interactions that shape the Universe, has provided an excellent description of the observational data [1], and has led to new insights into the problems of space and time, and their relation with the physical Universe. The recent detection [2] of the gravitational waves has again proved the existence of the excellent correspondence between experimental data and the theoretical predictions of GR in a range extending from weak to strong gravitational fields. The introduction of Riemann geometry into GR has provided a powerful mathematical framework to describe the properties of the gravitational field. However, despite its remarkable success, a number of recent observational results have raised some questions about the absolute validity of standard GR, which may still present some limitations especially on astrophysical scales exceeding the Solar System one. The two fundamental problems facing present day gravitational theories, and, in particular General Relativity, are the dark energy and the dark matter problems, respectively. For recent reviews on the dark energy and dark matter problems see [3][4][5][6][7][8][9][10][11][12][13]. To fix the theoretical imperfections of standard General Relativity, two important physical components, the dark matter and the dark energy, are introduced in the cosmological scenario in a rather ad hoc way, with the major goal of explaining and adapting the standard gravity model to describe realistic physical situations, related to the motion of massive particles around galaxies, and to solve the problem of the accelerated expansionary state of Universe. We may call the approach based on the introduction of two new physical components in the overall matter/energy balance of the Universe as the dark components model [14].
However, a second approach to gravitational phenomena is also possible, and it is called the dark gravity approach [14]. In this approach one assumes that both dark matter and dark energy can be explained by changing the nature of the gravitational force. Many modified theory of gravity, going beyond the standard GR model have been proposed in order to build a fundamental framework for the explanation mysterious dark matter and dark energy [15], and to give a solution of the present observational-theoretical contradictions and conflicts. In the framework of the Riemannian geometry one can naturally generalize the Einstein-Hilbert action by substituting the Ricci scalar R with an arbitrary functions f (R). This leads to the f (R) modified theory of gravity [16]- [23]. There are two approaches in f (R) gravity, namely, the metric formulation, in which the metric is considered as the only dynamical variable, and the Palatini formulation, in which the connection is considered as another independent variable, beside the metric tensor. One can find detailed discussions of the metric formulation of f (R) gravity in [24]- [29], and for the Palatini formulation in [30]- [32], respectively. The most obvious drawbacks of f (R) theory is that the scalar field in the Palatini formulation is not dynamical, which implies that no new degrees of freedom can be introduced, resulting into the existence of infinite tidal forces that generally are physically impossible [33]. On the other hand the extra freedom introduced from the metric formulation would lead to contradiction with the observational results obtained in the Solar System [34,35].
To allow for the generation of long-range forces and simultaneously passing the Solar system test, in [36]- [40] a new approach to gravitational effects was proposed. In this theory, called Hybrid Metric-Palatini Gravity, the Einstein-Hilbert action is supplemented with a correction term inspired by the Palatini formulation. Another interesting and important modification of gravity is the inclusion of a non-minimal coupling of geometry and matter into the action [41]- [48], by using arbitrary functions of the scalar curvature and Lagrangian density of matter (in the f (R, L m ) gravity theory [47]), or by considering a gravitational Lagrangian of the form f (R, T ) [48], where T is the trace of the matter energy-momentum tensor. In these classes of theories the covariant derivative of the energy-momentum tensor is always non-zero, which implies a non-geodesic motion of test particles, and the appearance of an extra force. For a recent review of some modified gravity theories in Riemann geometry one can refer to [14] and [49], respectively.
The standard GR theory is formulated in Riemann geometry. Hence, an alternative avenue for searching for a generalized description of gravity is the extensions of the geometrical framework on which GR is based. In an attempt to unify gravity and electromagnetism H. Weyl introduced in 1918 a generalization of the Riemann geometry [50]. In the Weyl geometry both the orientation and the length of vectors are allowed to vary under parallel transport, while in Riemann geometry only the variation of the orientation is allowed. Weyl geometry represents a completely consistent generalization of Riemannian geometry. In modern language, the vector field introduced by Weyl, which generates a new component in the connection, which in Weyl geometry is no longer metric-compatible, is actually the dilatation gauge vector. If the vector is the gradient of some function, a scale transformation of the formg µν = σ 2 g µν , where σ is the scale (conformal) factor, can be applied to cancel the Weyl vector. In this case the Weyl geometry is called integrable, and the length of vectors will be unchanged under a parallel transport along closed paths. For a detailed discussion of Weyl geometry see [51][52][53].
In an important mathematical and physical development E. Cartan introduced the anti-symmetric part of the connection, known as torsion, into the gravity theory, thus formulating an extension of GR [54], which is known as the Einstein-Cartan theory. The Weyl ge-ometry can be immediately generalized by including the torsion, which leads to the Weyl-Cartan geometry, and a corresponding geometric theory of gravitation [55][56][57][58][59][60][61][62][63][64][65][66]. For a review of geometrical properties and physical applications of Riemann-Cartan and Weyl-Cartan spacetimes one can refer to [67]. Another important mathematical development with important physical implications is related to the work by R. Weitzenböck [68], who developed a geometry with torsion and zero Riemann curvature. Since curvature is zero, Weizenböck spaces possess the interesting property of distant parallelism, known as teleparallelism. The teleparallel approach substitutes the metric tensor, which plays a central role as a basic physical variable in gravitational theories, with a set of tetrad vectors. In this approach, the torsion is generated by tetrad fields and it describes the gravitational field entirely, once the torsion is properly chosen to eliminate the curvature. This is the so-called teleparallel equivalent of general relativity (TEGR) [69][70][71], which is also known as f (T ) theory, whereT is the trace of torsion tensor. In f (T ) theory, torsion completely compensates the curvature, and the spacetime becomes flat. Weyl-Cartan-Weizenböck gravity (WCW) was introduced as an extension of the teleparallel gravity models in [72]. In this approach, the Weizenböck condition of the vanishing of the total curvature is implemented in a Weyl-Cartan spacetime. Moreover, in [73] the Weizenböck condition of the exact compensation of torsion and curvature was introduced in the action via the Lagrange multiplier approach, in Riemann-Cartan spacetime. For a review of teleparallel gravity theory see [74].
An important extension of f (Q) theory has been obtained in [97] by including a nonminimal curvaturematter coupling into the gravitational Lagrangian, with L = f 1 (Q) + f 2 (Q)L m , where f 1 and f 2 are arbitrary functions of the nonmetricity Q. Similarly to the nonminimal couplings between curvature and matter in Riemannian geometry [47,48], the coupling leads to the nonconservation of the matter energy-momentum tensor as well, which leads to the existence of a new term in the geodesic equation, which can be interpreted as an extraforce. A Bayesian statistical analysis using redshift space distortions data was performed to test a model of f (Q) gravity in [98]. The cosmological background evolution is similar to the ΛCDM one, but differences arise in the perturbations. The best fit parameters indicate that the σ8 tension between Planck and Large Scale Structure data can be relieved in this theory.
The f (Q) gravity model was extended to include a non-minimal coupling in the Lagrangian in [99], in which the gravitational action L is given by an arbitrary function f of the nonmetricity Q and of the trace of the matter energy-momentum tensor T , with L = f (Q, T ). Several cosmological applications of the theory were considered by chosing some simple functional forms of the function f (Q, T ), corresponding to additive expressions of f (Q, T ) of the form f (Q, T ) = αQ+βT , f (Q, T ) = αQ n +1+βT , and f (Q, T ) = −αQ−βT 2 , respectively, where α, β and n are constants. The Hubble function, the deceleration parameter, and the matter energy density were obtained in each case as a function of the redshift by using analytical and numerical techniques. For all considered cases the Universe experiences an accelerating expansion, ending with a de Sitter type evolution. Gravitational baryogenesis in f (Q, T ) gravity was considered in [100], and it was found that f (Q, T ) gravity can contribute significantly to this phenomenon. The various cosmological parameters in Friedmann-Lemaitre-Robertson-walker (FLRW) geometry have been obtained in [101] for different choices of the function f (Q, T ) in terms of the scale-factor and redshift z by constraining the energy-conservation law.
The observational constraints on the model have been obtained by fitting the model parameters using the available data sets like Hubble data sets H(z), Joint Light Curve Analysis (JLA) data sets and union 2.1 compilation of SNe Ia data sets. The various energy conditions for cosmological models in f (Q, T ) gravity were studied in [102]. The equation of state parameter w = −1 also supports the accelerating behavior of the Universe. In the considered f (Q, T ) models the null, weak, and dominant energy conditions are obeyed, while the strong energy conditions are violated during the present accelerated expansion. The late time cosmology in f (Q, T ) gravity was investigated in [103]. Constraints on the model parameters were imposed from the updated 57 points of Hubble data sets, and 580 points of union 2.1 compilation supernovae data sets. The performed analysis did show that f (Q, T ) gravity represents a promising approach for explaining the current cosmic acceleration, and it can provide a consistent solution to the dark energy problem.
A particular type of f (Q, T ) model was considered in [104], in which the scalar non-metricity Q µν of the spacetime was expressed in its standard Weyl form, and therefore it is fully determined by a vector field w µ . The field equations of the theory have been obtained under the assumption of the vanishing of the total scalar curvature, a condition which was added into the gravitational action via a Lagrange multiplier. The cosmological implications of the theory were also investigated for a flat, homogeneous and isotropic geometry, and the generalized Friedmann equations were obtained. Several cosmological models were investigated by adopting some simple functional forms of the function f (Q, T ), and the predictions of the theory have been compared with the standard ΛCDM model.
The main goal of the present paper is to investigate some fundamental properties of motion of test particles in the Weyl type f (Q, T ) gravity theory, introduced in [104]. We derive the geodesic deviation equation and the Raychaudhuri equation in the Weyl geometry and in the presence of the extra-force that describes the effects of the nonmetricity-matter coupling. As compared to standard general relativity a number of new terms do appear in both equations, indicating the existence of a complex dynamics resulting from the intricate interplay of the geometrical and matter factors. The weak field limit of the theory is also considered, and the Poisson equation, containing the corrections arriving from the Weyl geometry and nonmetricity-matter coupling are determined. As a physical application of the obtained results we consider the tidal force problem in Weyl type f (Q, T ) gravity, and a generalization of the Roche limit is obtained, under the simplifying assumption that the center of mass of the two-body system coincides with the geometrical center of the massive object with mass M .
The present paper is organized as follows. In Section II we review the basics of the Weyl geometry, and we introduce the Weyl type f (Q, T ) gravity theory, and its field equations. The Geodesic deviation equation is derived in Section III, while the Raychaudhuri equation in Weyl geometry and in the presence of the extra force induced by the nonmetricity-matter coupling is obtained in Section IV. In Section V we consider the weak field limit of the theory, and the generalized Poisson equation for the gravitational potential, containing the corrections due to the geometric and coupling effects, is obtained. The properties of the tidal forces, as well as the expression of the Roche limit in Weyl geometry and in the presence of extra force, are also investigated. We discuss and conclude our results in Section VI. Some mathematical results used in the derivation of the main relations of the paper are summarized and detailed in the Appendix A. In the present Section we briefly review the fundamentals of the Weyl geometry, and of the Weyl-type f (Q, T ) gravity.

A. Quick start for Weyl geometry
In the differential geometry of the Riemann spaces the connection describes the properties of the parallel transport, the fundamental concept used for the characterization of the various mathematical aspects associated to the geometrical objects in curved space-times. In order to keep the tensor properties of the differentiation operation acting on a vector, some new terms, containing the connection, must be introduced to assure the tensorial nature of the differentials of a vector. Hence, in Riemann geometry the covariant derivatives of the con-travariant and covariant vectors are given by [105][106][107] In Riemann spaces, the connection is of Christoffel type, it is compatible to the metric g µν , so that ∇ λ g µν = 0, and it is given by With this connection the parallel transportation in a Riemann space will preserve the length of the vectors, and only the orientation of the vector is changed. In 1918 H. Weyl proposed a new geometry [50] by introducing a connection with the property that under parallel transportation both the orientation and the magnitude of a vector change. The connection in Weyl geometry is no longer metric compatible. Moreover, a new vector field, known as Weyl vector field, is introduced, allowing to write the Weyl connection as [52,53], where in the following the tilde indicates the quantities defined in the Weyl geometry. With the use of the Weyl connection one can immediately obtain the fundamental result that the covariant derivative associated with the connectionΓ λ µν , when applied on the metric tensor, gives a non-zero value, In the standard representation of the Weyl geometry the variation of the length of a vector under parallel transport is given by where l 0 is the length of the vector before transportation. Thus the variation of the length of a vector under transport along a closed loop is, where δs µν is the area surrounded by the loop, and W µν is given by In Weyl geometry one can define a curvature tensor in the same way as in Riemann geometry, with the use of the intrinsic connection, The explicit expressions of the curvature tensor in Weyl geometry can be obtained as By contracting the first and the third indices of the curvature tensor, we find Finally, contracting the remaining two indices one can obtain the scalar curvature as In order to develop some physical applications we introduce two types of nonmetricities as follows, Moreover, the scalar nonmetricity Q is defined as where L λ µν is given by Hence we can rewrite the connection in Weyl geometry asΓ With the use of Eqs. (13) and (16), we obtain for the scalar nonmetricity the expression (see Appendix A 1 for the calculational details), The explicit expressions of the covariant derivatives for some vector and tensor quantities in Weyl geometry are given in Appendix A 2.
B. The Weyl-type f (Q, T ) gravity theory A particular model of the general f (Q, T ) theory was introduced in [104], and it is based on the following action, where we have denoted κ 2 = 1/16πG, and m is the mass of the particle associated to the vector field w µ . Moreover, with L m we have denoted the ordinary matter action. In the action (19), the second and third terms are the kinetic energy and the mass terms associated to the Weyl-type vector field, respectively. As for the first term, giving the gravitational Lagrangian, it is taken as an arbitrary function of the nonmetricity Q and of the trace T of the matter energy-momentum tensor.
Moreover, in the present approach, we assume the vanishing of the total scalar curvature, and we impose in the total action the Weizenböck condition with the help of the Lagrange multiplier λ. We impose this condition in order to follow the essential ideas of the teleparallel approach to gravity. The action (19) generalizes the GR equivalent symmetric teleparallel gravity theory, or f (Q) theory, by introducing a matter-geometry coupling, in a similar manner as in the f (R, T ) theory, and adding an extra matter distribution for the particles associated with the Weyl vector field. However, it should be mentioned that the present theory is different from f (R, T ) theory due to the presence of a boundary term, and hence it can be seen as a generalized equivalent theory of GR, in the sense discussed in [88].
In the case where f (Q, T ) = Q, and if we have a vanishing mass m = 0, the theory reduces to the symmetric telleparallel equivalent to GR, as is fully developed in [76].
By first varying the total action (19) with respect to the vector field w µ , we obtain the generalized Proca equation that describes the evolution of the vector field w µ in the Weyl geometry, where f Q = ∂f (Q, T )/∂Q. Next, we vary the action with respect to the metric tensor, and with the use of the Weizenböck condition, that is, after omitting the terms that explicitly containR, we obtain the generalized field equation of the Weyl-type f (Q, T ) gravity as follows, where in the field equations we have denoted by f T the partial derivative of f with respect to T , f T = ∂f /∂T , and where we have introduced the energy-momentum tensor S µν associated to the Weyl type vector field, defined according to As customary, the energy-momentum tensor T µν of the ordinary matter is defined as The tensor Θ µν is obtained from the variation of the energy-momentum tensor with respect to the metric. By adopting the assumption according to which the ordinary matter Lagrangian density L m depends only on the metric tensor components, and not on their derivatives, we obtain for Θ µν the simple expression, By contracting the indices µ and ν in Eq. (21), we obtain the scalar S, given by The trace of the matter energy-momentum tensor T can be obtained from the equation, By taking the covariant divergence of Eq. (21), and using the constraintR = 0, we obtain for the divergence of the ordinary matter energy-momentum tensor the expressioñ Hence, as one can see from the above equation, the ordinary matter energy-momentum tensor is not conserved in the Weyl-type f (Q, T ) theory. From a physical point of view the nonconservation of the matter energymomentum tensor can be interpreted as indicating the presence of an extra-force, acting on massive test particles, and making the motion nongeodesic. It should be noted that in the special case f T = 0, the energymomentum tensor becomes conserved, as one can be easily seem from Eq. (21). In this case the term in the action generated by the non-minimal coupling between matter and geometry does vanish, and the theory reduces to f (Q) coincident gravity [76].
There is a special case where the denominator of Eq. (27) vanishes identically, and therefore the field equations Eq. ( (21)) are not valid. This particular case corresponds to 1 + 2κ 2 f T = 0, giving immediately is an arbitrary function depending only on the scalar nometricity. For this form of f (Q, T ) the gravitational field equations Eqs. (21) become, where C Q = ∂C(Q)/∂Q. The conservation equation of the matter energy-momentum tensor reduces to the simple form As a first step in our further investigations of the geometric and physical properties of the Weyl type f (Q, T ) gravity, in the present Section we derive the geodesic deviation equation in this theory, and obtain its form by taking into account the presence of the extra-force generated by the nonconservation of the energy-momentum tensor.
A. The extra force in Weyl-type f (Q, T ) theory As we have already mentioned, due to the presence of nonmetricity, in Weyl geometry the length of a vector is no longer preserved under parallel transport. However, we define the four-velocity, as usual, according to the expression u µ = dx µ /dξ, where ξ is the affine parameter, which forms a tangent bundle of congruences of timelike curves. But in a Weyl geometry we need to normalize the four-velocity according to [108], where ℓ(x α ) is an arbitrary function of space and time coordinates. With this normalization, even if we set the proper time as dτ 2 = −g µν dx µ dx ν , ξ does not necessarily coincide with τ [108], and by chain-rule of differentiation we obtain the relation dτ /dξ = ±ℓ.
The presence of nonmetricity, as well as its specific properties related to the background spacetime strongly affects the nature of the hypersurfaces orthogonal to the timelike u µ velocity field. Therefore, in Weyl geometry, a generalized projection tensor operator needs to be introduced by the definition [108] The generalized projection operator has the usual properties h µν = h νµ , and h µν h µν = 3, respectively. The mixed tensor form of the projection operator can be obtained by raising one of its indices, and thus For mathematical convenience, here onwards we denote the temporal derivative by ′ , to indicate covariant differentiation with respect to ξ, and the spatial derivative D µ for a generalized tensor in Weyl geometry, respectively, as follows [108], Due to the presence of nonmetricity in Weyl geometry there exist two types of four-acceleration, denoted by A µ and a µ , respectively, and defined according to and respectively. The two accelerations are related by the important equation, From these two types of acceleration we can obtain some pure geometrical relations. Multiplying the two acceleration vectors by u µ , we obtain For the sake of clarity we mention that (ℓ 2 ) ′ = u µ∇ µ (ℓ 2 ). With the above definitions, the geodesic equation in Weyl geometry takes the following form, where generally the extra force f µ gives the supplementary acceleration induced by the nonminimal curvature-matter coupling, The energy-momentum tensor of a perfect fluid can be defined in a Weyl spacetime according to [108] T µν = (p + ρ) ℓ 2 u µ u ν + pg µν .
The four-velocity u µ of the perfect fluid must be ξparameterized according to the definition of velocity in Weyl geometry, and it obeys the normalization condition (30). The details of the derivation of the generalization of the perfect fluid model in Weyl geometry can be found in Appendix A 3. Now multiplying with the projection operator the covariant divergence of the energy-momentum tensor in the Weyl-type f (Q, T ) theory, we obtain the important correspondence between the matter energy-momentum tensor T µν and the four-acceleration A µ as With the use of Eqs. (35) and (40), from Eq. (42) we can obtain the expression of the extra force in Weyl-type f (Q, T ) gravity as follows (for the calculational details see Appendix A 4), By substituting the expression of the nonmetricity tensor in the above equation, we obtain for the extra force the final expression B. The geodesic deviation equation Consider now a one-parameter congruence of curves x µ (ξ; σ) satisfying the geodesic equation for the parameter ξ in Weyl geometry for each σ = σ 0 = constant. We introduce the four-vectors The second temporal derivative of U µ can be expressed as Using the definition of curvature tensor, we have Hence we obtain the geodesic deviation equation in Weyl geometry in the presence of an extra force as, By multiplying with δσ both sides of the above equation gives Eq. (49) represents the geodesic deviation equation in Weyl geometry in the presence of an extra force generated by the coupling between nonmetricity and matter, as follows from the mathematical and physical structure of f (Q, T ) gravity. In the special case when f T = 0, the energy-momentum tensor becomes conserved, and Eq. (49) reduces to the geodesic deviation equation in Weyl geometry, as considered in [108]. Note that in this case the present theory becomes equivalent to the coincident gravity theory [76]. In the more general case where the Weyl-vector also vanishes, the space-time becomes flat and the right hand side of Eq. (49) identically vanishes.

IV. GENERALIZED RAYCHAUDHURI EQUATION IN WEYL-TYPE f (Q, T ) THEORY
In the present Section we derive the Raychaudhuri equation in the Weyl-type f (Q, T ) theory, and in the presence of an extra force. Our approach basically follows, and generalizes, the similar analysis performed in [108], in which the presence of the extra-force has not been taken into account.
We begin our investigation by decomposing the covariant derivative of the four-velocity u µ into its temporal and spatial components, according to [108], where D ν u µ = h β ν h λ µ∇β u λ and ξ µ = u ν∇ µ u ν = u ν ∇ µ u ν − ℓ 2 w µ . We can see that by construction we have the relation ξ µ u µ = a µ u µ .
Similarly, the projected covariant derivative can be decomposed into where is the "volume" scalar in Weyl geometry. We can also construct the scalar D µ u µ = h µν∇ ν u µ . In a Weyl geometry and in the presence of nonmetricity the shear tensor is defined according to while the vorticity tensor is introduced based on the definition The shear tensor is symmetric and trace-free due to its construction. As for the vorticity tensor, it is naturally anti-symmetric. Thus, based on these mathematical properties, we find that the following relations are always satisfied, σ µ µ = 0 = ω µ µ , σ µν u ν = 0 = ω µν u ν , σ µν h µν = 0 = ω µν h µν and σ µν ω µν = 0.
From Eqs. (50) and (51) we obtain the explicit expression of the covariant derivative of u µ as (55) while for the covariant derivative of u µ we find Raychaudhuri's equation is a purely geometrical relation, and it immediately follows from a set of fundamental geometric relations, known as Ricci's identities. In a Weyl spacetime and in the presence of nonmetricity, the definition of the curvature tensor is, We multiply both sides of Eq.(57) with g λν u µ , thus obtaining The first term on the left hand side of Eq. (58) can be evaluated as while the second term on the left hand side of Eq. (58) evaluates to Hence in Weyl geometry the left hand side of Eq. (58) finally reads, The right hand side of the above equation can be evaluated as With the use of Eqs. (61) and (62) we finally obtain the Raychaudhuri equation in Weyl geometry, and in the presence of an extra force, as, It is worth mentioning again that the effect of the non-minimal matter-geometry coupling f (Q, T ) will enter into the Raychaudhury equation (63) through the expression of the extra force, given by Eq. (44). As we have already mentioned earlier, in the case of minimal matter-geometry coupling f T = 0, the above equation reduces to the generalized Raychaudhuri equation in the coincidence gravity theory [76]. As is well-known, the coincidence gravity is a generalization of the symmetric teleparallel gravity [76]. In this sense, Eq. (63) with f T = 0 could be considered as the Raychaudhuri equation of the generalized symmetric teleparallel equivalent to GR.
The special case with zero extra force of the Raychaudhuri equation has been obtained and discussed in [108]. The first three lines in Eq. (63) have analogue forms to the similar terms in modified gravity theories with geometry-matter coupling formulated in Riemann geometry [109]. The Raychaudhuri and optical equations for null geodesic congruences with torsion were investigated in [110].

V. WEAK FIELD APPROXIMATION, NEWTONIAN AND POST-NEWTONIAN LI,ITS, TIDAL FORCE AND ROCHE RADIUS IN WEYL-TYPE f (Q, T ) GRAVITY THEORY
In the present Section, based on the previous mathematical results, we will consider the weak field limit of the Weyl-type f (Q, T ) gravity theory. The generalized Poisson equation, as well as the expressions of the tidal force tensor are derived. As an astrophysical application of the obtained results we consider the modifications to the Roche limit that are induced by the assumption of the presence of nonmetricity in the space-time geometry.

A. The weak field approximation
If one considers the physical situation in which the motion of the test particles is slow, and that the gravitational field intensity created by the material particles is comparably weak, one could easily compare generalized metric theories of gravity with each other, with the experimental observations, as well as with Newtonian gravity. In this case, the first order approximation is adequately accurate to compare the theoretical predictions of the gravitational theories with past, present and future Solar System observations. Generally, this kind of approximation of gravitational theories, the so-called post-Newtonian limit, is valid within the near-zone of the system, corresponding to a spherical region with size smaller than one gravitational wavelength.
As a first step in our investigation of the physical properties of the f (Q, T ) gravity theory, in the following we investigate the linear approximation of the metric field, or the weak field approximation, by assuming the decomposition of the metric tensor as where η µν = diag(−1, 1, 1, 1) is the metric tensor in the Minkowski spacetime. In Weyl geometry, by keeping only the first order of H µν , the connection and curvature tensor can be respectively formed from their Riemann part [106], to which we add the non-Riemannian (Weyl) components of the connectioñ and of the curvature tensor, respectively. In the following the Latin letters denote the spatial components (i = 1, 2, 3) of the tensors, and we will use Greek letters to indicate both spatial and temporal components (µ = 0, 1, 2, 3). By assuming the metric is diagonal and timeindependent, indicating that H µν is also a static field, one can obtain the explicit weak-field expressions of the Christoffel part in curvature tensor as [106,107], From Eq. (66) we obtain the component of R γ 0β0 , and For the details of the calculations of the Riemann and Ricci tensors in the weak field approximation see Appendix A 7.

B. Post-Newtonian analysis
In this Section we will discuss in detail the Newtonian limit of f (Q, T ) gravity, by assuming that the perfect fluid filling the space-time is nonrelativistic. First of all, we will chose a reference frame where the motion of the fluid is static. Moreover, we assume that matter is in the form or pressureless dust with the property ρ ≫ p, and that the velocity of particles is small as compared to the speed of light.
Hence we can approximate the four-velocity of the fluid as u µ = (ℓ, u i ), and u µ = (−ℓ, u i ), respectively, and we keep only the first order terms in u i . In the framework of these approximations the matter energy-momentum tensor T µν reads, Assuming that λ = λ 0 + δλ and ℓ = 1 + δℓ, where δλ and δℓ are perturbation variables, from Eq.(5) and Eq. (20), we obtain the first-order perturbation of δλ and δℓ as Using the above equation together with equation (71), one finds that the only non-vanishing component of the energ-momentum tensor is T 0 0 = ρ. It should be noted that in the Newtonian limit, one decompose the tensor H µν as where Φ and Ψ are Newtonian potentials. Also we assume that the Weyl vector can be written as From these definitions, one can see that the tensor S µν in (22) becomes second order in perturbation variables and does not contribute to the Newtonian limit of the theory.
Considering the metric field equations (21) at background level, it follows that the background value of the Lagrange multiplier λ 0 should be constant, and also the condition f (0) = 0, where by (0) we mean the background value, must hold.
By expanding the field equation (21) up to first order in the perturbation parameters, and by using the background constraints derived above, one can obtain the (00) and (ii) components of the field equations as Also, from the off-diagonal components of the metric field equation (21), we obtain, In order to close the system of dynamical equations, one should also consider the constraint equationR = 0, and also the vector field equation (20) to first order in perturbations. The constraint equation becomes and the divergence of the vector field equation becomes Now, by substituting ∇ i ω i and ∆δλ from equations (77) and (76) in the (ii) component of the metric field equation (75), one can see that f (0) T ≈ 0. Now, solving the remaining equations (74) and (78) for the Newtonian potentials, we find where we have defined the generalized Newtonian gravitational constant G ef f as and the PPN γ parameter as One can see from its definition that the generalized Newton gravitational constant depends on the derivative of the function f (Q, T ) with respect to the trace of the energy-momentum tensor T , on the mass of the Weyl vector field, and on the background value of the Lagrange multiplier.
It is worth mentioning that the value γ = 1, which corresponds to the GR result, and which is confirmed by the observations at the Solar System level, occurs when f (0) Q = −λ 0 κ 2 − m 2 /12κ 2 . In this case we obtain λ 0 = 1/2G, where G is the Newtonian constant.
The present analysis does show that the effects of the matter geometry couplings f (Q, T ) do appear already in the first-order perturbation of the theory. This is in fact different from the generalized teleparallel theory, where the first-order perturbation analysis does not reveal the extra degrees of freedom, and one should take into account higher order perturbation analysis [117]. However, in order to explore the detailed dependence of the function f (Q, T ) on the trace of the energy-momentum tensor T , one should consider higher order perturbation analysis of the model.
The extra force (44), can be decomposed into an energy-momentum related component, and a geometry related component, where we have defined the tensor which can be written explicitly with the use of Eq. (27) as, where u µ F µ = 0 should be always fulfilled. Using the definition of the energy-momentum tensor and also expanding the above expression to first order in perturbation variables, one obtains where we have assumed that the energy density varies slowly in this limit. At first order in perturbations one obtains, F 0 = 0. Also note the second term in Eq. (86) is non-vanishing only in the case of ρ = 0. Since in the following calculations we will use only the i components, we will omit this term. In the particular case of 1 + 2κ 2 f T = 0, f T = −1/2κ 2 , we find lim fT →−1/2κ 2 ∇ ν f T / 1 + 2κ 2 f T = 0, and thus we still obtain for the extra-force the non-trivial expression F ρ = −ℓ 2 ω ρ + ℓu ρ ω 0 . The equation describing a world line in Weyl geometry with extra force reads, By adopting the linear approximation and the Newtonian limit, the spatial component in Eq. (87) becomes This equation represents the generalization of the Newtonian equation of motion in the Weyl-type f (Q, T ) gravity theory.
C. Tidal forces in Weyl type f (Q, T ) gravity theory In a Weyl geometries, the properties of parallel transport of a vector along a geodesic line are preserved with the important exception of the magnitude of a vector changing after transport. Similarly to Riemannian geometry, in the Weyl geometry one can always find a set of tangent spaces generated by the four velocity u µ for every point along a world line, where their axes (with the same index) remain parallel, under parallel transport along the world line. It should be noted that these four velocities are not necessarily normalized to one in every tangent space. We have already defined u µ u µ = −ℓ 2 , and in this case we will not take ℓ to one. Here onward we will consider geodesic reference frame in which all connection components vanish, indicating thatΓ λ µν = 0 point wisely.
We at first go back to Eq. (66) to evaluate the spatial total curvature tensorR i0j0 , which is of significant important in the following part. Hence we obtaiñ where Γ λ µν = −g µν w λ + δ λ µ w ν + δ λ ν w µ is obtained from the condition of zero total connection. Raising the index i, we haveR Hence the geodesic deviation equation Eq. (49) will become, If we consider the linear approximation, the Newtonian limit, and the zero total connection for this system, and we also use η 0 = 0 to indicate that the accelerations of the particles are compared at equal times, and f 0 = 0, respectively, indicating that the thermodynamic parameters of the matter do not depend on time, we obtain Now we write the explicit expression of the curvature tensor in Weyl geometry by using Eq. (90) as, and hence we can reformulate Eq. (93) by introducing a newly defined tidal force vector F i , and the tidal matrix K i j , which has been modified by the matter-curvature coupling and Weyl geometry, with the explicit expression of K i j given by where F i is the component of extra force, and has been defined in Eq. (86). The contraction of the tidal matrix gives, and by using Eq. (69), (70, and (73), we obtain and the scalar K reads, Substituting the expression (86) for F i , one obtains and respectively.

D. The Roche radius in Weyl type f (Q, T ) gravity
In Newtonian gravity, the spherical potential of a given particle with mass M is given by where we have assumed that 8πG = 1. In a frame of reference with the x axis passing through the particle's position, indicating that the particle is located at (x = r, y = 0, z = 0), the Newtonian tidal force tensor τ ij will be diagonal, and has only the following nonzero components [109], The components of the Newtonian tidal force F i can be written as respectively [107]. The Roche limit, an important astrophysical and astronomical concept, is defined as the closest distance r Roche that a cosmic object, having mass m, radius R m , and density ρ m , respectively, can come near a massive star of mass M , radius R M and density ρ M , respectively, without being torn apart by the tidal gravity of the star. In the following we consider a simplified case with M ≫ m, a condition which allows us to set the center of mass in the geometrical center of the mass M .
We consider a small object of mass ∆m located at the surface of the small body of mass m. The gravitation force from the small mass acting on ∆m is given by while the tidal force from the big massive body acting on ∆m is obtained as where r is the distance between the centers of the two celestial objects, and we have neglected the differences in the distances between ∆m and M , and ∆M and m, respectively. The Roche limit is reached if the two forces acting on ∆m are equal, F G = F . Thus we obtain the Roche limit in Newtonian gravity r Roche as With the use of Eq. (88) we obtain the modification of the gravitational force in Weyl-type f (Q, T ) gravity, which can be represented as where the index r indicates the radial components, F gravity is the Newtonian gravitational force, F geometry gives by modifications from geometry components beyond Riemann geometry, and F Extraf orce is the component generated by the geometry-matter coupling. Thus, in Weyl geometry, by using Eq. (97), the Roche limit r Roche is obtained as and the vectors containing index r (no summation upon r) must be evaluated in the coordinate system in which the Newtonian tidal tensor is diagonal. Considering that the gravitational effects due to the coupling between matter and curvature are small as compared to the Newtonian ones, we have Substituting F r from Eq. (86) one finally obtains

VI. DISCUSSIONS AND FINAL REMARKS
Abandoning the metricity conditions, and including nonminimal curvature-matter couplings are some promising ways to modify standard general relativity, and to explain the major challenges present day gravity theories face. A possible geometric avenue for the generalization of general relativity is represented by the so-called symmetric teleparallel gravity theory [75], and by its extensions [76,97,98]. In particular, the role of matter and of the geometry-matter couplings have been analyzed in [99] and [104], respectively. In the present paper we have extended the previous analyses of the Weyl type f (Q, T ) gravity, a particular version of the general f (Q, T ) type theories, by developing some basic theoretical tools that would allow not only to further investigations of the fundamental geometrical and physical properties of these gravity theories, but can also open the possibility of their observational testing.
More exactly, from the fundamental point of view of the analysis of the Weyl type f (Q, T ) theories, we have obtained two of the basic equations of the gravitational physics, namely, the geodesic deviation equation, and the Raychaudhuri equation, respectively. The geodesic deviation equation geodesic describes the way objects approach or recede from one another when moving under the influence of a spatially varying gravitational field. One of the important applications of the geodesic deviation equation is in the study of the tidal forces, which in modified theories of gravity acquire some extra terms due to the presence of the new terms that modify the gravitational interaction. Hence the geodesic equation can be used to observationally test the Weyl type f (Q, T ) gravity model through the observations of the effects of the tides produced by an extended mass distribution. Tidal effects play an important role in the eccentric inspiralling neutron star binaries [111]. The neutron stars can be modelled as a compressible ellipsoid, which can deform nonlinearly due to tidal forces, while the orbit evolution can usually be described with the post-Newtonian theory. In general, the tidal interaction can accelerate the inspiral, and cause orbital frequency and phase shifts. Tidal interactions have an essential effect on the star formation in galaxies, since tidal perturbations induced by close companions increase the gas accretion rates [112]. By using gravitational wave detector networks one can constrain the equation of state of binary neutron-stars, and extract their redshifts through the imprints of tidal effects in the gravitational waveforms [113]. The existence of light, fundamental bosonic fields is an attractive possibility that can be tested via black hole observations. The effect of a tidal field caused by a companion star or black hole on the evolution of superradiant scalar-field states around spinning black holes can test the existence of light bosonic fields [114]. For large tidal fields the scalar condensates are disrupted, and the impact of tides can be relevant for known black-hole systems such as the one at the center of our galaxy or the Cygnus X-1 system. The companion of Cygnus X-1 will disrupt possible scalar structures around the black hole for large gravitational couplings. Tidal effects in massless scalar-tensor theories were considered in [115], where a new class of scalar-type tidal Love numbers. It turns out that in a system dominated by dipolar emission, tidal effects may be detectable by LISA or third generation gravitational wave detectors. Another astrophysical situations in which the effects of the tides are of major importance are perturbations of the Oort cloud by the galactic field, globular clusters evolving under the influence of the galactic mass distribution, and galactic encounters [116]. As we have seen in our analysis of the F (q, t) gravity, the curvature-matter coupling significantly modifies the nature of the tidal forces, as well as the equation of motion in the Newtonian limit. Therefore, the comparison of the theoretical predictions of the Weyl type f (Q, T ) gravity about the modifications of the tidal forces with the observational evidences, coming from a large class of astrophysical phenomena could give, at least in principle, some insights into the fundamental aspects of the gravitational interaction, and its geometric description.
We have also obtained the generalization of the Poisson equation, describing the properties of the gravitational potential. The Poisson equation, and its solution, is an important tool in the investigation of many gravitational effects involving small velocities and low matter densities. The modifications of the gravitational potential, and the new terms appearing in the equation may provide a theoretical explanation for the observed dynamics of the particles moving on circular orbits around galaxies. These observations are usually explained by postulating the existence of dark matter, a mysterious major component of the Universe, which has not been detected yet. Hence the novel geometric effects induced by the Weyl-type f (Q, T ) gravity may provide a geometric explanation for the galactic dynamics of test particles without having to resort to the dark matter hypothesis. In present paper we have obtained the equation of motion of the particles in the Weyl-type f (Q, T ) gravity, and we have discussed it in detail. Note that the extra force has two components of different origins, coming from the matter distribution and from the geometrical properties of the space-time. The geometry provides an extra degree of freedom, with the nonminimal curvature-matter coupling also generating more degrees of freedom for the gravitational interaction. Hence these extra degrees of freedom contribute with new terms to the extra force, the tidal force, and the Roche limit. These extra terms may have observational (and even experimental) effects, which can be used to test the theoretical gravity model we have investigated in this paper.
The Raychaudhuri equation is of major importance in the investigation of the space-time singularities, and in construction of cosmological models. For the sake of completeness we briefly mention some cosmological applications of our results. Let's consider a flat Friedmann-Lemaitre-Robertson-Walker Universe, with metric given by ds 2 = −dt 2 + a 2 (t)δ ij dx i dx j , where a(t) is the scale factor. We take ℓ = 1, and we adopt a co-moving reference system, with u µ = (1, 0, 0, 0). In this case the prime operator is given by ′ = u µ∇ µ = u 0∇ 0 , and, when applied on a scalar, we have ′ = d/dτ = d/dt. We also introduce the Hubble function H = (1/a(t))da(t)/dt, describing the rate of change of a(t) with respect to time.
We consider a general model for the cosmological nonmetricity, which was introduced in [118], and according to which Q λµν = A(t)u λ h µν +B(t)h λ(µ u ν) +C(t)u λ u µ u ν , where A, B, C are time-dependent functions representing the behaviour of non-Riemannian degrees of freedom and h µν is the projection tensor previously defined. With the Weyl vector defined in Eq. (13), only one extra degrees of freedom is added by the nonmetricity, By assuming A = −C, and by taking B = 0, we obtain for the cosmological nonmetricity the expression Q λµν = Au λ g µν = −2w λ g µν , giving for the Weyl vector the simple expression By taking advantage of Eq.(83), we obtain for the extra force the expression, We also notice that under conditions that we took formerly in this approach, from a geometrical point of view we have, leading to F ρ = 0. This condition is a direct consequence of the use of the co-moving reference frame. In the laboratory reference system the extra-force does not vanish. If we keep the terms containing F ρ , the Raychaudhuri equation in the Weyl spacetime, given by Eq. (63), can now be written as (see Appendix A 8 for the calculational details), and it can also be reformulated into a equation of a(t), which represents the generalized cosmic acceleration equation. Considering the fact that F ρ = 0 in the adopted coordinate system, and after substituting the expression ofR µν u µ u ν , one can see that the above equation is identically satisfied. The most general form of the acceleration equation in the presence of torsion and nonmetricity has been obtained, in a Friedmann-lemaitre-Robertson-Walker geometry, with the help of the generalized Raychaudhuri equation with torsion and nonmetricity, in [118]. Cosmological Hyperfluids, representing fluids with intrinsic hypermomentum that induce spacetime torsion and non-metricity, were studied in [119], where the most general form of the Friedmann equations with torsion and non-metricity were also obtained.
To conclude, the present investigation of some fundamental aspects of the Weyl type f (Q, T ) gravity opens further possibilities for the theoretical, observational and even experimental study of the alternative purely geometrical theories of gravity in the presence of the coupling between geometry and matter. Moreover, the results obtained in the present paper may also be relevant for other classes of modified gravity theories.

The explicit expressions of the covariant derivative in Weyl geometry
We will present below the explicit expression of the covariant derivative in Weyl geometry of certain types of tensors. For vectors, For second order tensors we obtain,

The perfect fluid model in Weyl geometry
We should first notice that since the magnitude of the four-velocity in Weyl geometry is not preserved, and u µ u µ = −ℓ 2 , the perfect fluid model must be generalized to include the effects ℓ. Generally, a perfect fluid can be characterized by its four-velocity u, and the thermodynamic quantities -the proper density ρ, the isotropic pressure p, the temperature T , and the specific entropy s, or the specific enthalpy ω = (p + ρ)/n [106] (in this part ω is defined independently from the main body). Here n is the conserved baryon number density, and n does not change its magnitude during parallel transport. We also introduce the particle number density four-vector n µ , defined as Consequently, n = g µν n µ n ν ℓ 2 g , where g is the determinant of the matrix g µν . Next we need to introduce the matter Lagrangian, which we assume as depending only on the energy density scalar the in local rest frame of the fluid, The energy-momentum tensor of the fluid is given by T µν = −2∂L m /∂g µν + g µν L m under the constraints [106] δs = 0, δn µ = 0. (A11) From the thermodynamical relation, we immediately obtain δρ = ωδn. Using Eqs. (A8,A9,A12), we obtain δn = 1 2n n µ n ν ℓ 2 g δg µν − n µ n ν g µν ℓ 2 g 2 δg = n 2 − u µ u ν ℓ 2 δg µν + u µ u µ ℓ 2 g δg .
The derivative of the matter Lagrangian with respect to the metric tensor is given by and hence we finally obtain the energy-momentum tensor of the perfect fluid model in Weyl geometry as,