Astrophysical flows near \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f\,\,(T)$$\end{document}f(T) gravity black holes

In this paper, we study the accretion process for fluids flowing near a black hole in the context of f(T) teleparallel gravity. Specifically, by performing a dynamical analysis by a Hamiltonian system, we are able to find the sonic points. After that, we consider different isothermal test fluids in order to study the accretion process when they are falling onto the black hole. We find that these flows can be classified according to the equation of state and the black hole features. Results are compared in f(T) and f(R) gravity.


Introduction
One of the most important problems in modern cosmology is the dark energy issue, which is responsible for the accelerated expansion of the observed Universe. Over the last few decades, several studies have been focused on trying to tackle this problem. It is well known that this form of energy is acting as a repulsive gravitational force so that in General Relativity (GR) one needs to consider a further non-standard fluid with a negative pressure to justify this accelerated scenario. The simplest approach is to consider a cosmological constant in order to explain it. However, from quantum considerations, the necessary expected value of it must be extremely much larger than the observed value [1]. Another approach to the cosmic accelerated behavior comes from modified theories of gravities where, instead of searching for new material ingrea e-mail: ayyesha.kanwal@sns.nust.edu.pk b e-mail: azreg@baskent.edu.tr c e-mail: sebastian.beltran.14@ucl.ac.uk d e-mail: capozziello@na.infn.it e e-mail: mjamil@sns.nust.edu.pk dients, the philosophy is to address cosmic dynamics taking into account possible further degrees of freedom of the gravitational field. A very well-studied approach to modified gravity comes out from the "Teleparallel equivalent to General Relativity" (TEGR). This theory yields the same field equations as in General Relativity, so that TEGR is an alternative and equivalent theory. However, the geometrical interpretations of these theories are different. On the one hand, GR assumes a non-zero curvature and a vanishing torsion by choosing the symmetric Levi-Civita connection. On the other hand, TEGR considers an antisymmetric connection provided with a non-vanishing torsion and a zero curvature (Weitzenböck connection). In other words, one can say that GR uses the curvature to geometrize the space-time, meanwhile TEGR uses torsion to explain gravitational effects. In TEGR, we need to use tetrad fields as the dynamical variables in order to define the Weitzenböck connection (see [2][3][4][5][6][7][8][10][11][12][13], and also the review [14] for the basis in TEGR).
A natural generalization of TEGR is, instead of using the scalar torsion T , to consider an arbitrary and smooth function of the torsion f (T ) in the gravitational action [15][16][17][18]. This theory is the so-called " f (T ) gravity". The idea comes out naturally exactly as when GR is generalized to f (R) gravity [19][20][21][22]. An important problem related to f (T ) gravity is that it is no longer invariant under local Lorentz transformations so that different tetrads might give rise to different solutions. Therefore one needs to be very careful choosing the correct tetrad [23]. Although TEGR is equivalent to GR, it is important to mention that f (R) is no longer equivalent to f (T ) gravity [24]. One needs to consider a more general theory of gravity, the so-called " f (T, B) gravity" to obtain the teleparallel equivalent to f (R) gravity [25]. In addition, it is important to remark that f (T ) gravity contains only second order derivative terms; meanwhile f (R) gravity contains up to fourth order derivative terms in the metric formalism.
In the last few years, f (T ) gravity acquired a lot of interest in cosmology due to the possibility to explain by it the accelerated expansion of the cosmic Hubble fluid (see [26][27][28][29][30][31][32][33][34][35]). In addition, astrophysical studies related with compact objects as black holes has been considered among f (T ) gravity such as in [36][37][38][39][40]. However, it is worth noticing that this is not the only solution that can be achieved by the Noether symmetry approach. As shown in [41] for f (R) gravity, the symmetries select the form of the function and several Noether vectors can exist. In the specific case of f (T ) gravity, other solutions have been found as discussed in [42,43]. A very well-studied process, known as accretion, occurs when a fluid is situated in the vicinity of a black hole or a massive astrophysical object (see [44][45][46][47]). In this process, the compact object takes particles from the fluid and increases its mass. Accretion takes place regularly in the Universe, and it can be used to test gravitational theories using observational measurements [48][49][50]. The first study of accretion was performed using Newtonian gravity by Bondi [51]. He found transonic solutions for a gas accreting onto compact objects. Michel extended the later work considering GR for a Schwarzschild black hole [52]. An important work in this field has been pursued by Babichev et al., where they showed that the mass of the black hole decreases when a phantom fluid is in accretion onto it [53]. Later, Jamil and Qadir showed that primordial black holes decay earlier when the effect of accretion of phantom energy is considered [54]. In addition, Nayak and Jamil also found that primordial black holes accrete radiation, matter, and vacuum energy when they pass through radiation, matter, and vacuum dominated eras, respectively, with the result that they live longer during the radiation era [55]. After that, several works have been published on accretion onto compact objects (see [56][57][58][59][60]).
Recently, Ahmed et al. studied accretion for cyclic and heteroclinic flows near f (R) black holes [61]. In this paper, we will use a similar formalism in order to study the accretion process in a black hole in the context of f (T ) gravity. This paper is organized as follows: In Sect. 2, we briefly introduce the TEGR and f (T ) gravity. In Sect. 3, we discuss the metric representation of black holes in f (T ) gravity. Section 4 is devoted to finding the general equations for spherical accretion. In Sect. 5, we perform a dynamical system analysis using the Hamiltonian formalism and we study the system at the critical points (CPs). In Sect. 6, we obtain solutions for isothermal test fluids for different kind of fluids. In Sect. 7, we analyze the accretion process for a polytropic test fluid. Finally, in Sect. 8, we discuss our results and draw conclusions. Throughout the paper we will use the metric signature (−, +, +, +) and the geometric units G = c = 1.

Teleparallel equivalent of general relativity and f (T ) gravity
Let us briefly introduce TEGR and its generalization which is the so-called f (T ) gravity. We will adopt the notation used in [25]. In this theory, the dynamical variable is the tetrad field e μ a (or vierbein), where Latin and Greek index indicate tangents space and space-time index, respectively. The construction of this theory relies on the relationship between the tetrad field and the metric g μν in the following way: where g μν is the inverse of the metric, E μ a is the inverse tetrad, which satisfies the relation E μ a e a ν = δ μ ν , and η ab = (−1, 1, 1, 1) is the Minkowski metric. Therefore, at each point x μ of the manifold, the tetrad field forms an orthonormal basis for the tangent space.
As we discussed before, TEGR uses a specific connection (Weitzenböck connection) where the space-time is globally flat but is endowed with a non-zero torsion tensor. This connection is defined by Then we can construct the torsion tensor using the antisymmetric part of the Weitzenböck connection as follows: Using the torsion tensor, one can define the contorsion tensor In addition, it is useful to define where T μ = T λ λ μ is the contraction of the torsion tensor.
Using the above tensor, the torsion scalar T can be defined as The Riemann tensor can be expressed depending on the contorsion tensor as follows: Here ∇ μ represents the covariant metric derivative. Therefore, the Ricci scalar R and the torsion scalar T are related by where e = det(e a μ ). It is important to remark that B = 2 e ∂ μ (eT μ ) is a boundary term. Instead of using the Ricci scalar R as in GR, the TEGR Lagrangian density is described by the torsion scalar T Since B is a boundary term, from (9), one can see that the TEGR action will arise to the same field equations as the Einstein-Hilbert action, making these two theories equivalent.
One important and very well-studied generalization of TEGR is to consider an arbitrary smooth function of the scalar torsion to construct the action This theory is called " f (T ) gravity" and it has numerous and interesting applications, for example in cosmology (see [14] for a comprehensive review of those models). One important feature of this theory is that meanwhile TEGR is an equivalent theory to GR, f (T ) does not produce the same field equations as f (R) gravity (due to Eq. (9)) and therefore one needs to consider a generalization of (11) from f (T ) → f (T, B) to find the teleparallel equivalent to f (R) gravity as discussed in [25]. Starting from the action (11), the field equations read where the energy-momentum tensor is defined as follows: With these considerations in mind, let us start our discussion of black holes in f (T ) gravity.

Black hole in f (T ) gravity
The metric for a spherically symmetric black hole with mass M in f (T ) gravity is given by [40] where where and, X ≡ Here, all c 1 , c 2 , c 3 , c 4 , and c 5 are constants. The horizon is given by where we have introduced the new constants C 5 = c 5 /c 3 and X = c 1 /c 3 , which will turn out to be very useful in the study of the dynamical system. To ensure that r h > 0, C 5 , and X must have the same sign: C 5 / X > 0. Since A must be positive at spatial infinity, we must have X > 0 resulting in C 5 > 0. Upon performing the coordinate transformation we bring the metric (14) to the following form, where α(r ) = c 2 3 A(r ): This is precisely the general form of metric used in Ref. [61] where accretions of samples of f (R) black holes were investigated, among which we find the solution This will serve in Sect. 6.4 as a tool for comparing accretion onto the f (T ) black hole (14) with that onto the f (R) black hole (21). The metric (14) being equivalent to (20), all general equations expressed in terms of α, which were derived in Ref. [61], are thus applicable to our present investigation upon replacing α by c 2 3 A. However, because of their importance, we will outline their derivations below.

General equations for spherical accretion
Let n be the baryon number density in the fluid rest frame and u μ = dx μ /dτ be the four velocity of the fluid where τ is the proper time. We define the particle flux or current density by J μ = nu μ where n is the particle density. From the particle conservation law, we see that the divergence of current density is zero, i.e.
where ∇ μ is the covariant derivative. On the other hand, the energy-momentum tensor is explicitly given by where denotes the energy density and p is the pressure. We assume that the fluid is radially flowing in the equatorial plane (θ = π/2), therefore u θ = 0 and u φ = 0. For the sake of simplicity, we set u r = u. Using the normalization condition u μ u μ = −1 and (14), we obtain On the equatorial plane (θ = π/2), the continuity equation (22) yields or, upon integrating, where C 1 is a constant of integration. The thermodynamics of simple fluids is described by [62] where T is the temperature, s is the specific entropy, and is the specific enthalpy. On the other hand, a theorem of relativistic hydrodynamics [62] states that the scalar hu μ ξ μ is conserved along the trajectories of the fluid, where ξ μ is a Killing vector of space-time. Considering the timelike Killing vector ξ μ = (1, 0, 0, 0) of the metric (14), we obtain where C 2 is a constant of integration. It is easy to show that the specific entropy is conserved along the fluid lines: u μ ∇ μ s = 0. In fact, if we rewrite energy-momentum tensor μν (23) as nhu μ u ν +(nh−e)g μν [61], and then project the conservation formula of μν onto u μ , we obtain In the special case we are considering in this work where the fluid motion is radial and stationary (no dependence on time), and it conserves the spherical symmetry of the black hole, the latter equation reduces to ∂ r s = 0 everywhere, that is, s ≡ const. Thus, the motion of the fluid is isentropic and equations (27) reduce to Equations (26), (30), and (32) are the main equations that we will use to analyze the flow. Since s is constant, this reduces the canonical form of the equation of state (EOS) of a simple fluid e = e(n, s) to the barotropic form From the second equation (32), where the prime denotes differentiation with respect to n. Now, the first equation (32) which can be integrated by parts to derive We see that an EOS of the form p = G(n) is not independent of an EOS of the form = F(n). The relation between F and G can be derived upon integrating the differential equation The local three-dimensional speed of sound a is defined by Since the entropy s is constant, this reduces to a 2 = d p/d . Using (32), we derive Using (34), this reduces to Since the motion is radial in the plane θ = π/2, we have dθ = dφ = 0 and the metric (14) implies the decomposition where we have used u r = u = dr/dτ , u t = dt/dτ , u t = −Au t , and (24). This implies and (26) becomes These results will be used in the following Hamiltonian analysis.

Hamiltonian systems
We have derived two integrals of motion (C 1 , C 2 ) given in (26) and (30). Let H be the square of the left-hand side of (30): Using (41) the Hamiltonian (43) of the dynamical system reads as derived in Ref. [61] where f has been replaced by c 2 3 A. We can absorb the constant c 2 3 into a redefinition of the Hamiltonian, however, we will do that in a further step of our derivation.

Sonic points
With H given by (44), the dynamical system readṡ (here the dot denotes thet derivative). Evaluating the righthand sides we find Following the same approach as in Ref. [61], we arrive aṫ For the CP, the right-hand sides vanish if the conditions hold. Here A c ≡ A(r c ) and A c,r c ≡ A ,r | r =r c . They lead to If solutions to the system of equations (50) exist, we rewrite the constant C 2 1 in (42) as where we have used the second equation in (50). Using this in (42) we obtain the result If no solution to (50), we can keep (42) as it is or introduce any point (r 0 , v 0 ) from the phase portrait to obtain The above dynamical system allows one to perform the analysis of the fluids that we are considering.

Isothermal test fluids
Isothermal flow is often referred to the fluid flowing at a constant temperature. In this section we find the general solution of the isothermal EOS of the form p = k , that is, of the form p = k F(n) (33) with G(n) = k F(n) (37). Here k is the state parameter such that (0 < k ≤ 1) [63]. The differential equation (37) reads where we have chosen the constant of integration 1 so that (28) and (34) lead to the same expression for h: Using (53) or (54), we obtain and where all the constant factors have been absorbed into the redefinition of the timet and the Hamiltonian H. Now we will analyze the behavior of the fluid by taking different cases for the state parameter k. For instance, we have k = 1 (ultra-stiff fluid), k = 1/2 (ultra-relativistic fluid), k = 1/3 (radiation fluid) and k = 1/4 (sub-relativistic fluid). In the case of the metric (14), Eq. (51) reduces to and yields where r h is given by (18). It is easy to see that, in order to have r c > r h > 0, we must have C 5 / X > 0 and 1/3 < k < 1. This fixes the values of k that yield a critical flow with the presence of a CP given by (61) and v 2 c = k. In Ref. [61] we have shown that if the flow approaches the horizon with a vanishing three-dimensional speed, the pressure must diverge as if A(r ) = 0 has a single root.
6.1 Solution for ultra-stiff fluid (k k k = 1) The equation of state for the ultra-stiff fluids is p = k i.e. the value of state parameter is defined as k = 1. The Hamiltonian (59) reduces to   Fig. 1, which are the only existing solutions for k = 1, the speed v behaves asymptotically as v ∼ 1/r 2 . Using this and the fact that A ∼ r 2 in (42), we obtain n ∼ 1/r . 6.2 Solution for ultra-relativistic fluid (k k k = 1/2) Ultra-relativistic fluids are those fluids whose isotropic pressure is less than the energy density. In this case, the equation of state is defined as p = 2 , yielding k = 1/2. Using this expression in (61) reduces to r c = 5r h /2. Thus, we have two CPs given by The Hamiltonian (59) takes the simple form For some given value of H = H 0 , Eq. (65) can be solved for v 2 . Another way to represent the flow is to use contour plots, as shown in Fig. 2. For the global solutions depicted in the figure, the speed v has two different asymptotic behaviors.
Since H retains the same constant value and A ∼ r 2 , we have either (a) v → 0 as v ∼ cst/r or (b) v → 1 such that Recall that the analysis made in this paper considers the fluid elements as test particles not modifying the geometry of the f (T ) black hole. It is thus assumed that the accretion does not modify the mass of the black hole nor its other intrinsic properties. The flow, being non-geodesic, however, still obeys the simple rule that if r increases, v must be positive, and if r decreases, v must be negative. For instance, for v > 0, we see from Fig. 2 that the red plot has two branches. Consider the branch on the right of the vertical line r = r c . The flow along the segment of that branch along which v increases and r decreases is unphysical, for this is neither an accretion nor a flowout.
6.3 Solutions for radiation fluid (k k k = 1/3) and sub-relativistic fluid (k k k = 1/4) Radiation fluids (k = 1/3) are the fluids which absorb the radiation emitted by the black hole. It is the most interesting case in astrophysics and sub-relativistic fluids (k = 1/4) are those fluids whose energy density exceeds their isotropic pressure. The Hamiltonian (59) for these fluids takes the following expressions, respectively: As we concluded earlier in this section, there is no critical flow for these fluids and for all fluid cases where k ≤ 1/3; rather, simple fluid flow characterizes this class of fluids. Moreover, the fluid flow for this class of fluids is not global, in that it does not extend to spatial infinity except in the case k = 1/3 where the flow can be global and non-global. This conclusion can be derived from (59) as follows. If the flow is global, v behaves asymptotically as where α > 0, v 0 , and |v ∞ | ≤ 1 are constants. If we assume that the flow is global, that is, r may go to infinity, the Hamiltonian (59) behaves in the limit r → ∞ as Thus, in the case k < 1/3, the Hamiltonian diverges at spatial infinity. Since the Hamiltonian is a constant of motion, the assumption that r goes to infinity is not valid. For k = 1/3 global flow is possible, as we shall justify below; however, non-global flow is also realizable. Now, if k = 1/3 and 0 < |v ∞ | < 1, the Hamiltonian (69) has a finite limit as r → ∞, so global flow is possible. To achieve it, that is, to determine such global flow solutions, notice that the value of the Hamiltonian (66) in this case is For the global flow, we determine the particle density n as follows. Equation (70) with H given by the right-hand side of (66) yields Substituting this in (42) we obtain where all constants (X, c 3 , v ∞ ) have been grouped or absorbed into the new constant N 2 . Since asymptotically |v| → v ∞± , which is a non-zero constant, n ∼ r −3 .

Accretions in f (T ) and f (R) gravities
We draw a comparison between accretions in f (T ) and f (R) gravities. For that end we select from f (R) gravity black holes a similar solution (21) to the one considered here (14), that is, an anti-de Sitter-like f (R) black hole [61]. The following enumeration shows similarities and differences.
1. The accretion of an isothermal perfect fluid with k = 1 is characterized by the presence of global solutions, which are the only existing solutions with no CPs. The speed v and the particle density n behave asymptotically as v ∼ 0 and n ∼ 1/r for both gravities. This, however, is just a qualitative comparison. First of all notice that the black hole (21) of the f (R) gravity reduces to that of GR and the theory itself reduces to GR, f (R) = R + , if the f (R)-parameter β = 0. This is not the case with the black hole (14) of the f (T ) gravity which does not reduce to any of the known GR black holes no matter how the f (T )-parameters (X, C 5 ) are chosen.
A deeper investigation should focus on the evaluation of the rates of accretion and efficiencies of the outgoing spectra for different black holes and different gravity theories.
The efficiency of the conversion of gravitational (potential) energy into radiation is one of the open problems of radial accretion onto a black hole; this is if one assumes, as most workers concluded, that the infall velocity scales almost as the free fall velocity (the case of Fig. 1 or the case of the critical subsonic accretion followed by a supersonic accretion of Fig. 2). This efficiency problem becomes more involved if we consider the critical accretion of Fig. 2 along the branch where v vanishes as r → r h or accretions along the blue and magenta branches of the same figure. Here the three velocity has a deceleration phase from r c to r h and it does not scale as a free fall velocity. This is the main discovery in this work and in [61]. The deceleration of the fluid increases by far the conversion efficiency; moreover, the efficiency is roughly proportional to n 2 [64], which diverges by (54) as r → r h .
All that is out of the scope of this work and could be the aim and task of subsequent works. In a first step one may consider the simplest cases of the f (T ) = T [ f (R) = R or GR] gravity theory. We believe that, when all these tasks are performed (most likely numerically), the result that will be at hand will confirm the equivalence of these gravity theories.

Polytropic test fluids
The polytropic equation of state is where K and γ are constants. For ordinary matter, one generally works with the constraint γ > 1. Inserting (74) into the differential equation (37), it is easy to establish [61] the following expressions of the specific enthalpy: by integration, and the three-dimensional speed of sound from (39) where we have introduced the baryonic mass m. Since γ > 1, this implies a 2 < γ − 1 and, particularly, v 2 c < γ − 1. Using (53) or, preferably, the general expression (54), in (76) we arrive at where so we have a positive constant. If the CPs exist, Z takes the special form Inserting (77) into (44) we evaluate the Hamiltonian by where (c 3 m) 2 has been absorbed into a redefinition of (t, H).
The constraint X > 0, in (14), yields A ,r > 0 for all r , and this implies that the constant Z > 0 (recall that γ > 1).
Thus, the sum of the terms inside the square parentheses in (80) is positive, while the coefficient A/(1 − v 2 ) diverges as r → ∞ (0 ≤ 1 − v 2 < 1). So, the Hamiltonian too diverges as r approaches spatial infinity. Since the Hamiltonian has to remain constant on a solution curve, we conclude that there are no global solutions (solutions that extend to, or emanate from, spatial infinity). This conclusion is general and it extends to all anti-de Sitter-like solutions [61].
Since γ > 1, the solution curves do not cross the r axis at points where v = 0 and r = r h , for otherwise the Hamiltonian (80) would diverge there. The curves may cross the r axis at r = r h only. The horizon (18) being a single root to A(r ) = 0, if we assume v ∝ |r − r h | δ and δ > 0 near the horizon, it is easy to show that (81) being valid for δ > 0, we see that only physical solutions with 1 < γ < 2 may cross the r axis. For these values of γ , the pressure p = K n γ diverges at the horizon as Now, substituting into (76), we arrive at which along with Eq. (51) takes the form of the following expressions at the CPs: where we have used (14) to reduce the right-hand side of (51). The critical flow in the left plot of Fig. 5 is no difference of that of Fig. 2 (black plot). The only different feature is that the former flow is non-global while the latter flow is global. Similarly, the magenta and blue curves (corresponding to H > H c ) of the left and middle plots of Fig. 5 have branches which are subsonic for the whole process of accretion-flowout as is the case of the curves of Fig. 2 corresponding to H > H c . Another similarity emerges upon comparing the solutions with no CPs corresponding to Z < Z 0 (right plot of Fig. 5) with those of Fig. 3 where no CPs occur too.
A common conclusion we can draw upon comparing the solutions of this section with those of the previous one is that low pressure fluids (k and K small) do not develop critical flows (no CPs) and high pressure fluids develop critical flows but they may maintain purely subsonic, even non-relativistic, flows.
Barotropic fluids with γ > 2, if there are any, may have CPs but no critical flow and their accretion velocity never vanishes as depicted in Fig. 6. The accretion make take place along two different paths starting from rightmost point of the lower branch of Fig. 6. For large values of the Hamiltonian (this would be the case if Z is large, n, or K ), the accretion along one of these two paths is almost non-relativistic for r > r c , then the velocity jumps to supersonic and relativistic values as r approaches r h . For lower values of the Hamiltonian, the accretion takes place near the CP and the polytropic fluid never reaches the horizon.
As the title of this section indicates, the analysis made in this section and in the previous ones concern accretion of test fluids neglecting all back-reaction effects. This rules out any homoclinic flow and motion along closed paths, as those shown in Fig. 6, where v conserves the same sign but r increases and decreases.

Conclusions
In this paper, we discussed in detail the accretion process of a spherically symmetric black hole in the context of f (T ) gravity. In order to select the form of f (T ) model, we adopted the Noether symmetry approach, following [40]. In particular, we discussed spherically symmetric solutions coming from f (T ) = T m models (and, in general, analytic f (T ) models) that give rise to metrics of the form (20) and related gravitational potentials of the form (21); see [40] for details. We have analyzed the motion of isothermal relativistic and ultra-relativistic fluids by means of a Hamiltonian dynamical system capable of representing hydrodynamics around the black hole. The thermodynamical properties of the fluids have been discussed according to the suitable EOS. Furthermore, conserved quantities and CPs have been selected for any fluid. Roughly, the accretion mechanism can be classified as subsonic and supersonic according to the features of the black hole and the EOS. In particular, the three-dimensional velocity flow strictly depends on the EOS, the radius, and the CPs on the phase space. Finally, the results have been compared to the analog results in f (R) gravity putting in evidence similarities and differences.
Clearly, the accretion process of the fluids flowing the black holes strictly depends on the conserved quantities (Noether's symmetries) and the structure of CPs, as shown above. If conserved quantities are not identified, it could become extremely difficult to define the phase space structure of the dynamical problem and consequently the features of CPs. In conclusion, identifying the Noether symmetries allows one to fix the model (i.e. the form of f (T )), to derive the metric and the gravitational potentials, thanks to the reduction of the dynamical system, to define the form of the space phase. Models without these features are very difficult to handle.
From a very genuinely observational point of view, these studies could be related to the possible observable features of f (T ) black holes. In particular, the possibility to investigate f (T ) vs. f (R) black holes could be a powerful tool to discriminate between the curvature (GR) and torsional (TEGR) formulation of theories of gravity (see [14] for a detailed discussion). Specifically, the accretion process onto a black hole could be the feature capable of discriminating among competing models and, in general, between a curvature or a torsional formulation. A main role in this discussion is played by the stability conditions. For example, as discussed in [65] for the case of f (R) gravity, the stability conditions for any self gravitating object strictly depend on the theory. There it is demonstrated that the Jeans stability criterion is different if one considers f (R) instead of GR because effective mass, stability radius, Jeans wave length, and the other parameters characterizing any astrophysical object slightly change according to the underlying model. In general, if the accretor has a mass M and a radius R, the gravitational energy release is Clearly the accretion yield increases with the compactness M/R, that is, given a mass M, the yields depend on the accretor radius. Considering alternative theories of gravity, the above relation can be written as where the features of the given model can be summarized using the effective gravitational coupling G eff . This means that the effective potential (related, for example, to the g 00 component of the metric), determines the accretion process. For example, the potential (21) indicates that the extra terms with respect to the Newtonian one contribute to any accretion process by modifying the accretion yield. As discussed in Sect. 6D, differences and similarities between the f (T ) and f (R) pictures can be put in evidence by a detailed study of the accretion process. In particular, the number of CPs, the state parameter k and other features, besides the effective potential, can discriminate among competing models. From a genuine observational point of view, luminous phenomena powered by black holes could contain features capable of discriminating among theories as soon as the parameters G, M, and R are combined into a gravitational potential. For example, the accretion luminosity, is a feature directly related to these phenomena. HereṀ is the mass variation with time. If one considers a gamma ray burst, we have L ∼ 10 52 erg/s withṀ ∼ 0.1M /s. As shown in [66], this huge amount of energy can be addressed in a strong field regime by curvature corrections. In other words, the role of G eff for the adopted underlying model is crucial. Furthermore, other characterizing parameters, besides G eff , can be identified to discriminate observationally concurring accretion models: e.g. the Salpeter timescale [67], the blackbody temperature T b for thermalization, the Eddington limit [68], and so on. These arguments will be the topic of a forthcoming paper.