Solving Tolman-Oppenheimer-Volkoff equations in f(T) gravity: a novel approach

The torsion models have stood out among the proposals for an alternative description of gravity. The simplest of them, the Teleparallel theory, is equivalent to General Relativity and there are many studies that seek to study its extension to more general functions of the torsion $T$. The purpose of {our study } is to consider a family of $f(T)$ models and apply their corresponding Tolman-Oppenheimer-Volkof equations to compact objects such as neutron stars. Consequently, through a numerical analysis, calculate, among other things, the maximum mass allowed by the model for a neutron star, which would also allow us to evaluate which models are in accordance with observations. In the present paper, the first in the series, we show explicitly the set of equations that must be solved, and how to solve it, in order to model compact stars in $f(T)$ gravity without the need to adopt any particular form for the metric functions or consider any perturbative approach, as has been done in some works in the literature.