Acceleration in Friedmann cosmology with torsion

TL;DR

This paper explores a Friedmann cosmological model within Einstein-Cartan theory incorporating torsion linked to matter spin, demonstrating compatibility with observations and explaining the universe's accelerated expansion.

Contribution

It introduces a specific torsion function proportional to the Hubble parameter and matter density, fitting recent cosmological data and explaining acceleration.

Findings

Model reproduces accelerated expansion with transition redshift ~0.65

Torsion function form is compatible with observational data

Deceleration parameter indicates a transition from deceleration to acceleration

Abstract

A Friedmann like cosmological model in Einstein-Cartan framework is studied when the torsion function is assumed to be proportional to a single $\phi(t)$ function coming just from the spin vector contribution of ordinary matter. By analysing four different types of torsion function written in terms of one, two and three free parameters, we found that a model with $\phi(t)=- \alpha H(t) \big({\rho_{m}(t)}/{\rho_{0c}}\big)^n$ is totally compatible with recent cosmological data, where $\alpha$ and $n$ are free parameters to be constrained from observations, $\rho_m$ is the matter energy density and $\rho_{0c}$ the critical density. The recent accelerated phase of expansion of the universe is correctly reproduced by the contribution coming from torsion function, with a deceleration parameter indicating a transition redshift of about $0.65$.

Tables

Table 1. Table 1: Mean values of the free parameters and 95% c.l. constraints for Case I.

Parameter	95% limits
${𝛀}_{𝒎}$	${0.52}_{-0.20}^{+0.21}$
$𝜶$	${0.38}_{-0.11}^{+0.12}$
${𝑯}_{\mathrm{𝟎}}$	${69.6}_{-3.1}^{+3.1}$

Table 2. Table 2: Mean values of the free parameters and 95% c.l. constraints for Case II.

Parameter	95% limits
${𝛀}_{𝒎}$	${0.116}_{-0.057}^{+0.059}$
$𝜶$	$-{0.160}_{-0.076}^{+0.084}$
${𝑯}_{\mathrm{𝟎}}$	$59.2\pm 1.3$

Table 3. Table 3: Mean values of the free parameters and 95% c.l. constraints for Case III.

Parameter	95% limits
${𝛀}_{𝒎}$	${0.31}_{-0.12}^{+0.11}$
$𝜶$	${0.14}_{-0.12}^{+0.14}$
$𝒏$	$-{0.47}_{-0.36}^{+0.26}$
${𝑯}_{\mathrm{𝟎}}$	${68.8}_{-3.1}^{+3.0}$

Table 4. Table 4: Mean values of the free parameters and 95% c.l. constraints for Case IV.

Parameter	95% limits
${\Omega }_m$	${0.62}_{-0.62}^{+0.39}$
$𝜶$	${0.30}_{-0.31}^{+0.51}$
$𝒎$	$-{3.6}_{-7.6}^{+6.9}$
$𝒏$	$-{1.5}_{-2.7}^{+2.4}$
${𝑯}_{\mathrm{𝟎}}$	$68.7\pm 2.2$

Table 5. Table 5: Common cosmological parameters ( Ω m subscript Ω 𝑚 \Omega_{m} and H 0 subscript 𝐻 0 H_{0} ) from Cases I, II, III, IV and from Λ Λ \Lambda CDM, O Λ Λ \Lambda CDM models at 95% c.l. and respective statistical parameters. χ ν 2 = χ m ​ i ​ n 2 / ( n − p ) subscript superscript 𝜒 2 𝜈 subscript superscript 𝜒 2 𝑚 𝑖 𝑛 𝑛 𝑝 \chi^{2}_{\nu}=\chi^{2}_{min}/(n-p) , where n 𝑛 n is number of data and p 𝑝 p is number of free parameters. AIC = χ m ​ i ​ n 2 + 2 ​ p AIC subscript superscript 𝜒 2 𝑚 𝑖 𝑛 2 𝑝 \mathrm{AIC}=\chi^{2}_{min}+2p , BIC = χ m ​ i ​ n 2 + p ​ ln ⁡ N BIC subscript superscript 𝜒 2 𝑚 𝑖 𝑛 𝑝 𝑁 \mathrm{BIC}=\chi^{2}_{min}+p\ln N .

	Case I	Case II	Case III	Case IV	O$\Lambda$CDM	$\Lambda$CDM
${𝛀}_{𝒎}$	${0.52}_{-0.20}^{+0.21}$	${0.116}_{-0.057}^{+0.059}$	${0.31}_{-0.12}^{+0.11}$	${0.62}_{-0.62}^{+0.39}$	${0.238}_{-0.057}^{+0.056}$	${0.278}_{-0.028}^{+0.031}$
${𝑯}_{\mathrm{𝟎}}$	$69.6\pm 3.1$	$59.2\pm 1.3$	${68.8}_{-3.1}^{+3.0}$	$68.7\pm 2.2$	${68.0}_{-2.6}^{+2.7}$	${69.57}_{-0.98}^{+0.99}$
${𝝌}_{𝒎𝒊𝒏}^{\mathrm{𝟐}}$	1068.23	1097.90	1053.28	1054.20	1055.15	1057.74
${𝝌}_{𝝂}^{\mathrm{𝟐}}$	0.9747	1.0017	0.9619	0.9627	0.9627	0.9642
$𝒑$	3	3	4	4	3	2
AIC	1074.23	1103.90	1061.28	1062.20	1061.15	1061.74
BIC	1089.24	1118.91	1081.29	1082.21	1076.16	1071.74

Equations

$$R_{\mu \nu }-\frac{1}{2}R{g}_{\mu \nu }=\kappa T_{\mu \nu },$$

$$K_{\mu \nu }^{\alpha }=S_{\mu \nu }^{\alpha }+S_{\mu \nu }^{\alpha }+S_{\nu \mu }^{\alpha }$$

$$S_{\alpha \mu \nu }=-\frac{1}{4}\kappa (2s_{\mu \nu \alpha }+g_{\nu \alpha }{s}_{\mu }-g_{\alpha \mu }{s}_{\nu }),$$

$$S_{\alpha \mu \nu }=\varphi (h_{\alpha \mu }{u}_{\nu }-h_{\alpha \nu }{u}_{\mu })S_{\alpha }=-3\varphi u_{\alpha },$$

$$\bigg{(}\frac{\dot{a}}{a}\bigg{)}^{2}=\frac{8\pi G}{3}\rho-\frac{k}{a^{2}}-4\phi^{2}-4\bigg{(}\frac{\dot{a}}{a}\bigg{)}\phi\,,$$

$$\frac{\ddot{a}}{a}=-\frac{4\pi G}{3}(\rho+3p)-2\dot{\phi}-2\bigg{(}\frac{\dot{a}}{a}\bigg{)}\phi\,,$$

$$\dot{\rho }+3(1+\omega )H\rho +2(1+3\omega )\varphi \rho =0,$$

$$\rho=\rho_{0}\bigg{(}\frac{a_{0}}{a}\bigg{)}^{3(1+\omega)}\textrm{e}^{-2(1+3\omega)\int_{t_{i}}^{t}\phi(t)dt}$$

$$q=\frac{4\pi G}{3H^2}(\rho +3p)+2\frac{\dot{\varphi }}{H^2}+2\frac{\varphi }{H}.$$

$${\chi }_H^2=\sum _{i=1}^{51}\frac{{\left[H_{obs,i}-H(z_i,\mathbf{s})\right]}^2}{{\sigma }_{H_i,obs}^2},$$

$${\chi }_{SN}^2={\left[{\mathbf{m}}_{obs}-\mathbf{m}(z,\mathbf{s})\right]}^T{\mathbf{C}}^{-1}\left[{\mathbf{m}}_{obs}-\mathbf{m}(z,\mathbf{s})\right]$$

$$H^2=\frac{8\pi G}{3}{\rho }_m-\frac{k}{a^2}-4\alpha H_{0}H-4{\alpha }^2{H}_0^2,$$

$$\rho_{m}=\rho_{0m}\bigg{(}\frac{a}{a_{0}}\bigg{)}^{-3+2\alpha}$$

$$H^2=\frac{8\pi G}{3}{\rho }_m-\frac{k}{a^2}+4\alpha H^2-4{\alpha }^2{H}^2,$$

$$\frac{H}{H_0}=\sqrt{\frac{{\Omega }_m(1+z{)}^{3-2\alpha }+{\Omega }_k(1+z{)}^2}{1-4\alpha +4{\alpha }^2}},$$

$$\rho_{m}(a)=\rho_{0c}\frac{3^{1/n}}{\big{(}2\alpha+3C_{1}(a/a_{0})^{3n}\big{)}^{1/n}}\,,$$

$$H^{2}=\frac{8\pi G}{3}\rho_{m}-\frac{k}{a^{2}}+4\alpha H^{2}\bigg{(}\frac{\rho_{m}}{\rho_{0c}}\bigg{)}^{n}-4\alpha^{2}H^{2}\bigg{(}\frac{\rho_{m}}{\rho_{0c}}\bigg{)}^{2n}\,,$$

$$\frac{H}{H_0}=\frac{(3-2\alpha {\Omega }_m^n)+2\alpha {\Omega }_m^n(1+z{)}^{3n}}{(3-2\alpha {\Omega }_m^n)-4\alpha {\Omega }_m^n(1+z{)}^{3n}}\sqrt{\frac{3^{1/n}{\Omega }_m}{{\left[2\alpha {\Omega }_m^n+\frac{(3-2\alpha {\Omega }_m^n)}{(1+z{)}^{3n}}\right]}^{1/n}}+{\Omega }_k(1+z{)}^2},$$

$$H^{2}=\frac{8\pi G}{3}\rho_{m}+4\alpha H_{0}^{m+1}H^{-m+1}\bigg{(}\frac{\rho_{m}}{\rho_{0c}}\bigg{)}^{n}-4\alpha^{2}H_{0}^{m+2}H^{-2m}\bigg{(}\frac{\rho_{m}}{\rho_{0c}}\bigg{)}^{2n}\,,$$

$${\Omega }_m={\left(1+\frac{2{\varphi }_0}{H_0}\right)}^2=(1-2\alpha {\Omega }_m^n{)}^2$$

Full text

I Introduction

A more complete understanding of general relativity with the presence of matter can be obtained when one consider that the intrinsic angular momentum of fermionic particles (spin) promotes torsion effects in space-time. This can be achieved with the presence of asymmetric affine connection in the construction of a manifold, introducing the torsion of spacetime and therefore allowing emerge of new geometric degrees of freedom in the system. Thus, matter becomes responsible for being a source of torsion, enriching studies in cosmological scenarios, with more general prescriptions. An example is based on well-established studies of the Einstein-Cartan-Kibble-Sciama (ECKS) gravitational theory. This theory allows to describe in a more complete way the invariance of local gauge in relation to the group of Poincarè Sciama:1964wt ; Hehl:1971qi ; Hehl:1976kj ; Shapiro:2001rz , being very useful in studies of condensate of particles with half-integer spin and averaged as a spin fluid Hehl:1974cn ; Sabbata90 ; Sabbata94 besides scenarios with an effective ultraviolet cutoff in quantum field theory for fermions Poplawski:2009su . Even though there is no observational evidence to ponder the existence of torsion in spacetime, some suggestions for experimental tests involving spacetime studies with non-zero torsion for gravity can be found in Kostelecky:2007kx ; March:2011ry ; Hehl:2013qga . One of the major problems in finding this evidence is associated with the fact that effects of torsion become considerable mainly at high density and energies. There are other cosmological scenarios that have interesting consequences generated by torsion corrections. In Gonzalez-Espinoza:2019ajd the torsion effects generated by scalar fields contributes to explain inflation. In addition, non-minimal couplings with torsion effects have been studied to understand gravitational waves Valdivia:2017sat ; Izaurieta:2019dix .

However, Friedmann-Robertson-Walker (FRW) cosmological scenarios also can be addressed in presence of torsion. In particular, the very tiny value of the cosmological constant or dark energy needed to accelerate the universe could be mimicked due to contribution of the torsion. Moreover, the high symmetry of FRW spacetime preserves the symmetry associated to Ricci curvature tensor, which implies that the corresponding Einstein tensor and energy-momentum tensor also preserve a symmetric form. Such construction is very well motivated and discussed in Kranas2019 , which we recommend for further details. The whole effect of torsion due to spin of matter may be associated to a single scalar function, depending only on time. Such approach was also adopted in Olmo:2011uz ; Capozziello:2008kb ; Jimenez:2015fva ; Poplawski:2010kb . A dynamical system approach with weak torsion field was done recently in barrow2019 . A review on Friedmann cosmological models in Einstein-Cartan framework can be found in Medina2019 . The kinematics of cosmological spacetimes with nonzero torsion in the context of classical Einstein-Cartan gravity is given by Pasmat2017 and the first derivation of FRW equations with torsion was presented in tsamPRD .

The present paper aims to study torsion effects in FRW background for late time expansion of the universe, particularly the possibility to explain the recent accelerated phase of expansion as a consequence of torsion effects. It is assumed four different types of torsion function, parameterized by one, two and three free parameters. Constraints with observational data allows to fix the free parameters and compare the known parameters with the ones obtained from standard cosmological model, namely, the $\Lambda$CDM model parameters obtained from last Planck satellite observations Aghanim2018 .

The paper is organised as follows. Section 2 presents the main equations of Friedmann cosmology with torsion, based on Kranas2019 . In Section 3, the constraints from observational data are obtained for four different torsion functions. Section 4 analyses the torsion function and deceleration parameter evolution for the best function obtained in previous section. Conclusion is left to Section 5.

II Friedmann cosmology with torsion

We follow the same notation from Kranas2019 . The standard Einstein equation of gravitation maintain its original form in terms of Ricci tensor, Ricci scalar and energy momentum tensor,

[TABLE]

with $\kappa={8\pi G}$, however in a space-time with torsion the affine connection is endowed with an antisymmetric part, namely $\Gamma^{\alpha}_{~{}~{}~{}\mu\nu}=\tilde{\Gamma}^{\alpha}_{~{}~{}~{}\mu\nu}+K^{\alpha}_{~{}~{}~{}\mu\nu}$, where $\tilde{\Gamma}^{\alpha}_{~{}~{}~{}\mu\nu}$ defines the symmetric Christoffel symbols and $K^{\alpha}_{~{}~{}~{}\mu\nu}$ defines the contorsion tensor,

[TABLE]

written in terms of the torsion tensor $S^{\alpha}_{~{}~{}~{}\mu\nu}$, which is antisymmetric in its covariant indices, $S^{\alpha}_{~{}~{}~{}\mu\nu}=-S^{\alpha}_{~{}~{}~{}\nu\mu}$. In general case the energy momentum tensor is coupled to $S^{\alpha}_{~{}~{}~{}\mu\nu}$ by means of the Cartan field equations,

[TABLE]

where $s_{\alpha\mu\nu}$ and $s_{\alpha}=s^{\mu}_{~{}~{}\alpha\mu}$ are the tensor and vector spin of matter, respectively. Physically, torsion provide a link between the spacetime geometry and the intrinsic angular momentum of the matter Pasmat2017 . With the presence of torsion terms into Eq. (1), it is known as the Einstein-Cartan equation of gravitation.

In a homogeneous and isotropic Friedmann background, the torsion tensor and the associated vector are Kranas2019

[TABLE]

where $\phi=\phi(t)$ is an unique time dependent function representing torsion contribution due to homogeneity of space, $h_{\mu\nu}$ is a projection tensor, symmetric and orthogonal to the 4-vector velocity $u_{\mu}$.

In terms of the torsion field $\phi(t)$, the Friedmann equations are Kranas2019 :

[TABLE]

[TABLE]

where $k$ is the curvature parameter, $\rho$ and $p$ are the energy density and pressure of matter. Together the Friedmann equations, for a barotropic matter satisfying $p=\omega\rho$, the continuity equation reads Kranas2019

[TABLE]

with $H=\dot{a}/a$, whose solution with $\omega$ constant is

[TABLE]

where the subscript [math] denotes present values and $t_{i}$ some initial time. We see that torsion alters the energy density evolution of standard matter through the exponential term.

In order to better understand the influence of torsion function into recent accelerated phase of expansion of the universe, we look for the deceleration parameter, which can be written as

[TABLE]

For a constant and negative $\phi$ for instance, torsion tends to accelerate the expansion. For a flat and empty space ($k=\rho=0$), Eq. (5) leads to $\phi(t)=-H(t)/2$, which suggest a $H$ dependence to the torsion function. However, since the physical source of torsion is the spin of matter, a torsion function dependent on the matter density is also a much more realistic choice. The above discussions and dimensional arguments will guide us in the next section in order to build some torsion functions and use them to compare with observational constraints, constraining the free parameters.

III Constraints from observational data

In order to study the possibility of dark matter and dark energy being driven by torsion effects in cosmological evolution, let us analyse the constraints imposed by observational data in four different models of torsion fields together ordinary matter contribution. Cases I and II below are just phenomenological assumptions for the torsion function, the first a constant function and the second evolving with $H(t)$. Cases III and IV are more realistic once they are explicitly dependent on the matter energy density, the real sources of spin in the universe.

The data used here were 51 $H(z)$ data from Magaña et al. Magana2018 and 1048 SNe Ia data from Pantheon compilation Scolnic2018 .

In all analyses here, we have written a $\chi^{2}$ function for parameters, with the likelihood given by ${\mathcal{L}}\propto e^{-\chi^{2}/2}$. The $\chi^{2}$ function for $H(z)$ data is given by

[TABLE]

where $\mathbf{s}$ is the parameter vector. For Pantheon, instead, we included systematic errors, thus we had to deal with the full covariance matrix. In this case, the $\chi^{2}$ is given by

[TABLE]

where $\bm{C}$, $\bm{m}_{obs}$ and $\bm{m}$ are covariance matrix, observed apparent magnitude vector and model apparent magnitude, respectively. We have assumed flat priors for all parameters and have sampled the posteriors with the so called Affine Invariant Monte Carlo Markov Chain (MCMC) Ensemble Sampler by GoodWeare , which was implemented in Python language with the emcee software by ForemanMackey13 . In order to plot all the constraints on each model, we have used the freely available software getdist111getdist is part of the great MCMC sampler, COSMOMC cosmomc ., in its Python version.

III.1 Case I: $\phi(t)=\phi_{0}=-\alpha H_{0}$

For this simplest case of a constant torsion field, as already discussed by Kranas2019 , with $\alpha$ a dimensionless constant to be determined222The presence of $H_{0}$ warrants the correct dimension for the torsion term., we write the Friedmann equation (5) as:

[TABLE]

where $\rho_{m}$ is the matter density parameter obtained as a solution of (7) with $\omega=0$, namely \rho_{m}=\rho_{0m}\big{(}{a_{0}}/{a}\big{)}^{3}\textrm{e}^{2\alpha H_{0}t}, where $\rho_{0m}$ represents the present day matter energy density. Analytic solution of (12) exists just for spatially flat $(k=0)$ background, however a numeric treatment can be done in general case and the parameters $\alpha$, $\Omega_{m}$ and $H_{0}$ can be constrained with observational data333Here $\Omega_{m}=\frac{\rho_{0m}}{\rho_{0c}}$ as usual, and $\rho_{0c}=\frac{3H_{0}^{2}}{8\pi G}$ is the critical density..

Figure 1 shows the $1\sigma$ (68.3% c.l.) and $2\sigma$ (95.4% c.l.) contours for $\Omega_{m}$, $\alpha$ and $H_{0}$ parameters obtained with $H(z)$ and SNe Ia observational data. Table 1 presents the mean values of the parameters with $95\%$ c.l. constraints. For this model we see that $\Omega_{m}$ is just marginally compatible at $2\sigma$ with the last results for $\Lambda$CDM model from the Planck collaboration on the cosmological parameters444From Aghanim2018 , $\Omega_{m}=0.315\pm 0.007$ for matter density and $H_{0}=(67.4\pm 0.5)$ km/s/Mpc. Aghanim2018 , while $H_{0}$ is compatible at $1\sigma$.

III.2 Case II: $\phi(t)=-\alpha H(t)$

For this case, the analytic solution of (7) is

[TABLE]

and the Friedmann equation (5) turns:

[TABLE]

In terms of the density parameters, Eq. (14) is:

[TABLE]

where555$\Omega_{k}\equiv-\frac{k}{a_{0}^{2}H_{0}^{2}}$ is the curvature parameter $\Omega_{k}=1-\Omega_{m}-4\alpha+4\alpha^{2}$ and the redshift is introduced by $(1+z)=a_{0}/a$.

Figure 2 shows the constraints for $\Omega_{m}$, $\alpha$ and $H_{0}$ at $1\sigma$ and $2\sigma$ contours for $H(z)$ and SNe Ia observational data. Table 2 presents the mean values of the parameters with $95\%$ c.l. constraints. For this model we see that both $\Omega_{m}$ and $H_{0}$ are very small, not compatible with the $\Lambda$CDM model even at $2\sigma$.

III.3 Case III: \phi(t)=-\alpha H(t)\bigg{(}\frac{\rho_{m}(t)}{\rho_{0c}}\bigg{)}^{n}

This general case is much more interesting, since that the torsion function is proportional to matter density $\rho_{m}$ and it is expected that torsion contribution comes from spin of ordinary matter. Also, for this case it is easy to verify that a solution of Eq. (7) for the energy density with $\omega=0$ is

[TABLE]

where $C_{1}$ is a integration constant. In order to have $\rho_{m}(a_{0})=\rho_{m0}$, we set $C_{1}=-\frac{2}{3}\alpha+\Omega_{m}^{-n}$.

The Friedmann equation (5) is:

[TABLE]

In terms of the density parameters Eq. (17) is:

[TABLE]

where $\Omega_{k}=(1-2\alpha\Omega_{m}^{n})^{2}-\Omega_{m}$.

Figure 3 shows the constraints for $\Omega_{m}$, $n$, $\alpha$ and $H_{0}$ at $1\sigma$ and $2\sigma$ contours for $H(z)$ and SNe Ia observational data. Table 3 presents the mean values of the parameters with $95\%$ c.l.. We see that both $\Omega_{m}$ and $H_{0}$ are in very good agreement to the $\Lambda$CDM model at $1\sigma$, with a small positive $\alpha$ value and a negative $n$ value. With such parameters the model is totally compatible with the recent cosmic acceleration, with the dark energy component being represented by torsion function $\phi(t)$.

III.4 Case IV: \phi(t)=-\alpha H_{0}\bigg{(}\frac{H_{0}}{H(t)}\bigg{)}^{m}\bigg{(}\frac{\rho_{m}(t)}{\rho_{0c}}\bigg{)}^{n}

For this general case there is no analytic solution for the energy density $\rho_{m}$ and one must resort to numerical methods. Due to this model having many free parameters, we choose to work with the spatially flat case ($k=0$), which is favoured by inflation and recent CMB observations.

The Friedmann equation (5) for a spatially flat Universe is:

[TABLE]

Due to spatial flatness, there is a relation among the free parameters of this model, which can be obtained from the Friedmann equation, as:

[TABLE]

where the last equality holds only for Case IV. As we can see from this relation, due to nonlinearity, there is not a unique solution for parameters $\alpha$ or $n$. So, for this analysis, we choose to work with the free parameter $\varphi_{0}\equiv\frac{\phi_{0}}{H_{0}}$, we find the constraints over $\varphi_{0}$ and then we use it to obtain $\Omega_{m}$ by using (20). We choose a flat prior $\varphi_{0}\in[-1,0]$, which yields a flat prior $\Omega_{m}\in[0,1]$.

Figure 4 shows the constraints for $\Omega_{m}$, $n$, $m$, $\alpha$ and $H_{0}$ at $1\sigma$ and $2\sigma$ contours for $H(z)$ and SNe Ia observational data. Table 4 presents the mean values of the parameters with $95\%$ c.l.. We see that $\Omega_{m}$ and $H_{0}$ are compatible with the $\Lambda$CDM model values. Also, one can see that $\Omega_{m}$ is poorly constrained by this analysis, due to having too many free parameters.

Now we can compare the mean values for $\Omega_{m}$ and $H_{0}$ cosmological parameters obtained for the above four cases with the ones from $\Lambda$CDM model, together with the statistical parameters as $\chi^{2}_{min}$, AIC and BIC Akaike74 ; Schwarz78 ; JesusEtAl16 . We choose to compare the torsion models with both concordance models, namely spatially flat $\Lambda$CDM model, as well as with the O$\Lambda$CDM model, which allows for nonzero spatial curvature, once that only Case IV is spatially flat. Table 5 present the parameters at 95% c.l. for the combined analysis $H(z)$+Pantheon.

As we can see on this Table, AIC favours O$\Lambda$CDM and Case III over the concordance flat $\Lambda$CDM. But it is a known fact JesusEtAl16 that AIC does not penalize enough the excess of parameters, and BIC is in general most accepted, being an approximation of the Bayesian Evidence. According to BIC, the flat $\Lambda$CDM model and O$\Lambda$CDM are favoured by this analysis. Case III is the torsion model with the lowest BIC, although significantly above $\Lambda$CDM.

IV Torsion evolution and transition redshift

In order to better reproduce the standard model constraints and obtain a cosmic acceleration in agreement with the latest observational data (see footnote 3), Case III above is the better one, with both $\Omega_{m}$ and $H_{0}$ compatible within $1\sigma$ c.l. and with lower $\chi^{2}$ and BIC parameters.

For this case it is interesting to analyse the evolution of torsion function and the transition redshift. From (16) and (17), we have obtained the mean $\phi(z)$ from the parameters MCMC chains, jointly with its variance. The evolution of the torsion function is shown in Figure 5, for the mean $\phi(z)$ (blue line) and for $1\sigma$ c.l. (orange and green lines). The behaviour of the torsion function on the past is strongly dependent on the values of parameters, specifically on the $n$ parameter. At present, the behaviour is similar in all cases, showing an increase on the absolute value of torsion function just in recent times, which coincides with the late time acceleration phase of expansion of universe. In this sense, torsion function makes the role of a dark energy acting during the whole history of the universe. In the past the matter energy density dominates over the torsion contribution and today is the torsion function that dominates, driving the acceleration. A similar observation was made in a recent work within the context of a scalar-tensor theory with torsion Cid:2017wtf , where the authors showed that an effective torsional fluid plays an important role in recent acceleration phase of expansion of universe but becomes subdominant in the past, where a pressureless matter component dominates at redshift around 200.

The behaviour of the deceleration parameter is better to understand the recent evolution of the universe and is presented on Figure 6. For larger $z$ values the deceleration parameter seems to converge to $0.5$, a value characteristic of a matter dominated universe, as expected from standard model. As seen above, for $z\lesssim 1$ the torsion function start to increase and a transition to accelerated phase occurs, dominated by torsion term. The transition redshift $z_{t}$ occurs at about $z_{t}=0.65$, in good agreement to standard model.

V Conclusion

We have analysed a Friedmann like universe with the contribution of a torsion function in Einstein-Cartan cosmology. The torsion function is represented by $\phi(t)$, and for four different types of function written in terms of one, two and three free parameters we have studied the cosmic evolution and constrained the free parameters with observational data from $H(z)$ and SN Ia.

From the four different functions, Case III presents a very good agreement to observational data, in the sense that both matter energy density parameter and $H_{0}$ are completely compatible with the results for $\Lambda$CDM model obtained from last Planck mission observations. Statistical parameters as BIC and $\chi^{2}$ also favours Case III, as pointed out in Table 5, comparing the four models with standard $\Lambda$CDM model.

The effect of torsion function is to act as a dark energy fluid at late time, correctly explaining the present accelerated phase of expansion of the universe. A model-independent result concerning Einstein-Cartan gravity was obtained in ivanov2016 , showing that torsion can be responsible for the vacuum energy density or the cosmological constant, exactly as dark energy density in the universe.

Finally, the deceleration parameter obtained for the model furnish a desirable transition to accelerated phase at about $z_{t}=0.65$, coming from a matter dominated phase in the past, as occurs for standard model of cosmology.

Acknowledgements.

This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - Brasil (CAPES) - Finance Code 001. SHP would like to thank CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico, Brazilian research agency, for financial support, grants number 303583/2018-5 and 400924/2016-1. RCL would like to thank CAPES. JFJ is supported by Fundação de Amparo à Pesquisa do Estado de São Paulo - FAPESP (Process no. 2017/05859-0).