Welcome to visit Communications in Theoretical Physics,
Gravitation Theory, Astrophysics and Cosmology

f(T, B) gravity with statistically fitting of H(z)

  • S H Shekh , 1, * ,
  • N Myrzakulov , 2, 3 ,
  • A Bouali , 4 ,
  • A Pradhan , 5
Expand
  • 1 Department of Mathematics, S.P.M. Science and Gilani Arts, Commerce College, Ghatanji, Yavatmal, Maharashtra 445301, India
  • 2 L.N. Gumilyov Eurasian National University, Astana 010008, Kazakhstan
  • 3 Ratbay Myrzakulov Eurasian International Center for Theoretical Physics, Astana, 010009, Kazakhstan
  • 4Laboratory of Physics of Matter and Radiation, Mohammed I University, BP 717, Oujda, Morocco
  • 5Centre for Cosmology, Astrophysics and Space Science (CCASS), GLA University, Mathura-281 406, Uttar Pradesh, India

*Author to whom any correspondence should be addressed.

Received date: 2023-02-14

  Revised date: 2023-06-19

  Accepted date: 2023-07-04

  Online published: 2023-08-10

Copyright

© 2023 Institute of Theoretical Physics CAS, Chinese Physical Society and IOP Publishing

Abstract

Some recent developments (accelerated expansion) in the Universe cannot be explained by the conventional formulation of general relativity. We apply the recently proposed f(T, B) gravity to investigate the accelerated expansion of the Universe. By parametrizing the Hubble parameter and estimating the best fit values of the model parameters b0, b1, and b2 imposed from Supernovae type Ia, Cosmic Microwave Background, Baryon Acoustic Oscillation, and Hubble data using the Markov Chain Monte Carlo method, we propose a method to determine the precise solutions to the field equations. We then observe that the model appears to be in good agreement with the observations. A change from the deceleration to the acceleration phase of the Universe is shown by the evolution of the deceleration parameter. In addition, we investigate the behavior of the statefinder analysis, equation of state (EoS) parameters, along with the energy conditions. Furthermore, to discuss other cosmological parameters, we consider some well-known f(T, B) gravity models, specifically, f(T, B) = aTb + cBd. Lastly, we find that the considered f(T, B) gravity models predict that the present Universe is accelerating and the EoS parameter behaves like the ΛCDM model.

Cite this article

S H Shekh , N Myrzakulov , A Bouali , A Pradhan . f(T, B) gravity with statistically fitting of H(z)[J]. Communications in Theoretical Physics, 2023 , 75(9) : 095401 . DOI: 10.1088/1572-9494/ace3ae

1. Introduction

Currently, one of the most exciting problems in astrophysics and cosmology is finding the physical mechanism responsible for the cosmic acceleration of our Universe in late time [1, 2]. Various models have been proposed to describe the nature of this phenomenon. One of these is the cosmological constant, which is considered the phenomenologically simplest possibility when cold dark matter constitutes the cosmological standard model known as the ΛCDM. A number of dark energy cosmological models without the cosmological constant have been proposed to explain cosmic acceleration. There are a quintessence (canonical scalar field) [35], a phantom (non-canonical scalar field) [68], fermion fields [912], tachyon fields [13, 14], Chaplygin gas with a special equation of state (EoS) [15, 16] and so on. All these models have been successfully investigated in the framework of Riemannian gravitational theories with a Levi-Civita connection, where space-time is mediated by curvature in general relativity (GR). Additionally, the scientific community has stimulated interest in modifications of the Einstein GR action, in order to include a higher-order curvature invariant with respect to the curvature. Theories of modified gravity have attracted much attention in the explanation of both early-time and late-time acceleration [17, 18]. Specific models R2, RabRab and RabcdRabcd corrections are considered in the literature. These corrections to GR were found to be important and close to the Planck scale. Confrontation with observational data in the case of f(R) gravitational theory [19] was performed in [20, 21]. Additionally, comparisons with solar system data were performed in [22, 23]. Some other of the latest investigations in f(R) gravity are given in [2427].
An alternative gravitation theory that describes gravitational interactions in terms of torsion is the known teleparallel equivalent of general relativity (TEGR) introduced by Einstein [28, 29]. This theory is conceptually different from the GR, at least at the level of the gravitational equations of motion. Linear frames and tetrads are two basic objects that will be fundamental in the construction of TEGR. In this theory, the Levi-Civita connection is replaced with a so-called Weitzenbock connection. In recent years, extending torsional gravity, namely f(T) gravitational theory has been extensively proposed and investigated in the literature [30, 31]. The f(T) theory of gravity construction can provide a theoretical interpretation of the early-time and late-time acceleration. The Lagrangian of f(T) is taken to be a nonlinear arbitrary function of the TEGR. Note that local Lorentz invariance in the formulation of f(T) gravity is strongly restricted.
Extended and modified gravities have attracted the attention of many cosmologists, because they provide a geometric and systematic approach to the explanation of cosmological observations. Recently, some authors proposed an interesting extension of modified f(T) gravity namely $f(T,{ \mathcal T })$ gravity, where T is torsion scalar and is the ${ \mathcal T }$ trace of the matter energy-momentum tensor [32]. Inspired by the torsional formulation, the teleparallel equivalent of Gauss-Bonnet gravity has been constructed and is known as f(T, TG) [33]. These gravity theories are tested by observational tests at cosmological scales.
Boundaries/boundary terms are fundamental concepts in many areas of theoretical physics. Our interest is another extension of the f(T) theory, in which the Lagrangian takes the form of f(T, B), where T and B are scalars (T is the torsion scalar and B is the boundary term) that characterize the equivalency with GR [34]. Paper [35] investigated finding the exact solutions for spherically symmetric Lemaitre-Tolman-Bondi dust models for f(T, B) gravity. Paper [36] explored the cosmological evolution of the Universe with the Lagrange multiplier in teleparallel f(T, B) gravity. The authors performed dynamical analysis and investigated their stationary points. The f(T, B) teleparallel gravity in a five-dimensional brane scenario was studied and the gravitational perturbations were investigated in [37]. Work [38] discussed classes of exact and perturbative spherically symmetric solutions in modified teleparallel f(T, B) gravity, while in [39] the authors examined the second law of thermodynamics in f(T, B) theory and using cosmological reconstruction technique, showed that some models of f(T, B) can mimic the de-Sitter, power-law and ΛCDM models. Analysis of the Tsallis holographic dark energy and energy conditions in teleparallel f(T, B) gravity theory were considered in [40, 41].
In cosmology, Markov Chain Monte Carlo (MCMC) simulations [42, 43] are used to find the best-fitting distance modulus for each cosmological model and for the most probable free parameter value with given certain physical constraints. In this study, we explore the parametrization of the Hubble parameter to obtain the scenario of an accelerating universe. The best fit values of model parameters were obtained from recent observational data: the Hubble datasets H(z), consisting of a list of 57 measurements that were compiled from the differential age method and others; the Type Ia supernovae sample called Pantheon datasets, consisting of 1048 points covering the redshift range 0.01 < z < 2.26; and the Baryon Acoustic Oscillation fs(BAO) datasets, consisting of six points [44, 45]. Our analysis uses the combination of the H(z), Pantheon samples and BAO datasets to constrain the cosmological model. In order to recreate the shape of the f(T) alteration in a model-independent manner, Cai et al [46] developed the Gaussian processes analysis for the case of f(T) gravity, utilising as the sole input the observational data sets of the Hubble function measurements H(z). Santos et al [47] have investigated the observational constraints on f(T) gravity from model-independent data. A model-independent method with phantom dividing line crossing in Weyl-type f(Q, T) gravity has been recently investigated by Koussour [48]. Mu et al [49] discussed the most recent cosmic observations, including Pantheon + SNe Ia samples, to reconstruct the modified gravity, which is defined by the modified factor μ in linear matter density perturbation theory, in a completely data-driven and model-independent manner. Model-independent constraints on modified gravity from current data and from the Euclid and SKA future surveys are investigated by Taddei et al [50]. In the current study, model-independent constraints on the modified f(T, B) gravity have been examined using data from the H(z), Pantheon samples, and BAO datasets.
The structure of this article is as follows: We present a brief description of f(T, B) gravity in section 2. In section 3, a statistical fitting of H is given. Observational constraints are studied in section 4, while section 5 covers certain particular cosmographic parameters. Some specific cosmological models of f(T, B) gravity are investigated in section 6. Finally, we discuss our findings in section 7.

2. f(T, B) gravity with homogeneous space-time

Basically, the modification of the left side of the Einstein field equations (by some arbitrary function) are called modified theories of gravity, which are the probable access to describe the accelerating expansion of the Universe. A few inspirations to walk around the modified theories of gravity are f(R), f(T), f(R, T), f(G), f(R, G), f(T, B). Among these, one model which is established on the T and B or the combination of f(R) and f(T) gravity, namely f(T, B) gravity [34], is as follows:
$\begin{eqnarray}S=\int e\left(\displaystyle \frac{f(T,B)}{{\kappa }^{2}}+{L}_{m}\right){{\rm{d}}}^{4}x,\end{eqnarray}$
where κ2 = 8πG and f(T, B) is a function of the torsion scalar T and the boundary term $B=\tfrac{2}{e}{\partial }_{\mu }({{eT}}^{\mu })$, in which ${T}_{\mu }={T}_{v\mu }^{\mu }$. Lm and $e=\det ({e}_{\mu }^{i})$ are the matter action and determinant of the tetrad components, respectively.
By varying the action given in equation (1) with respect to the tetrad, the field equation is defined as:
$\begin{eqnarray}\begin{array}{l}2e{\delta }_{\nu }^{\lambda }{{\rm{\nabla }}}^{\mu }{{\rm{\nabla }}}_{\mu }{\partial }_{B}f-2e{{\rm{\nabla }}}^{\lambda }{{\rm{\nabla }}}_{\nu }{\partial }_{B}f\\ \quad +\,{eB}{\partial }_{B}f{\delta }_{\nu }^{\lambda }+4e\left({\partial }_{\mu }{\partial }_{B}f+{\partial }_{\mu }{\partial }_{T}f\right){S}_{\nu }^{\mu \lambda }\\ \quad +\,4{e}_{\nu }^{a}{\partial }_{\mu }\left({{eS}}_{a}^{\mu \lambda }\right){\partial }_{T}f-4e{\partial }_{{Tf}}{T}^{\sigma }{}_{\mu \nu }{S}_{a}^{\lambda \mu }\\ \quad -\,{ef}{\delta }_{\nu }^{\lambda }=16\pi {{eT}}_{\nu }^{\lambda }.\end{array}\end{eqnarray}$
In the standard cosmological principle, our universe is filled with perfect fluid. Consequently, the energy-momentum tensor for this perfect fluid is defined as the following form:
$\begin{eqnarray}{T}_{v}^{\lambda }=(\rho +p){u}_{v}{u}^{\lambda }-p{\delta }_{v}^{\lambda },\end{eqnarray}$
where p and ρ are the pressure and the energy density of the fluid, respectively. Note that ${u}^{\nu }=\left(0,0,0,1\right)$ is the four-velocity vector of the fluid with uνuν = 1. The non-vanishing elements of the energy-momentum tensor are
$\begin{eqnarray}{T}_{0}{}^{0}=\rho ,{T}_{1}{}^{1}={T}_{2x}{}^{2}={T}_{3}{}^{3}=-p.\end{eqnarray}$
Here, we explore the spatially, isotropic and homogeneous Friedmann−Robertson−Walker (FRW) line element as
$\begin{eqnarray}{{ds}}^{2}=-{{dt}}^{2}+{\delta }_{{ij}}{g}_{{ij}}{{dx}}^{i}{{dx}}^{j},\,\,\,\,i,j=1,2,3,\ldots ..N,\end{eqnarray}$
where gij are the function of (t, x1, x2, x3) and t refers to the cosmological time. In four dimensional FRW space-time, from above equation we have
$\begin{eqnarray}{\delta }_{{ij}}{g}_{{ij}}={a}^{2}(t,x).\end{eqnarray}$
The above relations show that, for the FRW universe, all three metrics are equal (i.e g11 = g22 = g33 = a2(t, x)). For the FRW metric (5), the (0 − 0) and the (ii) components of Einstein’s gravitational equation (2) become the following forms:
$\begin{eqnarray}-3{H}^{2}(3\partial {}_{B}f+2{\partial }_{T}f)+3H{\partial }_{B}\dot{f}-3\dot{H}{\partial }_{B}f+\displaystyle \frac{1}{2}f={\kappa }^{2}\rho ,\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}-(3{H}^{2}+\dot{H})(3{\partial }_{B}f+2{\partial }_{T}f)-2H{\partial }_{T}\dot{f}\\ \quad +\,{\partial }_{B}\ddot{f}+\displaystyle \frac{1}{2}f=-{\kappa }^{2}p.\end{array}\end{eqnarray}$
The overhead dot denotes the derivative with respect to cosmological time t.
The torsion scalar for the FRW line element (5) is
$\begin{eqnarray}T=-6{H}^{2}.\end{eqnarray}$
The Ricci curvature scalar and torsion scalar are related as
$\begin{eqnarray}R=-T+B.\end{eqnarray}$
The curvature scalar R is used in the creation of the entire conventional relativity theory. However, in f(T, B) gravity theory, it is made up of the torsion scalar T and the boundary term B. Be aware that the boundary term [51] makes these models unique. The boundary term B for the FRW metric (5) has
$\begin{eqnarray}B=-6\left(\dot{H}+3{H}^{2}\right).\end{eqnarray}$
The curvature scalar R from (10) is obtained as
$\begin{eqnarray*}R=-6\left(\dot{H}+2{H}^{2}\right).\end{eqnarray*}$
The standard first and second Friedman equations are
$\begin{eqnarray}3{H}^{2}={\kappa }^{2}{\rho }_{\mathrm{tot}},\end{eqnarray}$
$\begin{eqnarray}2\dot{H}+3{H}^{2}=-{\kappa }^{2}{\bar{p}}_{\mathrm{tot}}.\end{eqnarray}$
The parameters ρtot and ptot in f(T, B) gravity are found as
$\begin{eqnarray}{\rho }_{\mathrm{tot}}=\rho +{\rho }_{d},\end{eqnarray}$
$\begin{eqnarray}{p}_{\mathrm{tot}}=p+{p}_{d}.\end{eqnarray}$
Using equations (7), (8) and (12)–(15), we find the effective isotropic pressure pd and effective energy density ρd towards f(T, B) gravity as
$\begin{eqnarray}\begin{array}{l}{\kappa }^{2}{\rho }_{d}=3{H}^{2}(1+3{\partial }_{B}f+2{\partial }_{T}f)\\ \quad -\,3H{\partial }_{B}\dot{f}+3\dot{H}{\partial }_{B}f-\displaystyle \frac{1}{2}f,\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}{\kappa }^{2}{p}_{d}=-3{H}^{2}(1+3{\partial }_{B}f+2{\partial }_{T}f)\\ \quad -\,\dot{H}(2+3{\partial }_{B}f+2{\partial }_{T}f)-2H{\partial }_{T}\dot{f}+{\partial }_{B}\ddot{f}+\displaystyle \frac{1}{2}f,\end{array}\end{eqnarray}$
where the quantities pd and ρd are the parts of the pressure and energy density respectively that appear from f(T, B) gravity. The expressions of ρd and pd given in equations (16) and (17) slightly differ from those in the equations given in (7) and (8) in view of the standard Friedmann equations provided in (12) and (13). Hence, we consider that equations (16) and (17) are the pressure and density for f(T, B) gravity. Now, we have explored the f(T, B) gravity model of the form
$\begin{eqnarray}f=f(T,B)={{aT}}^{b}+{{cB}}^{d}.\end{eqnarray}$
For the above-said model of f(T, B) gravity, the set of field equations (16) and (17) become
$\begin{eqnarray}\begin{array}{l}{\rho }_{d}=\displaystyle \frac{1}{2{\kappa }^{2}}\left\{a{6}^{b}(2b+1){\left(-{H}^{2}\right)}^{b}\right.\\ \quad \left.+\,c(d+1){\left(-18{H}^{2}-6\dot{H}\right)}^{d}+2{f}^{{\prime} }\right\},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{l}{p}_{d}=\displaystyle \frac{1}{6{\kappa }^{2}}\left\{a{6}^{b}{\left(-{H}^{2}\right)}^{b-1}\left((6b+3){H}^{2}+2b\dot{H}\right)\right.\\ \quad \left.-\,3c(d+1){\left(-18{H}^{2}-6\dot{H}\right)}^{d}-6{f}^{{\prime\prime} }\right\}.\end{array}\end{eqnarray}$
The next important parameter of the Universe is the equation of state (EoS). As we know, the equation of state parameter is associated with energy density and pressure, which classifies the expansion of the Universe. If the value of the equation of state parameter is exactly 1 then it represents the static fluid era. If it is 0 then the Universe represents the matter-dominated era, while for $\omega =\tfrac{1}{3}$, the Universe represents the radiation-dominated era. Whereas if it lies in between −1 to 0, i.e. −1 < ω < 0, then the Universe shows the quintessence era, while ω = − 1 shows the cosmological constant, i.e., the ΛCDM era and the phantom era is observed when ω < − 1. From the above equations, we obtained the equation of state parameter ω as
$\begin{eqnarray}\omega =-\,\displaystyle \frac{a{6}^{b}{\left(-{H}^{2}\right)}^{b}\left((6b+3){H}^{2}+2b\dot{H}\right)+3c(d+1){H}^{2}{\left(-18{H}^{2}-6\dot{H}\right)}^{d}+6{f}^{{\prime\prime} }{H}^{2}}{3{H}^{2}\left(a{6}^{b}(2b+1){\left(-{H}^{2}\right)}^{b}+c(d+1){\left(-18{H}^{2}-6\dot{H}\right)}^{d}+2{f}^{{\prime} }\right)}\end{eqnarray}$

3. Statistical fitting of H

Note that in general, regression equations express the relationship between two variables. One of our goals is to explore linear and nonlinear regression equations. In the present study, we examine the relationship between the regression equation with the redshift parameter z as a predictor and the Hubble parameter H. To obtain a suitable relation between the redshift z and Hubble parameter H, we have to get the best possible regression relation for our model. To find the best fit we need to calculate R2, where ${R}^{2}=\tfrac{{SSR}}{{SST}}$. The term SST is the sum of squares total. This term mathematically indicates the sum of square deviations of the value of the dependent variables around their mean. And the SSR term stands for the regression sum of squares. This term mathematically indicates the sum of square differences between the sample mean of the prediction and the regression predictions. An application of the regression method presented in [52]. Now we will try to study three models in the following forms:
$\begin{eqnarray}\begin{array}{rcl}H(z) & = & {b}_{0}+{b}_{1}z+{b}_{2}{z}^{2}\,\,\,{H}_{1}(z)={b}_{0}+{b}_{1}z\,\,\,{\rm{and}}\\ {H}_{2}(z) & = & {b}_{0}{e}^{({b}_{1}z)}.\end{array}\end{eqnarray}$
Among the regression relations, whichever gives the highest value of R2 close to 1, will be assumed to be the best model. In equations (22), z and H(z), H1(z), H2(z) represent the redshift and Hubble parameter, respectively. Also, b0 represents the regression constant and b1, b2 are the regression coefficients of the given models. The three models of regression represented in equations (22) are now fitted to the Supernovae type Ia, Cosmic Microwave Background (CMB), Baryon Acoustic Oscillation and Hubble data sets. The three regression models and the corresponding R2 towards H(z), H1(z), H2(z) are given in table 1. All of the models have been proven to be almost equally efficient through the high value of R2. Four decimal places have been retained to understand the relative efficiency of each model. Although the three models have equal efficiency, the H(z) model is found to be more accurate than H1(z) and H2(z) (see table 1).
Table 1. The three regression models and their corresponding R2.
Regression model R2
H(z) = b0 + b1z + b2z2 0.9404
H1(z) = b0 + b1z 0.9394
${H}_{2}(z)={b}_{0}{e}^{({b}_{1}z)}$ 0.9246
Hence, considering the highest value of R2, the first model of the form is considered to be the best model to establish the functional relationship between z and H(z).
$\begin{eqnarray}H(z)={b}_{0}+{b}_{1}z+{b}_{2}{z}^{2}.\end{eqnarray}$
The expression of H(z) with its present value H = 0 is obtained as
$\begin{eqnarray}H(z)={H}_{0}+{b}_{1}z+{b}_{2}{z}^{2}.\end{eqnarray}$
Table 2. Summary of the best fit and the mean values of the cosmological parameters.
Model Par Prior Best fit Mean
ΛCDM Ωm [0.001, 1] ${0.312308}_{-0.00607812}^{+0.00607812}$ ${0.312583}_{-0.00607362}^{+0.00607362}$
h [0.4, 1] ${0.678603}_{-0.00447778}^{+0.00447778}$ ${0.678435}_{-0.00447417}^{+0.00447417}$
Ωbh2 [0.001, 0.08] ${0.0224102}_{-0.000134672}^{+0.000134672}$ ${0.0224002}_{-0.00013449}^{+0.00013449}$

Polynomial model b1 [0.001, 100] ${55.87}_{-1.4800}^{+1.4800}$ ${55.0}_{-1.500}^{+1.500}$
b2 [0.001, 100] ${0.961}_{-0.0171}^{+0.0171}$ ${0.960}_{-0.017}^{+0.017}$
H0 [40, 100] ${66.16}_{-1.2900}^{+1.2900}$ ${66.2}_{-1.3000}^{+1.3000}$
Table 3. Summary of the ${\chi }_{{red}}^{2}$ and ${\chi }_{\min }^{2}$.
Model ${{\chi }_{{tot}}^{2}}^{\min }$ ${\chi }_{{red}}^{2}$
ΛCDM 1073.9795 0.9816
Polynomial model 1125.2628 1.0195
The first derivative of H(z) is observed as
$\begin{eqnarray}\dot{H}(z)=-(1+z)({b}_{1}+2{b}_{2}z)H(z).\end{eqnarray}$

4. Observational constraints

We have briefly described the f(T, B) gravity and solved the field equation with a new parametrization of the Hubble parameter. In order to extract the best fit values, the considered form of H(z) was constrained by SNIa from Pantheon, CMB from Planck 2018, BAO and 36 data points from Hubble datasets using the MCMC approach. In what follows, we describe in detail the methodology adopted and data used in our analysis. Figures 1 and 2 represent the behavior of the distance modulus and Hubble parameter in terms of the redshift z. The results of our study are shown in the contour plots (two-dimensional) with 1 − σ and 2 − σ errors (figure 3). In figure 3 the chains are run sufficiently to get convergent results. In addition, from the MCMC contour plots, one can notice that the posteriors are smooth with only one maximum. We have again constrained our model parameters by running 10 Markov chains. Furthermore, we have performed a Gelman-Rubin convergence test (please see the table below). The R-1 is smaller than 0.05 for all of the parameters, indicating that our chains have completely reached the convergence region.
R-1 b1 b2 H0
0.017025 0.01672 0.01692
Figure 1. The distance modulus behavior in terms of the redshift z.
Figure 2. The Hubble parameter behavior in terms of the redshift z.
Figure 3. The 1 − σ and 2 − σ confidence contours obtained from SNIa+BAO+H(z) data obtained for the statistical fitting model.

4.1. Supernovae type Ia

We use supernovae from the Pantheon compilation made of 1048 spectroscopically confirmed Type Ia Supernovae, distributed in the redshift range 0.01 < z < 2.26 [53]. So far, the Pantheon compilation is the largest and it contains measurements from different supernovae surveys such as SDSS, SNLS and HST. For the purpose of estimating the best fit parameters, we compute the chi-square
$\begin{eqnarray}{\chi }_{{SN}}^{2}={\left({\mu }_{{obs}}-{\mu }_{{th}}\right)}^{T}.{C}_{{Pantheon}}^{-1}.({\mu }_{{obs}}-{\mu }_{{th}}),\end{eqnarray}$
where μobs is the observed distance modulus and CPantheon is the covariance matrix of the Pantheon data.

4.2. Cosmic microwave background

The χ2 for CMB is expressed as follows:
$\begin{eqnarray}{\chi }_{\mathrm{CMB}}^{2}({b}_{1},{b}_{2},h)={{\bf{X}}}_{\mathrm{CMB}}^{T}({b}_{1},{b}_{2},h).{{\bf{C}}}_{\mathrm{CMB}}^{-1}.{{\bf{X}}}_{\mathrm{CMB}}({b}_{1},{b}_{2},h),\end{eqnarray}$
where the CMB covariance matrix CCMB [54].

4.3. Baryon acoustic oscillation

The Baryon Acoustic Oscillation seen by galaxy surveys plays a crucial role in the determination of the evolution of the Universe. From BAO datasets we can measure the angular diameter distance DA(z) using clustering perpendicular to the line of sight. In addition, we can measure the expansion rate of the Universe H(z), through clustering along the line of sight. The angular diameter distance DA(z) and the spherical averaged scale DV are related to H(z) as follows
$\begin{eqnarray}{D}_{V}(z)={\left[{\left(1+z\right)}^{2}{D}_{A}^{2}(z)\displaystyle \frac{{cz}}{H(z)}\right]}^{1/3}.\end{eqnarray}$
The peak positions of BAO are in general given in terms of DV(z)/rs(z), DA(z)/rs(z) and H(z)rs(z) measured at the drag epoch zdrag, i.e., where baryons were released from photons. In this paper, we use correlated BAO data and uncorrelated ones [5557]. Hence, the total chi-square for BAO, ${\chi }_{{BAO}}^{2}$, is expressed as
$\begin{eqnarray}\begin{array}{l}{\chi }_{{BAO}}^{2}({b}_{1},{b}_{2},h)=\ {\chi }_{6{dFGS}}^{2}+{\chi }_{{SDSS}}^{2}+{\chi }_{{BOSS}-{LOWZ}}^{2}\\ \quad +\,{\chi }_{{BOSS}-{CMASS}}^{2}+{\chi }_{{WiggleZ}}^{2}+{\chi }_{{BOSS}-{DR}12}^{2}.\end{array}\end{eqnarray}$

4.4. Hubble data

For tighter constraints, we also make use of the Hubble measurements H(z). In general, the Hubble rate can be inferred either from the clustering of galaxies/quasars by measuring the BAO in the radial direction [58], or from the differential age method. Both methods lead to a compilation of 36 data points of the Hubble parameter H(z). The chi-square of the Hubble measurements is given by
$\begin{eqnarray}{\chi }_{H(z)}^{2}({b}_{1},{b}_{2},h)=\displaystyle \sum _{i=1}^{36}{\left[\displaystyle \frac{{H}_{{obs},i}-H({z}_{i},{b}_{1},{b}_{2},h)}{{\sigma }_{H,i}}\right]}^{2}.\end{eqnarray}$
Finally, the total chi-square, ${\chi }_{{tot}}^{2}$, is given by the sum of all the χ2 previously defined:
$\begin{eqnarray}{\chi }_{{tot}}^{2}({b}_{1},{b}_{2},h)={\chi }_{{SN}}^{2}+{\chi }_{\mathrm{CMB}}^{2}+{\chi }_{{BAO}}^{2}+{\chi }_{H(z)}^{2}.\end{eqnarray}$
In the previous sections, we have briefly described the f(T, B) gravity and solved the field equation with a new parametrization of the Hubble parameter. The considered form of H(z) contains two model parameters b1 and b2, which have been constrained through some observational data for further analysis. We have used some external datasets, such as an observational Hubble dataset of a recent compilation of 57 data points, the Pantheon compilation of SNeIa data with 1048 data points, and also the Baryonic Acoustic Oscillation dataset with six data points, to obtain the best fit values for these model parameters (see tables 2 and 3), in order to validate our technique. 

5. Some cosmographic parameters

Geometrical parameters played a significant role in analyzing the model for any gravity theory. In this section we discuss some cosmographic dimensionless parameters such as the deceleration q(z) and jerk j(z) parameter. The deceleration parameter q(z) is a quantity that shows how the expansion rate changes over time. The jerk parameter j(z) defines the rate of change of acceleration.
The deceleration parameter q(z) in H(z) form is observed as
$\begin{eqnarray}q(z)=-1-\frac{\dot{H}}{{H}^{2}}.\end{eqnarray}$
According to astrophysical data, our universe is in the stage of cosmic acceleration. To understand the entire cosmic history of the Universe, a cosmological model must include both the decelerated and accelerated phases of the expansion. Consequently, it is necessary to investigate the behavior of the deceleration parameter q(z). The behavior of q(z) for the associated values of the model parameters constrained are shown in figure 4. It is clear that the parameter q(z) shows a transition from a decelerated to an accelerated phase. Furthermore, the range at which a transition takes place resembles that of recent observations [59].
Figure 4. Deceleration parameter $\left(q(z)\right)$ versus redshift $\left(z\right)$.
In particular, the jerk parameter, a dimensionless third derivative of the scale factor a(t) with respect to cosmic time t, can provide us with the simplest approach to search for departures from the ΛCDM model. It is defined as $j(t)=\tfrac{\dddot{a}}{{{aH}}^{3}}$ [60]. In terms of redshift z and the deceleration parameter q(z), the jerk parameter j(z) can be observed as
$\begin{eqnarray}j(z)=(1+z)\frac{{dq}}{{dz}}+q(1+2q).\end{eqnarray}$
Blandford et al [61] described how the jerk parameterization provides an alternative and a convenient method to describe cosmological models close to the ΛCDM model. A powerful feature of j(z) is that for the ΛCDM model j(z) = 1 (constant) always. It should be noted here that Sahni et al [62, 63] drew attention to the importance of j(z) for discriminating different dark energy models, because any deviation from the value of j(z) = 1 would favor a non-ΛCDM model. The behavior of j(z) for the associated values of model parameters constrained is shown in figure 5. It is clear that the parameter j(z) shows that the jerk parameter is positive throughout the evolution and finally acquires the value that tends to one, which describes the ΛCDM model.
Figure 5. Jerk parameter $\left(j(z)\right)$ versus redshift.

6. Some specific models of f(T, B) gravity

One method to solve the equations is by assuming a form of f(T, B) gravity. This way seems more physically interesting because one can propose a specific form of theory f(T, B) gravity. A similar analysis can be done in f(T, B) for the different b and d provided in (18). In order to illustrate this, let us consider some specific models. From the standard considered model (18), we will study the behaviors of energy density and the equation of state parameters, as well the energy conditions for specific models.

6.1. Model-I

In this subsection, we define the f(T, B) gravity model towards b = 2 and d = 2. With this model the expressions of ρd and ωd from equations (19) and (21) are of the form
$\begin{eqnarray}{\rho }_{d}=18\left(5{{aH}}^{4}+3c{\left(3{H}^{2}+\dot{H}\right)}^{2}\right)+{f}_{1}.\end{eqnarray}$
The behavior of energy density of the statistical f(T, B) gravity model versus redshift is presented in figure 6. and it is observed that the energy density is consistently non-negative and increases with the passage of redshift.
Figure 6. Energy density parameter of statistical f(T, B) = aT2 + cB2 gravity model $\left({\rho }_{d}\right)$ versus redshift $\left(z\right)$.
The equation of state parameter is obtained as
$\begin{eqnarray}{\omega }_{d}=-1+\frac{-24{{aH}}^{2}\dot{H}+{f}_{1}-{f}_{2}}{18\left(5{{aH}}^{4}+3c{\left(3{H}^{2}+\dot{H}\right)}^{2}\right)+{f}_{1}}.\end{eqnarray}$
Equation (36) represents the expression for the equation of state parameter of the statistical f(T, B) gravity model and its behavior is clearly shown in figure 7 with redshift. One can see in figure 7 that at the late Universe (z > 0) towards the value of the equation of state parameter of statistical f(T, B) gravity model is less than −1, which represents the Universe involving phantom field dark energy era, while at z = 0 and −1, the equation of state parameter of statistical f(T, B) gravity model is −1. Hence for both z = 0 and −1 the Universe involves the ΛCDM era, whereas in the early Universe (z > 0) it consists of both barotropic as well as of a quintessence field dark energy era.
Figure 7. Equation of state parameter of statistical f(T, B) = aT2 + cB2 gravity model $\left({\omega }_{d}\right)$ versus redshift $\left(z\right)$.
Energy conditions The most famous energy conditions NEC, DEC and SEC are observed as
NEC:
$\begin{eqnarray}{p}_{d}+{\rho }_{d}=-24{{aH}}^{2}\dot{H}+{f}_{1}-{f}_{2},\end{eqnarray}$
DEC:
$\begin{eqnarray}\begin{array}{l}{\rho }_{d}-{p}_{d}=12\left(a\left(15{H}^{4}+2{H}^{2}\dot{H}\right)\right.\\ \quad \left.\,+\,9c{\left(3{H}^{2}+\dot{H}\right)}^{2}\right)+{f}_{1}+{f}_{2},\end{array}\end{eqnarray}$
SEC
$\begin{eqnarray}\begin{array}{l}{\rho }_{d}+3{p}_{d}={f}_{1}-3\left(12a\left(5{H}^{4}+2{H}^{2}\dot{H}\right)\right.\\ \quad \left.\,+\,36c{\left(3{H}^{2}+\dot{H}\right)}^{2}+{f}_{2}\right).\end{array}\end{eqnarray}$
Equations (37), (38) and (39) represent the expression for energy conditions of the statistical f(T, B) gravity model. The behavior of the energy conditions is clearly shown in figure 8 with redshift. In the present universe, it is observed that NEC and DEC are well satisfied throughout cosmic evolution. However, the SEC is violated for our model. The violation of the SEC is due to the accelerated expansion of the Universe.
Figure 8. Null, Dominant and Strong energy conditions of f(T, B) = aT2 + cB2 gravity models versus redshift $\left(z\right)$.

6.2. Model-II

In this subsection, we define the f(T, B) gravity model towards b = 2 and d = 3. With this model the expressions of ρd and ωd from equations (19) and (21) are of the form
$\begin{eqnarray}{\rho }_{d}=18\left(5{{aH}}^{4}-24c{\left(3{H}^{2}+\dot{H}\right)}^{3}\right)+{f}_{3}.\end{eqnarray}$
The behavior of the energy density of the statistical f(T, B) gravity model versus redshift is presented in figure 9. and it is observed that the energy density is consistently non-negative and increases with the passage of redshift, which is the same as that of model-I.
Figure 9. Energy density parameter of statistical f(T, B) = aT2 + cB3 gravity model $\left({\rho }_{d}\right)$ versus redshift $\left(z\right)$.
The equation of state parameter is obtained as
$\begin{eqnarray}{\omega }_{d}=-\frac{6a\left(15{H}^{4}+4{H}^{2}\dot{H}\right)-432c{\left(3{H}^{2}+\dot{H}\right)}^{3}+{f}_{4}}{18\left(5{{aH}}^{4}-24c{\left(3{H}^{2}+\dot{H}\right)}^{3}\right)+{f}_{3}}.\end{eqnarray}$
Equation (41) represents the expression for the EoS parameter of the statistical f(T, B) gravity model and its behavior is clearly shown in figure 10 with redshift. One can see in figure 10 that its behavior also is the same as that of model-I, i.e., at the late Universe (z > 0) toward the value of the equation of state parameter of the statistical f(T, B) gravity model is less than −1, which represents the Universe involving phantom field dark energy era, while at z = 0 and −1, the equation of state parameter of the statistical f(T, B) gravity model is -1. Hence for both z = 0 and −1 the Universe involves the ΛCDM era, whereas in the early Universe (z > 0) it consists of both barotropic as well as of a quintessence field dark energy era which resembles model-I.
Figure 10. Equation of state parameter of statistical f(T, B) = aT2 + cB3 gravity model $\left({\rho }_{d}\right)$ versus redshift $\left(z\right)$.
Energy conditions
The most famous energy conditions NEC, DEC and SEC are observed as NEC:
$\begin{eqnarray}{p}_{d}+{\rho }_{d}=-24{{aH}}^{2}\dot{H}+{f}_{3}-{f}_{4},\end{eqnarray}$
DEC:
$\begin{eqnarray}\begin{array}{l}{\rho }_{d}-{p}_{d}=12a\left(15{H}^{4}+2{H}^{2}\dot{H}\right)\\ \,-\,864c{\left(3{H}^{2}+\dot{H}\right)}^{3}+{f}_{3}+{f}_{4},\end{array}\end{eqnarray}$
SEC
$\begin{eqnarray}\begin{array}{l}{\rho }_{d}+3{p}_{d}=-36a\left(5{H}^{4}+2{H}^{2}\dot{H}\right)\\ \,+\,864c{\left(3{H}^{2}+\dot{H}\right)}^{3}+{f}_{3}-3{f}_{4},\end{array}\end{eqnarray}$
Equations (41), (42) and (43) represent the expression for energy conditions of the statistical f(T, B) gravity model. The behavior of the energy conditions is clearly shown in figure 11 with redshift. In the present universe, it is observed that NEC and DEC are well satisfied throughout cosmic evolution. However, the SEC is violated for our model. The violation of the SEC is due to the accelerated expansion of the Universe.
Figure 11. Null, Dominant and Strong energy conditions of f(T, B) = aT2 + cB3 gravity models versus redshift $\left(z\right)$.

6.3. Model-III

In this subsection, we define the f(T, B) gravity model towards b = 3 and d = 2. With this model the expressions of ρd and ωd from the equations (19) and (21) are of the form
$\begin{eqnarray}{\rho }_{d}=54\left(c{\left(3{H}^{2}+\dot{H}\right)}^{2}-14{{aH}}^{6}\right)+{f}_{5}.\end{eqnarray}$
The behavior of energy density of the statistical f(T, B) gravity model versus redshift is presented in figure 12. and it is observed that the energy density is consistently non-negative and increases with the passage of redshift, which is the same as that of model-I.
Figure 12. Energy density parameter of statistical f(T, B) = aT3 + cB2 gravity model $\left({\rho }_{d}\right)$ versus redshift $\left(z\right)$.
The equation of state parameter is obtained as
$\begin{eqnarray}{\omega }_{d}=-1+\frac{216{{aH}}^{4}\dot{H}+{f}_{5}-{f}_{6}}{54\left(c{\left(3{H}^{2}+\dot{H}\right)}^{2}-14{{aH}}^{6}\right)+{f}_{5}}.\end{eqnarray}$
Equation (45) represents the expression for the EoS parameter of the statistical f(T, B) gravity model and its behavior is clearly shown in figure 10 with redshift. One can see in figure 13 that its behavior is also the same as model-I, i.e., at the late Universe (z > 0) toward the value of the equation of state parameter of the statistical f(T, B) gravity model is less than −1, which represents the Universe involving phantom field dark energy era, while at z = 0 and −1, the equation of state parameter of statistical f(T, B) gravity model is -1. Hence for both z = 0 and −1 the Universe involves the ΛCDM era, whereas in the early Universe (z > 0) it consists of both barotropic as well as of a quintessence field dark energy era.
Figure 13. Equation of state parameter of the statistical f(T, B) = aT3 + cB2 gravity model $\left({\rho }_{d}\right)$ versus redshift $\left(z\right)$.
Energy conditions
The most famous energy conditions NEC, DEC and SEC are observed as NEC:
$\begin{eqnarray}{p}_{d}+{\rho }_{d}=216{{aH}}^{4}\dot{H}+{f}_{5}-{f}_{6},\end{eqnarray}$
DEC:
$\begin{eqnarray}\begin{array}{l}{\rho }_{d}-{p}_{d}=108\left(c{\left(3{H}^{2}+\dot{H}\right)}^{2}-2{{aH}}^{4}\left(7{H}^{2}+\dot{H}\right)\right)\\ \,+\,{f}_{5}+{f}_{6},\end{array}\end{eqnarray}$
SEC
$\begin{eqnarray}\begin{array}{l}{\rho }_{d}+3{p}_{d}={f}_{5}-3\left(-72{{aH}}^{4}\left(7{H}^{2}+3\dot{H}\right)\right.\\ \,\left.+\,36c{\left(3{H}^{2}+\dot{H}\right)}^{2}+{f}_{6}\right),\end{array}\end{eqnarray}$
Equations (46), (47) and (48) represent the expression for energy conditions of the statistical f(T, B) gravity model. The behavior of the energy conditions is clearly shown in figure 14 with redshift. In the present universe, it is observed that NEC and DEC are well satisfied throughout cosmic evolution. However, the SEC is violated for our model. The violation of the SEC is due to the accelerated expansion of the Universe.
Figure 14. Null, Dominant and Strong energy conditions of f(T, B) = aT3 + cB2 gravity models versus redshift $\left(z\right)$.

7. Conclusions

Teleparallel theories of gravity and their adaptations have received a lot of attention recently, in order to answer numerous cosmological problems. These theories are located in a torsionally supported globally flat manifold. GR has an equivalent teleparallel representation (TEGR) based on torsion (and tetrads) rather than curvature, as is well known (and metrics). Numerous modified teleparallel hypotheses have been put forth from this angle. The first, known as f(T) gravity, is a natural generalization of the TEGR action that is achieved by altering the torsion scalar T in the action. This method is comparable to the metric counterpart of f(R) gravity. The cosmological behavior of the cosmos has been remarkably well described by these two ideas. A modified teleparallel theory of gravity known as the f(T, B) theory, which, within certain bounds, can recover either f(T) or f(R) gravity, was developed with the purpose of unifying both f(R) and f(T) gravity and examining how these theories are related. In this manuscript, we have presented a cosmological analysis for a teleparallel f(T, B) theory of gravity for flat FRW space-time.
The main results for these cosmological models are:

We have considered a well-known f(T, B) gravity model, namely f(T, B) = aTb + cBd, to examine other cosmological parameters.

By parametrizing the Hubble parameter and estimating the best fit values of the model parameters b0, b1, and b2 imposed from Supernovae type Ia, Cosmic Microwave Background, Baryon Acoustic Oscillation, and Hubble data using the MCMC method, we propose a method to determine the precise solutions to the field equations. We then observed that the model appears to be in good agreement with the observations.

The Hubble parameter H and the redshift parameter z are related in the current study’s regression equation. We have used the best regression relation feasible to construct an appropriate functional relation between the Hubble parameter H and redshift z. We have calculated R2, where ${R}^{2}=\tfrac{{SSR}}{{SST}}$, to determine the best fit.

We employed supernovae from the Pantheon compilation, which consists of 1048 spectroscopically verified Type Ia Supernovae spread in the redshift range 0.01 < z < 2.26 [53], in our investigation for various dark energy models.

We calculated the angular diameter distance DA(z) from BAO datasets by clustering perpendicular to the line of sight. Additionally, by seeing clustering along the line of sight, we were able to determine the Universe’s expansion rate, H(z).

Two model parameters, b1 and b2, are present in the considered form of H(z), and they have been restricted for further investigation using certain observational data. For the purpose of validating our methodology, we have used a number of external datasets, including a recently compiled observational Hubble dataset with 57 data points, the Pantheon compilation of SNeIa data with 1048 data points, and a Baryonic Acoustic Oscillation dataset with six data points, to determine the best fit values for these model parameters.

We have also obtained three specific models for f(T, B) gravity for different values of d and b. In Model I, b = d = 2; in Model II, b = 2, d = 3, and in Model III, b = 3, d = 2. In all three of these models, we have seen that the cosmos has an ΛCDM era for z = 0 and −1, but the early universe for z > 0 has both a barotropic and a quintessence field dark energy era. Thus, all three models of f(T, B) gravity are in good agreement with the ΛCDM era. Using information from the H(z), Pantheon samples, and BAO datasets, the model-independent constraints on the modified f(T, B) gravity have been investigated in the current study. This type of model-independent constraint on modified gravity from current data has been recently investigated [4650] and is already mentioned in the introduction.

Also of note is that in the current epoch, all models under study anticipate a violation of the strong energy condition. It is a well-known fact that applying the energy conditions to modified theories of gravity is an open investigation that ultimately leads to a conflict between theory and observation.

In recent years, f(T, B) gravity has been used to study a number of cosmic aspects. [39] demonstrates how f(T, B) gravity can mimic the de Sitter universe, power law, and CDM models using cosmological reconstruction methods. [64] used this method to present the bouncing solution in this extended teleparallel gravity and to explore the singularity and little rip cosmology, and [65] carried out a stability analysis on the cosmological models that were modelled in conflict with the observational data for the accelerating universe. Four cosmological models that have shown promise in meeting late-time cosmic acceleration measurements that can generate quintessence behavior and experience transition along the phantom-divide line are presented in [66] as the observational side of this class of models.

In 2023, Kadam et al [67] investigated a different type of cosmological model in f(T,B) gravity, where the Witzenb $\ddot{o}$ ck connection has been used instead of the usual levi-Civita connection by using f(T, B) = αT + βBn. A number of writers have looked into two different f(T,B) gravity model types so far: the first using the Levi-Civita connection and the second using the Witzenb $\ddot{o}$ ck connection. Kadam et al [67] and other authors obtained the cosmological models in f(T,B) gravity by using f(T, B) = αT + βBn. But in the present paper, we have investigated the cosmological models in f(T,B) gravity by using f(T, B) = aTb + cBd. It follows that our model is superior and more general than the prior one.

All of the aforementioned findings give us reason to assume that an extension in the form of f(T, B) gravity could result in intriguing situations where we can think about how torsion and boundary terms combine with astrophysical evidence to shed insight on the study of the late-time accelerating universe.

Acknowledgments

The authors (S H Shekh and A Pradhan) are appreciative of the help and resources given by the IUCAA in Pune, India. Additionally, the Science Committee of the Ministry of Education and Science of the Republic of Kazakhstan provided funding for this study (Grant No. AP09058240). The authors are thankful to the anonymous reviewers and editor for their constructive comments, which helped to improve the paper’s quality in its present form.

Declaration of competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
1
Riess A G 1998 Observational evidence from supernovae for an accelerating universe and a cosmological constant Astron. J. 116 1009

DOI

2
Perlmutter S 1999 Measurements of Ω and Λ from 42 high-redshift supernovae Astrophys. J. 517 565

DOI

3
Koussour M Myrzakulov N Shekh S H Bennai M 2022 Quintessence Universe and cosmic acceleration in f(Q, T) gravity Int. J. Mod. Phys. D. 31 2250115

DOI

4
Tsujikawa S 2013 Quintessence: a review Class. Quant. Grav. 30 214003

DOI

5
Liu J C H 2022 A quintessence dynamical dark energy model from ratio gravity Eur. Phys. J. C 82 165

DOI

6
Sakti M F A R Sulaksono A 2021 Dark energy stars with phantom field Phys. Rev. D 103 084042

DOI

7
Karakasis T Papantonopoulos E Vlachos C 2022 f(R) gravity wormholes sourced by a phantom scalar field Phys. Rev. D 105 024006

DOI

8
Tripathy S K Mishra B 2020 Phantom cosmology in an extended theory of gravity Chin. J. Phys. 63 448 458

DOI

9
Ribas M O Devecchi F P Kremer G M 2005 Fermions as sources of accelerated regimes in cosmology Phys. Rev. D 72 123502

DOI

10
Ribas M O Devecchi F P Kremer G M 2008 Cosmological model with non-minimally coupled fermionic field EPL 81 19001

DOI

11
Myrzakulov N Bekov S Myrzakulova S Myrzakulov R 2019 Cosmological model of (T) gravity with fermion fields via Noether symmetry J. Phys. Conf. Ser. 1391 012165

DOI

12
Myrzakulov N Tursumbayeva G Myrzakulova S 2021 The (2+ 1) dimensional metric (R) gravity non-minimally coupled with fermion field J. Phys. Conf. Ser. 2090 012065

DOI

13
Xiao K 2020 Tachyon field in loop cosmology Phys. Lett. B 811 135859

DOI

14
Tsyba P Razina O Suikimbayeva N 2021 Analysis cosmological tachyon and fermion model and observation data constraints Int. J. Mod. Phys. D 30 2150114

DOI

15
Barreiro T Sen A A 2004 Generalized Chaplygin gas in a modified gravity approach Phys. Rev. D 70 124013

DOI

16
Elmardi M Abebe A Tekola A 2016 Int. J. Geom. Meth. Mod. Phys. 13 1650120

DOI

17
Nojiri S Odintsov S D 2003 Modified gravity with negative and positive powers of the curvature: unification of the inflation and of the cosmic acceleration Phys. Rev. D 68 123512

DOI

18
Nojiri S Odintsov S D 2011 Unified cosmic history in modified gravity: from F(R) theory to Lorentz non-invariant models Phys. Rep. 505 59 144

DOI

19
De Felice A Tsujikawa S 2010 f(R) theories Living Rev. Rel. 13 3

DOI

20
Arapoglu A S Deliduman C Eksi K Y 2011 Constraints on perturbative f(R) gravity via neutron stars J. Cosmol. Astropart. Phys. JCAP07(2011)020

DOI

21
Santos J Alcaniz J S Carvalho F C Pires N 2008 Latest supernovae constraints on f(R) cosmologies Phys. Lett. B 669 14

DOI

22
Ruggiero M L Iorio L 2007 Solar System planetary orbital motions and f(R) theories of gravity J. Cosmol. Astropart. Phys. JCAP01(2007)010

DOI

23
Chiba T Smith T L Erickcek A L 2007 Solar System constraints to general f(R) gravity Phys. Rev. D 75 124014

DOI

24
Godani N 2022 Thin-shell wormholes in non-linear f(R) gravity with variable scalar curvature New Astron. 96 101835

DOI

25
Godani N Samanta G C 2022 Cosmologically viable f(R) model and wormhole solutions Int. J. Geom. Methods Mod. Phys. 19 2250224

DOI

26
Godani N Samanta G C 2021 Charged traversable wormholes in f(R) gravity Int. J. Geom. Methods Mod. Phys. 18 2150098

DOI

27
Godani N Samanta G C 2020 Traversable wormholes in f(R) gravity with constant and variable redshift functions New Astron. 80 101399

DOI

28
Einstein A 1928 Riemann-geometrie mit aufrechterhaltung des begriffes des fernparallelismus. preussische akademie der wissenschaften Phys.-math. Klasse Sitzungsberichte 217 221

29
Hayashi K Shirafuji T 1979 New general relativity Phys. Rev. D 19 3524 3553

DOI

30
Bengochea G Ferraro R 2009 Dark torsion as the cosmic speed-up Phys. Rev. D 79 124019

DOI

31
Cai Y F Capozziello S De Laurentis M Saridakis E N 2016 f(T) teleparallel gravity and cosmology Rept. Prog. Phys. 79 106901

DOI

32
Harko T Lobo F S N Otalora G Saridakis E N 2014 $f(T,{ \mathcal T })$ gravity and cosmology J. Cosmol. Astropart. Phys. JCAP12(2014)021

DOI

33
Kofinas G Saridakis E N 2014 Teleparallel equivalent of Gauss–Bonnet gravity and its modifications Phys. Rev. D 90 084044

DOI

34
Bahamonde S Boehmer C G Wright M 2015 Modified teleparallel theories of gravity Phys. Rev. D 92 104042

DOI

35
Najera S Aguilar A Rave-Franco G A Escamilla-Rivera C Sussman R A 2022 Inhomogeneous solutions in f(T, B) gravity Int. J. Geom. Methods Mod. Phys. 19 2240003

DOI

36
Paliathanasis A Leon G 2021 Cosmological evolution in f(TB) gravity Eur. Phys. J. Plus 136 1092

DOI

37
Moreira A R P Silva J E G Lima F C E Almeida C A S 2021 Thick brane in f(TB) gravity Phys. Rev. D 103 064046

DOI

38
Bahamonde S Golovnev A Guzman M-J Levi Said J Pfeifer C 2022 Black Holes in f(TB) gravity: exact and perturbed solutions J. Cosmol. Astropart. Phys. JCAP01(2022)037

DOI

39
Bahamonde S Zubair M Abbas G 2018 Thermodynamics and cosmological reconstruction in f (T, B) gravity Phys.Dark Univ. 19 78 90

DOI

40
Shekh S H Chirde V R Sahoo P K 2020 Energy conditions of the f(T, B) gravity dark energy model with the validity of thermodynamics Commun.Theor.Phys. 72 085402

DOI

41
Shekh S H Myrzakulov N Pradhan A Mussatayera A 2023 Observational constraints on F(T, TG) gravity with Hubble’s parametrization Symmetry 15 321

DOI

42
Dunkley J Bucher M Ferreira P G Moodley K Skordis C 2005 Fast and reliable Markov chain Monte Carlo technique for cosmological parameter estimation Month. Not. Roy. Astron. Soc. 356 925 936

DOI

43
Hogg D W Foreman-Mackey D 2018 Data analysis recipes: Using Markov Chain Monte Carlo* Astrophys. J. Suppl. Ser. 236 11

DOI

44
Scolnic D M 2018 The complete Light-curve sample of spectroscopically confirmed Type Ia supernovae from Pan-STARRS1 and cosmological constraints from the combined Pantheon sample ApJ 859 101

DOI

45
Peng P The Pantheon Sample analysis of cosmological constraints under new models arXiv:2303.10095

46
Cai Y-F Khurshudyan M Saridakis E N 2020 Model-independent reconstruction of f(T) gravity from Gaussian processes ApJ 882 62

DOI

47
dos Santos F B M Gonzalez J E Silva R 2022 Observational constraints on f(T) gravity from model-independent data Eur. Phys. J. C 82 823

DOI

48
Koussour M 2023 A model-independent method with phantom divide line crossing in Weyl-type f(Q,T) gravity Chin. J. Phys. 83 454

DOI

49
Mu Y Li E-K Xu L 2023 Data-driven and model-independent reconstruction of modified gravity J. Cosmol. Astropart. Phys. 06 022

DOI

50
Taddei L Martinelli M Amendola L 2016 Model-independent constraints on modified gravity from current data and from the Euclid and SKA future surveys J. Cosmol. Astropart. Phys. JCAP12(2016)032

DOI

51
Bahamonde S Capozziello S 2017 Noether symmetry approach in f(TB) teleparallel cosmology Eur. Phys. J. C 77 107

DOI

52
Sarkar A Chattopadhyay S Gudekli E 2021 A statistical analysis of observational Hubble parameter data to discuss the cosmology of holographic Chaplygin gas Symmetry 13 701

DOI

53
Scolnic D M 2018 The complete light-curve sample of spectroscopically confirmed SNe Ia from Pan- STARRS1 and Cosmological Constraints from the combined pantheon sample Astrophys. J. 859 101

DOI

54
Komatsu E (WMAP Collaboration) 2018 Five-year Wilkinson microwave anisotropy probe (WMAP) observations: cosmological interpretation J. Suppl. 180 330

DOI

55
Anderson L (BOSS Collaboration) 2014 The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: baryon acoustic oscillations in the Data Releases 10 and 11 Galaxy samples Mon. Not. Roy. Astron. Soc. 441 24

DOI

56
Alam U Bag S Sahni V 2017 Constraining the cosmology of the phantom brane using distance measures Phys. Rev. 95 023524 D

DOI

57
Beutler F 2011 The 6dF galaxy survey: baryon acoustic oscillations and the local Hubble constant Mon. Not. Roy. Astron. Soc. 416 3017

DOI

58
Gaztanaga E Cabre A Hui L 2009 Clustering of luminous red galaxies IV: baryon acoustic peak in the line-of-sight direction and a direct measurement of H(z) Mon. Not. Roy. Astron. Soc. 399 1663

DOI

59
Mehrabi A Rezaei M 2021 Cosmographic parameters in model-independent approaches Astrophys. J. 923 274 284

DOI

60
Rapetti D Allen S W Amin M A Blandford R D 2007 A kinematical approach to dark energy studies Mon. Not. Roy. Astron. Soc. 375 1510

DOI

61
Blandford R D 2004 ASP Conf. Ser. 339 27

62
Sahni V Saini T D Starobinsky A A Alam U 2003 Statefinder—A new geometrical diagnostic of dark energy J. Exp. Theor. Phys. Lett. 77 201

DOI

63
Alam U Sahni V Saini T D Starobinsky A A 2003 Exploring the expanding Universe and dark energy using the statefinder diagnostic Mon. Notices Royal Astron. Soc. 344 1057

DOI

64
Caruana M Farrugia G Said J L 2020 Cosmological bouncing solutions in f(T, B) gravity Eur. Phys. J. C 80 640

DOI

65
Franco G A R Rivera C E Said J L 2020 Stability analysis for cosmological models in f(T, B) gravity Eur. Phys. J. C 80 677

DOI

66
Rivera C E Said J L 2020 Cosmological viable models in f(T, B) theory as solutions to the H0 tension Class. Quantum Grav. 37 165002

DOI

67
Kadam S A Said J L Mishra B 2023 Accelerating models in f(T, B) gravitational theory Int. J. Geom. Methods Mod. Phys. 20 2350083

DOI

Outlines

/