f(Q, T) gravity

We propose an extension of the symmetric teleparallel gravity, in which the gravitational action L is given by an arbitrary function f of the non-metricity Q and of the trace of the matter-energy-momentum tensor T, so that L=f(Q,T)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L=f(Q,T)$$\end{document}. The field equations of the theory are obtained by varying the gravitational action with respect to both metric and connection. The covariant divergence of the field equations is obtained, with the geometry–matter coupling leading to the nonconservation of the energy-momentum tensor. We investigate the cosmological implications of the theory, and we obtain the cosmological evolution equations for a flat, homogeneous and isotropic geometry, which generalize the Friedmann equations of general relativity. We consider several cosmological models by imposing 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\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f(Q,T)=\alpha Q+\beta T$$\end{document}, f(Q,T)=αQn+1+βT\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f(Q,T)=\alpha Q^{n+1}+\beta T$$\end{document}, and f(Q,T)=-αQ-βT2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f(Q,T)=-\alpha Q-\beta T^2$$\end{document}. The Hubble function, the deceleration parameter, and the matter-energy density are obtained 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. The theoretical predictions are also compared with the results of the standard Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda $$\end{document}CDM model.


I. INTRODUCTION
The development of gravitational theories closely followed the advances in differential geometry. In all geometric descriptions of gravity it is assumed, following [1], that the space-time is endowed with a metric structure in a general space based on the element of arc ds = F x 1 , ..., x n ; dx 1 , ..., dx n , where F (x; y) is a positive (for y = 0) function defined on the tangent bundle T M . Moreover, it is generally assumed that F is homogeneous of degree one in y [2]. An important special case is represented by the choice F 2 = g µν dx µ dx ν , with the corresponding geometry called generally Riemannian geometry. Riemannian geometry lays at the foundations of general relativity [3][4][5], a geometric theory of gravity, which has become, together with quantum mechanics, one of the cornerstones of present day physics. General relativity is presently considered to be the most successful theory of gravity ever proposed. Its remarkable predictions on the perihelion advance of Mercury, on the deflection of light by the Sun, gravitational redshift [6], or radar echo delay [7,8] have been confirmed observationally at an unparalleled level of accuracy. Moreover, predictions such as the orbital decay of the Hulse-Taylor binary pulsar, due to gravitational -wave damping, have also fully confirmed the obser-vationally weak-field validity of the theory [9]. For a review on the experimental and observational tests of general relativity see [10]. The detection of the gravitational waves [11] did give the opportunity to evaluate the predictions of general relativity in the final stages of binary black hole coalescence, corresponding to the limiting case of the strong gravitational fields.
On the other hand recent observational advances in cosmology have provided strong evidence that recently our Universe did enter in an accelerated expansion phase [12][13][14][15][16][17][18][19]. Moreover, the same observations indicate the surprising result that around 95 -96% of the content of the Universe is in the form of two mysterious components, called dark energy and dark matter, respectively, with only about 4 -5% of the total composition represented by baryonic matter [20,21]. These observations have shown the limitations of standard general relativity, which despite its important achievements, and its remarkable success at the Solar System scale, may not be adequate to fully explain gravitational phenomena on galactic and cosmological ranges. Hence standard general relativity may not be the ultimate theory of the gravitational force, since it cannot give satisfactory explanations to the two fundamental problems present day cosmology is confronted with: the dark matter problem and the dark energy problem, respectively. Moreover, since Einstein's standard theory predicts the existence of space-time singularities in the Big Bang and inside black holes, general relativity represents an incomplete physical model. To solve the singularity problem a consistent prolongation of general relativity into the quantum domain is probably needed.
To explain the observational results of cosmology many different approaches at the classical level have been proposed recently. However, a satisfactory theory of gravity has yet to be found. One possibility to build new gravitational theories is to assume that at large scales the Einstein gravity model of general relativity breaks down, and a more general action than the standard Hilbert one, given by S = R/2κ 2 + L m √ −gd 4 x, where R is the Ricci scalar, κ is the gravitational coupling constant, and √ −g is the determinant of the metric tensor, respectively, describes the gravitational field. An important theoretical direction of study is represented by approaches in which, by keeping the geometrical background as strictly Riemannian, the standard Hilbert-Einstein action is replaced by a more general action. One of the simplest possibilities of extending Einstein's gravity is to introduce an arbitrary function f of the Ricci scalar R into the gravitational action [22,23], which thus becomes S = f (R)/2κ 2 + L m √ −gd 4 x. In this framework, a geometric solution to the dark matter problem can also be obtained [24]. A second approach to extend the Hilbert-Einstein action is to assume the existence of a nonminimal coupling between geometry and matter. This direction of research leads to distinct classes of gravitational theories, called f (R, L m ) gravity [25][26][27][28], with action given by S = f (R, L m ) √ −gd 4 x, and to the f (R, T ) gravity theory [29], with action given by S = f (R, T ) √ −gd 4 x, where T is the trace of the energy-momentum tensor, respectively. Another theoretical approach, called hybrid metric-Palatini gravity, and which combines both the metric and Palatini formalisms of modified gravity theories was proposed in [30,31] to construct a new type of gravitational Lagrangian. For extensive reviews and discussions of the modified gravity theories and of their implications see [32][33][34][35][36][37][38][39][40][41][42][43][44].
The properties as well as the astrophysical and cosmological implications of the f (R, T ) gravity theory have been investigated in detail . An interesting feature of the theory is its possible interpretation as an effective description of some quantum gravity phenomena. y. As suggested in [45], by adopting a nonperturbative approach for the quantization of the gravitational metric, proposed [46][47][48], a particular type of f (R, T ) gravity naturally comes up due to the quantum fluctuations of the metric tensor, with the action given by This interesting theoretical result may imply the existence of a deep connection between the quantum field theoretical description of the gravitational interaction in curved backgrounds, which automatically involves particle creation in the gravitational field, and the corresponding effective classical description within the framework of the f (R, T ) gravity theory [48]. Since general relativity is basically a geometric theory, formulated in the Riemann metrical space, a second promising approach for obtaining generalized theories of gravity consists in looking for more general geometric structures that could describe the gravitational field. Hence more general geometries than the Riemannian one, which may be valid at the Solar System level only, may provide an explanation of the behavior at large cosmological scales of the matter in the Universe.
The first attempt to create a more general geometry than the Riemannian one is due to Weyl [78], which is a classic example of the fruitful interplay between mathematics and physics. The main goal of the study by Weyl was to obtain a geometrical unification of electromagnetism and gravitation. The fundamental concept in Riemann geometry is the metriccompatible Levi-Civita connection, which allows the comparison of lengths. Weyl did replace the metric field by the class of all conformally equivalent metrics, and he did introduce a connection that does not contain any information about the length of a vector in the parallel transport. In order to obtain information on the vector length, Weyl introduced an extra connection, the length connection, which does not contain any knowledge about the direction of a vector on parallel transport. The only role of the length connection is to fix, or gauge, the conformal factor. The covariant divergence of the metric tensor is non-zero in Weyl's theory, and this property can be expressed mathematically in terms of a new geometric quantity, called non-metricity. In the physical applications of this geometry the length connection was identified with the electromagnetic potential. Dirac [79] proposed a generalization of Weyl's theory, which is based on the idea of the existence of two metrics, the physically undetectable metric ds E , altered by the transformations in the standards of length, while the second metric, a measurable one, is given by the conformally invariant atomic metric ds A . Weyl's theory has a remarkable intrinsic mathematical beauty, associated with a rich physical structure. However, it was largely ignored by physicists, and it did not become a mainstream research topics mainly because of Einstein's very early criticism [80] that "...in Weyl's theory the frequency of spectral lines would depend on the history of the atom, in complete contradiction to known experimental facts." However, another important development in geometry, which led to a new class of generalized geometric theories of gravity, took place due to the work of Cartan, who, based on his geometric work [81], proposed an extension of general relativity [82][83][84], known today as the Einstein-Cartan theory [85]. The torsion field, representing the new geometric element of the theory, is usually interpreted, from a physical point of view, as the spin density [85]. The Weyl geometry can be naturally extended to include the torsion. The resulting geometry is called the Weyl-Cartan geometry, and it was widely studied from both mathematical and physical points of view [86][87][88][89][90][91][92][93][94]. Torsion was included in the geometric framework of the Weyl-Dirac theory in [95][96][97], leading to an action integral from which one can construct a general relativistic massive electrodynamics, gauge covariant in the sense of Weyl. For a review of the of the physical applications and geometric properties of the Riemann-Cartan and Weyl-Cartan geometries see [98].
A third independent mathematical development that quickly did find important physical applications took place through the work of Weitzenböck [99], who introduced what are presently known as the Weitzenböck spaces. A Weitzenböck manifold is characterized by the properties ∇ µ g σλ = 0, T µ σλ = 0, and R µ νσλ = 0, where g σλ , T µ σλ and R µ νσλ are the metric tensor, the torsion tensor, and the curvature tensor of the manifold, respectively. When T µ σλ = 0, the Weitzenböck manifold is reduced to a Euclidean manifold. The torsion tensor has different values on different regions of the Weitzenböck manifold. Since the Riemann curvature tensor of a Weitzenböck space is zero, these geometries have the important property of distant parallelism, a property also known as absolute parallelism, or teleparallelism. Weitzenböck type space-times were first applied in physics by Einstein, who proposed a unified teleparallel theory of electromagnetism and gravity [100].
In the teleparallel approach to gravity the basic idea is to replace the metric g µν of the space-time, the basic physical variable describing the gravitational properties, by a set of tetrad vectors e i µ . The torsion, generated by the tetrad fields, can then be used to entirely describe gravitational effects, with the curvature replaced by the torsion. Thus we arrive to the so-called teleparallel equivalent of General Relativity (TEGR), which was initially introduced in [101][102][103], and is also known presently as the f (T ) gravity theory. Hence, in teleparallel, or f (T ) type theories, torsion exactly compensates curvature, with the important consequence that the space-time becomes flat. An important advantage of the f (T ) gravity theory is that the field equations are of second order, unlike in f (R) gravity, which in the metric approach is a fourth order theory. For a detailed discussion of teleparallel theories see [104]. f (T ) gravity theories have been widely applied to the study of astrophysical processes, and to cosmology, and in particular they are extensively used to explain the late-time accelerating expansion of the Universe, without the need of introducing dark energy [105][106][107][108][109][110][111][112][113][114][115][116][117][118][119][120][121][122][123].
In [124] an extension of the teleparallel gravity models, called WCW gravity, was proposed. In this theory, the Weitzenböock condition of the vanishing of the sum of the curvature and torsion scalar is imposed in a background Weyl-Cartan type space-time. A basic difference with the standard teleparallel theories is that this the model is formulated in a four-dimensional curved space-time, and not in a flat Euclidian geometry. WCW gravity leads to a purely geometrical description of dark energy, with the late time acceleration of the Universe fully determined by the intrinsic properties of the space-time. An extension of the Weyl-Cartan-Weitzenböck (WCW) and teleparallel gravity in which the Weitzenböck condition of the exact cancellation of curvature and torsion in a Weyl-Cartan geometry is inserted into the gravitational action via a Lagrange multiplier was considered in [125]. As a particular model the case of the Riemann-Cartan space-times with zero nonmetricity, which mimics the teleparallel theory, was considered. Several classes of exact cosmological models were also investigated.
From the above presentation it turns out that general relativity can be represented in (at least) two equivalent geometric representations: the curvature representation (in which the torsion and the nonmetricity vanish), and the teleparallel representation (in which the curvature and the nonmetricity vanish), respectively. However, a third equivalent representation is also possible, in which the basic geometric variable describing the properties of the gravitational interaction is represented by the nonmetricity Q of the metric, which geometrically describes the variation of the length of a vector in the parallel transport. Such an approach, called symmetric teleparallel gravity, was initially introduced in [126], and it has the advantages of covariantizing the usual coordinate calculations in general relativity. It turns out that in symmetric teleparallel gravity the associated energy-momentum density is essentially the Einstein pseudotensor, which becomes a true tensor in this geometric representation. Symmetric teleparellel gravity was further developed into the f (Q) gravity theory (or coincident general relativity) in [127], and it is also known as nonmetric gravity. Different geometrical and physical aspects of symmetric teleparallel gravity have been investigated in the past two decades in a number of studies, with the interest for this theory increasing rapidly recently [128][129][130][131][132][133][134][135][136][137][138][139][140][141][142][143][144][145][146][147]. For a review of teleparallel gravity see [148].
The propagation of gravitational waves in various extensions of symmetric teleparallel gravity was investigated in [136], with a particular focus on their speed and polarization. For the simple symmetric teleparallel gravity, and for theories that arise from the generalized irreducible decomposition of symmetric teleparallel gravity, as well as for f (Q) gravity, the same speed and polarizations of the gravitational waves were obtained as in general relativity. A derivation of the exact propagator for the most general infinite-derivative, even-parity and generally covariant theory in the symmetric teleparallel spacetimes was presented in [137]. In this approach the action made up of the non-metricity tensor and its contractions was decomposed into terms involving the metric and a gauge vector field. The propagation velocity of the gravitational waves around Minkowski spacetime and their potential polarizations in a general class of symmetric teleparallel gravity theories, called "newer general relativity" class, was investigated in [141]. The theory is defined in terms of the most general Lagrangian that is quadratic in the nonmetricity tensor, does not contain its derivatives and is determined by five free parameters. As a result of this investigation it was found that all gravitational waves propagate with the speed of light. The Noether Symmetry Approach was used to classify all possible quadratic, first-order derivative terms of the nonmetricity tensor in the framework of Symmetric Teleparallel Geometry in [142]. The considered models were invariant under point transformations in a cosmological background. The symmetries of these models were used to reduce the dynamics of the system in order to find analytical solutions. The cosmology of the f (Q) theory and its observational constraints were investigated in [143] and [144], and it was shown that in this theory the accelerating expansion is an intrinsic property of the geometry of the Universe, without need of either exotic dark energy or extra fields. The dynamical system method was used to investigate the general properties of the cosmological evolution. The behaviour of the cosmological perturbations in f (Q) gravity was investigated in [146]. Tensor perturbations feature a re-scaling of the corresponding Newton's constant, while vector perturbations do not contribute in the absence of vector sources. In the scalar sector two additional propagating modes were found, indicating that f (Q) theories introduce, at least, two additional degrees of freedom.
An extension of symmetric teleparallel gravity was considered in [140] by introducing, in the framework of the metricaffine formalism, a new class of theories where the nonmetricity Q is nonminimally coupled to the matter Lagrangian. A Lagrangian of the form L = f 1 (Q) + f 2 (Q)L m was considered, where f 1 and f 2 are generic functions of Q, and L m is the matter Lagrangian. This nonminimal coupling leads to the nonconservation of the energy-momentum tensor, and consequently the appearance of an extra force in the geodesic equation of motion. Several cosmological applications were considered for some specific functional forms of the functions f 1 (Q) and f 2 (Q), such as power-law and exponential dependencies of the nonminimal couplings. The cosmological solutions lead to accelerating evolutions at late times.
It is the main goal of our present investigation to consider another extension of f (Q) gravity, which is based on the nonminimal coupling between the nonmetricity Q and the trace T of the matter energy-momentum tensor. More exactly, we assume that the Lagrangian density of the gravitational field is given by a general function of both Q and T , so that L = f (Q, T ). From this gravitational Lagrangian the geometric action can be constructed in the usual way. By varying the action with respect to the metric tensor we obtain the general field equations describing gravitational phenomena in the presence of geometry-matter coupling. By considering the covariant derivative of the field equations we obtain the basic result that the divergence of the matter energy-momentum tensor does not vanish in the present approach to the gravitational interaction. The cosmological implications of the f (Q, T ) theory are investigated for three classes of specific models. The obtained solutions describe both accelerating and decelerating evolutionary phases of the Universe, and they indicate that f (Q, T ) gravity can provide useful insights for the description of the early and late phases of cosmological evolution.
The present paper is organized as follows. The geometric background, the gravitational action, the field equations and the divergence of the matter energy-momentum tensor are presented in Section II. The cosmological formalism of f (Q, T ) gravity is investigated, for a homogeneous and isotropic flat geometry in Section III. Three specific cosmological models, corresponding to different choices of the function f (Q, T ), are analyzed in detail in Section IV. We discuss and conclude our results in Section V. The explicit calculations of the geometric and physical quantities necessary to obtain the field equations and the divergence of the matter energy-momentum tensor (the general expression of Q, the variation δQ, the variation of the gravitational action with respect to the connection, the divergence of the field equations, and the expression of Q for the cosmological case) are presented in detail in Appendices A-E.

II. FIELD EQUATIONS OF f (Q, T ) THEORY
In the present Section we briefly review the geometrical foundations of the gravitational theories based on the assumption of the existence of a general line element in the spacetime. Then we will introduce the variational principle of the f (Q, T ) gravitational theory, and we obtain the gravitational field equations of this geometric approach to the gravitational phenomena. The divergence of the matter energy-momentum tensor is also considered, and we show that due to the coupling between matter and geometry this tensor is not conserved.

A. Geometrical preliminaries
Weyl introduced an important generalization of the Riemannian geometry, representing the mathematical basis of general relativity, by assuming that during the parallel transport around a closed path, an arbitrary vector will not only be subjected to a change of its direction, but it will also experience a modification of its length [78]. To describe mathematically these two simultaneous changes, Weyl proposed the introduction of a new vector field w µ , which, together with the metric tensor g µν , represent the fundamental fields of the Weyl geometry. The Weyl geometric theory has the important characteristic that the mathematical properties of the vector w µ exactly coincide with those of the electromagnetic potentials. This suggests that the electromagnetic and gravitational forces, both long-range forces, may have a common geometric origin [79].
If in a Weyl space a vector of length l is carried along an infinitesimal path δx µ by parallel transport, the variation in its length δl is given by the expression δl = lw µ δx µ [79]. After the parallel transport of a vector around a small closed loop of area δs µν , the variation of the length of the vector is given by the expression δl = lW µν δs µν , where we have denoted and where the covariant derivative ∇ ν is defined with respect to the metric g µν . By performing a local scaling of lengths of the forml = σ(x)l, the field w µ changes asw µ = w µ + (ln σ) ,µ , while the metric tensor coefficients are modified according to the conformal transformationsg µν = σ 2 g µν andg µν = σ −2 g µν , respectively [98]. Another important property of the Weyl geometry is the existence of the semi-metric connection where Γ λ µν denotes the usual Christoffel symbol, obtained with the help of the metric g µν . In the Weyl geometryΓ λ µν is assumed to be symmetric in its lower indices, and with its help one can construct a gauge covariant derivative in the standard way [98]. By using the covariant derivative one can obtain the Weyl curvature tensor, which can be written as where we have defined the quantities respectively, with the square brackets denoting antisymmetrization. For the first contraction of the Weyl curvature tensor we find where by R µ ν we have denoted the Ricci tensor constructed from the metric. Finally, for the Weyl scalar we obtain the expressionR The Weyl geometry can be generalized by taking into account the torsion of the space-time, thus obtaining the Weyl-Cartan spaces with torsion. In a Weyl-Cartan space-time we can introduce a symmetric metric tensor g µν , which defines the length of a vector, and an asymmetric connection Γ λ µν , which determines the law of the parallel transport as dv µ = −v σΓµ σν dx ν [85,98]. In the case of the Weyl-Cartan geometry the connection can be decomposed into three irreducible parts as follows: the Christoffel symbol Γ λ µν , the contortion tensor C λ µν , and the disformation tensor L λ µν , respectively, so that generally one can write [85] The first term in the above equation, the Levi-Civita connection of the metric g µν , is given by its standard definition The contorsion tensor C λ µν in Eq. (8) can be obtained from the torsion tensorΓ λ [µν] , defined aŝ according to the following relation As one can see immediately from the above equation, the contorsion tensor is antisymmetric with respect to its first two indices. The disformation tensor is obtained from the nonmetricity as As for the non-metricity tensor Q λµν , it is defined as (minus) the covariant derivative of the metric tensor with respect to the Weyl-Cartan connectionΓ λ µν , ∇ σ g µν = Q σµν , and it can be obtained as [85], The comparison of Eqs. (2) and (8) immediately show that the Weyl geometry is a particular case of the Weyl-Cartan geometry, in which the torsion is zero, and the non-metricity is represented by the expression Q λµν = −2g µν w λ . Therefore in a Weyl-Cartan geometry the connection can be written in the form where is the contortion, while the Weyl-Cartan torsion T λ µν is defined according to With the use of the connection, one can define the Weyl-Cartan curvature tensor aŝ With the use of Eq. (14), one can find the curvature tensor R λ µνσ in the terms of the standard Riemann tensor, plus some new terms containing the Weyl vector, the torsion and the contortion. By contracting the resulting curvature tensor, one can obtain the Weyl-Cartan scalar of the geometry as followŝ where we have defined T µ = T ν µν , and all covariant derivatives are considered with respect to the metric.
The symmetric teleparallel gravity is a geometric description of gravity, which is fully eq1uivalent to general relativity. This equivalence can be easily proven in the so-called coincident gauge, for whichΓ λ µν ≡ 0. Now, by imposing the condition that the connection is symmetric, the torsion tensor identically vanishes, and the Levi-Civita connection can be expressed in terms of the disformation tensor as On the other hand, as it is well known ro standard general relativity, after eliminating the boundary terms in the expression of the Ricci scalar, the gravitational action can be reformulated in a (noncovariant) form as [6] By taking into account the relation (19), it turns out that in the coincident gauge the gravitational action can be reformulated in terms of the disformation tensor as The action given by Eq. (21), called the action of the symmetric teleparallel gravity, is thus equivalent with the standard Hilbert-Einstein action of general relativity. However, there are some fundamental differences between the two gravitational models. In the symmetric teleparallel gravity the overall geometry of the space-time is flat, due to the vanishing of he curvature tensor (17). Hence the global geometry is of Weitzenböck type. Moreover, the gravitational effects are carried out not because of the rotation of the angle between two vectors in the parallel transport, but because of the variation of the length of the vector itself. In the following we will consider an extension of the Lagrangian (21) of the symmetric teleparallel gravity, given by where g ≡ det (g µν ), and we have defined and respectively. By T we have denoted the trace of the energymomentum tensor. We define the trace of nonmetricity tensor as We also introduce the superpotential of our model, defined as Then, as explicitly shown in Appendix A, we obtain for Q the relation Next, we vary the action in Eq. (22) with respect to the components of the metric tensor. Hence, as a first step, we obtain, The explicit form of the variation of δQ is presented in Appendix B. Moreover, as usual, we define which means that δT = δ(T µν g µν ) = (T µν + Θ µν ) δg µν . Then we can easily find for the variation of the action the expression As for the term 2f Q √ −gP αµν ∇ α δg µν , after integration and with the use of the boundary conditions it turns out that it takes the form −2∇ α (f Q √ −gP αµν ) δg µν . Finally, after equating the variation of the gravitational action to zero, we obtain the field equations of the f (Q, T ) gravity theory as, Ref. [140] also has similar terms like Eq. (31), even that the considered basic physical model and action are somehow different from the present approach.
By using the Lagrangian Multiplier Method with two constrains T α βγ = 0 and R α βµν = 0, we can find the variation with respect to the connection. The explicit calculations are presented in Appendix C. Moreover, we define the hypermomentum tensor density as By taking into account the anti-symmetry property of µ and ν in the Lagrangian multiplier coefficients λ µν α and ξ βµν α , we can eliminate them by introducing ∇ µ ∇ ν into the original part of action variation. Hence, after taking the variation of the gravitational action with respect to the connection we obtain the field equations C. The energy-momentum tensor balance equation For a (1,1)-form tensor v µ ν we define its covariant derivative as Here we haveΓ α µν = Γ α µν + L α µν , while by Γ α µν we have denoted the Levi-Civita connection associated to the metric. D µ denotes the covariant derivative with respect to the Levi-Civita connection. From Eq. (24) one can easily check that L µ µρ = −1/2 Q ρ . The field equations in the (1,1)-form are given by The metric divergence of the field equations (35) is explicitly calculated in Appendix D, and it is given by To simplify the above equation, we first solve Eq. (33) by introducing the tensor A ν α , so that where we have the additional constrain that We can always add an anti-symmetrical tensor ∇ µ M [µν] α to √ −gA ν α without adding extra terms to Eq. (38). Now we simply combine Eq. (36) and Eq. (37), and we find another form of the energy-momentum balance equation, given by or, equivalently, Hence in the f (Q, T ) gravity theory the matter energymomentum tensor is not conserved, D µ T µ ν = B ν = 0, with the nonconservation vector a function of Q, T , and of the thermodynamics quantities of the system. For a perfect fluid, described by its energy density ρ and its pressure p, respectively, the energy-momentum tensor takes the form T µ ν = (ρ + p) u ν u µ +pδ µ ν , u µ is the four-velocity of the fluid, normalized as u µ u µ = −1. Then, as shown in [140], from the divergence of the energy-momentum tensor we obtain the energy balance and the momentum conservation equations aṡ and respectively, where we have denoted by an overdot the quantityḟ = u µ D µ f , while we have defined H = (1/3) D µ u µ . h µν is the projection operator, given by h µν = g µν + u µ u ν . Eq. (41) describes the energy balance in f (Q, T ) gravity. From a physical point of view it gives the amount of energy that enters or goes out in a specified volume of a physical system. The source term B µ u µ corresponds to the energy creation/annihilation. The total energy of the gravitating system is conserved only if the condition B µ u µ = 0 is satisfied in all points of the spacetime. If B µ u µ = 0, then energy transfer processes or particle production takes place in the given system. Eq. (42) represents the equation of motion of massive particles in f (Q, T ) gravity. As it can be seen immediately from the equation of motion, the dynamical evolution of the massive particles is not geodesic, and an extra-force with components F µ = h µν (B ν − D ν p) / (ρ + p) does appear, due to the coupling between Q and T . Hence in f (Q, T ) gravity a supplementary force is exerted on any particle, besides the usual gravitational force. F µ is orthogonal to the matter four-velocity u µ , since from the properties of the projection operator it follows that we always have F µ u µ = 0, which is the standard requirement for a physical force, for which only the components that are orthogonal to the four-velocity of the particle can contribute to its equation of motion.

III. COSMOLOGICAL EVOLUTION OF THE FRIEDMANN-LEMAITRE-ROBERTSON-WALKER UNIVERSE IN f (Q, T ) GRAVITY
We are going now to consider the cosmological applications of the f (Q, T ) theory, by assuming that the Universe is described by the isotropic, homogeneous and spatially flat Friedmann-Lemaitre-Robertson-Walker (FLRW) metric, given by where a(t) is the scale factor, and the lapse function N (t) = 1 is for the standard case. The expansion and dilation rates are defined as follows In cosmology H(t) is called the Hubble function. By adopting the coincident gauge, in the covariant derivatives reduce to ordinary derivatives, after straightforward calculations presented in the Appendix E, we find

A. The generalized Friedmann equations
To derive the two generalized Friedmann equations describing the cosmological evolution, we assume that the matter content of the Universe consists of perfect fluid, whose energy-momentum tensor is given by T µ ν = diag(−ρ, p, p, p). Then for the tensor Θ µ ν we obtain the expression To simplify the mathematical formalism we introduce the notations and respectively. By using the FLRW metric, from the field equations we can easily find By solving Eq. (49) and Eq. (50) we obtain where we have denoted By explicitly includingρ andṗ in the expression ofḟ = FQ + 8πGṪ , we obtain the generalized energy balance equation in f (Q, T ) gravity aṡ We can easily see from the above equation that when f has no T dependence, which means G = 0, the continuity equation is always valid.
Next, we consider the case when N = 1, which is the case of the standard FRW geometry. Thus we have Q = 6H 2 , M = f /2−6F H 2 , S = 2 Ḟ H + FḢ , and the generalized Friedmann equations reduce to Combining the above two equations, we obtain the evolution equation for the Hubble function H aṡ We can bring the cosmological evolution equations to a form similar to the standard general relativity Friedmann's equations, by defining an effective energy density ρ ef f and an effective pressure p ef f so that, Then it follows that the effective thermodynamic quantities satisfy the conservation equatioṅ An important cosmological quantity is the deceleration parameter q, which is an indicator of the accelerating/decelerating nature of the evolution of the Universe. The deceleration parameter is defined as where w = p ef f /ρ ef f is the parameter of the equation of state of the dark energy. Negative values of the deceleration parameter indicate an accelerating evolution, while positive values indicate decelerating expansion. Explicitly, the deceleration parameter can be expressed as, To obtain cosmological results that can allow a direct comparison of the model predictions with the astronomical observations, we introduce, instead of the time variable t, as independent variable the redshift z, defined according to where we have normalized the scale factor so that its present day value is one, a(0) = 1. Therefore for the time operator we obtain The deceleration parameter q can be obtained as a function of the cosmological redshift z as We will also compare the behavior of the cosmological parameters in the f (Q, T ) gravity with the standard ΛCDM model. We assume that the late Universe is filled with dust matter only, having negligible pressure. Then form the standard general relativistic energy conservation equationρ + 3Hρ = 0 we find for the variation of matter energy density the expression ρ ∼ 1/a 3 ∼ (1 + z) 3 . The evolution of the Hubble function is given by [17] where Ω DM , Ω b and Ω Λ are the density parameters of the cold dark matter, baryonic matter, and dark energy (interpreted as a cosmological constant), respectively. The density parameters satisfy the important constraint Ω DM + Ω b + Ω Λ = 1. The Hubble function H(z) = H 0 h(z) can be written as a function of the redshift in a dimensionless form as The redshift dependence of the deceleration parameter is obtained as In the following for the density parameters we adopt the numerical values Ω DM = 0.2589, Ω b = 0.0486, and Ω Λ = 0.6911 [17], obtained from the Planck data, giving for the total matter density parameter Ω m = Ω DM + Ω b the numerical value Ω m = 0.3089. From these numerical values of the cosmological parameters it follows that the present day value of the deceleration parameter as q(0) = −0.5381. As for the variation of the dimensionless matter density with respect to the redshift, we obtain the expression r(z) = Ω m (1 + z) 3 = 0.3089(1 + z) 3 .

B. The de Sitter solution
Before considering specific cosmological models of f (Q, T ) gravity, we would to find the vacuum solution for our field equations, and check if the theory admits a de Sitter type solution, which corresponds to the constrains ρ = p = 0 and H = H 0 =constant, respectively. For a vacuum Universe Eq. (51) suggests M = S = 0, where S = 0 gives a constant F = F 0 , which implies f = F 0 Q + Λ, with Λ also an arbitrary constant of integration. The condition M = 0 reduces to which simply gives H 0 = Λ/6F 0 . This result is similar to the one in Ref. [140], which also gives a result equivalent to the general relativistic case when F 0 = 1. Hence the f (Q, T ) theory admits the de Sitter type evolution in the limiting case of a vacuum Universe. As one can easily calculate, for the de Sitter solution we have q = −1 and w = −1, respectively.

IV. SPECIFIC COSMOLOGICAL MODELS
In the present Section we will investigate some specific cosmological models in the f (Q, T ) gravity theory, corresponding to different choices of the functional form of f (Q, T ). For the sake of generality we will assume that the cosmological matter satisfies an equation of state of the form p = (γ − 1)ρ, where γ is a constant, and 1 ≤ γ ≤ 2. Such a linear barotropic equation of state can describe the baryonic matter behavior in both the high density limit (corresponding to the early Universe), and in the low density limit, appropriate for the description of the present day Universe.
With the use of the barotropic equation of state, from Eqs. (54) and (56) we obtain for the matter density the general expression As a first example of the cosmological evolution in f (Q, T ) gravity we will consider the case in which the function f (Q, T ) has the simple form f (Q, T ) = αQ + βT , where α and β are constants. Then we immediately obtain F = F Q = α, and 8πG = f T = β. Hence Eq. (56) becomeṡ Eqs. (57) and (58) take the form and respectively, which leads to the consistency condition For ρ = 0, the above condition implies 1 + β/8π = 0, which in turn leads, with the use of Eq. (70), to the equationḢ = −, or H = H 0 = constant, and a(t) = e H0t . The cosmological evolution is of de Sitter type, in the presence of a nonvanishing matter energy density. The evolution of ρ can be obtained from the conservation equation (53), which taking into account that in the present model S = 0, becomeṡ giving ρ(t) = ρ 0 e −3γH0t , where ρ 0 is an arbitrary constant of integration. Hence, the exponential expansion of the Universe is associated, in this model of the f (Q, T ) gravity theory, with an exponential decrease of the matter content.
Then from Eq. (69) we obtain the expression of the matter density as By using this expression of the density it follows that the evolution equation for H, Eq. (56), takes the simple forṁ and it has the general solution where we have used the initial condition H (t 0 ) = H 0 . The evolution of the scale factor is given by where a 0 is an arbitrary constant of integration, and we have denoted a 1 = (n + 1) [4(β + 4π) − βγ]. The deceleration parameter is constant, and it is given by If the model parameters satisfy the constraint 3(β+8π)γ/(n+ 1) [16π − β(γ − 4)] < 1, the deceleration parameter takes negative values, and the expansion of the Universe is accelerating.
Finally, as a simple example of a cosmological model in f (Q, T ) gravity, we will consider the case when the function f (Q, T ) has the form f (Q, T ) = −αQ − βT 2 , where α > 0 and β > 0 are constants. Moreover, for the sake of simplicity, we will fix the equation of state of the cosmological matter from the beginning as dust, that is, we choose γ = 1, giving p = 0. Then we immediately obtain Eq. (69) gives for the matter density the simple algebraic equation which has the physical solution If the condition 3(1 + 32π)αβH 2 (t)/32π 2 << 1, is satisfied, by power expanding the square root in the above equation gives ρ(t) ∝ H 2 (t). Thus in this limit we recover the standard general relativistic result. The evolution equation for the Hubble function, Eq (69), takes for this model the forṁ We rescale now the Hubble function according to and we introduce the model parameter λ, defined as Then, in the new variables, Eq. (84) becomes (87) In terms of the redshift z Eq. (87) takes the form By introducing the new variable u(z) = h 2 (z), the above equation can be written as (89) In the limit h 2 (t) << 1, Eq. (87) can be approximated as with the general solution given by where h 0 = h (t 0 ). From H(t) = H 0 h(t) =ȧ/a we obtain the scale factor as a(t) = [2 + h 0 λ (t − t 0 )] 2H0/λ . The deceleration parameter is given by q = λ/2H 0 − 1 = √ 3/2 − 1 < 0, while the matter energy density varies as The variation of the Hubble function as a function of the redshift, obtained by numerically integrating Eq. (89), is presented in Fig. 1. The evolution equation for the Hubble function was integrated with the initial condition u(0) = 1, and we have considered the redshift range z ∈ [0, 1].
As one can see from Fig. 1, the Hubble function is a monotonically increasing function of the redshift (monotonically decreasing function of time) for all considered values of the model parameter λ. For small values of z the cosmological evolution is practically independent on the numerical values of λ, but at higher redshifts there is a significant effect of the parameter value on the cosmic expansion. For the sake of comparison we have also presented the variation of the Hubble function in the standard ΛCDM model, as given by Eq. (66). Despite the existence of some quantitative differences in the Hubble functions of the two models, at least on a qualitative level the two descriptions give relatively similar results. The For all adopted numerical values of the parameter λ the energy density is a monotonically increasing function of the redshift (a monotonically decreasing function of the cosmological time). The increase is almost linear, and for small redshifts it is almost independent on the numerical values of λ. However, a dependence on the model parameter can be seen at higher redshifts. The comparison with the matter energy density in the ΛCDM model shows that, if in the range z ∈ [0, 0.4] there is an approximate concordance between the predictions of the two models, for higher redshifts the differences in the matter densities are high. While in the ΛCDM model the matter energy density increases rapidly as (1 + z) 3 , the almost linear increase of ρ in this particular f (Q, T ) model predicts a much lower matter density at higher redshifts.
The variation of the deceleration parameter q is represented, as a function of the redshift, in Fig. 3, for the same values of the parameter λ as considered in the previous figures.
The deceleration parameter is a monotonically increasing function of z. The evolution of the Universe begins, at redshift z = 1, from a decelerating phase, with q > 0. The expansion of the Universe accelerates, and at a finite value of z it reaches the value q = 0, corresponding to the transition to the accelerated phase. The evolution of q is strongly dependent on the numerical values of the model parameter λ. Depending on these values a large range of present day values of the deceleration parameter can be obtained. The comparison with the deceleration parameter variation in the ΛCDM model show that there is a qualitative similarity between the two models. However, the present choice of the function f (Q, T ) cannot fully reproduce the standard cosmological evolution, in a quantitative way. However, it provides similar qualitative results. Hopefully, by using an advanced fitting procedure, based on the direct application of observational data, and more general functional forms of f (Q, T ), this gravity theory may provide an alternative to the standard ΛCDM model.

V. DISCUSSIONS AND FINAL REMARKS
After more than one hundred years since Einstein did propose the first geometric theory of gravity, general theory, we are presently witnessing the very interesting situation that at least three geometric descriptions of gravity are possible, based on the three basic quantities introduced in Riemannian geometry and its extensions (curvature, torsion and nonmetricity, respectively). These findings raise the fundamental question of the possibility of a unique geometric description of gravity. Are these three descriptions completely equivalent, or they are perhaps just some particular cases of a more general geometric theory, which is still needed to be found?
In the present paper we have investigated some theoretical aspects of the third geometric description of gravity, the symmetric teleparallel gravity, or f (Q) gravity, by introducing a new class of theories where the nonmetricity Q is coupled nonminimally to the trace of the matter energy-momentum tensor. From a mathematical point of view we have performed our analysis in the framework of the metric-affine formalism. Our theory is constructed in a similar way like the f (R, T ) theory [44], but with the geometric part of the action being replaced by the symmetric teleparallel formulation. Similarly to the in the standard curvature -trace of energy-momentum tensor couplings, in the f (Q, T ) theory the coupling between Q and T leads to the nonconservation of the energy-momentum tensor. This nonconservation has important physical implications, implying significant changes in the thermodynamics of the Universe, similarly to those in the f (R, T ) theory [44], and, due to the nongeodesic motion of test particles, to the appearance of an extra force. Our approach may also lead to an improvement of the geometrical formulation of gravity theories with geometry-matter coupling. Implemented in both matter and geometry sectors, our approach allows a consistent and workable representation of the nonminimal curvature-matter coupling theories. In this context we have derived the gravitational field equations of the f (Q, T ) gravity theory from a variational principle that generalizes the variational principle of the f (Q, T ) theory, and we have obtained the general relation describing the nonconservation of the matter energy-momentum tensor.
As a theoretical test of our theory we have analysed its cosmological applications. As a first result in this respect we have obtained the generalized Friedmann equations of the f (Q, T ) theory describing the cosmological evolution in a flat, homogeneous and isotropic Friedmann-Lemaitre-Robertson-Walker type geometry. The generalized Friedmann equations can be reformulated as the standard equations of general relativity, but with the ordinary matter energy density and pressure replaced by some effective quantities. The effective quantities depend on the Lagrangian f of the theory, and on its derivatives with respect to Q and T . Interestingly enough, both effective thermodynamic quantities contain linear combinations of the ordinary matter energy density and pressure. In fact the coupling between the trace of the energy-momentum tensor and the Q field introduces two types of corrections. The first is the presence of an additive term of the form f /4F that independently appears in both Friedmann equations. Secondly, we have the term 4π/F , multiplying the linear combination constructed from the components of the energy-momentum tensor (energy density and pressure) in both Friedmann equations. The coefficients of the linear combinations of the energy density and pressure are constructed from the functionG ∝ f T . Consequently, the basic equations describing the cosmological evolution in f (Q, T ) gravity can be formulated in terms of an effective energy density and pressure, which both depend on the energy and pressure components of the energy-momentum tensor, and on the functions f (Q, T ), f Q (Q, T ) and f T (Q, T ), respectively. An important indicator of the cosmological evolution, the deceleration parameter, given by Eq. (61), has a complicated dependence on the Lagrangian f and its derivatives, indicating that, depending on the functional form of f (Q, T ), a large number of cosmological models can be obtained. We have also shown explicitly that for the vacuum case, when the f (Q, T ) theory reduces to f (Q) theory, for late times, the Universe enters into an exponentially accelerating de Sitter type phase.
We have also considered three explicit classes of cosmological models, obtained by imposing some specific simple mathematical forms for the function f (Q, T ). In all three example we have considered that Q and T enter in an additive form in the Lagrangian, neglecting the possible existence of some cross terms of the type QT , or functions of it. In the case f (Q, T ) = αQ + βT , the cosmological evolution is of de Sitter type, with the Universe expanding exponentially. The model f (Q, T ) = αQ n + βT leads to a power law type form of the scale factor, and to a constant deceleration parameter. However, by an appropriate choice of the model parameters α, β and γ accelerating expansions can be obtained easily. The third model with f (Q, T ) = −αQ − betaT 2 leads to a complex cosmological dynamics, involving the transition from a decelerating to an accelerating state. The results can be obtained only by numerically integrating the generalized Friedmann equations. The nature of the cosmological evolution is strongly dependent on the numerical values of the model parameters. We have also compared the theoretical predictions of the f (Q, T ) theory with the corresponding results in the standard ΛCDM cosmology. For the specific range of cosmological parameters we have considered it did follow that the Universe began its recent evolution in a decelerating phase, and in the large time limit it can reach a de Sitter phase. Depending on the model parameters, a large spectrum of present day values of the deceleration parameter can be obtained. The theoretical predictions of the Hubble parameter are similar to those of the standard general relativistic cosmological model in the presence of the cosmological constant. However, some significant deviations appear for the behavior of the matter energy density and of the deceleration parameter. But if investigated for a larger range of parameters and functional forms of f this model may represent an alternative to the ΛCDM cosmology, with the late time de Sitter phase induced by the coupling between nonmetricity and matter.
The f (Q, T ) gravity theory is also valid when instead of ordinary matter one includes scalar fields in the action. Another possible application of the f (Q, T ) theory is to consider inflation in the presence of scalar fields, an approach that may provide a completely new perspective on the geometrical, gravitational, and cosmological processes that did play a major role in the very early dynamics of the Universe. The analysis of structure formation in f (Q, T ) theory is also a major topics of research that could be investigated, with the use of a background metric. For different nonmetricity-trace of the energy-momentum tensor coupling models, the SNIa, BAO, CMB shift parameter data can be used to obtain constraints for the respective models. Moreover, such an approach may allow the detailed exploration and analysis of structure formation from a different perspective. An interesting issue is to obtain the Newtonian and the post-Newtonian limits of the f (Q, T ) gravity, and to investigate the constraints the local gravity at the Solar System level impose on the theory. The Newtonian limit can also help in finding constraints arising from other astrophysical observations.
To conclude, in the present investigation we have introduced a new version of the symmetric teleparallel theory, and we have proven its theoretical consistency. This approach also motivates and encourages the study of further extensions of the f (Q) type family of theories. We have shown that the presented approach predicts de Sitter type expansions of the Universe, and thus it may represent a geometric alternatives to dark energy. Hence this study offers some basic theoretical tools for the in depth investigation of the geometric aspects of gravity, and of its cosmological implications.