Study on Anisotropic Strange Stars in f ( T , T ) Gravity

: In this work, we study the existence of strange stars in the background of f ( T , T ) gravity in the Einstein spacetime geometry, where T is the torsion tensor and T is the trace of the energy-momentum tensor. The equations of motion are derived for anisotropic pressure within the spherically symmetric strange star. We explore the physical features like energy conditions, mass-radius relations, modiﬁed Tolman–Oppenheimer–Volkoff (TOV) equations, principal of causality, adiabatic index, redshift and stability analysis of our model. These features are realistic and appealing to further investigation of properties of compact objects in f ( T , T ) gravity as well as their observational signatures.


Introduction
In relativistic astrophysics, compact stars, which are more compact (i.e., possess larger mass and smaller radius) than ordinary stars, have received much attention from researchers to study their ages, structures and evolutions. A neutron star is the final stage of a gravitationally collapsed star which, after exhausting all its thermo-nuclear fuel, gets stabilized by degenerate neutron pressure. A massive neutron star may again collapse into a black hole, but a lower mass neutron star may be converted into a quark star. The strange quark matter [1][2][3] consists of u, d, s quarks and some electrons, which may be more stable than the ordinary nuclear matter. Actually, a neutron star consists of neutrons [4] whereas a strange star consists of quarks or strange matters. In this context, it should be noted that the study of neutron stars (NS) has provided us important information about the equation of state (EOS) of neutron-rich nuclear matter [5,6]. Theoretically, NS may be converted to (strange) quark stars, which are made purely of deconfined u, d and s quark matter, i.e., strange quark matter (SQM) [7][8][9]. However, there may have some electron-type leptons due to charge neutrality and β-equilibrium.
In 1916, Schwarzschild et al. [10,11] first gave the interior stellar solution for isotropic fluid (having equal radial pressure (p r ) and transversal pressure (p t )). However, much later in 1972, Ruderman [12] first observed that the interior geometry of the nuclear matter of the steller system with density of order 10 15 gm/cc possesses anisotropic behavior (i.e., p r = p t ). After that Herrera et al. [13] have analyzed the local anisotropic nature for self-gravitating systems. It is interesting to note that in the experimental area, Hewish et al. [14] explained the observed compactness of many astrophysical objects (rotating neutron stars) like Her X − 1, 4U 1820 − 30, RX J 1856 − 37 and SAX J 1808.4 − 3658. Thus, theoretical modeling gradually became complementary to the observational results.
The study of anisotropic stars in the background of General Relativity (GR) is important in the current research. Hossein et al. [15] studied the stable structure of stellar objects for anisotropic systems. In the context of the present Standard cosmology, leading beyond the standard model of particle theory, such analysis inevitably involves additional elements in the standard model of gravity. Consequently, Kalam et al. [16] have studied the anisotropic neutron star model by considering quintessence dark energy. To study compact stellar models, the researchers have used well known Krori-Barua (KB) metric [17,18]. Paul et al. [19] have analyzed a new exact solution for compact stars. Several astrophysical natures of neutron stars as well as quark stars have been investigated [3,20]. Bhar [21] has proposed a stable anisotropic quintessence strange star model whereas Abbas et al. [22] have studied anisotropic strange stars for quintessence dark energy model. Using the MIT bag model, Rahaman et al. [23] have observed the existence of strange stars. By considering the MIT bag model, the stability of the strange stars with anisotropic fluid has been analyzed by Abrañil and Malheiro [24]. Furthermore, for the MIT bag model, Murad [25] studied anisotropic charged strange stars.
On the other hand, the study of anisotropic stars in the background of modified gravity theory is of great interest to many researchers [26][27][28][29][30][31]. The compact stars in f (R), f (G), f (T) and f (R, G) gravity theories have been studied by several authors [22,[32][33][34][35][36][37], where R, G and T are Ricci scalar, Gauss-Bonnet scalar and the torsion scalar, respectively. Strange stars in f (T) gravity with an MIT bag model have been studied in [38]. Compact stars, in the form of a neutron stars as well as quark stars, in the Rastall theory of gravity, have been studied by several researches [39][40][41][42][43][44], whereas strange stars with the MIT bag model in the Rastall gravity have been studied in Reference [45]. Charged anisotropic collapsing stars with heat flux in f (R) gravity has been investigated by Nazar and Abbas [46]. Studies on anisotropic stars in different types of modified gravity theories are available in literature [41,[47][48][49][50][51]. Saha and Debnath [52] have investigated anisotropic stars in f (T) gravity with modified Chaplygin gas. The anisotropic compact stars with a quintessence field and modified Chaplygin gas in the framework of f (T) gravity model has been studied by Saha and Debnath [53]. On the other hand, strange stars as well as gravastars in f (R, T ) gravity have been studied in [54][55][56][57][58][59][60][61], where T is the trace of the energy-momentum tensor.
We note from literature survey that almost all the works under the background of f (T, T ) gravity have been performed in different cosmological context [62][63][64][65][66][67][68][69][70][71][72][73]. However, there are a few astrophysical applications of f (T, T ) gravity as can be found in the works of Pace and Said [74,75]. Pace and Said have studied quark star incorporating the MIT bag model [74] and neutron star models using a perturbative approach [75] in f (T, T ) gravity theory whereas gravastars have been investigated by Ghosh et al. [73] in the same alternative theory of gravity.
Our motivation to study anisotropic strange stars in f (T, T ) gravity lies in the above background in the presence of Einstein spacetime within the spherically symmetric geometry. Thus we investigate in the present work strange quark star with MIT bag model [74] assuming anisotropy in the matter distribution and would like to note any appreciable effects due to anisotropy in addition to GR. The organization of the paper is as follows: in Section 2, the mathematical formulation of f (T, T ) gravity in Einstein spacetime is given. In Section 3, we investigate the basic stellar equations in the framework of f (T, T ) gravity. In Section 4, we find out the solutions to the Einstein field equations due to anisotropic fluid source. In Section 5, we explore the physical features like energy conditions, mass-radius relations, modified Tolman-Oppenheimer-Volkoff (TOV) equations, principal of causality, adiabatic index, redshift and stability analysis of our model whereas in Section 6 we perform a comparison with the standard results of GR, in order to show the effect of our model, and why it justifies our analysis and give a relevant discussion on the results. It is to note that for a comparative study we have considered, on the arbitrary ground only, the strange star candidate LMC X − 4 for plotting and exhibiting physical features which has specific total mass and radius as M = 1.29 M Sun [76] and R = 9.711 km [77] respectively. In Section 7, we deliver the conclusions of the work.

Basic Mathematical Formalism of f (T, T ) Gravity and the Einstein Field Equations
The modified theories of teleparallel gravity are those for which the scalar torsion of teleparallel action is substituted by an arbitrary function of this latter. As it is done in teleparallel, the modified versions of this theory are also described by the orthonormal tetrads and its' components are defined on the tangent space of each point of the manifold. The line element is written as with the following definitions Note that η ij = diag(1, −1, −1, −1) is the Minkowskian metric and the {e i µ } are the components of the tetrad which satisfy the following identity: In General Relativity, one use the following Levi-Civita connection which preserves the curvature whereas the torsion vanishes, such as In the teleparallel theory and its modified version, one keeps the scalar torsion by using Weizenbock's connection defined as: From this connection, one obtains the geometric objects. The first is the torsion as defined by from which we define the contorsion as where the expression • Γ λ µν designs the above defined connection. Then we can write The two previous geometric objects (the torsion and the contorsion) are used to define another tensor by The torsion scalar is usually constructed from torsion and contorsion as follows: In the modified versions of teleparallel gravity, one can use a general algebraic function of scalar torsion instead of the scalar torsion only as it is done in the initial theory. So, the modified action can be given by where T µν is the energy-momentum tensor of the strange quark matter (SQM) distribution and L m represents the Lagrangian for the matter distribution. We define T µν as follows and also define the trace of T µν as T = g µν T µν . As L m depends only on the metric tensor components g µν and not on their derivatives, so we find Varying the action with respect to the tetrad, one obtains the equations of motion [62] as with /∂T 2 and T σ a represents the stress-energy tensor for the anisotropic fluid distribution defined as where p r and p t represents the radial and tangential pressures of the SQM distribution, u µ is the four-velocity which satisfies the conditions u µ u µ = 1 and u µ ∇ ν u µ = 0. By using some transformations, we can establish the following relations: Hence from the combination of Equations (16) and (17), the field equation (14) can be written as (under the relativistic geometrized units where Newton's gravitational constant G and the speed of light c are fixed to be the unity, i.e., G = c = 1): where ( The covariant derivative of Equation (19) reads as In the current work, we specifically focus on the existence of anisotropic strange stars in the extended teleparallel gravity and for this purpose, we use an algebraic function following the prescription by Harko et al. [62] to comply with the observational data, such that f (T, T ) = T n T − 2Λ, where , n and Λ are arbitrary constants, specifically represents the matter-geometry coupling constant, n is a pure number (for simplicity assumed to be unity here) and Λ can be recognized as the Cosmological constant as has been treated in their above-mentioned work by Harko et al. [62]. For the above-mentioned functional form of f (T, T ) the Equation (19) reduces to where and

Explicit Stellar Equations in f (T, T ) Gravity under the Einsteinian Spacetime
To describe the interior spacetime of the spherically symmetric static stellar system, we take the metric as follows: where the metric potentials ν and λ are the functions of the radial coordinate r only.
In order to preserve the invariance under local Lorentz transformations of line element and to get the form (2), the choice of tetrad matrices is not unique. Following the References [78][79][80][81][82] one can get e = det e i µ = e (ν+λ)/2 r 2 sin (θ) and with Equations (3)-(10), the torsion scalar in terms of r becomes where the prime ( ) denotes the derivative with respect to the radial coordinate r. Now the nonzero components of the Einstein tensors can be provided as follows: where primes stand for derivatives with respect to the radial coordinate r only. Substituting Equation (15) into Equation (22), we find the explicit form of the Einstein field equations for the interior metric (25) as follows: Here ρ eff , p eff r and p eff t represent the effective energy density, radial pressure and tangential pressure for our system and given as In the present model we have introduced p r in a very simple bag model and p t by hand, treating only gravity in a serious way. However, the physical origin of p t could be either for electrically charged strange stars [25] or for solid quark stars with rigidity [83]. Astrophysical origin as well as consequence of p t are also available along with the latter one in the following references [12,13,18,[84][85][86][87][88][89][90][91][92][93][94][95][96][97].
We assume that the SQM distribution inside the strange stars is governed by the simple phenomenological MIT Bag model equation of state (EOS) [98]. In the Bag model, by introducing ad hoc Bag functions, all the corrections of the energy and pressure functions of SQM have been maintained. We also consider that the quarks are non-interacting and massless in a simplified Bag model. The quark pressure therefore can be defined as where p f is the individual pressure of the up (u), down (d) and strange (s) quark flavors and B is the vacuum energy density (also well known as the 'Bag' constant) which is a constant quantity within a numerical range. In the present paper, we consider the value of the Bag constant as B = 83 MeV/ f m 3 [23]. Now the individual quark pressure (p f ) can be defined as p f = 1 3 ρ f , where ρ f is the energy density of the individual quark flavor. Hence, the energy density of the de-confined quarks inside the Bag is given by Using Equations (36) and (37), we have the EOS for SQM given as It is observed that ignoring critical aspects of the quantum particle physics in the framework of GR several authors [59,60,[99][100][101][102][103][104][105][106] have successfully introduced this simplified form of the MIT Bag EOS to study stellar systems made of SQM.
To have non-singular monotonically decreasing matter density inside the spherically symmetric stellar system, following Mak and Harko [107], we assume a simplified form of ρ given as where ρ c and ρ 0 are two specific constants and denote the maximum and minimum values of ρ at the center and on the surface, respectively. Now following Moraes et al. [108] we consider p t is related to ρ by a relation given as where c 1 and c 2 are purely constants. We define the mass function of the spherically symmetric stellar system as [55] At this juncture, we consider the Schwarzschild metric to represent the exterior spacetime of our system given by where M is the total mass of the stellar system. Now, substituting Equation (41) into Equation (30) we find where one can note that the metric potential e λ is of Schwarzschild type. In this case, the conservation equation (24) in f (T, T ) gravity takes the following form The essential stellar structure equations required to describe static spherically symmetric sphere in f (R, T ) gravity theory are given by dm dr = 4πr 2 ρ eff , dp r dr = 1 −1 + 16π+2 1 − dρ dp r + 4 dp t dp − (ρ + p r )

Solution to the Einstein Field Equations in f (T, T ) Gravity for Stellar Modeling
Substituting Equation (38) in Equations (39) and (40), we have the following pressures Again, substituting Equations (38), (46) and (47) into Equations (33), (34) and (35), we get the following physical parameters: Using Equation (48), we integrate Equation (45) which yields the mass function as where d is the integrating constant.

Physical Features of Strange (Quark) Stars in f (T, T ) Gravity
In this section, we study some physical features of the compact star, in order to examine the physical validity and stability of the system in the f (T, T ) gravity.

Energy Conditions
According to relativistic classical field theories of gravitation, anisotropic fluid model must obey the energy conditions. There is often a linear relationship between the energy density and pressure of the matter by the name (i) Null Energy Condition (NEC), (ii) Weak Energy Condition (WEC), (iii) Dominant Energy Condition (DEC) and (iv) Strong Energy Condition (SEC) [95,109]. In the f (T, T ) gravity theory these are as follows: NEC : ρ e f f + p e f f r ≥ 0, ρ e f f + p e f f t ≥ 0, WEC : ρ e f f + p e f f r ≥ 0, ρ e f f ≥ 0, ρ e f f + p e f f t ≥ 0,

SEC
: ρ e f f + p e f f r ≥ 0, ρ e f f + p e f f t ≥ 0, ρ e f f + p e f f r + 2p e f f t ≥ 0, DEC :

Mass-Radius Relation
From Equation (51), we can rewrite the mass function for our stellar system in f (T, T ) gravity as follows: where where m GR represents the mass function of the SQM distribution in the Einstein frame whereas in the expression of m(r) the second term, i.e., m MG , represents the mass of the matter distribution and originates due to coupling between the matter and geometry in the f (T, T ) gravity. In order to plot the Mass (M/M ) versus Radius (R) curve for the strange stars due to different values of , we can rewrite (57), the mass function for our stellar system in f (T, T ) gravity, as follows: where It is already mentioned that the strange star candidate LMC X − 4 has specific total mass and radius as M = 1.29 M [76] and R = 9.711 km [77] respectively. Based on this data set it can be shown that the central density of the star is ∼ 1.961 × 10 15 gm/cm 3 whereas the surface density is ∼ 1.045 × 10 15 gm/cm 3 [35]. It is argued that according to quantum chromodynamics at sufficiently high energy densities, hadronic nuclear matter undergoes a deconfinement transition to a new phase of quarks and gluons [110]. Based on astrophysical observations and theoretical calculations in a model-independent way, Annala et al. [111] have inferred that the matter in the interior of maximally massive stable neutron stars with mass M ≈ 2 M exhibits characteristics of the deconfined phase, which can be interpreted as evidence for the presence of quark-matter cores.

Stability of the Stellar Model
In order to study the stability of the stellar model in a more formal manner in f (T, T ) gravity, the following three issues will be discussed: (i) Modified form of the TOV equation, (ii) Herrera cracking concept and (iii) adiabatic index in the following sub-subsections.

Modified TOV Equation in f (T, T ) Gravity Theory
To check the hydrostatic equilibrium of our proposed model, the TOV [112,113] condition is very useful which is followed by the Equation (44). There exist four different forces, such as the gravitational force (F g ), hydrostatic force (F h ), anisotropic force (F a ) and an extra force in connection to f (T, T ) gravity (F e ) provided in Equation (44) and can be rewritten as

Principle of Causality
Anisotropic fluid stellar model will be physically admissible if we pay special attention to making the boundaries of the radial and the transversal speeds of sound in (0, 1) within the matter distribution, i.e., they are less than the speed of light c (in relativistic geometrized units, the speed of light c becomes 1). This is known as the Causality condition which is based on the concept of the Cracking condition as provided by Herrera [114] which is generally expressible as 0 ≤ v 2 s = dp dρ ≤ 1. However, instead of physical p and ρ, we have here employed p e f f and ρ e f f which are involved with the coupling factor . This allows us to see the contribution of the considered f (T, T ) theory of gravity, in particular, the effect of the coupling parameter. Thus, the radial and transversal speeds of sound of our model are v 2 sr = dp eff r

Cracking Condition
There is an another way to check the stability of the system through the Herrera cracking condition [114,115]. Here, we check whether the square of radial speed of sound exceeds the square of transversal speed of sound or not for this anisotropic fluid star model. So, this criterion sets up Using Equation (66), we calculate the difference of two speeds and they are

Adiabatic Index
The adiabatic index is a crucial feature to investigate regarding the stable status of an anisotropic fluid star model. It has been proposed that the radial adiabatic index should be greater than 4 3 , which represents the rigidity of the EOS parameter of the model [116][117][118][119]. The radial adiabatic index is defined as  Figure 9. Plot of Γ r and Γ t versus r for the strange star candidate LMC X − 4. The red, purple, magenta, blue and black colors represent, respectively, cases = 0, = 0.5, = 1, = 1.5 and = 2.

Compactification Factor and Redshift
According to Biswas et al. [61], the compactness of the anisotropic fluid model is given by the compactification factor (71) and the twice of the compactification factor [? ] is followed by The surface redshift can be defined as

A Comparison Study Between f (T, T ) gravity and Standard GR
In this section, we perform a comparative study with the standard results of GR, in order to show the effect of our model in the framework of f (T, T ) gravity and why it justifies our analysis, and thereafter give a relevant discussion on the results. These aspects can be discussed as follows: (i) We have plotted in Figure 1 the behavior of the metric potentials e ν and e λ with respect to the fractional radial coordinate r/R respectively both of which exhibit normal features as are expected, i.e., they should exhibit regular behavior at the center as well as on the surface of the stellar system. On the other hand, in Figure 2 the variation of the effective energy density ρ eff and the effective pressures p eff r and p eff t are shown with respect to r/R respectively whereas the anisotropic feature is shown in Figure 3. We observe from all the plots of Figure 2 that a kind of monotony is evolving from the center to a minimum value and hence indicates the physically viable character of the f (T, T ) gravity model under consideration. One can note that the standard solutions to the Einstein field equations for stellar modeling can be recovered for = 0 and Λ = 0.
(ii) From Figure 4, we can conclude that the conditions NEC, WEC, DEC and SEC are satisfied for different values of coupling matter geometry parameter for our anisotropic fluid model absence of which GR features can be retrieved.
(iii) Since r → 0, m(r) → 0 in Figure 5 so one can notice that the mass function is regular at the origin. From this mass-radius relationship, we also note a striking feature which depicts that with increasing value of radius the mass of the physical system also increases and beyond R ∼ 12 km mass starts gradually decreasing. Thus, we observe a stable configuration of strange stars around 12 km, beyond which it suffers from instability. Moreover, we observe that in contrast to the GR case the curve allows for a massive quark star at the same radius.
(iv) Based on Equation (63), i.e., the modified TOV equation, we have drawn Figure 6 which indicates that our propounded anisotropic fluid model for different values of the coupling matter-geometry parameter is in hydrostatic equilibrium state under the combination of considered forces.
(v) From Equations (64) and (65), we observe that the Causality condition depends on the coupling matter-geometry parameter and the constant c 1 . The related Figure 7 exhibits that the Principle of Causality is preserved inside the anisotropic strange star in f (T, T ) as here 0 ≤ v sr ≤ 1 and 0 ≤ v 2 st ≤ 1 if 0 ≤ c 1 ≤ 1.
(vi) In a similar way, Equations (67) and (68) indicate that the Cracking condition depends on the coupling matter-geometry parameter and c 1 . If 0 ≤ c 1 ≤ 1 3 then we have −1 ≤ v 2 st − v 2 sr ≤ 0 and 0 ≤ |v 2 st − v 2 sr | ≤ 1 which is evident from Figure 8 and inconformity with the criteria imposed on the system [114,120]. It means that we have a potentially stable anisotropic star fluid model.
(vii) From Figure 9, we note that the variation of both Γ r and Γ t is a monotonically increasing function at all the interior points of our model with Γ r > 4 3 and hence provides confirmation on the stability of our model as predicted for the inside of an anisotropic, relativistic and dynamically stable stellar system [117,121,122].
(viii) Figure 10 (left panel) indicates that u(r) is lying between the range 1 4 and 1 2 . The 2u(r) follows the maximum allowed value 8 9 for our proposed model. On the other hand, we have plotted a figure for the surface redshift with the help of Equation (73) in Figure 10 (right panel) which shows that it is a decreasing function starting from the core with higher value and on the surface of the star acquires a finite non-zero value. In this context, it is to be noted that Barraco and Hamity [123] proved that for an isotropic star and in the absence of the cosmological constant Z s < 2. According to Böhmer and Harko [124] the surface redshift for an anisotropic star can reach a maximum higher value Z s ≤ 5, in the presence of the cosmological constant. Though, according to Ivanov [89], the maximum acceptable value of the surface redshift will be restricted up to 5.211. As our result for the surface redshift maintains the limit, i.e., z s ≤ 5 always, so the proposed model is quite realistic.
In a nutshell, in all the above figures, one can note that we get back the GR result for = 0 whereas all the other non-zero values of correspond to f (T, T ) gravity contribution, e.g., in Figure 6, the null effect due to modified gravity exhibits through the flat curve which coincides with the x-axis and hence shows the standard GR result. However, the effect of becomes gradually prominent as we go on increasing its value towards a higher magnitude.

Discussions and Conclusion
In this work, we present mathematical modeling of strange stars in the background of f (T, T ) gravity in Einstein-Maxwell spacetime in order to obtain viable physical properties. We develop the equations of motion using anisotropic property within the spherically-symmetric strange star and thereafter explores the physical features, like the energy conditions, mass-radius relations, modified TOV equations, principle of causality, adiabatic index, redshift and stability analysis of our model.
The salient features of the presented model can be put forward as follows: (1) For the matter-geometry coupling parameter = 0 the usual results as available in GR can be retrieved from the results due to modified gravity f (T, T ). The above statement obviously also follows from Figures 9 and 10, e.g., for = 0 the usual form of TOV equation as in GR can easily be retrieved because, in this situation, the force due to modified gravity vanishes, i.e., the effect of coupling also becomes ineffective. In order to verify model parameters, we have used the data for the strange star candidate LMC X − 4 and have drawn all the Figures 9 and 10 for different values of the matter-geometry coupling parameter .
(2) In the present work, we have considered tetrad formalism to correspond to the metric (25), however, it is known that in spherically-symmetric solutions the diagonal tetrad is not the consistent choice. In this connection it is of note that in literature this type of approach is not at all unavailable even in the f (T, T ) gravity theory [62,[73][74][75]. The tetrad formalism has been used by Harko et al. [62] in cosmology whereas Pace and Said [74,75] as well as Ghosh et al. [73] have employed the technique in astrophysics.
(3) In this work, as a first step, we considered a very specific ansatz, viz. f (T, T ) = ωT n T − 2Λ, to present a simple mathematical model following the prescription by Harko et al. [62] to comply with the observational data. However, after successful execution of this specification, there are ample scopes to opt for other specified functional forms of f (T, T ) gravity in the future projects following other ansatz, e.g., f (T, T ) = αT + γT 2 which has been also proposed by Harko et al. [62] or linear function in the forms (i) f = αT(r) + βT (r) + ϕ as proposed by Pace and Said [74] and (ii) f = T(r) + T (r) + h(T, T ) + O( 2 ) as proposed by Pace and Said [75]. In the above ansatz, which has been employed in the present work, the value of n is assumed to be unity, i.e., n = 1, however n = 1 case also can be considered for further study.
(4) We would also like to point out here that in the present investigation, following Rahaman et al. [23], we have considered the value for Bag constant as B = 83 MeV/ f m 3 . However, it seems that this value is quite arbitrary and needs a specific range with lower and upper values [125]. This means that further extension of the present work will provide a more realistic stellar model to validate on the physical ground.
(5) Theoretically under specific conditions some of the the up (u) and down (d) quarks get transformed into strange (s) quarks. It is argued that the cold strange matter is the true ground state of nuclear matter [1,3,[126][127][128] and hence the u and d quarks once converted into strange matter. As a result, the entire quark matter gets converted into strange matter and in turn, the neutron star totally gets converted into a strange quark star [129]. It is pointed out through some recent simulations [130,131] that the merger process of two strange stars is possible. Consequently, the gravitational wave signal detection of these stars may enhance the probability of the strange quark matter (SQM) hypothesis.
(6) Now the curious issue may be whether the f (T, T ) gravity passes the tests of binary compact star merger (BH-BH, BH-NS, NS-NS), especially clear the binary BH merging? The answer seems positive as far as the simulations based work by Bauswein et al. [130,131] on the merger process of two strange stars indicate. However, eventually, that depends on the puzzle of alternative gravity, i.e., the degeneracy in understanding both the equation of state and modified gravity (MG) [132]. On constraining recent gravitational wave observations of binary BH mergers, Motohashi and Minamitsuji [133] have classified GR solutions in MG based on their retrieval status from MG → GR. According to this classification scheme, f (T, T ) gravity where GR solutions are realized by the cancellation of the coupling constant, can pass the tests of binary merger (for further study following References [134,135] can be consulted).
So, based on the features as shown by Figures 9 and 10 and other aspects, we can overall conclude that our propounded model is quite realistic and appeals to further study of theoretically observable features as well as observational signatures of its predicted compact objects.