Reconstruction of energy conditions from observations and implications for extended theories of gravity

The attempt to describe the recent accelerated expansion of the Universe includes different propositions for dark energy models and modified gravity theories. Establish their features in order to discriminate and even rule out part of these models using observational data is a fundamental issue of cosmology. In the present work we consider a class of extended theories of gravity (ETG) that are minimally coupled to the ordinary matter fields. In this context, and assuming a homogeneous and isotropic spacetime, we derive the energy conditions for this ETG class. We then put constraints on these conditions using a model-independent approach to reconstruct the deceleration function along with the Joint Light-curve Analysis (JLA) supernova sample, 11 baryon acoustic oscillation and 22 cosmic-chronometer measurements. We also consider an additional condition imposing the strong energy condition only on the ordinary matter. This is to guarantee the presence of only attractive matter in the energy-momentum tensor, at least in the redshift range of the observations, i.e., the recent accelerated expansion of the Universe is due solely to the modifications in the gravity theory. The main result of this work is a general reconstruction of the energy conditions valid for every considered ETG.


I. INTRODUCTION
In the last decades, a great amount of cosmological observational data have been accumulated, endowing the modern cosmology with the capability of quantitatively reproducing the details of many observed cosmic phenomena, including the late time accelerating stage of the Universe. The evidence of such an acceleration comes, for example, from the measurements of the distance modulus of type Ia supernovae (SNe Ia), however it still lacks a satisfactory physical explanation. Most of the current studies relates it to some unknown field (dark energy) or to some extension of Einstein theory of General Relativity (GR).
In this context the energy conditions (ECs) are useful tools to evaluate some features of the Universe evolution, since they can be derived considering a minimum set of assumptions on the gravitation theory. The ECs have been used to obtain more information from data assuming GR and homogeneous and isotropic metric [1][2][3][4][5][6]. Among the interpretations of the ECs in GR, we have the positivity of the energy density (weak EC) and the attractiveness of gravity -focusing theorem -(strong EC) [7]. In Refs. [5,6] the authors found a violation of the strong EC, i.e., an indication of the accelerated expansion of the universe, with more than 99.73% confidence interval in the redshift range [∼ 0.1, ∼ 0.2] using SNe Ia data.
The ECs have also been addressed in the context of the Extended Theories of Gravity (ETGs) [8,9], such as f (R) gravity, which modifies the Einstein-Hilbert Lagrangian by the introduction of an arbitrary function of the curvature scalar R. As shown by Santos et al. [10], the ECs requirements may constrain the parameter space of each specific f (R) functional form. On the other hand, Bertolami and Sequeira [11] and Wang et al. [12] have generalized the ECs for f (R) theories with non-minimal coupling to matter. They found that the ECs are strongly dependent on the geometry. The conditions to keep the effective gravitational coupling positive as well as gravity attractive were also obtained. In a similar fashion, Wang and Liao [13] have also obtained the ECs requirements for such a theory. Furthermore, imposing the fulfillment of the ECs at the present time they have constrained the parameter space for a particular Lagrangian using the current measurements of the Hubble parameter H 0 , the deceleration parameter q 0 and the jerk parameter j 0 . Recently, Santos et al. [14] studied the attractive/non-attractive character of f (R) gravity considering the strong energy condition.
As a consequence of the ECs derived for a generic f (R) theory and of the equivalence between f (R) gravity and scalar-tensor theories, Atazadeh et al. [15] have studied the ECs in Brans-Dicke theory and put constraints on the parameters involved.
Using the parameters H 0 , q 0 and j 0 together with the appropriate ECs, García et al. [16], Banijamali et al. [17] and Atazadeh and Darabi [18] have shown the viability of some formulations of the modified Gauss-Bonnet gravity. The ECs inequalities were also used to constraint f (R, T ) theories [19], where T is the trace of the energymomentum tensor. An extension of such a theory which also takes into account an arbitrary function of the quantity Q = R ab T ab in the Lagrangian was considered in the scope of ECs by Sharif and Zubair [20] and bounds on the parameters were obtained. Modified theories of gravity with non-null torsion also have been constrained with ECs. In this respect, Sharif et al. [21] and Azizi and Gorjizadeh [22] were able to put bounds on the parameters of some particular formulations of the theory.
It is worth mentioning that differently from GR, the interpretation of the ECs in modified gravity theories is not necessarily direct as pointed out by Capozziello et al. [9]. For instance, for theories with non-minimal coupling between matter and geometry, the energy-momentum tensor is not conserved and test particles follow non-geodesic paths [8,23].
In the majority of the work mentioned above, the constraints on the parameter space of each specific theory were obtained by imposing that the ECs are satisfied at the present time. However, such a statement does not imply the fulfillment (or not) of the ECs in the whole cosmic history. Consequently, the parameters could be further constrained if the ECs were extended for redshifts beyond z = 0. The aim of the present work is to give a general treatment for such an issue in the context of ETGs.
In order to achieve this goal we consider a class of metric torsion-free ETGs that presents a minimal curvaturematter coupling and the energy-momentum tensor is conserved. We also assume a homogeneous and isotropic spacetime and conformal transformation such that the extra degrees of freedom of the theory can be written as an effective energy-momentum tensor. This comprises the f (R) theories, f (R, T ) theories for which the particular case ∇ b T ab = 0 is assumed, scalar-tensor theories such as Brans-Dicke theory and several other possible formulations [8]. We write the ECs in terms of the modified gravity functions. In addition, we derive another condition by assuming that gravity is always attractive when considering only ordinary matter on scales where the ETGs reduces to GR. We then apply the modelindependent reconstruction method of the deceleration function introduced by Vitenti and Penna-Lima [24] to obtain observational constraints on the functions of the modified term, given the EC and the additional inequations, for a redshift range and not only for the present time. In particular, we use the Joint Light-curve Analysis (JLA) supernova sample [25], 11 baryon acoustic oscillation (BAO) data points from 6dF Galaxy Survey and the Sloan Digital Sky Survey (SDSS) [26][27][28][29][30] and 22 cosmicchronometer [H(z)] measurements [31][32][33][34].
The layout of the paper is as follows. In Sec. II we introduce the ECs in the context of a class of ETGs, and then we derive the EC inequations considering a homogeneous and isotropic spacetime. In Sec. III we recall the main steps of the model-independent reconstruction approach [24], and obtain the deceleration and Hubble function estimates using SNe Ia, BAO and H(z) data. The observational bounds on the ECs and the additional condition and their respective implications in the context of GR and on the ETG functions are discussed in Sec. IV. Finally, we present our conclusions in Sec. V. Throughout the article we adopt the metric signature (−, +, +, +).

II. ENERGY CONDITIONS
In this section, we examine the ECs in the context of a class of ETGs for which the field equations can be written in the following generic form [9] where G ab ≡ R ab − 1 2 g ab R, R ab is the Ricci tensor, R = R a a is the Ricci scalar, G is the gravitational constant, and T ab is the energy-momentum tensor of the matter fields. The tensor H ab encapsulates the additional geometrical information of the modified theory. For instance, it may depends on scalar and (or) vector fields, scalars made out of Riemann and Ricci tensors, and derivatives of these quantities (see [8,35] and references therein). Finally, Ψ i refers to these mentioned fields and geometric quantities, and the modified coupling with the matter fields is driven by g 1 (Ψ i ) and g 2 (Ψ i ), where the latter includes explicit curvature-matter couplings [9,36].
In this work we consider a class of ETGs that presents a minimal curvature-matter coupling 1 , i.e., g 2 (Ψ i ) = 1, and the matter action is invariant under diffeomorphisms, then ∇ a T ab = 0 [7]. Given the conservation of the energymomentum tensor and the twice-contracted Bianchi identity, ∇ a G ab = 0, we have that [9] ∇ a H ab = − 8πG Many authors have been discussing the ECs in the context of different ETGs such as f (R) [10][11][12][13][14], scalar-tensor gravity theories [21] and massive gravity [37].
A common procedure, though, is to consider the modified gravity term H ab as a source of the effective energymomentum tensor, However, as also pointed out by Capozziello et al. [9], these fictitious fluids can be related to scalar fields or other further degrees of freedom. It is, thus, somewhat misleading to apply the standard ECs obtained in GR to the resulting effective energy-momentum tensors derived in such theories. In this sense, one should perform suitable conformal transformations in order to better define the energy conditions in terms of T eff ab [9,38,39]. Nevertheless, this misleading aspect rises from the fact that two of the energy conditions (strong and null) result from the convergence conditions applied in a GR setting, relating them to the energy-momentum tensor [40]. For this reason, here we use the convergence conditions directly in these two cases, avoiding the issue mentioned above.

A. Energy conditions for ETGs with minimal coupling
Consider a timelike geodesic congruence with a tangent vector field t a and let τ be a parameter of these timelike curves. As matter is minimally coupled to geometry, test particles follow geodesics, i.e., we are assuming the modified gravity does not introduce any additional force. In this case the Raychaudhuri equation is given by where θ = ∇ a t a is the expansion of the congruence. Analogously, for a congruence of null curves parametrized by λ and with tangent vector k a , we have In the above equations, σ ab (σ ab ) is the shear tensor, ω ab (ω ab ) is the vorticity tensor, and the hat means that these quantities are projected onto the subspace normal to the null vectors. Note that the necessary and sufficient conditions for a congruence be hypersurface-orthogonal is ω ab =ω ab = 0 [7]. Considering zero vorticity and given that the second term of Eq. (4) is nonpositive, the convergence of timelike geodesics (i.e., focusing) occurs if this inequation is also know as timelike convergence condition [40].
Equation (1) with g 2 = 1 can be written in the following form Thus, in GR (g 1 = 1 and H ab = H = 0), Eq. (6) implies the so called strong energy condition (SEC), which expresses the attractive character of gravity. Now, for the considered ETGs, Eq. (6) reads In this case, as pointed out by Capozziello et al. [9], even if the matter fields do not contribute positively (e.g., a matter field with negative pressure), Eq. (9) can still be fulfilled given the geometric term. In other words, the attractiveness of gravity can remain in the presence of a dark energy (DE) like fluid depending on the modified gravity theory. Similarly, from Eq. (5) we have that the condition for the convergence of null geodesics is R ab k a k b ≥ 0, the null convergence condition [40]. Thus, the null energy condition (NEC) in GR is and in the context of the ETGs considered here, On the other hand, the weak and dominant energy conditions, WEC and DEC, respectively, are direct restrictions on the energy-momentum tensor T ab (they do not originate from the Raychaudhuri equation). As we will see in more details in Sec. II B, the WEC states that the matter-energy density is positive for every time-like vector, i.e., and the DEC states that the speed of the energy flow of matter is less than the speed of light. This condition can be written in the form In this work we compare the bounds above with the reconstruction of the geometry. This means that we estimate R ab directly from the data and then apply the bounds to it. For this reason, the ECs stemming directly from T ab (WEC and DEC) will depend explicitly on the ETG functions, while SEC and NEC will not depend on them.
Inasmuch as in all known gravitational tests made in non-cosmological scales gravity is always attractive, we also consider an additional condition. The attractiveness of gravity in scales where it is well described by GR requires Eq. (8). Then our new condition will be the imposition of this inequality for T ab containing only ordinary matter. 2 Therefore, the way we defined SEC [Eq. (9)] is such that it measures if the ETG is attractive regardless of modification. As indicated by observational data, that it is not the case and SEC is violated, for this reason this bound serves to measure with what statistical confidence the violation occur.
On the other hand, the condition imposed by Eq. (8) has a more interesting interpretation. Consider for a moment the ETG equations of motion written as GR with an effective energy-momentum tensor [as in Eq. (3)]. Then imposing Eq. (8) is equivalent to requiring that all DE like behavior comes solely from H ab . Thus, this condition must be satisfied if one hopes to solve the acceleration problem using only an ETG, and not a mixture of ETG and a DE like matter content.

B. Homogeneous and isotropic spacetime
In this section we derive the ECs for a homogeneous and isotropic universe. This is described by the Friedmann-Lemaître-Robertson-Walker (FLRW) metric, where a(t) is the scale factor and k = 0, 1 or −1, whose flat, spherical and hyperbolic functions are S k (r) = (r , sin(r), sinh(r)), respectively. In this case, R ab and H ab are diagonal tensors and their components are functions of the time t only, namely, and, analogously, As stated above, to be compatible with a Friedmann metric, the tensor H ab must have the form of Eq. (16). This means that all information about the ETG is encoded in two time dependent functions h t (t) and h s (t).
In turn, the energy-momentum tensor for the matter fields can be written as where ρ is the matter-energy density, p is the pressure, and the four-velocity of the fluid is U a (where U a U a = −1). Hence, the Friedmann equations for the considered ETGs acquire the following form, Finally, giving a normalized timelike vector t a = N U a + N a , where t a t a = −1 and U a N a = 0, and a null vector k a , we rewrite the ECs by substituting Eqs. (16) and (17) into Eqs. (9), (11) (12) and (13), thus 8πG Note that, imposing that the conditions are fulfilled for any timelike/null vector field (parameterized as the timelike vector field t a defined above), we get two conditions for SEC, WEC and DEC when we apply it to a energymomentum tensor compatible with a Friedmann metric. In order to confront the ECs with observational data, such that one can infer local (in redshift) information about the fulfillment of these conditions, Lima et al. [5,6] showed that it is convenient to write the ECs in terms of the Hubble function, H(z) ≡ȧ/a = H 0 E(z), and the deceleration function, q(z) ≡ −äa/ȧ 2 , where 1 + z = a 0 /a. Therefore, by using the Friedmann equations, Eqs. (18) and (19), the energy conditions are rewritten as where Ω 0 k = −k/(a 0 H 0 ) 2 and the subscript (superscript) 0 stands for the present-day quantities. Conditions WEC1 and WEC2 refers respectively to the first and second inequalities in (22). Except for WEC, the above expressions refer only to the first inequation of each condition in Eqs. (21-23) (ρ − p ≥ 0 for DEC), and provide a complete unambiguous set of conditions.
The recent accelerated expansion of the universe evinced by SNe Ia [41,42], large scale structure (LSS) [43,44] and cosmic microwave background radiation (CMB) [45,46] data represents the violation of SEC [5,6], i.e., q(z) < 0. As a consequence, this fact requires a modified theory of gravity and/or the existence of an exotic fluid. On the other hand, a variety of tests in the solar system, for example, indicates that the GR is so far the most successful theory of gravitation [47].
Giving this, we now introduce a fifth condition imposing that ordinary matter, characterized by T ab , implies that GR is always attractive, as discussed in the end of Sec. II A. Thereby, from Eq. (18), we state that That is, any contribution for a negative value of the deceleration function, and consequently a DE like behavior, originates exclusively from the modified gravitational term H ab , thus One could also think that a similar new condition could be obtained from NEC. However, that is not the case, this condition coincides exactly with the second one obtained from WEC, i.e., WEC2. It is easy to observe this by comparing the second condition in Eq. (22) with Eq. (20).
In this work we aim to put observational constraints on combinations of h t (z) and h s (z) by requiring the fulfillment of the ECs and Eq. (30). For this, it is clear that Eqs. (24)(25)(26)(27)(28)(29)(30) require estimates of q(z) and E(z). Hence in the following section we will describe the methodology and observational data sets to obtain these estimates.
Here we apply the model-independent method to reconstruct the deceleration function introduced by Vitenti and Penna-Lima [24] (hereafter VPL), using SNe Ia, BAO and H(z) data (see Sec. III A). The approach is model-independent because q(z) is reconstructed without specifying the matter content of the universe or the theory of gravitation. Particularly, we just assume that the universe is homogeneous and isotropic.
In short, the VPL approach consists in approximating the q(z) function by a cubic spline over the redshift range [z min , z max ], where the minimum and maximum redshifts are defined by the observational data. Then, choosing the number of knots, n + 1, we write the reconstructed curvê q(z) in terms of the parameters {q i }, i = 0, ..., n. As discussed in Ref. [24], the complexity of the reconstructed function depends on n. However, instead of varying the number of knots, VPL addressed this question by including a penalty function, which is parametrized by σ rel , such that small values of σ rel (e.g., σ rel = 0.05) forcê q(z) to be a linear function, whereas large values (e.g., σ rel = 1.5) provide a high-complexity function. By construction, the errors of the reconstructed curves are dominated by biases (small σ rel ) and over-fitting (large σ rel ).
Considering different values of σ rel and fiducial models for q(z), VPL validated the reconstruction method via the Monte Carlo approach, using SNe Ia, BAO and H(z) mock catalogs. Then, evaluating the bias-variance trade-off for each case, which is characterized by a fiducial model and a σ rel value, the best reconstruction method was determined. That is, the best σ rel value is the one that minimizes the mean squared error, requiring the bias to be at most 10% of this error.
Finally, here we use the VPL approach considering 12 knots, σ rel = 0.3 and the redshift interval [0.0, 2.33] to reconstruct the deceleration function (for further discussions, see VPL [24]). The observable quantities such as E(z) and the transverse comoving distance D M (z) are written, respectively, in terms of q(z) as and for Ω k = 0, K −1 sinh (Kχ(z)) for Ω k > 0, (32) where K = H0 √ |Ω 0 k | c and the comoving distance is A. Data In the present work, we use some current available observational data for small redshifts (z ≤ 2.33) and their likelihoods, namely, the Sloan Digital Sky Survey-II and Supernova Legacy Survey 3 years (SDSS-II/SNLS3) combined with Joint Light-curve Analysis (JLA) SNe Ia sample [25], BAO data [26][27][28][29][30] and H(z) measurements [31][32][33][34]: where d and θ comprehend the observational data sets and the parameters to be fitted, respectively, and the last term is the penalization factor withq i = (q i−1 +q i+1 )/2.
The JLA SNe Ia likelihood is given by where m B is a vector of the 740 measured rest-frame peak B-band magnitudes, C −1 SNIa and |C SNIa | are the inverse and the determinant of the covariance matrix, respectively, and The SNe Ia astrophysical parameters α, β, M 1 and M 2 are related to the stretch-luminosity, colour-luminosity and the absolute magnitudes, respectively. The luminosity distance is D L (z hel , z cmb ) = (1 + z hel )D M (z cmb ), where z hel and z cmb are the heliocentric and CMB frame redshits. The BAO likelihood is where b and C BAO represent, respectively, 8 BAO measurements and the respective covariance matrix [26,28,30] as displayed in Table I and r d , i.e., the sound horizon at the drag redshift, where c s (z) is the sound wave speed in the photon-baryon fluid. Since the VPL reconstruction method is defined in a small redshift interval, we cannot calculate the integral above. Furthermore, this integral is model dependent, and for this reason in this analysis we treat r d as a free parameter. The last two terms of Eq. (37) are computed using the likelihood distribution from Refs. [27,29], respectively. We use 22 H(z) measurements obtained from Refs. [31][32][33][34] through the cosmic chronometers formalism. The H(z) likelihood is where the data points H obs i , errors σ i and the respective references are listed in Table II.

B. Analysis
We now apply the VPL methodology to reconstruct q(z) (and, consequently, H(z) and the cosmological distances) along with the observational data, described in Secs. III and III A, respectively, to put constraints on h t (z) and h s (z) from the EC bounds (see Eqs. (26), (27), (28) and (30)). We also evaluate the violation/fulfillment of the ECs in GR. The numerical as well as the postprocessing analyses carried out in this work made use of the Numerical Cosmology library (NumCosmo) [48].
First, following the same procedure as in VPL, we use the NcmFitESMCMC function to perform a Markov Chain Monte Carlo (MCMC) analysis given an ensemble sampler with affine invariance approach [49]. Thus, from Eq. (34) and to avoid any further assumptions on the astrophysical and cosmological dependencies of the SNe Ia, BAO and H(z) likelihoods, we reconstruct q(z) by fitting its coefficients along with Ω 0 k , the SNe Ia parame- ters, the drag scale (present in the BAO likelihood) and the Hubble parameter H 0 , i.e., where i = 0, ..., 11.
Here we consider three different cases regarding the spatial curvature. In the first two we fit Ω 0 k by assuming zero-mean Gaussian priors with standard deviation equal to 0.05 and 0.10, respectively. This choice is consistent with the Planck results [46]. The third case refers to a flat universe, i.e., we fix Ω 0 k = 0.0. As described in Sec. III, the "prior" on {q i } corresponds to the penalization factor in Eq. (34). By construction, the Riess et al. [34]  Thus, for each case, we ran the NcmFitESMCMC algorithm computing 5 × 10 6 sampling points in the 18-and 19-dimensional parameter spaces distributed among 50 chains. The convergence is attained as indicated by the multivariate potential scale reduction factor (MPSRF) of about 1.015 and the effective sample size (see [50] for details of the convergence tests and criteria). The variance of −2 ln L and all 19 parameters (18 in the flat case) also converged, e.g., Var(−2 ln L) 36.08 (flat case) which is consistent with a chi-squared distribution with 18 degrees of freedom (χ 2 18 ).

IV. RESULTS
We first analyze the ECs in the context of GR. From Eqs. (24)(25)(26)(27)(28), we have that the ECs are fulfilled at any redshift value if they are equal or greater than zero. Note that WEC2 correspond to the NEC in GR. Figures 1 and  2 show the MCMC results for the case where we consider a Gaussian prior on Ω 0 k with zero mean and scatter 0. 10 Fig. 1) is fulfilled in the entire redshift interval, and the DEC (lower panel) presents a violation within 99.7% CI just for roughly z 2.17 in the three study cases, which is a highly degenerated region due to the few number of data points.
In turn, the upper panel of Fig. 1 shows that NEC is violated within 99.7% CI for z 0.06, 1.3 z 1.8 and z 2.08, where the last is in the degenerated region of the q(z) [and H(z)] reconstruction. We note that these violation intervals are narrower than those obtained in Refs. [5,6] using only SNe Ia data and z ∈ [0, 1.0], where the violations were present for z 0.1 and z 0.8. Despite of using a larger SNe Ia sample and BAO and H(z) data, it is worth mentioning that the reconstruction method used in the present work is less restrictive than that in [5,6], 3 and it was calibrated such that the bias contributes only with 10% of the total error budget. The NEC, WEC and DEC results obtained for the flat and Ω 0 k = 0.0±0.05 cases are pretty similar to those presented in Fig. 1.
The most interesting result concerns SEC, since this is the only EC we expect to be violated, which is tightly linked with the recent accelerated expansion of the universe. Figure 2 shows the reconstructed deceleration function, namely the mean q(z) curve along with the 68.3%, 95.5% and 99.7% CIs. We note (upper panel) that the evidence for SEC violation takes place over the entire redshift interval. In fact, we obtain SEC's fulfillment within 68.3% CI just for 0.97 two cases, we have z ∈ [∼ 0.9, ∼ 1.1]). In VPL [24] this fulfillment was observed in the range 1.84 z 2.13. Despite the reconstruction methodology be the same used here, the differences result from some distinct and new BAO and H(z) data points (e.g., [29,30,33]) used in the analysis. Contrarily to the new constraints for NEC, SEC violation is stronger than those obtained in previous work. The lower panel of Fig. 2 [5,6], respectively, using only SNe Ia data. In short, the most current SNe Ia, BAO and H(z) data strengthens the evidence of an accelerated expansion of the universe. For instance, calculating the posterior for q min = min(q(z)) for z ∈ (0, 0.5) we found that q min < 0 for all points in our sample of the posterior. This means that the probability of finding q min ≥ 0 is smaller than 1/(5 × 10 6 ) (one in the number of posterior sampled points), which translates to at least 5.22σ confidence level.
As shown above, the fulfillment/violation of the ECs in the GR case can be directly tested since the bounds need just to be compared to constant values (dashed lines in Figs. 1 and 2). Regarding the class of ETGs considered in this work, this is just valid for NEC and SEC as they have the same form in GR as in these ETGs and, therefore, the results and analyses presented above (upper panel of Fig. 1 and Fig. 2) are also valid for these ETGs. This is to be expected, since these conditions were obtained from the convergence conditions imposed directly on R ab and our reconstruction method output is exactly R ab .
On the other hand, the remaining ECs, i.e., WEC1, WEC2 and DEC (Eqs. (26), (27) and (28), respectively), and also the fifth condition [Eq. (30)] involve not only q(z), E(z) and Ω 0 k , but also the arbitrary functions h t (z) and h s (z) of the modified gravity tensor H ab . Consequently, instead of checking whether a condition is satisfied or not, our methodology allows one to put constraints on these functions [and their combinations, say f (h t , h s )].
From the reconstructed q(z) and E(z) curves, we obtain now the observational constraints on the functions f (h t , h s ) by requiring that WEC1, WEC2, DEC and the fifth condition are fulfilled (see Eqs. (26), (27), (28) and (30)). Figure 3 shows the 68.3% and 99.7% CIs of the observational bounds obtained for the flat and Ω 0 k = 0 ± 0.1 cases. Note that WEC1 (left upper panel), WEC2 (right upper) and DEC (left lower) provide upper bounds on their respective functions f (h t , h s ) whereas the fifth condition (right lower panel) gives a lower bound. The shaded areas in all four panels of Fig. 3 indicates the values per redshift for which the respective conditions are violated. These forbidden areas are wider, i.e., more restrictive as the redshift values decrease due to the larger number of data points. Figure 4 shows the same results for the better constrained region z ∈ [0, 0.5]. We see that the GR thresholds, which correspond to the case where H ab ≡ 0, i.e., h t (z) = h s (z) = 0, of WEC1, WEC2, and DEC are fulfilled in the entire redshift interval.
An interesting result is verified for the fifth condition. The right lower panel of Fig. 4 shows that this condition is violated in GR theory within 99.7% for 0.01 z 0.26. Consequently, GR is excluded as a solution at least at this redshift interval. Recall that the fifth condition states the attractiveness of gravity for all ordinary matter. Since for GR the fifth condition is equivalent to SEC, the present result is consistent with that discussed previously (Fig. 2), where SEC violation in GR indicates that the universe presents an accelerated expansion, which must be driven by a cosmological constant or an exotic fluid with equation of state, for example, of the type p = wρ, with w < −1/3. To go further without choosing a specific ETG, we need to make some assumptions about h t and h s . For this it is useful to use the language of an effective fluid. First we define ρ ETG ≡ h t /(8πG) and p ETG ≡ −h s /(8πG). Rewriting DEC and the fifth condition using these variables, we can impose upper and lower bounds on Ω ETG ω ETG , i.e., where ω ETG ≡ p ETG /ρ ETG and Ω ETG ≡ 8πGρ ETG /(3H 2 0 ). If the DE was explained by a dark fluid instead of an ETG, the same bounds above would be obtained. Nonetheless, in this case the fifth condition requires again that all fluids but the DE fluid satisfy the SEC. In addition, the upper bound (derived from DEC) requires a similar assumption, all fluids but the DE one must satisfy the DEC. It is interesting to note that the last mentioned condition is not necessary in the ETG case since, by definition, DEC already imposes conditions only on the ordinary matter.
Assuming that Ω ETG > 0, we divide the bounds above by Ω ETG to obtain bounds on ω ETG . In Fig. 5 we plot these bounds considering Ω ETG constant. Even though we do not expect, a priori, that Ω ETG is constant, these bounds plotted with different values of Ω ETG show approximately the bounds one would obtain if Ω ETG smoothly evolved inside the chosen redshift values. This amounts to show how these bounds can impose constraints on the behavior of the ETG. For instance, in Fig. 5 we note that the lower the Ω ETG value, the higher the evidence to obtain a phantom-like behavior. Naturally, the value of Ω ETG today must be closer to the one estimated by current data (Ref. [46], for example, gives Ω ETG 0 ≈ 0.7). Nevertheless, for z > 0 the evolution of Ω ETG can take it to different values depending on the specific dynamic of the ETG and, therefore, these bounds show the restrictions depending on how this function evolves in time.

V. CONCLUSION
Numerous propositions of modified gravity theories and also of the DE equation of state in the context of GR have been introduced and discussed in the last two decades to explain the recent accelerated expansion of the universe. Therefore, one major task in cosmology is to constrain these models and also scrutinize their viability using observational data.
In this work we introduced a methodology to obtain observational bounds on a class of ETGs and also on the parameter space of a specific theory. This approach consists in requiring the fulfillment of the ECs, from which we obtain the theoretical bounds for the ETGs' functions and/or paramaters.
We derived the ECs for ETGs for which matter is minimally coupled to geometry [Eq. (7)]. Then, assuming a homogeneous and isotropic metric, we wrote these ECs in terms of the functions q(z) and E(z), and the parameter Ω 0 k . We also considered a fifth condition stating that gravity is attractive for all ordinary matter. Using the VPL model-independent reconstruction method [24] and SNe Ia, BAO, H(z) data, we obtained observational bounds on combinations of h t (z), h s (z) as given in Eqs. (26), (27), (28) and (30).
We first studied the ECs in GR. We verified that WEC and DEC are fulfilled in the redshift interval [0, 2.33]. NEC is violated at very low redshifts z 0.06 (and also at higher z values where the variance of the estimated curve is big, see Sec. IV). It is worth emphasizing that NEC violation is weaker than those obtained in [5,6]. On the other hand, the evidence of SEC violation obtained in the present work is stronger than the previous ones [5,6,24]. In particular, there is an indication bigger than 5.22σ CI of the recent accelerated expansion of the universe.
We also obtained the allowed/forbidden values for the WEC1, WEC2, DEC and fifth condition functions of h t (z) and h s (z), as showed in Figs.3 and 4. Particularly, we have that GR is excluded within 99.7% CI for z ∈ [0.01, 0.26] since it violates the fifth condition (equivalent to SEC).
The present study considered a general class of ETGs and, therefore, we obtained just lower and upper bounds for different functions of h t (z) and h s (z). However, applying this approach for a specific model, one will be able to obtain upper and lower bounds for a given function or parameters. At last, in this regard we provided an example of possible results. By assuming a positive Ω ETG , we obtained bounds on the effective ETG equation of state as showed in Fig. 5. This proves to be potentially useful to determine the behavior of the ETG in the reconstructed redshift interval, from a phantom like comportment ω ETG < −1 or the opposite direction ω ETG > −1 depending on how Ω ETG evolves in time. Therefore, we believe this method is a useful tool to constrain the parameter spaces of different ETGs. For instance, in Ref. [37] we considered an ETG whose modified gravity term, i.e., the tensor H ab , acts as a cosmological constant. In this context, we also studied two bimetric massive gravity theories putting constraints on their parameter space.