Model-independent reconstruction of f (T) teleparallel cosmology

We propose a model-independent formalism to numerically solve the modified Friedmann equations in the framework of f (T) teleparallel cosmology. Our strategy is to expand the Hubble parameter around the redshift z = 0 up to a given order and to adopt cosmographic bounds as initial settings to determine the corresponding f (z) f (T (H(z))) function. In this perspective, we distinguish two cases: the first expansion is up to the jerk parameter, the second expansion is up to the snap parameter. We show that inside the observed redshift domain z <= 1, only the net strength of f (z) is modified passing from jerk to snap, whereas its functional behavior and shape turn out to be identical. As first step, we set the cosmographic parameters by means of the most recent observations. Afterwards, we calibrate our numerical solutions with the concordance Lambda CDM model. In both cases, there is a good agreement with the cosmological standard model around z <= 1, with severe discrepancies outer of this limit. We demonstrate that the effective dark energy term evolves following the test-function: f (z) = A + Bz(2)e(Cz). Bounds over the set {A, B, C} are also fixed by statistical considerations, comparing discrepancies between f (z) with data. The approach opens the possibility to get a wide class of test-functions able to frame the dynamics of f (T) without postulating any model a priori. We thus re-obtain the f (T) function through a back-scattering procedure once f (z) is known. We figure out the properties of our f (T) function at the level of background cosmology, to check the goodness of our numerical results. Finally, a comparison with previous cosmographic approaches is carried out giving results compatible with theoretical expectations.