Anisotropic fluid spheres in the framework of $f(R,\mathcal{T})$ gravity theory
S. K. Maurya Francisco Tello-Ortiz


APS/123-QED 



 sunil@unizwa.edu.om
  Department of Mathematics and Physical Science,
College of Arts and Science, University of Nizwa, Nizwa, Sultanate of Oman

 francisco.tello@ua.cl
Departamento de Física, 
Facultad de ciencias básicas, Universidad de Antofagasta, Casilla 170, 
Antofagasta, Chile

  





The main aim of this paper is to obtain analytic relativistic anisotropic spherical solutions in f(R,𝒯) scenario. To do so we use modified Durgapal-Fuloria metric potential and the isotropic condition is imposed in order to obtain the effective anisotropic factor Δ̃. Besides, a notable and viable election on f(R,𝒯) gravity formulation is taken. Specifically f(R,𝒯)=R+2χ𝒯, where R is the Ricci scalar, 𝒯 the trace of the energy-momentum tensor and χ a dimensionless parameter. This choice of f(R,𝒯) function modifies the matter sector only, including new ingredients to the physical parameters that characterize the model such as density, radial, and tangential pressure. Moreover, other important quantities are affected such as subliminal speeds of the pressure waves in both radial and transverse direction, observational parameters, for example, the surface redshift which is related with the total mass M and the radius r_s of the compact object. Also, a transcendent mechanism like equilibrium through generalized Tolman-Oppenheimer-Volkoff equation and stability of the system are upset. We analyze all the physical and mathematical general requirements of the configuration taking M=1.04 M_⊙ and varying χ from -0.1 to 0.1. It is shown by the graphical procedure that χ<0 yields to a more compact object in comparison when χ≥0 (where χ=0.0 corresponds to general relativity theory) and increases the value of the surface redshift. However, negative values of χ introduce in the system an attractive anisotropic force (inward) and the configuration is completely unstable (corroborated employing Abreu's criterion). Furthermore, the model in Einstein gravity theory presents cracking while for χ>0 the system is fully stable. The relationship between effective pressures and effective density  ρ̃ is discussed and obtained. This is achieved by establishing the corresponding equation of state.



                              
 Anisotropic fluid spheres in the framework of f(R,𝒯) gravity theory
    Francisco Tello-Ortiz
    December 30, 2023
====================================================================






§ INTRODUCTION

Put forward by Harko and his collaborators <cit.>, f(R,𝒯) gravity theory was designed to face the late-time acceleration of the Universe and the existence of dark matter. All these issues provided by recent observational data <cit.>. At present f(R,𝒯) theory is an active research field in the cosmological context <cit.>. Other interesting works available in the literature are for example the study by Sharif et.al about the non-static line element for collapsing of a spherical body having anisotropic fluid <cit.>, the static spherical wormhole solutions found in <cit.>. Moreover, perturbation techniques were used by Bhatti et.al in the study of spherical stars <cit.>. The effects on gravitational lensing due to f(R,𝒯) gravity were discussed by Houndjo in <cit.>. Furthermore, Baffou et al. <cit.> employed perturbation on de-Sitter space-time and power-law models in order to explore some cosmic viability bounds.

Even though the study of the Universe as a whole is a very intriguing current problem of dealing with, studies and investigations of structures within it such as neutron stars, white dwarfs and black holes among others, which constitute real laboratories to analyze in a fragmented way the most hidden secrets of the Universe, it is expected once these secrets are revealed, they provide us with the expected response that will finally put the pieces of this great puzzle in the corresponding place. In this direction, many authors have investigated the existence of collapsed structures within the framework of  f(R,𝒯) theory,
exploring how different models of f(R, 𝒯) affect on the principal properties of these kind of objects, besides contrasting the reported results with general relativity theory (GR hereinafter)<cit.>. Of course, is not an easy task solve f(R, 𝒯) field equations (as in GR case), for this reason one needs to prescribe additional information such as a suitable metric potential, an electric field (in the corresponding case), an adequate anisotropy factor or an equation of state (EoS from now on). Concerning the latter, the obtaining and not imposing of the EoS leads to a better understanding of how the matter confined in the stellar interior behaves. Moreover, it is possible to determine what type of material constitutes the distribution, for example, ordinary matter such as neutrons or strange matter such as quarks <cit.> (and references contained therein). Additionally, the EoS provides the relation between the macro observables parameters of the star such as mass and radius.
Following the same spirit, in this work we study the existence of compact objects, specifically neutron stars (or the possibility of quark stars also known as strange stars) in the framework of f(R,𝒯) gravity. The main goal is to build up the full geometrical description of the interior space-time taking as departure point Durgapal-Fuloria <cit.> metric potential ξ and introduce an anisotropic behavior of the matter distribution. The latter is achieved imposing the isotropic condition at the level of the field equations considering p_r≠ p_t once a viable form of the f(R,𝒯) function is chosen. In this opportunity we have selected the modified gravity model to be f(R,𝒯)=R+2χ𝒯 <cit.>, being R the Ricci scalar, 𝒯=g^μνT_μν the trace of the energy-momentum tensor (As was pointed out by Harko et.al <cit.> the dependence from 𝒯 may be induced by exotic imperfect fluids or quantum effects) and χ a coupling dimensionless constant. This coupling constant in some sense quantifies the effect on the matter sector modifications and in the geometrical one also (it is included in the Durgapal-Fuloria ansatz).   

The study of compact objects driven by anisotropic matter distributions must fulfill some general and basic requirements in order to be a physical and mathematical admissible model from the astrophysical point of view. Taking into account such conditions, with the support of graphic analysis we have studied and corroborated the fulfillment of each of them, which have been established throughout history from the pioneering work by Bowers and Liang <cit.>, in research developed by Herrera, Ponce de León, Cosenza, Di Prisco and Ivanov, to name a few <cit.>. Although these demands have been widely developed and applied to the study of collapsed structures within the framework of GR <cit.> (and references contained therein). The same in general must be satisfied in the scenario of modified gravity theories, since these theories under some limit must reproduce GR and its results, which are known to be very precise and have been verified several times. Of course, the chosen f(R,𝒯) model in this research reproduces the original Durpgapal-Fuloria solution developed in the context of GR in the limit χ=0.0 and α=1.
 
So, the obtained model has interesting properties that can be compared with the obtained ones in the GR frame. In considering -0.1≤χ≤0.1 and M=1.04M_⊙ (regarding also the GR limit ı.e χ=0.0) we have performed a full study and checked all the necessary and sufficient conditions in order to describe an acceptable compact object. So, in section <ref> we derive the complete set of equations corresponding to f(R,𝒯) theory, and the matter content is fixed ı.e the energy-momentum tensor T_μν and the Lagrangian matter ℒ_m. In Sec. <ref> the full model has specified ı.e the inner space-time geometry and the material content. In Sec. <ref> the internal geometry is smoothly joining with the exterior Schwarzschild geometry. In Secs. <ref> and <ref> we analyze the full main salient features of the obtained model. Studying the physical and mathematical behavior of the space-time geometry and all the thermodynamic variables such as the effective radial and tangential pressure, the effective energy density. In Sec. <ref> by means of energy conditions we check if the energy-momentum tensor is driven the material content if well behaved at every point within the star. In Sec. <ref> causality condition and its implications are discussed.
Furthermore, we discuss the influence of anisotropies fluids distributions and f(R,𝒯) model on the surface redshift. Sec. <ref> is devoted to the analysis of the Tolman-Oppenheimer-Volkoff equilibrium equation in the f(R,𝒯) context. In order to explore how the whole system is affected by the different forces. We provide in Sec. <ref> the corresponding equation of state, examining its properties and discussing its importance in the study of collapsed configurations such as neutron stars. In Sec. <ref> the stability of the system is analyzed using Abreu's criteria. Finally, in Sec. <ref> we reinforce all the good properties obtained in this work and comparisons between f(R,𝒯) gravity theory and GR are given. 




§ A FORMULATION FOR MODIFIED F(R, 𝒯) GRAVITY THEORY 

 Let us consider the integral action S for modified f(R,𝒯) gravity theory as

    S=1/16π∫ f(R,𝒯)√(-g) d^4x+∫ℒ_m√(-g)d^4x,
 
where relativistic geometrized units were employed i.e c=G=1. Here, f(R,𝒯) is an arbitrary function of the Ricci scalar R
and the trace 𝒯 of the energy-momentum tensor T_μ ν, while ℒ_m denotes the matter Lagrangian density.   
  
By variation of action S with respect to the metric tensor g_μν yields the following field equation 

    ( R_μν- ∇_μ∇_ν)f_R (R,𝒯) + f_R (R,𝒯)g_μν - 1/2 f(R,𝒯)g_μν
     = 8π T_μν  - f_𝒯(R,𝒯) (T_μν  +Θ_μν)
 

where f_R (R,𝒯) denote the partial derivative of f (R,𝒯) with respect to R and f_𝒯 (R,𝒯) is the partial derivative of f (R,𝒯) with respect to 𝒯, while R_μν is the Ricci tensor. The box operator ≡∂_μ(√(-g) g^μν∂_ν)/√(-g) is called the D'Alambert operator, and ∇_μ represents the
covariant derivative associated with the Levi-Civita connection of metric tensor g_μν. 
The stress energy tensor T_μν and Θ_μν are defined as follows, 

    T_μν   =   g_μνℒ_m-2∂ℒ_m/∂ g^μν,
    Θ_μν   =    g^αβδ T_αβ / δ g^μν.


Using Eq.(<ref>), the Einstein tensor G_μν can be written as,

    0.98!G_μν=1/f_R(R,𝒯)[8π T_μν+1/2(f(R,𝒯)-Rf_R(R,𝒯))g_μν
    .93!-T_μν f_𝒯(R,𝒯)-(g_μν-∇_ν∇_μ)f_R(R,𝒯) ].
 
After taking the covariant derivative of the Eq.(<ref>), we get


    ∇^μT_μν=f_𝒯(R,𝒯)/8π -f_𝒯(R,𝒯)[(T_μν+Θ_μν)∇^μln f_𝒯(R,𝒯)  
    +∇^μΘ_μν-1/2g_μν∇^μ𝒯]


The equation (<ref>) shows that the covariant derivative of stress-energy momentum tensor T_μν not vanishes in modified f(R,𝒯) gravity as in other theories of gravity. Throughout in our study,
we consider the Lagrangian matter  ℒ_m=-𝒫 where 𝒫=1/3(p_r+2p_t)<cit.>. Then using Eq.(4) we obtain Θ_μν=-2T_μν-𝒫 g_μν.

In order to obtain the effective stress-energy momentum tensor for modified theory of gravity  we consider the simplest linear functional form of f(R,𝒯) (proposed by Harko et al. <cit.>) as follows 

    f(R,𝒯)=R+2χ𝒯
  
where χ is a coupling constant. The above f(R,𝒯) function has been used widely used to develop the different f(R,𝒯) gravity compact objects. By inserting the value f(R,𝒯) from Eq. (<ref>) in Eq. (<ref>) we obtain

    G_μν=8π T_μν+χ𝒯g_μν+2χ(T_μν+𝒫 g_μν)
    =8πT̃_μν,


Here, we consider the energy momentum tensor T_μν corresponding to anisotropic fluid distribution which can be defined as,

    T_μν=(ρ+p_t)u_μ u_ν-p_tg_μν+(p_r-p_t)v_μ v_ν,

where u_ν is the four velocity, satisfying u_μu^μ=-1 and u_ν∇^μu_μ=0. Here, ρ is matter density,  while p_r and p_t are radial pressure and tangential pressure respectively. 
Now the modified energy-momentum tensor T̃_μν can be written as  

    T̃_μν= T_μν (1+χ/4π)+χ/8π(𝒯+2𝒫)g_μν
 
 
By inserting the value of f(R,𝒯)=R+2χ𝒯 in Eq.(<ref>) we get 

    ∇^μT_μν=-1/2 (4π+χ)χ[g_μν∇^μ𝒯+2 ∇^μ(𝒫  g_μν)].

By using the Eqs.(<ref>) and (<ref>) we can write,

    ∇^μT̃_μν=0

     



§ THE FIELD EQUATION FOR ANISOTROPIC MATTER DISTRIBUTIONS IN MODIFIED F(R,𝒯) GRAVITY

Let us consider the spacetime being static and spherically symmetric, which describes the interior of the object can be written in the following form

    ds^2 = -e^ν(r)   dt^2+ξ^-1 dr^2 +r^2(dθ ^2 +sin ^2θ  dϕ ^2),


Since the modified energy tensor T̃_μν is sum θ_μν which will clearly generates the anisotropic pressure within the effective matter distribution. Hence, using Using Eqs. (<ref>) and (<ref>) together with line element (<ref>) The field equations for the spherically symmetric anisotropic stellar system can be written as,


    8π ρ̃   =   [1/r^2-ξ/r^2-ξ^'/r]   
     
    8π p̃_r    =   [ξ(1/r^2+ν^'/r)-1/r^2]
    
    8π p̃_t   =   [ξ/4(2ν^''+ν^'2+2ν^'/r)+ξ^'/4(ν^'+2/r)].
 
The primes denote differentiation with respect to the radial coordinate r. Using the Eqs.(<ref>) and (<ref>) the effective quantities like effective density (ρ̃), effective radial pressure (p̃_r) and effective tangential pressure (p̃_r) can be written in terms of anisotropic pressures p_r and p_t, energy density ρ for anisotropic matter distribution as,  

    8 π ρ+χ /3( 9 ρ-p_r-2 p_t)   =   8 π ρ̃
    
    8 π p_r-χ /3( 3 ρ-7 p_r-2 p_t)   =   8 π p̃_r,
    
    8 π p_t-χ /3 ( 3 ρ-p_r-8 p_t)    =   8 π p̃_t,

Eqs.(<ref>) defines the gravitational mass inside the star of radius r. Using the Eqs.(<ref>) and (<ref>) with ξ=1-2 m(r)/r we get,

    m̃ =4π∫_0^r r^2 ρ̃ dr

For this modified anisotropic matter distribution, it is also necessary that the anisotropic fluid must satisfy another additional equation (known as conservation equation) as,

    p^'_r+ν^'/2 (ρ+p_r)-2/r (p_t-p_r)=χ (3ρ^'-p^'_r-2p^'_t)/6 (4 π+χ).
 
The above conservation Eq.(<ref>) is also known as modified Tolman-Oppenheimer-Volkoff (TOV) equation for modified theory of gravity. It is note that this modified conservation equation reduce into conservation equation for General relativity when χ=0.0.  On the other hand, the energy density (ρ), radial pressure (p_r) and tangential pressure (p_t) for anisotropic stellar model in modified f(R,𝒯) gravity theory can be written as, 

    ρ   =   8 π ρ̃/(8π+4χ)+8 π (3ρ̃+p̃_̃r̃+2p̃_t) χ/3 (8π+2χ) (8π+4χ)
    
    p_r   =   8 π p̃_̃r̃/(8π+2χ)+8 π (3ρ̃-p̃_̃r̃-2p̃_t) χ/3 (8π+2χ) (8π+4χ)
    
    p_t   =   8 π p̃_̃t̃/(8π+2χ)+8 π (3ρ̃-p̃_̃r̃-2p̃_t) χ/3 (8π+2χ) (8π+4χ)



 §.§ Solution of field Equations (<ref>)- (<ref>) in f(R, 𝒯) gravity theory:
  

To solve the equations (<ref>), (<ref>) and (<ref>), we use the isotropy condition in Eqs. (<ref>) and (<ref>) which leads a second order differential equation of the form as, 

    ξ(ν^''/2+ν^'2/4-ν^'/2 r-1/r^2)   +   ξ^'ν^'/4+2 ξ^'/r-1/r^2
     = 8π(p̃_t-p̃_r)   =   (8π+2 χ) (p_t-p_r),

The above pressure isotropy equations in f(R,𝒯) gravity will be same as in GR if χ=0. Since  Eq.(<ref>) contain three unknowns ν, ξ and Δ=p_t-p_r. Then to solve this equation, we choose a modified anstaz of gravitational potential ξ, proposed by Durgapal-Floria <cit.>, of the form as,

    ξ   =   1-8 x (3+x)/7 (1+χ+x)^2,


The choice of above metric potential (<ref>) is well motivated, because it is free from physical and mathematical singularities everywhere within the compact structure. Furthermore, yields to a finite, well defined and decreasing outward energy density at all points inside the star. The inclusion of χ in potential ξ will cause the effect in modified energy density (ρ̃).  Now it can be noted that if χ=0 and p_t=p_r, then ν=4 ln(1+A r^2) will satisfy the isotropy Eq.(<ref>). By keeping this point in our mind, we construct the expression for Δ, by using isotropy condition (<ref>) and potential ξ, of the form as,   

    Δ   =   8 x [5 - χ + x + 2 α Δ_1(x) + 
       α^2 (Δ_2(x) + χ Δ_3(x))]/7  (1 + χ + x)^3 (1 + α  x)^2 (8π+2 χ).

with   Δ_1(x)=[(-6 + 7  x + x^2 - χ  (6 + 5  x))],

Δ_2(x)=[ 7 + 7 χ^3 - 15  x - 2  x^2 + 21 χ^2  (1 + x) ],

   Δ_3(x)=[ (21 + 6  x + 4  x^2) ]. 
  

It is observe from Eq. (<ref>), the anisotropy Δ is zero at centre, and then p_t=p_r at centre. However, the other details of physical features for Δ has been discussed in Sec. (V). Now by substituting the gravitational potential ξ and anisotropic factor Δ from Eqs.(<ref>) and (<ref>) into Eq.(<ref>), and by using the transformation ν=2lnΨ and x=Ar^2  we obtain,


    d^2Ψ/dx^2+ 4  [3 - x + χ (3 + 2 x)]/(1 + χ + x) [7 + 7 χ^2 - 10 x - x^2 + 14 χ (1 + x)] dΨ/dx-2 (1 + α x)^-2 [-4 α  (3 - x + 3 χ+ 2 χ x) + α^2 f (x)] /(1 + χ + x) [7 + 7 χ^2 - 10 x - x^2 + 14 χ (1 + x)] Ψ=0.


where,   f(x)=[7 + 7 χ^3 - 15 x - 7 x^2 - x^3 + 21 χ^2 (1 + x) + χ (21 + 6 x + 5 x^2)]. 
It is note that the value Ψ=(1+α x)^2 satisfy the above differential Eq.(<ref>) which implies that this value of Ψ leads a particular solution of Eq.(<ref>). Then most general solution of Eq.(<ref>) is given (using the change of dependent variable method) as,


    Ψ(x)=Ψ_1(x)[C+D (√(7) √(Ψ_ 2 (x))[f_ 1 (x) + f_ 2 (x)]/3 Ψ_ 3 (x) Ψ_1(x) √(Ψ_1(x)) +√(7) f_3 (x) /(Ψ _ 3 (x))^7/2 ln[-a (Ψ_3 (x))^5/2  f_ 5 (x)/ 2   
         f_ 4 (x)  (1 + a x)])],
      
where C and D are arbitrary constant of integration and expression of the used coefficients are as follows:  


Ψ_ 1 (x) = (1 + α  x)^2,    Ψ_ 2 (x) = [7 + 7 χ^2 - 10 x - x^2 + 14 χ (1 + x)],  Ψ_ 3 (x) = [-1 + 7 α^2 (1 + χ)^2 - 2 α (-5 + 7 χ)],


f_ 1 (x) = (1 - α - α χ) Ψ^2_ 3 (x) + Ψ_ 3 (x) [-1 + 6 αχ + α^2 (-23 - 16 χ + 7 χ^2)] (1 + α x),


f_ 2 (x) = [1 + α (7 - 23 χ) + α^2 (63 + 174 χ - 105 χ^2) + α^3 (-359 + 201 χ + 399 χ^2 - 161 χ^3)] (1 + α x)^2,


f_ 3 (x) = 4 [-1 + 2 χ + α (-1 - 23 χ + 14 χ^2) + α^2 (37 + 36 χ - 147 χ^2 + 70 χ^3) + α^3 (-163 + 257 χ + 21 χ^2 - 301 χ^3 + 98 χ^4)],


f_ 4 (x) = [-1 + 2 χ + α (-1 - 23 χ + 14 χ^2) + α^2 (37 + 36 χ - 147 χ^2 + 70 χ^3) + 
    α^3 (-163 + 257 χ + 21 χ^2 - 301 χ^3 + 98 χ^4)],


f_ 5 (x) = 5 - 7 χ + x + √(Ψ_ 3 (x)) √(Ψ_ 2 (x)) + α[7 + 7 χ^2 - 5 x + 7 χ (2 + x)]. 
 


By plugging the value of ξ, ν=2lnΨ from Eqs. (<ref>), and (<ref>) into Eqs.(<ref>)- (<ref>) we find the effective energy density ρ̃, effective radial pressure (p̃_r) and effective tangential pressure (p̃_t) in f(R,𝒯) gravity theory,


    ρ̃   =   8 [9 + 2 x + x^2 + χ (9 + 5 x)]/56 π (1 + χ + x)^3,
    p̃_r   =   A [-8 (3+x)+4 Ψ_2(x) Ψ_4(x)]/56 π (1+χ+x)^2,
 

    p̃_t   =   A [-8 (3+x)+4 Ψ_2(x) Ψ_4(x)]/56 π (1+χ+x)^2
           +
    8 A x [ Ψ_5(x)+ Ψ_6(x)+χ (21 + 6  x + 4  x^2)]/56 π  (1 + χ + x)^3 (1 + α x)^2.
 
where, 

    Ψ_4(x)   =   √(7) D  (1 +χ+x)/Ψ_1(x) Ψ(x) √(Ψ_2(x)) + 2 α/√(Ψ_1(x)), 
    Ψ_5(x)   =   α^2 [5 - χ + x+2 α (-6 + 7  x + x^2 - χ  (6 + 5  x))], 
    Ψ_6(x)   =   α^2 [ 7 + 7 χ^3 - 15  x - 2  x^2 + 21 χ^2  (1 + x) ].






§ MATCHING CONDITION

Since all compact structures are bounded objects, to ensure a well behaved stellar interior i.e finite material content and smooth geometry at the surface Σ≡ r=r_s (where r_s is the radius of the sphere) of the configuration, one needs to join the inner space-time ℳ^- with the corresponding outer space-time ℳ^+. In this case we are trying with uncharged anisotropic fluid sphere described by (<ref>) and (<ref>). Moreover, due to the election given by (<ref>) on the f(R,𝒯) function, the appropriated exterior space-time ℳ^+ corresponds to Schwarzschild geometry given by   
 
    ds^2 = -(1-2M/r)^-1 dr^2-r^2 (dθ ^2+sin ^2θ  dϕ ^2 )
    +(1-2M/r) dt^2,

this is so because the modification introduced in matter sector represented by 𝒯 is vanishing beyond Σ. So, to join the interior geometry with Schwarzschild outer space-time one requires to impose the so called first and second fundamental forms across Σ. The first fundamental form refers the continuity of the intrinsic metric g_μν induced by both metrics ℳ^- and ℳ^+ on Σ. Explicitly it reads 

    [ds^2]_Σ=0 ⇒ e^λ^-(r_s)=e^λ^+(r_s)     e^ν^-(r_s)=e^ν^+(r_s).

So, for the present model we have
 
    1- 8 x_s (3+x_s)/7 (1+χ+x_s)^2   =   1-2M/r_s,
    Ψ^2(r_s)   =   ( 1-2M/r_s).

The second fundamental form is related with the continuity of the extrinsic curvature K_μν induced by ℳ^- and ℳ^+ on Σ. The continuity of K_rr component across Σ yields to 

    p̃_r(r_s)=0.

The above requirement determines the size of the object ı.e the radius r_s which means that the material content is confined within the region 0≤ r≤ r_s. The continuity of the remaining components K_θθ and K_ϕϕ leads to

    m̃(r_s)=M.

Equation (<ref>) is the total effective mass contained in the sphere which is expressed by Eq. (<ref>). 
After solving Eqs.(<ref>)-(<ref>) we obtain the parameters C, D and A as,

    D/C   =   2 (1 + α x_s)^3 √(Ψ_2(x_s)) [3 + x_s-α Ψ_7(x_s) ]/[F_ 1 (x_s) + x^2_s  F_ 2 (x_s) + x^3_s   F_ 3 (x_s) + F_11(x_s)],     
    
    C   =   1/Ψ_1(x_s) [1+D/C F(x_s)]√(7-10 x_s-x^2_s/7(1+χ+x_s)^2)
    
    A   =   6r_s-7 (1 + χ) M ± 2 √(7r_s (χ^2- χ-2) M + 9 r^2_s)/(7 M - 4 r_s)  r^2_s,
 
where, x_s=Ar^2_s,
Ψ_7(x_s)=(7+7  χ^2-13   x_s-2 x^2_s+14 χ (1 + x_s)),
F_11(x_s)=x _s  F_ 4 (x_s) - 7  χ^2   F_ 6 (x_s) - χ [F_ 7 (x_s) + x^2_s   F_ 5 (x_s)]. The expressions for the coefficients used here are mentioned in the Appendix <ref>. 
Then, Eqs. (<ref>)-(<ref>) obtained from Israel-Darmois <cit.> junction conditions (first and second fundamental forms) are the necessary conditions to get the compete constants parameter that characterize the model.



§ GEOMETRIC CHARACTERIZATION OF THE MODEL

The studied model is described by a spherically symmetric static manifold whose temporal and radial metric tensor components are given by Eqs. (<ref>) and (<ref>), respectively. These metric potentials are free from physical and mathematical singularities throughout the compact object. This fact guarantees a well behaved space-time region.
Once the parameters A, C and D have been determined, the behavior of the metric potentials is studied through a graphical analysis. As shown in Fig. <ref> both ξ^-1 and e^ν are regular function with increasing radial coordinate everywhere inside the star for all values of the parameter χ. As usual ξ^-1 takes the value 1 at r=0 while e^ν>0 at the same point. The choice of Durgapal-Fuloria type potential ξ has well-founded physical reasons, it is clearly free from physical and mathematical discontinuities.  The extension carried out in this work to the context of f(R,𝒯) gravity theory maintains the same spirit. However, we have made a small modification in order to include the effects of this modified gravitational theory in the full obtained model. From Fig. <ref> the dashed red line (χ=0.0) represents the anisotropic Durgapal-Fuloria model in the GR framework. Moreover, the effect of χ on both metric potentials is clear. Taking χ=-0.1 the potential e^ν at the center is smaller than when χ=0.1 but for ξ^-1 we have the reverse situation ı.e for χ=-0.1 the curve takes larger values towards the boundary Σ than χ=0.1. On the other hand, GR curve behaviour is an intermediate value between them.



§ THERMODYNAMIC OBSERVABLES 

In this section we study and analyze the behaviour of the main salient features of the model ı.e matter density ρ̃, radial p̃_r and tangential p̃_t pressure respectively. Also we examine the role played by the effective anisotropy factor Δ̃ inside the stellar structure. 

It is well known that the main physical parameters of any compact object describing stellar interiors should be free from physical and mathematical drawbacks. Furthermore, they should be  monotonic decreasing functions of the radial coordinate towards the surface, with their maximum values attached at the center of the configuration. These general requirements ensure in principle a well behaved model which can serves to describe some natural objects like white dwarf, neutron stars even quark stars. Moreover, in the study of compact structures there are other ingredients as essential as the aforementioned, which provide a more realistic view of the behavior of celestial bodies. For example, the inclusion of anisotropies in the material content contained in the fluid sphere. Anisotropy in this context means that the pressure in the radial direction is different from the pressure in the angular directions ı.e p_r≠ p_t. So, the effective anisotropy factor is defined by Δ̃=p̃_t-p̃_r. The inclusion of anisotropies within the stellar content introduces improvements in stability and balance mechanisms and increases the value of the surface redshift. However, regarding the equilibrium mechanism the contribution that it will give depends on the sign, that is, whether it is positive Δ̃>0⇒p̃_t>p̃_r or negative Δ̃<0⇒p̃_t < p̃_r. In the first case the system experiences a repulsive force that helps to counteract the gravitational gradient and in the second case, the force due to the anisotropy helps the gravitational force compress the object. If the pressure exerted by the nuclear force fails to overcome the gravitational attraction, the structure eventually will continue to collapse until its Schwarzschild radius. At this point, the object forms a black hole with many unusual properties.  This means that the presence of an attractive force due to anisotropies damages the balance and stability of the configuration. It is clear that the collapse of the structure towards a singularity depends on the gradient pressure (hydrostatic force) exerted by the matter inside the star.
Figure <ref> shows the behaviour of all thermodynamic observables and anisotropy factor. From the upper panels we can see the behaviour of both effective radial and tangential pressure (left and right respectively). These physical quantities have their maximum values at the center of the configuration and are monotonic decreasing functions with increasing radial coordinate. It is observed that for negative values of χ the maximum value is greater than the values obtained considering 0.0 (GR limit) and 0.1. Respect to the effective density (lower right panel), it has its maximum value attained at the center corresponding with χ=-0.1, is monotonic decreasing function towards the surface and positive defined everywhere within the star. Then, all the thermodynamic observables increases at the center of the star when χ moves from -0.1 to 0.1. The behavior of the effective anisotropy factor Δ̃ (lower left panel) is strongly dependent on the value that χ takes. For χ=0.1 its behaviour is positive at all points within the star, vanishing at the center and increasing function with increasing radius. As we explained above, this conduct introduces in the system a repulsive force (outward). On the other hand, for χ=-0.1 the system is subject to an attractive force (inward). The effect on the system caused by this attractive force will be analyzed in the dynamical equilibrium section <ref>. Finally, the GR limit corresponding to χ=0.0 shows a positive anisotropy factor throughout the object and attains its maximum value within the configuration. Furthermore, comparing GR (χ=0.0) with χ<0 it is observed that in f(R,𝒯) gravity the object are more compact than in GR. Table <ref> displays the effect of χ on various physical parameters of the star such as the radius r_s, central and surface effective energy density, central effective pressure, surface redshift Z_s and compactness factor u=M/R for the same Mass=1.04M_⊙. As we can see, moving χ from -0.1 to 0.1 all these quantities (except r_s and Z_s) are decreasing in magnitude.
 



§ ENERGY CONDITIONS

It is well known that the matter distribution that makes up celestial bodies can be composed of a large number of material fields. Despite knowing the components that describe this material content inside the compact structure, it could be very complex to describe exactly the shape of the energy-momentum tensor. In fact, one has some ideas on the behaviour of the matter under extreme conditions of density and pressure.  

On the other hand, there are certain inequalities which are physically reasonable to assume for the energy-momentum tensor. So, in this section we are willing to verify these inequalities at all points in the interior of the star. In the literature these inequalities are known as energy conditions. Then we have the null energy condition (NEC), dominant energy condition (DEC), strong energy condition (SEC) and weak energy condition(WEC). Explicitly, these are given by 

    WEC   :    T_μνl^μ l^ν≥ 0  ρ≥  0, ρ+p_i ≥ 0 
    NEC   :    T_μνt^μ t^ν≥ 0  ρ+p_i ≥  0
    DEC   :    T_μνl^μ l^ν≥ 0   ρ≥ |p_i|
            T_μνl^μ∈
    SEC   :    T_μνl^μ l^ν - 1  2 T^λ_λ l^σ l_σ≥ 0   ρ+∑_i p_i ≥ 0.

where i≡ (radial r, transverse  t), l^μ and t^μ are time-like vector and null vector respectively. 
 
To verify a well defined energy-momentum tensor everywhere within the compact configuration the above inequalities must satisfy simultaneously.
We will check the energy conditions with the help of graphical representation. In Fig. <ref>, we have plotted the L.H.S of the above inequalities which verifies that  all the energy conditions are satisfied at the stellar interior.

These energy conditions beyond capturing the idea that the energy must be positive defined, have a clear physical and geometric interpretation <cit.>. From the physical point of view NEC means that an observer traversing a null curve will measure the ambient (ordinary) energy density to be positive. WEC implies that the energy density  measured by an observer crossing a timelike curve is never negative. SEC purports that the trace of the tidal tensor measured by the corresponding observers is always non-negative and finally DEC stand for mass-energy can never be observed to be flowing faster than light. Furthermore, violations of the energy conditions have sometimes been presented as only being produced by unphysical stress energy tensors. Usually SEC is used as a fundamental guide will be extremely idealistic. Nevertheless, SEC is violated in many cases, e.g. minimally coupled scalar field and curvature-coupled scalar field theories. It may or may not imply the
violation of the more basic energy conditions i.e. NEC and WEC. It is worth mentioning that both SEC and DEC imply NEC and DEC implies also WEC. Additionally, the fulfillment of SEC and DEC conditions imposes strong restrictions on the maximum plausible bound of the surface redshift Z_s of the compact structure when there are anisotropies in the stellar interior. These implications will be discussed in more details in the next section. 





§ CAUSALITY AND SURFACE REDSHIFT

Among the modifications introduced by the presence of anisotropies in the stellar interior. We have the velocities of propagation associated with the pressure waves in the main directions of the sphere, that is, in the radial and transverse directions and the modification of the upper bound of the surface redshift Z_s. First of all, the subliminal speeds corresponding to each direction are defined by

    v^2_r=dp̃_r/dρ̃     v^2_t=dp̃_t/dρ̃.

In order to obtain a physically admissible model, both speeds v_r and v_t must be bounded  by the speed of light (c=1 in relativistic geometrized units). This tells us that the pressure (sound) waves in the fluid do not propagate at arbitrary speeds. This is known as causality condition. This condition is peremptory regardless if the material content of the star is isotropic or anisotropic. The only difference between the mentioned cases is that for the anisotropic case there is propagation in the two main directions of the sphere ı.e radial and transverse directions. Moreover, in the isotropic case the subliminal sound speed, should be a decreasing function. However, this is not true in the case where there is anisotropy, since the  speed behaviour depends on the rigidity of the material. So, causality condition reads

    0≤ v_r≤ 1       0≤ v_t≤ 1.

Causality condition (<ref>) has strong implications on the behavior of the matter distribution within the object. One of them is related to energy-momentum tensor that describes the material content. If causality is preserved then the energy-momentum tensor is well defined. Secondly, imposing this important condition and once the radial pressure p̃_r is obtained, a relation between it and the density ρ̃ can be established ı.e an equation of state p̃_r=p̃_r(ρ̃). This last statement is very important because usually in searching  solutions to Einstein field equations, the equation of state is imposed. However, it often results in the violation of causality. Additionally, the fact of having different speeds in the directions mentioned above, influences the stability of the system (this subject will be discussed in section A). From Fig. <ref> it is appreciated that v^2_r and v^2_t are satisfying causality condition for all χ. For χ<0 both subliminal sound speeds are decreasing in nature. Nevertheless, v^2_t is greater than v^2_r at all points in the star. In the case χ>0 the radial speed is always greater than the tangential one, whilst in GR scenario the radial sound velocity is greater than the tangential velocity for 0<r<8.802 and then has an opposite behavior for r>8.802.

The surface redshift Z_s a significant observational parameter that relates the mass m̃(r_s)=M and the radius r_s of the star, is affected when anisotropies are introduced into the system, regardless the mechanism that originated them. For isotropic fluid spheres its maximum value is Z_s=2. This value is determined by Buchdahl constraint on the compactness factor u=2M/r_s≤ 8/9 <cit.>. The explicit relation between Z_s, M and r_s is given by

    Z_s=(1-2M/r_s)^-1/2-1.

So, the effect of anisotropies on Zs has a long history. For example, Bowers and Liang <cit.> considered an hypothetical model containing a constant density ρ=ρ_0 (incompressible fluid) and a specific form of the anisotropy factor Δ. They concluded that when the anisotropy factor is null i.e Δ=0 ⇒ p_r=p_t the maximum value for the surface redshift corresponds to Z_s=4.77, and in the case of a positive anisotropy factor Δ>0 ⇒ p_t>p_r the above value can be exceed (otherwise if Δ<0). Moreover, if the anisotropy factor is extremely large then the surface redshift will be too.
Moreover, Ivanov studies shown that for realistic anisotropic star models obeying SEC the maximum surface redshift is Z_s=3.842 (this value corresponds to a model without cosmological constant) meanwhile for models satisfying DEC is given by Z_s=5.211 <cit.>. These values correspond to the following mass-radius relation 0.957 and 0.974, respectively. 
In Fig. <ref> the surface redshift Z_s has a monotonic increasing behaviour towards the boundary with its maximum value attained at the boundary of the object. Besides, for negative values of χ the surface redshift takes larger values in comparison with χ≥0. This is also appreciated in table <ref>. Although Ivanov's research reveals that the value of Z_s in the presence of anisotropic fluids exceeds the Buchdahl bound, the values reported for Z_s in this work are below the maximum values reported for both isotropic and anisotropic distributions, but are in complete agreement with what has already been established, in the sense that they do not exceed such values.







§ DYNAMICAL EQUILIBRIUM CONDITION OF STELLAR MODEL IN F(R,𝒯) THEORY GRAVITY
 
 In this section we discuss the dynamical equilibrium condition of the stellar model by using Tolman-Oppenheimer-Volkoff  (TOV) <cit.> equation in the framework of f(R,𝒯) theory gravity. This modified TOV equation for f(R,𝒯) theory, as already mentioned by Eq.(<ref>), is given by 
  
    -dp_r/dr-ν^'/2( ρ+      p_r )+2/r(p_t-p_r)
    +      χ/6(4π+χ)(3 dρ/dr-dp_r/dr-2 dp_t/dr)=0.

where we denote first term -dp_r/dr=F_h, second term -ν^'/2( ρ+p_r )=F_g, third term 2/r(p_t-p_r)=F_a and fourth term χ/6(4π+χ)(3 dρ/dr-dp_r/dr-2 dp_t/dr)=F_χ. These term describe the hydrostatic force (F_h), gravitational force (F_g), anisotropic force (F_a) and coupling force (F_χ), respectively. 

In the case of isotropic fluid spheres (p̃_r=p̃_t) and regarding χ=0.0 (GR limit), this equation drives the equilibrium of relativistic compact structures such as neutron stars, white dwarfs, etc. In considering the inclusion of anisotropies and the effect of relativistic modified gravity theories, this equation still drives the balance of the system. However, its form change a little bit as shown Eq. (<ref>). As we can see this equation relates the effective thermodynamic quantities with the metric potential e^ν. In order to keep the system in equilibrium and prevent it from collapsing below its Schwarzschild radius, it is necessary to have a relationship between the thermodynamic variables p̃_r and ρ̃, that is, an equation of state (EoS) p̃_r=p̃_r(ρ̃) that links them. Nonetheless, in this case where there are contributions from the anisotropies and from the theory considered, for certain conditions it may not be enough to withstand the gravitational attraction. Thus
the structure equations (<ref>) and (<ref>) imply that there is a maximum mass that a star can have. As was pointed out before, the present model is under four forces. The impact that these have on the system is shown in fig. <ref>. We can observe that the system is in equilibrium, the gravitational gradient is counterbalances by the hydrostatic F_h, anisotropic F_a and coupling F_χ forces (although its contribution is very small) when χ>0. On the other hand, when χ<0 the anisotropic force takes negative values. It means that the system is under an attractive force. This fact his can damage the balance of the system if the hydrostatic gradient is not strong enough to counteract the force exerted by the gravitational attraction and the anisotropic force. Thus the system can collapse towards a singularity. However, as noted, the negative anisotropic force is very small in magnitude and the pressure gradient overcomes the action F_g+F_a.

Another interesting point to note is that for positive values of χ the hydrostatic gradient is lower in f(R,𝒯) theory than in the corresponding one in GR (the inverse situation is presented for χ<0). The same happens with the gravitational attraction.



 



§ THE EQUATION OF STATE (EOS)

In the study of compact structures such as neutron stars it is very important to know how the principal thermodynamic variables are connected. This relation known as the equation of state (EoS) drives a relationship between the effective radial pressure p̃_r and the effective energy density ρ̃. The microphysics, as described by the EoS, is linked to the macroscopic properties of the neutron star, in particular, their masses and radii, via the Tolman-Oppenheimer-Volkoff equations, which provide the direct relation that is necessary to use astrophysical observations to constrain nuclear physics at very high densities. However, the composition of a neutron star chiefly depends on the nature of strong interactions. Depending on the type of interaction, the models can be grouped into three broad categories: non-relativistic potential models,
relativistic field theoretical models, and relativistic Dirac-Brueckner-Hartree-Fock models. In addition, in each of these perspectives,  the presence of softening components such as hyperons, Bose condensates or quark matter, can be incorporated <cit.>. On the other hand, one can classify the EoS in two classes:  First, normal equations of state have a
pressure which vanishes as the density tends to zero. Second, self-bound equations of state
have a pressure which vanishes at a significant finite density. Respect to the self-bound EoS the most famous example is the MIT bag model EoS. It was pointed out by Witten <cit.> that strange quark matter is the ultimate ground state of matter. This leads to the fact that the internal and external vacuum densities of the hadrons are completely different and that the vacuum pressure of the bag wall balances the pressure of the quarks, stabilizing the whole system <cit.>. So, the MIT bag EoS model reads

    p=1/3(ρ-4B),

where B is the so called bag constant and represents the difference between the energy density of the perturbative and non-perturbative QCD vacuum. In this model the interactions of quarks and gluons are enough small, neglecting quark masses and supposing that quarks are confined to the bag volume. Concerning normal matter the EoS describes an interacting nucleon gas above a transition
density 1/3ρ_s to 1/2ρ_s (being ρ_s the surface density). Below this density, the ground state of matter consists of
heavy nuclei in equilibrium with a neutron-rich, low-density gas of nucleons. Nonetheless, the equilibrium of the system exists below the transition density <cit.>. So, in order to explain the structural properties of compact stars model at high densities, several authors have proposed the EoS P=P(ρ) should be well approximated by a linear function of the energy density ρ
<cit.>. Some authors have also expressed more convincing approximated forms of the EoS P =P(ρ)  as linear function of energy density ρ <cit.>. Furthermore, a linear relation between pressure P and energy density ρ ensure the preservation of causality condition.

For the present model the EoS is not a linear relation between pressure and energy density. The functional relation is more complicated. The effective radial pressure p̃_r in terms of the surface ρ̃_s and central ρ̃_c effective energy density the EoS has the following form

    p̃_r   =   ρ̃_c [-8 (3+ρ_1)+4 Ψ_2(ρ_1) Ψ_4(ρ_1)] (1+χ)^2/72  (1+χ+ρ_1)^2,     
    p̃_t   =   ρ̃_c (1+χ)^2/72(-8 (3+ρ_1)+4 Ψ_2(ρ_1) Ψ_4(ρ_1)/(1+χ+ρ_1)^2
           +
    8 ρ_1 [ Ψ_5(ρ_1)+ Ψ_6(ρ_1)+χ (21 + 6 ρ_1 + 4 ρ_1^2)]/(1 + χ + ρ_1)^3 (1 + α ρ_1)^2).
  
However the expressions for other used coefficients are mentioned in the Appendix <ref>. In Fig. <ref> we can graphically appreciate the shape of the EoS of the model under study. In spite of the complex relation given by <ref>, the behavior that appears from the surface to the core of the object is linear (the curve grows from regions of low to high densities). This curve can be described approximately by the following polynomial linear interpolation (keeping only the first order in ρ̃_s) as follows

    p̃_r=α(ρ̃-ρ̃_s),

where α is a non-negative constant. It is clear from Eq. (<ref>) that when r=r_s then p̃_r=0. This is so because ρ̃(r_s)=ρ̃_s ı.e the energy density at zero pressure (surface energy density). Additionally, we can observe from Fig. <ref> the effect of χ on the EoS. For χ=-0.1 the effective central pressure p̃_r, the effective central energy density ρ̃ and the effective surface energy density ρ̃_s are greater than when χ≥ 0. It can also be checked in table <ref>.





§ STABILITY

 In this section we analyze the stability of the model by means of Abreu's criterion <cit.>. Basically, the method consists in contrasting the speed of pressure waves in the two principal directions of the spherically symmetric star: the subliminal radial sound speed with the subliminal tangential sound speed, and then based on those values at particular points in the object, one could potentially conclude whether the system is stable or unstable under cracking instability. Put forward by Herrera <cit.>, cracking involves the possibility of smashing the fluid sphere in view of the appearance of total radial forces of different signs, and therefore in different directions, at different points within the configuration. It should be emphasized that this effect has never been observed, however under appropriate physical assumptions, it is a likely scenario.

Cracking process is mechanism to study  instability when anisotropy matter distributions are present. Nevertheless, this mechanism can be characterized most easily through the subliminal speed of pressure waves. Following we have

    δΔ̃/δρ̃∼δ(p̃_t-p̃_r)/δρ̃∼δp̃_t/δρ̃-δp̃_r/δρ̃∼ v^2_t-v^2r.

Moreover, from causality condition one has 0≤ v^2_r≤ 1 and 0≤ v^2_t≤ 1 which implies 0≤ |v^2_t-v^2_r|≤ 1. Explicitly it reads

    -1≤ v^2_t-v^2_r≤ 1  =
         {[   -1≤ v^2_t-v^2_r≤ 0   Potentially stable;    0< v^2_t-v^2_r≤ 1 Potentially unstable ]}.

Therefore, the main idea behind Abreu's criterion is that if the subliminal tangential speed v^2_t is larger than the subliminal radial speed v^2_r; then this could potentially result in cracking instabilities to occur in the object, rendering the latter an unstable configuration. So, with the help of graphical analysis one can determines the potentially stable/unstable regions within the star and then conclude whether the system is stable or not. From Fig. <ref> (upper plot) it is observed that the system presents all the regions completely stable, that is, the whole system is stable (for χ>0). Completely unstable behavior (for χ<0), it is also observed that the region 0<r<8.802 is stable while the region r>8.802 is unstable (for χ=0). Then, there is cracking (change in sign of v^2_t-v^2_r ) in GR theory for the present model. This is so because from r>8.802 the tangential velocity is greater than the radial velocity of the pressure waves. Instability of the system for negative values of χ could be anticipated from Fig. <ref>, since v^2_t is greater than v^2_r at every point inside the star. In distinction for positive values of χ where v^2_t is always less than v^2_r everywhere within the object. Although |v^2_t-v^2_r| lies between 0 and 1 (lower panel Fig. <ref>), it does not means that the system is stable. Notwithstanding, it warns us of the presence of cracking in the sphere, like in the case of GR. As seen in the lower panel of Fig. <ref>, the red curve (GR) is decreasing up to certain value of the radial coordinate r and then suddenly changes its behavior to a growing one. Therefore, in this opportunity for the specific value χ=0.0, the system has stable and unstable regions.





§ SUMMARY AND OUTLOOK

In the present paper we have obtained an analytic relativistic anisotropic spherical model in the framework of f(R,𝒯) gravity theory. This was achieved through the imposition of three ingredients, which form the fundamental pillars of the obtained model. The first of these was to consider a simple, notable and viable modified gravity model given by f(R,𝒯)=R+2χ𝒯 <cit.>, being R the usual Ricci scalar, 𝒯 the trace of the energy momentum-tensor and χ a constant coupling. In addition we have taken the Lagrangian density matter to be ℒ_m=1/3(p_r+2p_t). The second one is the imposition of the metric potential ξ, which corresponds to a modification of the original potential proposed by Durgapal-Fuloria <cit.> in the context of GR. The modification on ξ considers the inclusion of χ (see Eq. (<ref>)). Finally, we have imposed the isotropic condition under the restriction p_r≠ p_t in order to obtain the effective anisotropy factor Δ̃. With these ingredients in hand we arrive at the differential equation given by Eq. (<ref>) in order to get the e^ν metric potential. Once this differential equation is solved the inner geometry of the whole system is completely specified. The election of ξ is well motivated because it is free from physical and mathematical singularities everywhere inside the compact object. Moreover, the complete internal manifold (<ref>)-(<ref>) is well behaved at all point within the star. After that we proceeded to obtain the constant parameters of the solution. For this purpose we have made the joint between the model and the external Schwarzschild solution on the surface Σ of the compact structure. Thus, the first and second fundamental form provide the corresponding parameter space that characterizes the model. The junction between the collapsed configuration and the Schwarzschild space-time was possible due to the modifications introduced by 𝒯 in the matter sector remain finite and bounded by the object. It means that beyond Σ we have an empty space-time. 

It is of our interest to check the outcomes and the predictions of one of the extended gravity, i.e., f(R,𝒯) theory regarding the existence, stability and equilibrium of spherical stars. Therefore, we have explored the behaviour of the main salient features such as the effective radial and tangential pressure, the effective energy density and the effective anisotropy factor. The behaviour of all these quantities is influenced by χ. Throughout the study we have taken -0.1≤χ≤ 0.1, M=1.04M_⊙ and  α=1.12. For these values the effective thermodynamic observables p̃_r, p̃_t and ρ̃ remain positives and their maximum values are attained at the core of the compact star. It means that they have a monotonic decreasing behaviour from the center to the boundary. It is remarkable to note that χ=-0.1 induces greater central values than χ=0.1 and χ=0.0 (GR limit). Respect to the effective anisotropy factor, this quantity has a more intricate behavior, that is, for χ=-0.1 it is  negative decreasing function at every point inside the star, meanwhile for χ≥0 it is  positive increasing function everywhere within the object. A negative anisotropy factor means that the system is subjected to an attractive force (inward). This fact can damage the equilibrium and stability mechanism. In distinction with the case χ≥0 where the induced force is completely repulsive (outward) and helps to counteract the gravitational attraction in order to avoid point singularities or the formation of a black hole (the latter occurs when the object is compacted below its Schwarzschild radius). The proper behavior of the thermodynamic variables leads to a well-defined energy-momentum tensor. It is checked employing the so called energy conditions. In this study energy conditions holds for all values of χ. The realization of these conditions has strong implications for the preservation of causality (the material content inside the object does not travel faster than light) and the maximum possible value of the surface redshift, which can not be arbitrarily large when there is anisotropy. Fortunately, for this model the causality is completely preserved and the maximum value that the surface redshift takes for all χ values considered, does not exceed the maximum level reported in the literature <cit.>. The analysis of balance mechanism is studied using the modified version of the TOV
equation in the realm of couple of f(R,𝒯) model. For the particular choice of f(R𝒯) model, it is shown that the system is under four forces, namely the hydrostatic force F_h, gravitational force F_g,  anisotropic force F_a and the coupling force F_χ,
attaining the equilibrium condition by keeping all the forces sum to be zero. It is worth mentioning that the usual GR forces are
being modified due to f(R,𝒯) model, thus producing some
extra effects in the forces F_h, F_g and F_a. Furthermore, we have arrived at a more compact configuration for χ<0. This is so because, the attractive force due to negative anisotropy helps compact the object in conjunction with the gravitational force. Moreover, we have obtained the corresponding equation of state (EoS) by analytic procedure. This important relation between the effective radial p̃_r and the effective energy density ρ̃ tell us (in some sense) how is the material composition of the star. In addition, the EoS allows to determine the mass-radius relation (it is also related with the surface redshift), ı.e from the microphysic one can get the macro observables of the compact structure. Despite of the complex expression describing the EoS of the model given by Eq. (<ref>), its behaviour is approximately linear. This indicates that the interior structure of
the relativistic compact stellar objects is composed of normal ordinary matter. 
The stability of our
compact star depends upon the choices of the  parameter involved in f(R,𝒯) model. The difference of the squares subliminal sound speeds, i.e, |v^2_t-v^2_r| has been found to be within
[0, 1], for χ=0.0 the system presents cracking, meanwhile for χ=-0.1 and χ=0.1 we can not conclude anything. To check if the system is stable/unstable we need to study v^2_t-v^2_r. It was found that the system is completely unstable when χ=-0.1 and completely stable for χ=0.1. Furthermore, for the GR case (χ=0.0) the system has stable and unstable regions (cracking). All the above discussion is supported by an extensive graphic analysis which is states in Figs. <ref>, <ref>, <ref>, <ref>, <ref>, <ref>, <ref> and <ref>. Finally, we wish to comment that by taking χ=0.0 and α=1 the corresponding results of the Durgapal-Fuloria model reported in <cit.> are recovered and that in the case χ>0 the obtained model satisfies all the general requirements in comparison with GR. Then, the existences of compact structures within the framework of f(R,𝒯) is a good opportunity to understand many phenomena in the strong gravitational field regime.

 
S. K. Maurya acknowledge continuous support
and encouragement from the administration of University of Nizwa. F. Tello-Ortiz thanks 
the financial support by the CONICYT PFCHA/DOCTORADO-NACIONAL/2019-21190856, grant Fondecyt No. 1161192, Chile and project ANT-1855 at the Universidad de Antofagasta, Chile.


99

r1 T. Harko, F. S. N. Lobo, S. Nojiri and S. D. Odintsov, Phys. Rev. D, 84, 024020 (2011)

r2 A. G. Riess et al., Astron. J. 116, 1009 (1998)

r3 S. Perlmutter et al., Astrophys. J. 517, 565 (1999)

r4P. de Bernardis et al.,
Nature 404, 955 (2000)

r5 S. Hanany et al., Astrophys. J. 545, L5 (2000)

r6 P. J. E. Peebles and B. Ratra, Rev. Mod. Phys. 75, 559 (2003)

r7T. Padmanabhan, Phys. Repts. 380, 235 (2003)

r8P. H. R. S. Moraes, Astrophys. Space Sci., 352, 273 (2014)

r9P. H. R. S. Moraes, Eur. Phys. J. C 75, 168 (2015)

r10P. H. R. S. Moraes, Int. J. Theor. Phys., 55, 1307 (2016)

r11C. P. Singh and P. Kumar, Eur. Phys. J. C, 74, 11 (2014)

r12 H. Shabani and M. Farhoudi, Phys. Rev., D, 90, 044031 (2014)

r13 D. R. K. Reddy and R. S. Kumar, Astrophys. Space Sci., 344, 253 (2013)

r14 P. Kumar and C. P. Singh, Astrophys. Space Sci., 357, 120 (2015)

r15M. Sharif and Z. Yousaf, Astrophys. Space Sci., 354, 471 (2014)

r16Z. Yousaf, M. Ilyas and M. Z. Bhatti, Eur. Phys. J. Plus, 132, 268 (2017)

r17Z. Yousaf, M. Ilyas and M. Z. Bhatti, Mod. Phys. lett., A 32, 1750163 (2017)

r18M. Z. Bhatti and Z. Yousaf, Ann. Phys., 387, 283 (2017)

r19M. J. S. Houndjo, Int. J. Mod. Phys. D, 21, 1250003 (2012)

r20E. H. Baffou, A. V. Kpadonou, M. E. Rodrigues, M. J. S. Houndjo and J. Tossa,
Astrophys. Space Sci., 356, 173 (2014)

r21S. S. Yazadjiev, D. D. Doneva and K. D. Kokkotas, Phys. Rev. D, 91, 084018 (2015)

r22 M. Sharif and Z. Yousaf, Gen. Relativ. Gravit., 47, 48 (2015)

r23 M. Z. Bhatti, Z. Yousaf and M. Ilyas, Eur. Phys. J. C, 77, 690 (2017)

r24 Z. Yousaf, M. Sharif, M. Ilyas and M. Z. Bhatti, Eur. Phys. J. C, 77, 691 (2017)

r25P. H. R. S. Moraes, J. D. V. Arbañil and M. Malheiro, J. Cosmol. Astropart. Phys., 06,
005 (2016)

r26D. Deb, Rahaman, Farook, S. Ray and B. K. Guha, Phys. Rev. D, 97, 084026 (2018)

r27D. Deb, Rahaman, Farook, S. Ray and B. K. Guha,  J. Cosmol. Astropart. Phys., 1803, 044 (2018)

r28S. Hansraj and A. Banerjee Phys. Rev. D, 97, 104020 (2018)

harko1 K S. Cheng and T. Harko, A&A, 385, 947 (2002)

Durgapal M.C. Durgapal, R.S. Fuloria, Gen. Relativ. Gravit. 17, 671 (1985)


r29 R. L. Bowers and E. P. T. Liang, Astrophys. J., 188, 657, (1974)

r30 M. Cosenza, L. Herrera, M. Esculpi and L. Witten, J. Math.
Phys., 22, 118 (1981)

r31 M. Cosenza, L. Herrera, M. Esculpi and L. Witten, Phys. Rev.
D, 25, 2527 (1982)

r32 L. Herrera and J. Ponce de León, J. Math. Phys., 26, 2302 (1985)

r33 J. Ponce de León, Gen. Relativ. Gravit., 19, 797 (1987)

r34 J. Ponce de León, J. Math. Phys., 28, 1114 (1987)

r35 R. Chan, S. Kichenassamy, G. Le Denmat and  N.O. Santos, Mon. Not. R. Astron. Soc., 239, 91 (1989)

r36 R. Chan, L. Herrera and N.O. Santos, Class. Quantum Grav., 9, 133 (1992)

r37 L. Herrera, Phys. Lett. A, 165, 206 (1992)

r38 R. Chan, L. Herrera and N.O. Santos, Mon. Not. R. Astron. Soc., 265, 533 (1993)

r39 A. Di Prisco,  E. Fuenmayor, L. Herrera and V. Varela, Phys. Lett. A, 195, 23 (1994).

r40 L. Herrera and  N. O. Santos,  Phys. Rep., 286, 53 (1997)

r41 A. Di Prisco, L. Herrera and V. Varela, Gen. Relativ. Gravit., 29, 1239 (1997)

r42 L. Herrera, A.D. Prisco, J. Ospino and E. Fuenmayor, J. Math. Phys., 42, 2129 (2001)

r43 L. Herrera, J. Ospino and A.D. Prisco, Phys. Rev. D, 77, 027502 (2008)

r44 M. K. Gokhroo and A. L. Mehra,  Gen. Rel. Grav., 26, 75 (1994)

r45 B.V. Ivanov, Phys. Rev. D, 65, 104011 (2002)

r46 S.K. Maurya, A. Banerjee, and S.Hansraj, Phys. Rev. D, 97, 044022 (2018)

DebMNRS D. Deb , S. V. Ketov, S. K. Maurya, M. Khlopov, P. H. R. S. Moraes and  S. Ray, Mon. Not. R. Astr. Soc. 485, 5652 (2019) 

r47 W. Israel, Nuovo Cim. B, 44, 1 (1966)

r48 G. Darmois, Mémorial des Sciences Mathematiques (Gauthier-Villars, Paris, 1927), Fasc. 25 (1927)

curiel E. Curiel, Einstein Stud., 13, 43 (2017) 

r51 H. A. Buchdahl, Phys. Rev. D, 116, 1027 (1959)

r49 R. C. Tolman,  Phys. Rev., 55, 364 (1939)

r50 J.R. Oppenheimer and G.M. Volkoff, Phys. Rev., 55, 374 (1939)

lattimer J. M. Lattimer, M. Prakash, D. Masak and A. Yahil, ApJ., 355, 241 (1990)

prakash M. Prakash, I. Bombaci, Manju. Prakash, J. M. Lattimer, P. J. Ellis and R. Knorren, Phys. Rep., 280, 1 (1997)

witten E. Witten,  Phys. Rev., D 30, 272 (1984)

alcock C. Alcock, E. Farhi and A. Olinto, AJ, 310, 261 (1986)

haensel P. Haensel and J. L. Zdunik, Nature, 340, 617 (1989)

baym G. Baym, C. J. Pethick and P. Sutherland, ApJ, 170, 299 (1971)

lattimer1 J. M. Lattimer, C. J.  Pethick, D. G. Ravenhall and D. Q. Lamb,  Nucl. Phys., A432, 646 (1985)

dey M. Dey, I. Bombacci, J. Dey, S. Ray and B. C. Samanta,
Phys. Lett. B, 438, 123 (1998)

harko K. S. Cheng and T. Harko, Phys. Rev. D, 62, 083001 (2000)

gondek D. Gondek-Rosinska, T. Bulik, L. Zdunik, et al., A&A, 363, 1005 (2000)
mauryaprd2 S.K. Maurya, A. Banerjee, M K. Jasim, J. Kumar, A K. Prasad, and A. Pradhan, Phys. Rev. D 99, 044029 (2019)

r52H. Abreu, H. Hernández and L. A. Núñez,Calss. Quan-tum. Grav., 24, 4631 (2007)






§ THE EXPRESSIONS FOR USED COEFFICIENTS IN EQS.(<REF>)  AND (<REF>) ARE AS FOLLOWS:
  

    Ψ_ 1 (x_s)    =    (1 + α  x_s)^2,     Ψ_ 2 (x_s) = [ 7 + 7χ^2 - 10 x_s - x_s^2 + 14χ (1 + x_s)],   Ψ_ 3 (x_s) = [-1 + 7α^2 (1 + χ)^2 - 2α (-5 + 7χ)], 
    
    F(x_s)   =   √(7) √(Ψ_ 2 (x_s))[f_ 1 (x_s) + f_ 2 (x_s)]/3 Ψ_ 3 (x_s) Ψ_1(x_s) √(Ψ_1(x_s)) +√(7) f_3 (x_s) /(Ψ _ 3 (x_s))^7/2 ln[-a (Ψ_3 (x_s))^5/2  f_ 5 (x_s)/ 2   
         f_ 4 (x_s)  (1 + α x_s)],
    
    f_ 1 (x_s)     =     (1 - α - α χ)Ψ^2 _ 3 (x_s) + Ψ_ 3(x_s) [-1 + 
        6 α χ + α^2  (-23 - 16χ + 7χ^2)] (1 + α x_s), 
    
    f_ 2 (x_s)    =    [1 + α (7 - 23χ) + α^2 (63 + 174χ - 105χ^2) + α^3 (-359 + 201χ + 399χ^2 - 161χ^3)]  (1 + α  x_s)^2, 
    
    f_ 3 (x_s)     =     4[-1 + 2χ + α (-1 - 23χ + 14χ^2) + α^2 (37 + 36χ - 
           147 χ^2 + 70 χ^3) + α^3  (-163 + 257 χ + 21 χ^2 - 301 χ^3 + 
            98 χ^4)],
    
    f_ 4 (x_s)     =    [-1 + 2χ + α (-1 - 23χ + 14χ^2) + α^2 (37 + 36χ - 
           147χ^2 + 70χ^3) + α^3 (-163 + 257χ + 21χ^2 - 301χ^3 + 
            98χ^4)], 
    
    f_ 5 (x_s)     =     5 - 7χ +  x_s + √(Ψ_ 3 (x_s)) √(Ψ_ 2 (x_s)) + α[
     7 + 7χ^2 - 5 x_s + 7χ (2 + x_s)], 
    
     F_ 1 (x_s)     =     - 7 √(7) - 7 √(7) χ^3 + (6 - 14  α + 2  α^3  (7 + 13  α)  x^4_s  + 4  α^4   x^5_s  )  F (x_s) √(Ψ_ 2 (x_s));, 
    
    F_ 2 (x_s)     =    11 √(7) + (10  α + 96  α^2 - 42  α^3 ) F (x_s)√(Ψ_ 2 (x_s)),   F_ 3 (x_s)  = √(7)+ (18 α^2 + 84  α^3 -  14  α^4 )  F (x_s)  √(Ψ_ 2 (x_s)), 
    
    F_ 4(x_s)     =    3 √(7) + 2 (1 + 22  α - 21 α^2) F (x_s)  √(Ψ_ 2 (x_s)),   F_ 5 (x_s) = 13 √(7) + (84  α^2 + 84  α^3 ) F (x_s)  √(Ψ_ 2 (x_s)),
    
    F_ 6 (x_s)    =    3 √(7) + 3 √(7) x_s + (2 α + 6  α^2   x_s  + 6  α^3   x^2 _s  + 2  α^4   x^3 _s ) F (x_s)  √(Ψ_ 2 (x_s))
    
    F_ 7 (x_s)     =     18 x_s  √(7) +21 √(7) + [28 α^3 (3 + α)  x^3_s + 28 α^4  x^4_s  + 28 α + 
    2 x_s (14  α + 42  α^2 )] F (x_s) √(Ψ_ 2 (x_s)).




§ THE EXPRESSIONS FOR USED COEFFICIENTS IN EQS.(<REF>)  AND (<REF>) ARE AS FOLLOWS:
 
Here ρ̃, ρ̃_c and  ρ̃_s are effective energy density, effective central density and effective surface density respectively. 



    ρ_1   =   [32 - 84 ρ̃_1   (1 + χ)/84 ρ̃_1 +
     8  2^2/3 [(8 + 63 ρ̃_1  χ)+ (ψ_ 1(ρ_ 1))^2/3]/84 ρ̃_1 (ψ_ 1 (ρ_ 1))^1/3],  ρ̃_1=72 ρ̃/7 (1+χ)^2 ρ̃_c
    ρ_1s   =   [32 - 84 ρ̃_1s   (1 + χ)/84 ρ̃_1s +
     8  2^2/3 [(8 + 63 ρ̃_1s χ)+ (ψ_ 1(ρ_ 1s))^2/3]/84 ρ̃_1s (ψ_ 1 (ρ_ 1s))^1/3],  ρ̃_1s=72 ρ̃_̃s̃/7 (1+χ)^2 ρ̃_c
    ψ_ 1 (ρ_ 1)    =    32+378 ρ̃_1 χ+(2646 +  1323  χ - 1323 χ^2) ρ̃^2_1 + √([32 + 378 ρ̃_ 1 χ +1323 ρ̃^2_ 1  (2 + χ - χ^2)]^2-2  (8 + 63  ρ̃_ 1 χ)^3)
    ψ_ 1 (ρ_ 1s)    =    32+378 ρ̃_1s χ+(2646 +  1323  χ - 1323 χ^2) ρ̃^2_1s + √([32 + 378 ρ̃_ 1s χ +1323 ρ̃^2_ 1s  (2 + χ - χ^2)]^2-2  (8 + 63  ρ̃_ 1s χ)^3)
    
    G(ρ_1s)   =   2 (1 + α ρ_1s)^3 √(Ψ_2(ρ_1s)) [3 + ρ_1s-α Ψ_7(ρ_1s) ]/[F_ 1 (ρ_1s) + ρ^2_1s  F_ 2 (ρ_1s) + ρ^3_1s  F_ 3 (ρ_1s) + F_11(ρ_1s)],   H(ρ_1s)=1/Ψ_1(ρ_1s) [1+G(ρ_1s)  F(ρ_1s)]√(7-10 ρ_1s-ρ^2_1s/7(1+χ+ρ_1s)^2) ,
    Ψ(ρ_1)   =   Ψ_1(ρ_1)[H(ρ_1s)+L(ρ_1s)  F(ρ_1)],   F(ρ_1)=√(7) √(Ψ_ 2 (ρ_1))[f_ 1 (ρ_1) + f_ 2 (ρ_1)]/3 Ψ_ 3 (ρ_1) Ψ_1(ρ_1) √(Ψ_1(ρ_1)) +√(7) f_3 (ρ_1) /(Ψ _ 3 (ρ_1))^7/2 ln[-a (Ψ_3 (ρ_1))^5/2  f_ 5 (ρ_1)/ 2   f_ 4 (ρ_1)  (1 +α  ρ_1)],
    
    F(ρ_1s)   =   √(7) √(Ψ_ 2 (ρ_1s))[f_ 1 (ρ_1s) + f_ 2 (ρ_1s)]/3 Ψ_ 3 (ρ_1s) Ψ_1(ρ_1s) √(Ψ_1(ρ_1s)) +√(7) f_3 (ρ_1s) /(Ψ _ 3 (ρ_1s))^7/2 ln[-a (Ψ_3 (ρ_1s))^5/2  f_ 5 (ρ_1s)/ 2  f_ 4 (ρ_1s)  (1 + α ρ_1s)],
    
    L(ρ_1s)   =   G(ρ_1s) H(ρ_1s),   
    F_11(ρ_1s)=ρ_1s  F_ 4 (ρ_1s) - 7  χ^2   F_ 6 (ρ_1s) - χ [F_ 7 (ρ_1s) + ρ^2_1s  F_ 5 (ρ_1s)],
    Ψ_ 1 (ρ_1)    =   (1 + α ρ_1)^2,    Ψ_ 2 (ρ_1) = [7 + 7 χ^2 - 10 ρ_1 - ρ_1^2 + 14 χ (1 + ρ_1)],  Ψ_ 3 (ρ_1) = [ 7 α^2 (1 + χ)^2 - 2 α (-5 + 7 χ)-1],
    Ψ_ 1 (ρ_1s)    =    (1 + α ρ_1s)^2,   Ψ_ 2 (ρ_1s) = [ 7 + 7χ^2 - 10 ρ_1s - ρ_1s^2 + 14χ (1 + ρ_1s)],   Ψ_ 3 (ρ_1s) = [7α^2 (1 + χ)^2 - 2α (-5 + 7χ)-1], 
    Ψ_ 4 (ρ_ 1)     =    √(7)  L(ρ_1s)  (1 + χ + ρ_ 1)/Ψ_ 1 (ρ_ 1) Ψ (ρ_1) √(Ψ_ 2 (ρ_ 1)) + 2 α/√(Ψ_ 1(ρ_ 1)),    Ψ_ 5 (ρ_ 1) = α^2  [5 - χ + ρ_ 1 + 2α  (-6 + 7 ρ_ 1 + ρ_ 1^2 - χ  (6 + 5 ρ_ 1))], 
    Ψ_ 6 (ρ_ 1)     =    α^2  [  7 + 7χ^3 - 15 ρ_ 1 - 2 ρ_ 1^2 +  21 χ^2  (1 + ρ_ 1) ],  Ψ_7(ρ_1s)=[7+7  χ^2-13  ρ_1s-2 ρ^2_1s+14 χ (1 + ρ_1s)],   
     
    f_ 1 (ρ_1)    =    (1 - α - α χ) Ψ^2_ 3 (ρ_1) + Ψ_ 3 (ρ_1) [-1 + 6 αχ + α^2 (-23 - 16 χ + 7 χ^2)] (1 + αρ_1),
    
    f_ 2 (ρ_1)    =   [1 + α (7 - 23 χ) + α^2 (63 + 174 χ - 105 χ^2) + α^3 (-359 + 201 χ + 399 χ^2 - 161 χ^3)] (1 + α ρ_1)^2,
    
    f_ 3 (ρ_1)    =    4 [-1 + 2 χ + α (-1 - 23 χ + 14 χ^2) + α^2 (37 + 36 χ - 147 χ^2 + 70 χ^3) + α^3 (-163 + 257 χ + 21 χ^2 - 301 χ^3 + 98 χ^4)],
    
    f_ 4 (ρ_1)    =   [-1 + 2 χ + α (-1 - 23 χ + 14 χ^2) + α^2 (37 + 36 χ - 147 χ^2 + 70 χ^3) +
        α^3 (-163 + 257 χ + 21 χ^2 - 301 χ^3 + 98 χ^4)],
  

    f_ 5 (ρ_1)    =    5 - 7 χ + ρ_1 + √(Ψ_ 3 (ρ_1)) √(Ψ_ 2 (ρ_1)) + α[7 + 7 χ^2 - 5 ρ_1 + 7 χ (2 + ρ_1)],
    
    f_ 1 (ρ_1s)     =     (1 - α - α χ)Ψ^2 _ 3 (ρ_1s) + Ψ_ 3(ρ_1s) [-1 +
        6 α χ + α^2  (-23 - 16χ + 7χ^2)] (1 + αρ_1s), 
    
    f_ 2 (ρ_1s)    =    [1 + α (7 - 23χ) + α^2 (63 + 174χ - 105χ^2) + α^3 (-359 + 201χ + 399χ^2 - 161χ^3)]  (1 + α ρ_1s)^2, 
    
    f_ 3 (ρ_1s)     =     4[-1 + 2χ + α (-1 - 23χ + 14χ^2) + α^2 (37 + 36χ -
           147 χ^2 + 70 χ^3) + α^3  (-163 + 257 χ + 21 χ^2 - 301 χ^3 +
            98 χ^4)],
    
    f_ 4 (ρ_1s)     =    [-1 + 2χ + α (-1 - 23χ + 14χ^2) + α^2 (37 + 36χ -
           147χ^2 + 70χ^3) + α^3 (-163 + 257χ + 21χ^2 - 301χ^3 +
            98χ^4)], 
    
    f_ 5 (ρ_1s)     =     5 - 7χ +  ρ_1s + √(Ψ_ 3 (ρ_1s)) √(Ψ_ 2 (ρ_1s)) + α[
     7 + 7χ^2 - 5 ρ_1s + 7χ (2 + ρ_1s)], 
    
     F_ 1 (ρ_1s)     =     - 7 √(7) - 7 √(7) χ^3 + (6 - 14  α + 2  α^3  (7 + 13  α) ρ^4_1s  + 4  α^4  ρ^5_1s  )  F (ρ_1s) √(Ψ_ 2 (ρ_1s));, 
    
    F_ 2 (ρ_1s)     =    11 √(7) + (10  α + 96  α^2 - 42  α^3 ) F (ρ_1s)√(Ψ_ 2 (ρ_1s)),   F_ 3 (ρ_1s)  = √(7)+ (18 α^2 + 84  α^3 -  14  α^4 )  F (ρ_1s)  √(Ψ_ 2 (ρ_1s)), 
    
    F_ 4(ρ_1s)     =    3 √(7) + 2 (1 + 22  α - 21 α^2) F (ρ_1s)  √(Ψ_ 2 (ρ_1s)),   F_ 5 (ρ_1s) = 13 √(7) + (84  α^2 + 84  α^3 ) F (ρ_1s)  √(Ψ_ 2 (ρ_1s)),
    
    F_ 6 (ρ_1s)    =    3 √(7) + 3 √(7)ρ_1s + (2 α + 6  α^2 ρ_1s  + 6  α^3 ρ^2_1s  + 2  α^4 ρ^3_1s ) F (ρ_1s)  √(Ψ_ 2 (ρ_1s))
    
    F_ 7 (ρ_1s)     =     18 ρ_1s√(7) +21 √(7) + [28 α^3 (3 + α) ρ^3_1s + 28 α^4  ρ^4_1s  + 28 α +
    2 ρ_1s (14  α + 42  α^2 )] F (ρ_1s) √(Ψ_ 2 (ρ_1s)).
  





