Quark Stars in $f(T, \mathcal{T})-$Gravity

We derive a working model for the Tolman-Oppenheimer-Volkoff equation for quark star systems within the modified $f(T, \mathcal{T})$-gravity class of models. We consider $f(T, \mathcal{T})$-gravity for a static spherically symmetric space-time. In this instance the metric is built from a more fundamental tetrad vierbein from which the metric tensor can be derived. We impose a linear $f(T)$ parameter parameter, namely taking $f=\alpha T(r) + \beta \mathcal{T}(r) + \varphi$ and investigate the behavior of a linear energy-momentum tensor trace, $\mathcal{T}$. We also outline the restrictions which modified $f(T, \mathcal{T})$-gravity imposes upon the coupling parameters. Finally we incorporate the MIT bag model in order to derive the mass-radius and mass-central density relations of the quark star within $f(T, \mathcal{T})$-gravity.


I. INTRODUCTION
In recent years it has been shown that the Universe is accelerating in its expansion [1,2]. In order to explain this one can introduce the concept of the cosmological constant [3,4]. Together with the inclusion of dark matter we get the ΛCDM model which explains a whole host of phenomena within the universe. [5][6][7]. Another approach to explaining this acceleration is to modify the gravitational theory itself with alternative theories of gravity an example of which is f (R)-gravity [8][9][10][11].
f (T )-gravity uses a "teleparallel" equivalent of GR (TEGR) [12] approach, in which instead of the torsion-less Levi-Civita connection, the Weitzenböck connection is used, with the dynamical objects being four linearly independent vierbeins [13,14]. The Weitzenböck connection is curvature-free and describes the torsion of a manifold. In the current case we consider a pure tetrad [15], meaning that the torsion tensor is formed by a multiple of the tetrad and its first derivative only. The Lagrangian density can be constructed from this torsion tensor under the assumption of invariance under general coordinate transformations, global Lorentz transformations, and the parity operation [11,12,14,15]. Also the Lagrangian density is second order in the torsion tensor [12,14]. Thus f (T )-gravity generalises the above TEGR formalism, making the gravitational lagrangian * mark.pace.10@um.edu.mt † jackson.said@um.edu.mt a function of T [10][11][12].
Our study involves deriving a working model for the TOV equation within a new modification of f (T ) class gravity, namely f (T, T )-gravity. There is no theoretical reason against couplings between the gravitational sector and the standard matter one [10]. f (T, T )-gravity takes inspiration from f (R, T )-gravity [10,16] where instead of having the Ricci scalar coupled with the trace of the energy momentum tensor T , one couples the torsion scalar T with the trace of the matter energy-momentum tensor T [10,11,16].
Recently a modification to this theory has been propose, that of allowing for a general functional dependence on the energy momentum trace scalar, T µ µ = T .
Our interest is in studying the behaviour of spherically symmetric compact objects in this theory with a specific linear function being considered, namely f (T, T ) = αT (r) + βT (r) + ϕ where α, β are arbitrary constants, and ϕ we take as the cosmological constant. We consider the linear modification since it is the natural first functional form to consider, and the right place to start to understand how the trace of the stress-energy tensor might effect f (T, T ) gravity. In particular, our focus is on quark stars in f (T, T ) gravity. Besides the possibility of the existence of these exotic stars, this is also a good place to study the behavior of modified theories of gravity in terms of constraints. Moreover, this also opens the door to considerations of stiff matter in early phase transitions [17].
The plan of this paper is as follows; In sec-arXiv:1701.04761v1 [gr-qc] 17 Jan 2017 tion 2 we outline the theoretical background of the model. In section 3 we consider the rotated tetrad and use this to derive the TOV equation in f (T, T )-gravity in section 4. Section 5 will then present the contrasting mass-radius relations derived using the MIT bag model, which we derive numerically. Finally in section 6 we discuss the results.
The non-vanishing torsion [18][19][20] is given by In TEGR one uses the teleparallel spin connection, which by construction gives vanishing curvature, thus all the information of the gravita-tional field is embedded in the torsion tensor, while the gravitational Lagrangian is the torsion scalar [20]. The contorsion tensor is then defined as while the superpotential of teleparallel gravity is defined by [18,19] S µν The torsion scalar [18,19] is then given as As in the analogous f (R, T ) theories [21], the gravitational lagrangian is generalized to f (T, T ) giving [22,23] where T = δ ν µ T µ ν and is the trace of the energymomentum tensor while L m is the matter Lagrangian density [22]. In this instance f is an arbitrary function of the torsion scalar T and the trace of the energy-momentum tensor T [22]. The variation of the action defined in Eq.(5) with respect to the tetrad leads to the field equations In our case we take the spin connection as ω i λν = 0 from the start [20,[24][25][26][27].
We take a spherically symmetric metric for our system which has a diagonal structure [28] ds 2 = −e A(r) dt 2 +e B(r) dr 2 +r 2 dθ 2 +r 2 sin 2 θdφ 2 , and consider the fluid inside the star to be that of a perfect fluid which yields a diagonal energymomentum tensor of e-m T ν λ = diag(−ρ(r), p(r), p(r), p(r)), where ρ(r) and p(r) are the energy density and pressure of the fluid respectively and the time dependence will be suppressed for brevity [28]. These also make up the matter functions which, along with the metric functions, A(r) and B(r), are also taken to be independent of time. Thus the system is taken to be in equilibrium [7,28]. The equation of conservation of energy is given by dp(r) dr = −(ρ(r) + p(r)) dA(r) dr .
As is in [15] we use the following rotated tetrad We take this form of vierbein because it gives us more degrees of freedom [29] and it allows us to obtain a static and spherically symmetric wormhole solution in our standard formulation of f (T, T )-gravity [29,30].
Moreover, this particular form of the tetrad is what's called a pure tetrad [20]. This means that the spin connection elements of this tetrad vanish and the ensuing field equations do not need to consider spin connection terms [20].
Inserting this vierbein into the field equa-tions, from Eq.(4) we get the resulting torsion scalar For the only non-vanishing non-diagonal element Together these equations govern the behaviour of the compact star.

IV. TOV EQUATIONS IN f (T, T )-GRAVITY
We take f = αT (r)+βT (r)+ϕ as our lagrangian function where T (r) = ρ(r)−3p(r) [23]. Consider Eq. (12), and solving for A (r) we find Substituting this into Eq.(11) and reducing yields We invoke the density equation for an inhomogeneous body We invoke the MIT bag model [31] since it represents the EoS of quark stars We thus get Substituting this into Eq. (14) and then invoking the conservation Eq.(9) we get a relation between pressure, p(r), and radius, r, in this form dp(r) dr = (ρ(r) + p(r)) α 8πp(r)r + M (r) The mass-radius relation is also derived [32,33], using our modified Schwarzschild solution found in Eq.(18) and we get the following result Taking α = 1, β = 0, and ϕ = 0 we recover the GR TOV equations.

V. NUMERICAL MODELLING AND TESTING
In order to obtain a graphical relations of the TOV equations, we numerically integrate our derived TOV equations of the MIT bag model to this f (T, T )-gravity model for quark stars.
We use the MIT bag model because it is the simplest equation of state for quark matter [31,34]. This is obtained because a quark star is a self-gravitating system consisting of deconfined u, d, and s quarks and electrons [35]. These deconfined quarks are the fundamental elements of the colour superconductor system [31]. In the comparison with the standard hadron matter, they lead to a softer equation of state, the MIT bag model which is given in Eq. (17).
The value of ω in Eq. (17) is dependent on the mass m s of the strange quark [31]. In the case of radiation, we have m s = 0 and the parameter is ω = 0 [31]. In the case of a more relativistic model having m s = 250 MeV, the parameter would be ω = 0.28 [34,36]. The parameter γ 0 lies within the intervals 58.8 < γ 0 < 91.2 which has units MeV/fm 3 [37].

A. Mass Profile Curve
In Fig.(1) we show the mass profile curve of a quark star by setting the values of α = 1, ϕ = 2.036 × 10 −35 (cosmological constant) [38] then varying the value of β.
We take three values of β in this case to contrast between the GR case where β = 0, the case where the function for T (r) = ρ(r) − 3p(r) [23] is included i.e. β = −1. Finally we include the case where this function is magnified by including β = −10 so as to see the behaviour at various levels.
As we decrease the value of β we allow for a smaller quark star structure. With the inclusion of the T (r) element some variations to arise. The quark star's maximum mass has however increased showing that having β at lower orders of magnitude allows for a much denser quark star structure.
To show these variations properly we plot the curve for β = −10 in Fig.(1). Here we note that for β = −10 a more massive quark star is allowed in such a gravity framework however smaller in size. In Fig.(2) we also plot the central density-radius graph where we again set the values of α = 1, ϕ = 2.036 × 10 −35 (cosmological constant) [38] then vary the value of β.
Again we contrast with the GR case when taking β = 0. When we decrease the value of β we may note that the central density figure of the quark star is more reluctant to drop however when reaching a certain radius it then decreases at a more rapid rate.
To further magnify this effect we again plot the results which are given by taking β = −10. In contrast to the GR case we see that the curve allows for a slightly denser quark star at a certain radius.

VI. CONCLUSION
In this study we study the TOV equation and its derivative behaviour for the rotated, pure, spherically symmetric tetrad. We then contrasted this result to the GR case. Our model has responded well when the MIT bag model is considered.
Our main goal throughout this work was to keep our terms as general as possible, with the possibility to revert back to the GR case whenever we needed to. This fact was very useful in checking our results throughout the derivation.
Numerical techniques were required to solve the TOV equation where reasonable boundary conditions were used. We apply an equation of state so that we may eliminate one of the four variables i.e. make one of the variables dependent on another variable.
For future work we hope to be able to apply a lagrangian which is not linear however thus far we have been unable to yield working TOV equations.

VII. ACKNOWLEDGEMENTS
The research work disclosed in this publication is funded by the ENDEAVOUR Scholarship Scheme (Malta). The scholarship may be part-financed by the European Union -European Social Fund (ESF) under Operational Programme II -Cohesion Policy 2014-2020, "Investing in human capital to create more opportunities and promote the well being of society".