Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Anisotropic spheres via embedding approach in $\mathcal{R}+\beta\mathcal{R}^{2}$ gravity with matter coupling

Anisotropic spheres via embedding approach in $\mathcal{R}+\beta\mathcal{R}^{2}$ gravity with... Anisotropic spheres via embedding approach in R + βR gravity with matter coupling 1, ∗ 2, † 3, ‡ 4, § G. Mustafa, Xia Tie-Cheng, Mushtaq Ahmad, and M. Farasat Shamir Department of Mathematics, Shanghai University, Shanghai, 200444, Shanghai, People’s Republic of China Department of Mathematics, Shanghai University, Shanghai, 200444, Shanghai, People’s Republic of China. National University of Computer and Emerging Sciences, Chiniot-Faisalabad Campus, Pakistan. National University of Computer and Emerging Sciences, Lahore Campus, Pakistan. The manifesto of the current article is to investigate the compact anisotropic matter profiles in the context of one of the modified gravitational theories, known as f(R, T ) gravity, where R is a Ricci Scalar and T is the trace of the energy-momentum tensor. To achieve the desired goal, we capitalized on the spherical symmetric space–time and utilized the embedding class-1 solution via Karmarkar’s condition in modeling the matter profiles. To calculate the unidentified constraints, Schwarzschild exterior solution along with experimental statistics of three different stars LMC X-4, Cen X-3, and EXO 1785-248 are taken under consideration. For the evaluation of the dynamical equations, a unique model f(R, T ) = R + βR + λT has been considered, with β and λ being the real constants. Different physical aspects have been exploited with the help of modified dynamical equations. Conclusively, all the stars under observations are realistic, stable, and are free from all singularities. Keywords: Anisotropic spheres; f(R, T ) gravity; Compact stars; Embedding Class I. I. INTRODUCTION Late time evolution of stellar configurations, triggered by an immense gravitational pull has been anticipated to a great extent in the field of astrophysics and the modified gravitational theories. It expedites the examination of diverse attributes regarding the gravitating source by physical phenomena. Baade and Zwicky [1] forecast the inception of highly dense stellar objects inaugurating the debate that a supernova might be revolutionized into a highly dense star. This reality came into existence when exceptionally magnetized as well as rotating neutrons stars were detected. Therefore, a fundamental shift regarding normal stars to compact stars came into existence. By the newly discovered concept, the normal stars shifted into an extensive range, such as quark stars, neutron stars, gravastars, dark stars, and finally black holes. The actuality of the extensive range of these stars led the researchers to curiosity, regarding the formation of these stars. The stellar death of a normal star occurs, that is when the nuclear fusion reactions cease to act and burn all of their nuclear fuel results in the formation of new compact stars. The newly formed compact stars are primarily distinguished from the normal stars in two ways. Since all the fuel has been utilized by the star, hence the star cannot sustain against the gravitational collapse due to thermal pressure. Analogous to that, the white dwarf is stabilized due to strong degenerate electron pressure, while the neutron star is stabilized due to degenerate neutron pressure. Whereas black holes are entirely the collapsed remnants, therefore there is neither a thermal pressure nor a degenerate pressure sufficient enough to repress the centripetal pull of gravity; as a result, it leads towards the gravitational singularities and the event horizon. The formed compact stellar remnants consist of huge density and relatively small radii in contrast to the normal stars. The intention to investigate the physically stable models, leads us to an analytical approach regarding the Einstein field equations. One of the essential tools is to adopt the embedding class I space-time which transforms a four dimensional manifold into a Euclidean space of higher dimension. The conversion of curved embedding class space–time into higher dimensional space–time is substantial to develop exact new models in the field of astrophysics. The class I embedding condition leads towards a differential equation in the framework of spherically symmetric Electronic address: gmustafa3828@gmail.com Electronic address: xiatc@shu.edu.cn Electronic address: mushtaq.sial@nu.edu.pk Electronic address: farasat.shamir@nu.edu.pk; farasat.shamir@gmail.com arXiv:2101.00208v1 [gr-qc] 1 Jan 2021 2 space–time which connects the gravitational potentials i.e., g and g , the condition is also recognized as the rr tt Karmarkar condition [2]. The Karmarkar’s condition appears to be very influential in exploring new solutions for the astrophysical models. Schlia [3] was the pioneer in developing the Karmarkar condition for a spherically symmetric space–time. The embedding theorem based on the isometrics has been presented by Nash [4]. Maurya et al. [5]-[12] were the first explorers in the aspect of applying the embedding approach to the anisotropic matter configurations. After the new dawn of general relativity (GR), f(R) theory is considered to be quite a fascinating tool for the amplification of GR. Further, many researchers presented different versions of this theory, which were also very prosperous in diverse fields. The recent extension of this theory is regarded as f(R,T ) gravity, which was presented by Harko et al. [13]. The f(R,T ) theory has been the center of attention by many analysts and consequently many intriguing cosmological aspects have been unraveled [14]-[17]. The analysis regarding isotropic matter profile of the self-gravitating system and its stability has been done by Sharif et al. [18]. Alhamzawi and Alhamzawi [19] construed the occurrence of lensing of gravitation in the context of a modified f(R,T ) theory. Moraes et al. [20] numerically investigated the stability of the gravitational lensing by utilizing the Tolman–Oppenheimer–Volkov (TOV) equations in f(R,T ) gravity. Das et al. [21] formulated a family of solutions by characterizing the interior geometry of compact stars, permitting conformal motion under the influence of f(R,T ) gravity. Moraes et al. [22] investigated the configurations consisting of hydrostatic equilibrium along with fluids whose pressure was computed from equation of state (EoS) in the light of f(R,T ) gravitational theory. Yousaf et al. [23] investigated the formation of relativistic stellar profiles in the regime of f(R,T ) gravity by utilizing the Krori and Barura model. The study of dense anisotropic profiles consisting of charge has been investigated by Maurya and Aurtiz [24]. In this regard, they utilized the Durgapal–Fuloria model in the context of f(R,T ) gravity and applied gravitational decoupling utilizing geometric deformation. Waheed et al. [25] analyzed the existence of highly dense stellar configurations by utilizing Karmarkar along with the Pandey–Sharma condition. To do so, they used spherically symmetric space–time in the context of f(R,T ) gravity. Mustafa et al. [26] analyzed the Class 1 embedding condition in the presence of anisotropy matter profile and utilized the interior geometry of Schwarzschild along with Kohler–Chao solutions in modified gravity. The matter configuration consisting of nuclear density of 10 gm/cc exhibits the behavior of anisotropy i.e. p − p 6= 0 which exists due to certain factors involving magnetic flux, viscosity, phase transition, etc. t r In this regard, Ruderman [27] is the pioneer who argued about the anisotropy existing at the interior of the stars. Modified gravitational theories have provided an overwhelming approach in analyzing the anisotropic stellar con- figurations inheriting high matter profiles [28]-[33]. This work aims to analyze the modified f(R,T ) gravity to devise a realistic configuration which in nature is anisotropic. For this purpose, we take into account three different matter profiles i.e. LMC X-4, Cen X - 3, and EXO 1785–248, and apply a well-known embedding class 1 approach. Moreover, the structural aspect of anisotropic profiles has been examined by making use of spherically symmetric space–time along with the categorical f(R,T ) gravity model. The layout of this article is as follows: In Section 2, modified field equations have been formulated by utilizing the Karmarkar condition. Section 3 is to provide the matching condi- tions by considering Schwarzschild’s solution. The physical investigation has been done comprehensively in Section 4. Conclusive remarks have been provided in the last Section. II. f(R, T ) THEORY OF GRAVITY The modified form of Einstein-Hilbert action for extended f(R,T ) theory of gravity is defined as follows: S = [L + f(R,T )] −gd x, (1) where L and f(R,T ) denote the matter Lagrangian density and a modified function, respectively. Here, R and T are known as scalar curvature and trace of the energy-momentum tensor, respectively. Now, by varying Eq.(1), we get the following modified set of equations (1 − f (R,T ))T − f (R,T )Θ = − f(R,T )g + (R + (g ✷ − ∇ ∇ ))f (R,T ), (2) T νµ T νµ νµ νµ νµ µ ν R where, αβ 2 g δT ∂ L νµ m αβ Θ = = −2g − 2T + g L , νµ νµ νµ m νµ νµ αβ δg ∂g ∂g νµ ∂ ( −gg ∂ ) ∂f(R,T ) ∂f(R,T ) µ ν ✷ = √ , f (R,T ) = , f (R,T ) = , R T −g ∂R ∂T 3 with ∇, representing the covariant derivative. The energy-momentum tensor with the anisotropic matter source is defined as T = ρU U + p V V + p (U U − g − V V ), (3) νµ µ ν r µ ν t µ ν νµ µ ν where U represents the vector for 4-velocity and V is a vector in the direction of radial pressure. Further, the µ µ expressions, i.e., ρ, p and p are used to define define energy density, tangential and radial components of pressure, t r respectively. Using Eq.(3) in Eq.(2), we get the following set of equation: 1 1 G = (1 + f (R,T ))T + (∇ ∇ − g ✷)f (R,T ) + (f(R,T ) − Rf (R,T ))g T νµ µ ν νµ R R νµ νµ f (R,T ) 2 − ρg f (R,T )) . (4) νµ T We assume a static and spherically symmetric line element, which is defined as: 2 a(r) 2 b(r) 2 2 2 ds = −e dt + e dr + r dΩ , (5) 2 2 2 2 2 b(r) a(r) where the expression dΩ defines the g = r dθ and g = r sin θdφ components, e and e denote the θθ φφ gravitational components of stellar geometry. Further, we fix a quadratic model of f(R,T ), which is defined as: f(R,T ) = R + β × R + λ × T . (6) The considered model f(R,T ) = R + βR + λT involves a particular case of well-known Starobinsky model [35] with matter coupling. It is an interesting point that, in the Starobinsky model R + βR , a maximum value of M/M or beyond is reached when the value of the parameter β is selected to be negative. But, this leads to an issue; specifically, the Ricci scalar performs a damped oscillation. On the other hand, the Ricci scalar smoothly decreases to zero as we approach towards infinity for positive values of parameter β, for which the star can support a maximum mass lower than 2M/M . Now, we elaborate an eminent Karmarkar condition concisely which is the integral tool for current study. The infrastructure connecting the Karmarkar condition is established on the class I space of Riemannian geometry. A sufficient condition comprises of a second order symmetric tensor and the Riemann Christoffel tensor, given as Σ(Λ Λ − Λ Λ ) = R , ηµ υγ γµ νη υµ ηγ Λ ; n − Λ ; ν = 0. νµ νη Here ; stands for covariant derivative whereas Σ = ±1. These values signify a space-like or time-like manifold relying on the sign considered as − or +. Now, the Karmarkar condition is defined as R R = R R + R R . (7) 1414 2323 1224 1334 1212 3434 These Riemann tensor components are given below as follows. a(r) ′′ ′ 2 ′ ′ 2 2 b(r) e (2a (r) + a (r) − a (r)b (r)) r sin θ(e − 1) R = , R = , 1414 2323 b(r) 4 e ′ 2 ′ a(r) rb (r) rsin θb (r)e − b(r) R = , R = , 1212 3434 2 2 R = R sin θ, R = 0, 1334 1224 1224 where, R 6= 0. A differential equation can be achieved by utilizing the Karmarkar condition using Eq. (7) as ′ ′ a (r)b (r) ′ ′ ′ 2 ′′ ′ 2 b(r) − a (r)b (r) + a (r) − 2 a (r) + a (r) = 0, e 6= 1. (8) b(r) 1 − e Integration of Eq. (8) provides a connection between two main gravitational components of the space-time as follows b(r) a(r) ′ 2 e = e × a (r) + 1 + K, (9) 4 where K is a constant of integration. We choose a specific model for a g component which is expressed as tt a(r) 2 e = ψ r ψ + 1 , (10) 1 2 where ψ , ψ are assumed as constants, n is an integer. By plugging Eq. (10) in Eq. (9), we get the g component, 1 2 rr which is calculated as b(r) 2 2 n−2 e = r ψ ψ r ψ + 1 + 1, (11) 2 3 2 where ψ = 4n × ψ × ψ × K. It is mentioned here that we get realistic results for n > 2. Now, we are able to 3 2 1 calculate the following set of modified field equations for the anisotropic stellar configuration. 2 2 2 1 ψ 12β(λ + 1)n r ψ (−Υ ) 2 2 1 ρ = × − + Υ 2 3−n 4 n−2 2 3 Υ Υ + r ψ ψ Υ 2 3 1 2(λ + 1)(2λ + 1) r ψ ψ Υ + r 1 1 2 3 4 6 3 n+2 4β(3λ + 1)n r ψ 8β(3λ + 2) ψ Υ Υ 3 16 2 2 4 2 2 1 2 + − × − + ψ r 4nΥ r ψ − 3 − 3n r ψ − 6r ψ + 1 2 1 2 2 4 4 2 Υ Υ Υ 1 1 5 n+1 2 2ψ Υ 2nr (−Υ ) Υ 3 1 10 n+2 1 2 2 2n 4 − × r ψ ψ (Υ Υ ) Υ − ψ Υ Υ + (n − 2)Υ Υ − − Υ − Υ + Υ 2 1 8 3 7 6 9 11 12 3 1 1 1 3 2 Υ Υ 5 1 2 2 4 n r ψ Υ n 2r ψ ψ Υ 2 13 3 2 15 n−2 2 + + 2ψ Υ 2β(9λ + 5) + Υ (2λ + 1)r − 2(3βλ + β) − × λΥ 3 14 14 2 1 3−n Υ Υ 1 Υ + r ψ ψ Υ 1 2 3 1 2 2 − 4 6β + r + β 24λ + 2r ψ + Υ − Υ , (12) 2 17 18 3−n Υ + r ψ ψ Υ 2 3 1 −1 2ψ Υ Υ 4r ψ Υ 3 15 24 2 23 p = − 4βrψ + + Υ + Υ − 2Υ Υ + 4rψ r 2 21 22 14 20 3 2 3−n 4 4(λ + 1)(2λ + 1)r Υ Υ Υ + r ψ ψ Υ 14 2 3 1 1 2 3 2 2n r ψ Υ 2nrψ 25 2 n−2 2 2 2 × ψ Υ + + Υ (2λ + 1)r − 2(3βλ + β) − 2β(15λ + 7) β 8(5λ + 4) + 2r ψ 14 2 1 2 Υ Υ n−2 8(λ + 2)ψ Υ Υ 4(5λ + 6)n (−Υ ) ψ Υ Υ 3 16 1 3 15 28 × + + + Υ − Υ − Υ 26 27 14 2 2 3−n Υ Υ Υ + r ψ ψ Υ 5 1 2 3 1 2λr ψ ψ Υ 2 3 15 × + 4 (λ + 1)r − 6βλ , (13) 3−n Υ + r ψ ψ Υ 2 3 1 2 2 ψ 4βΥ n r ψ Υ 2ψ Υ Υ 1 2 31 2 33 3 15 32 p = + + − × n (β (−8λ 2 4 2 3−n n−2 3 Υ Υ Υ Υ + r ψ ψ Υ 2 3 1 1 2(λ + 1)(2λ + 1) r ψ ψ Υ + r 1 1 2 3 1 n−2 2 2 2 2 2 4(8λ + 7)ψ Υ Υ 4(3λ + 5)n (−Υ ) 4(7λ + 6)nr ψ r ψ − 3 4(24λ + 19)r ψ ψ Υ 2 2 3 16 1 2 2 1 3 15 + 2r ψ + − − 2 2 3 3−n 2 2 Υ Υ Υ Υ + r ψ ψ Υ 5 1 1 2 3 1 ! ! ! 4 2 3 3 Υ 12(3λ + 2)r ψ ψ Υ Υ 36 35 2 3 15 + + + + 8 + Υ Υ + Υ − Υ , (14) 14 34 29 30 3−n 3−n 3 2 3 Υ + r ψ ψ Υ Υ + r ψ ψ Υ Υ 2 3 1 2 3 1 5 1 1 where Υ , {i = 1, ..., 36}, are given in the Appendix (I). III. COMPARISON OF EXTERIOR AND INTERIOR SOLUTION By considering the Jebsen-Birkhoff’s theorem, the spherically symmetric vacuum solution of GR field equations must be asymptotically flat. In particular, the spacetime is of the form 2 2 −1 2 2 2 2 2 ds = εdt − ε dr − r (dθ + sin θdφ ), (15) 2M where ε = 1 − . Here, M denotes the stellar mass of the star. Now, considering the constraint p (r = R ) = 0 at r ǫ the boundary r = R and using metric coefficients g and g from Eq. (5) and Eq. (15), we calculate the following ǫ tt rr 5 TABLE I: Predicted values of ψ , ψ , ψ , and λ with β = 2 and and (Radii=9.711 km, & Mass =1.29M/M ). 1 2 3 ⊙ LMC X-4 n ψ ψ ψ λ 1 2 3 3 0.4329 0.001273 4.7820 0.066676 5 0.4363 0.000729 7.6640 0.040511 10 0.4387 0.000352 14.900 0.021340 20 0.4399 0.000173 29.392 0.011919 50 0.4406 0.000068 72.882 0.006322 100 0.4408 0.000034 145.373 0.004466 −6 500 0.4410 6.827742×10 725.285 0.002985 TABLE II: Predicted values of ψ , ψ , ψ , and λ with β = 2 and (Radii= 10.136 km, & Mass =1.49M/M ). 1 2 3 ⊙ Cen X - 3 n ψ ψ ψ λ 1 2 3 3 0.3765 0.001421 4.5683 0.055841 5 0.3807 0.000805 7.2751 0.032774 10 0.3838 0.000386 14.082 0.015565 20 0.3852 0.000186 28.267 0.007027 50 0.3861 0.000074 68.656 0.007027 100 0.3864 0.000036 138.251 0.000239 −6 500 0.3866 7.447003×10 682.743 0.000113 expressions 2M ψ ψ R + 1 = 1 − , (16) 1 2 −1 n−2 2M 2 2 ψ ψ R ψ R + 1 + 1 = 1 − , (17) 3 2 2 ǫ ǫ n−1 M nψ ψ ψ R + 1 = , (18) 1 2 2 p (r = R ) = 0. (19) r ǫ Utilizing these boundary conditions from Eqs. (16-19), we get the following relations −n (R − 2M) 1 − 2nM−nR +M ψ = , (20) ψ = , (21) R (nR − (2n + 1)M) 1−n ψ = 2n 1 − , (22) 2nM − nR + M λ + λ + λ + λ 1 2 3 4 λ = , (23) λ + λ + λ + λ 5 6 7 8 where λ , {i = 1, ..., 8} are given in the Appendix (II). The estimated values of the above parameters, i.e., ψ , ψ , ψ , & λ are given in Table I, Table II and Table III. 1 2 3 6 TABLE III: Predicted values of ψ , ψ , ψ , and λ with β = 2 and (Radii=8.849 km, & Mass =1.30M/M ). 1 2 3 ⊙ EXO 1785-248 n ψ ψ ψ λ 1 2 3 3 0.3768 0.001862 4.5698 0.076448 5 0.3811 0.001055 7.2777 0.044225 10 0.3842 0.000496 14.361 0.020799 20 0.3856 0.000248 27.733 0.009361 50 0.3864 0.000098 68.685 0.002592 100 0.3867 0.000048 136.944 0.000352 −6 500 0.3869 9.759977×10 683.029 0.000134 IV. PHYSICAL ANALYSIS In this section, we briefly present the results by analyzing the physical attributes along with different aspects of the stellar configurations under the acquired f(R,T ) gravity model. In order to fulfill the purpose, experimental data of distinct stars i.e. LMC X-4, Cen X – 3 and EXO 1785–248 is used. All the attributes of the stellar configurations are depicted in tabular form as well as graphically. A. Gravitational Metric Potential The existence of anomalies within the sphere such as geometric singularities are contemplated to be an essential peculiarity in the investigation of stellar spheres. In order to unravel the existence of singularities, we examine the a b nature of gravitational potential g = e and g = e at the core r = 0 of the sphere. Physical essence and endurance tt rr of the models rely upon the gravitational metric potentials and it should be decreasing on regular intervals within the spherical structures. It can be observed from Fig. 1 that the metric potential with in the interior of the sphere b(r=0) a(r=0) exhibits the behavior of e = 1 and e 6= 0 which is consistent and physically valid. It can also be observed that both of the metric potentials exhibit minimum values at the center and show non-linear increasing behavior towards the boundary. FIG. 1: Visual representation of gravitational potentials with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) B. Energy Density and Pressure Evolutions Prior to analysis of the anisotropy, we investigate the evolutional change of the matter profiles in connection to the energy density ρ along with anisotropic stresses such as p and p . The energy density along with pr and pt exhibits r t the exceptional behavior of high density of matter configuration. The phenomenal high density is due to the strong forces of attraction which are regarded as dipole interactions and intermolecular forces. Numerical values of density and the components of the pressure for the three compact spheres are provided in Tables IV-VI. All of the physical attributes remain positive and appear to be finite at the core. It confirms that the current system is independent 7 of all singularities. From the Figs. 2-4, it is evident that the matter configuration under consideration attains the maximum mass at the core and tends to zero at the boundary of the star, which depicts the high compactness of the stellar spheres. These graphical plots establish the presence of anisotropy of the compact sphere under the influence of our f(R,T ) model. FIG. 2: Visual representation of ρ with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) FIG. 3: Visual representation of p with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) FIG. 4: Visual representation of p with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) t 8 C. Anisotropy and Gradients In order to model the interior geometry of the relativistic stellar configuration under current circumstances, the role of anisotropy is crucial for the compact sphere modeling and it is represented as Δ = p − p . (24) t r It depicts the information regarding the anisotropic nature of the stellar configuration. If p > p then the anisotropy t r is considered to be non negative and is drawn outwards and depicted as Δ < 0. Whereas if p > p then the anisotropy r t turns out to be negative and this shows that anisotropy is drawn inwards. From 5 it is observed that for our ongoing study anisotropy remains positive, hence, directed outwards. The deviation of radial derivatives of the energy density dρ dp dp r t and pressure components , i.e., , and , are shown in Figs. 6-8 such that dr dr dr FIG. 5: Visual representation of Δ with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) dρ dp dp r t < 0, < 0, < 0. (25) dr dr dr It can be observed from the second order derivatives that the pressure components and energy density show the maximum value at the core r = 0, i. e., 2 2 2 d ρ d p d p r t > 0, > 0, > 0. (26) 2 2 2 dr dr dr FIG. 6: Visual representation of gradient of ρ with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) 9 FIG. 7: Visual representation of gradient of p with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) FIG. 8: Visual representation of gradient of p with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) D. Energy Conditions Energy conditions appear to be quite helpful in analyzing the realistic distribution of matter. These attributes play a decisive role to classify the exotic and normal mater distribution within the stellar model. The energy conditions have been crucially important in debating the issues related to cosmology and astrophysics. The energy condition are classified as NEC : ρ > 0, ρ + p ≥ 0, ρ + p ≥ 0, DEC : ρ ≥ |p |, ρ ≥ |p |, r t r t WEC : ρ − p ≥ 0, ρ − p ≥ 0, r t SEC : ρ − p ≥ 0, ρ − p ≥ 0, ρ − p − 2p ≥ 0. (27) r t r t TABLE IV: Predicted values of physical parameters at center and boundary. LMC X-4 a(r=0) b(r=0) 3 2 3 n e e ρ (g/cm ) p = p (dyne/cm ) ρ (g/cm ) p /ρ = p /ρ R r t 0 r 0 t 0 0 0 0 0 14 34 14 3 0.431 1.0 0.7546 × 10 13.5220 × 10 1.8378 × 10 0.074 14 34 14 5 0.432 1.0 0.8479 × 10 14.1351 × 10 1.7467 × 10 0.081 14 34 14 10 0.433 1.0 0.9735 × 10 15.9780 × 10 1.6629 × 10 0.094 14 34 14 20 0.434 1.0 1.0004 × 10 16.1932 × 10 1.5492 × 10 0.110 14 34 14 50 0.440 1.0 1.0875 × 10 17.1260 × 10 1.4207 × 10 0.124 14 34 14 100 0.441 1.0 1.1787 × 10 18.1910 × 10 1.3328 × 10 0.135 14 34 14 500 0.442 1.0 1.2433 × 10 18.9280 × 10 1.2558 × 10 0.136 10 TABLE V: Predicted values of physical parameters at center and boundary. Cen X - 3 (mass =1.49M/M & radii= 10.136 km) a(r=0) b(r=0) 3 2 3 n e e ρ (g/cm ) p = p (dyne/cm ) ρ (g/cm ) p /ρ = p /ρ R r t 0 r 0 t 0 0 0 0 0 14 34 14 3 0.376 1.0 1.7236 × 10 16.4320 × 10 1.9578 × 10 0.095 14 34 14 5 0.377 1.0 1.8429 × 10 17.1451 × 10 1.8769 × 10 0.103 14 34 14 10 0.378 1.0 1.9325 × 10 18.9380 × 10 1.7629 × 10 0.125 14 34 14 20 0.379 1.0 2.0137 × 10 19.1632 × 10 1.5092 × 10 0.140 14 34 14 50 0.380 1.0 2.0945 × 10 20.8460 × 10 1.4007 × 10 0.164 14 34 14 100 0.381 1.0 2.3527 × 10 21.1220 × 10 1.4428 × 10 0.169 14 34 14 500 0.382 1.0 2.5723 × 10 21.9180 × 10 1.4958 × 10 0.172 TABLE VI: Predicted values of physical parameters at center and boundary. EXO 1785-248 (Mass =1.30M/M & Radii=8.849 km) a(r=0) b(r=0) 3 2 3 n e e ρ (g/cm ) p = p (dyne/cm ) ρ (g/cm ) p /ρ = p /ρ R r t 0 r 0 t 0 0 0 0 0 14 34 14 3 0.376 1.0 0.8876 × 10 15.4634 × 10 1.9684 × 10 0.082 14 34 14 5 0.381 1.0 0.9574 × 10 16.7635 × 10 1.8073 × 10 0.099 14 34 14 10 0.384 1.0 1.0165 × 10 17.3650 × 10 1.7871 × 10 0.114 14 34 14 20 0.385 1.0 1.0969 × 10 18.7532 × 10 1.6063 × 10 0.131 14 34 14 50 0.386 1.0 1.1607 × 10 19.8591 × 10 1.4534 × 10 0.146 14 34 14 100 0.387 1.0 1.2485 × 10 20.1813 × 10 1.3115 × 10 0.153 14 34 14 500 0.388 1.0 1.4524 × 10 20.9994 × 10 1.3089 × 10 0.156 Here, NEC stands for null energy condition, SEC for strong energy condition , DEC for dominant energy condition and WEC for week energy condition. It can be seen from Fig. 9 that all the energy bounds exhibit decreasing behavior with the increase in radii of the compact stellar sphere. FIG. 9: Evolution of energy bounds with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) E. Analysis of Stability The stability of the stellar configuration plays a decisive role in analyzing the consistency of the acquired model. Many analytical discussions have been done in order to find the stability of the matter configuration but Herrera’s cracking conception emerged to be very effective [34]. The radial and tangential speed of sound are defined as dp dp r t 2 2 v = , and v = . (28) sr st dρ dρ For the conservation of causality condition, components of the speed of sound must be with the bounds of the interval 2 2 2 2 [0, 1] i.e 0 ≤ v and v ≤ 1. It can be observed from the Figs. 10-11 that the condition i.e 0 ≤ v and v ≤ 1 is r t r t 2 2 satisfied by both velocity components. Apart from that Abreu condition i.e −1 ≤ |v − v | ≤ 0 has also been satisfied t r 11 and can be observed from the Fig. 11. The validity of both of the aspects affirm the viability and the effectiveness of 2 2 our model. Further, it can also be observed that inverse Abreu condition i.e 0 ≤ |v − v | ≤ 1 is also satisfied. t r 2 2 FIG. 10: Visual representation of v and v with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and r t n = 500(⋆) 2 2 2 2 FIG. 11: Visual representation of |v − v | and |v − v | with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), t r r t n = 100(⋆), and n = 500(⋆) F. Equilibrium Analysis for Modified f(R, T ) Gravity In this section, we will analyze the equilibrium condition by considering the stability of the acquired solution of the three different stellar configuration. For the purpose, we make use of the TOV equation [36]-[41] dp dp dρ r t ′ λ − − 2 + 3 2Δ dp a dr dr dr − − (ρ + p ) + = 0. (29) r dr 2 3(2λ + 1) The above equation characterizes the necessary and sufficient condition for the hydrostatic-equilibrium. It comprises of four different forces dp dp dρ r t ′ λ − − 2 + 3 dr dr dr a dp 2Δ F = (ρ + p ), F = , F = , F = . (30) g r h a e 2 dr r 3(2λ + 1) • F represents the anisotropy force. • F represents the hydrostatic force. • F represents the gravitational force. g 12 • F represents the extra force. Consequently, the TOV equation can also be written as F + F + F + F = 0. From the attained graph as shown g h a e in Fig. 12, it is deduced that all of the forces sum up to neutralize the total effect, and this confirms the existence of the stable stellar structures. FIG. 12: Visual representation of F , F , F , and F with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), g h a e n = 100(⋆), and n = 500(⋆) G. Evolution of Adiabatic Index For the energy density, EoS stiffness can be better described by the adiabatic index. The stability of relativistic as well as the non-relativistic stellar structures can be explained through the adiabatic index. The concept of the dynamical stability via the radial adiabatic index was presented by Chandrasekhar [42]. This was further utilized by many authors [43]-[48]. For the system to be dynamically stable, the adiabatic index must go beyond 4/3. Adiabatic index corresponding to the radial stress is given as ρ + p dp ρ + p r r r Γ = ( ) = v . (31) p dρ p r r One of the quite fascinating fact of the above equation is that the stability of the Newtonian matter configuration is achieved when Γ > 4/3. While if Γ = 4/3 then a neutral equilibrium is achieved whereas if Γ < 4/3 then an r r r unstable matter configuration consisting of anisotropy is achieved. Anisotropic matter profile via adiabatic index can be elaborated as 4 ρ p 4 p − p i ri ti ri Γ = + ( r + ). (32) ′ ′ 3 2 | p | 3 | p | r ri ri From the Fig. 13, the graphical behavior of Γ with respect to increasing radii can be observed. It is noted that Γ r r shows the monotonically increasing conduct for all the stellar spheres and Γ is always greater than 4/3. Hence, Γ r r is consistent for the stability of our model in f(R,T ) gravity. H. Equation of state The evolution of the emergence of compact stars can be determined by EoS of matter. Moreover, EoS has a strong impact on the conditions of nucleosynthesis. Therefore, EoS is a vital tool in many astrophysical simulations. The EoS is considered to be a ratio of the pressure terms p and p with density. The components of EoS for the study r t of stellar configuration are i.e., ω and ω and are mathematically connected as r t p p r t ω = , ω = . (33) r t ρ ρ From Fig. 14, it can be observed that with the increase in radii, the components of the EoS shows monotonically decreasing behavior and are always less than 1. Moreover, the positive nature is observed for both of the components of EoS i.e. ω and ω with in the matter configuration. The accomplishment of the condition i.e. 0 ≤ ω and ω < 1 r t r t unveils that the our obtained solutions are valid and legitimate. 13 FIG. 13: Evolution of Γ with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) FIG. 14: Evolution of energy EoS with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) I. Compactness factor and Surface Redshift For the existence of any matter configuration mass function, compactness factor along with surface redshift function are contemplated to be an essential constituent. The fundamental relation for the mass function is given as m(r) = 4πρr dr. (34) FIG. 15: Evolution of m(r) with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) 14 FIG. 16: Evolution of μ(r) with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) FIG. 17: Evolution of z with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) The generalize compactness factor i.e. μ(r) is represented as m 1 μ(r) = = 4πρr dr. (35) r r The strong intermolecular interaction forces with in the stellar matter configuration and its corresponding EoS can be characterized by the term surface redshift i.e., Z . The generalized relation is Z = − 1. (36) 2 n 1/2 (ψ (ψ r + 1) ) 1 2 Figs. 15-17 represent the evolution of redshift function along with compactness factor corresponding to the increasing radii. It can be observed that the surface redshift is always less than 5 i.e., 5 i.e. Z ≤ 5 and the compactness factor remains less than 0.30 i.e., μ(r) ≤ 0.30. All of the mentioned functions are positive throughout the configuration. Hence, our models are stable. V. CONCLUSION The manifesto of the current study is to identify the realistic and stable configuration for the stellar sphere in the modified f(R,T ) gravitational theory. For the analysis of the stellar matter configuration, a viable model f(R,T ) = R + βR + λT is considered along with spherically symmetric space–time. In order to achieve the current objective, observational data of three compact spheres i.e. LMC X-4, Cen X-3 and EXO 1785–248 inheriting an anisotropic matter distribution has been utilized. This realistic range for masses of stars under this study is 1.29 to 1.5 solar mass. In this context, the considered models LMC X-4 (mass 1.29 M/M ), Cen X-3 (mass 1.49 M/M ) and EXO 1785-248 ⊙ ⊙ (mass 1.29 M/M ) are with in suitable range under this study. Moreover, the radii are also with in prescribed range 8.849 km to 10.136 km. In general, the study is valid for other models of stars and those within the given ranges of 15 mass and radius under this study. Embedding class 1 condition is used to find the potential i.e. g by considering rr a(r) 2 the primary potential as i.e g = e = χ r χ + 1 . To find the unknown constraints, Schwarzschild’s exterior tt 1 2 solution has been utilized. All the obtained results can be summarized as: • Metric potentials: The presence of singularities within the stellar configuration is an essential topic, worthy of debate. The stability of the stellar matter configuration depends on it. Therefore, the aspects of metric potential b(r=0) play a decisive role. The graphical behavior of Fig. 1 depicts that the fundamental condition i.e. e = 1 a(r=0) and e 6= 0 has been encompassed by the gravitational potential. Increasing attribute of the potentials is also observed throughout the configuration. Therefore, the potentials are free from any singularity and so is our model. • Energy density and stress constraints: From the Figs. 2-4 the evolution of density and components of stresses i.e., p and p can be observed. The density along with stress components show decreasing behavior and are non- r t negative throughout the configuration. The peak value is accomplished at the core while decreasing evolution is observed with the increase in radii and it tends to 0 towards the boundary. • Anisotropy and gradients: Fig. 5 depicts anisotropy behavior for our current configuration. It is noted that p 6= p and p > p , therefore, △ > 0 so the anisotropy is positive and is directed outwards. From Fig. 6-8 t r t r gradient of density and stress components are reviewed and it is noticed that all the gradients are negative and dρ dp dp r t exhibit decreasing behavior i.e., < 0, < 0, < 0. Since the non-positive behavior of all gradients along dr dr dr with their vanishing attribute at r = 0 is observed, therefore, our stellar configuration is stable. • Energy bounds: From the Figs. 2-4 and Fig. 9 the characteristic behavior of functions ρ, p , p , ρ − p , ρ − p , r t r t ρ−p − 2p are depicted. Decreasing characteristic behavior is observed with the increase in radii which further r t fulfills the bounds of energy such as (NEC), (SEC), (DEC) and (WEC). Hence the matter profile is realistic and viable. • Causality analysis The behavior of the constraints of sound speed is depicted in Figs. 10 and Fig. 11. The decreasing attribute is observed for both the components of the speed and it is shown that they are always 2 2 within the limits i.e., 0 ≤ v ≤ 1 and 0 ≤ v ≤ 1. The fulfillment of Aberu condition is also observed i.e., r t 2 2 2 2 −1 ≤ |v − v | ≤ 0. For the current model v > v . Therefore, fulfillment of all the conditions confirms the t r t r viability of the stellar sphere. • Equilibrium and EoS analysis: The balancing nature of all the forces, i.e., F , F and F is depicted in Fig. a h g 12. As all these forces add up to 0 and balance the effect of each other, therefore, the equilibrium condition is satisfied. From Fig. 14, the attributes of the parameters of EoS are observed. It is concluded that constraints i.e., ω and ω of EoS are positive in the interior of stellar profile and are in the stability bounds of 0 ≤ ω and r t r ω < 1. • Adiabatic index stability analysis: The behavior of adiabatic index Γ can be seen in Fig. 13 showing that Γ > 4/3. It also depicts the positive and decreasing nature, justifying the effectiveness of our system in the framework of f(R,T ) theory. • Redshift, mass function and compactness factor: Fig. 15-17 exhibit the compactness factor, mass function along with gravitational redshift. It can be seen that both μ(r) ≤ 0.30 and m(r) show increasing behavior. Apart from this, Z shows the decreasing attribute and Z ≤ 5 which is in alliance with the stability of the configuration. s s It is worth mentioning here that our obtained solutions in current study represent more dense stellar structures as compared to past related works on compact objects in f(R,T ) gravity [23, 25, 26, 49, 50]. 16 Appendix (I) 4 3 2 2 3 48β(3λ + 2)r ψ ψ (n − 1)r ψ + 1 3 2 2 2 Υ = r ψ + 1, Υ = , Υ = Υ (n − 1)r ψ − n + 4 , 1 2 2 3 1 2 3−n 2 3 Υ + r ψ ψ Υ 2 3 1 2 2 n 2 2 2 Υ = Υ (n − 1)(2n − 3)r ψ + 5n − 8 , Υ = r ψ ψ Υ + Υ , Υ = r ψ (n − 1)(2n − 3)r ψ + 9n − 15 + 6, 4 1 2 5 2 3 6 2 2 1 1 2 4 2 2 Υ = r ψ (n − 1)(n + 1)(2n − 3)r ψ + 3(n(n + 4) − 9)r ψ + 6n − 3 + 3, 7 2 2 22(3λ+2)nr ψ (−Υ ) 2 2 2 2 2 1 2βr ψ ψ (n − 1)r ψ + 1 48λ − + 24 2 2 2 2 2 1 Υ = r ψ (n − 1)r ψ − 3n + 9 + 3, Υ =  , 8 2 2 9 3−n 2 2 Υ + r ψ ψ Υ 2 3 1 n−2 2 2 n 8β(3λ + 2)ψ ψ Υ Υ Υ − r ψ ψ Υ Υ 3 2 4 3 2 3 1 1 1 2 n−2 Υ = λ + + λr ψ ψ Υ , 10 2 3 2rψ ψ (6λr+r) (n−1)r ψ +1 3 2 ( 2 ) 3 4 2 4βn r ψ − 4(3λ + 1) 3−n 2 2 Υ +r ψ ψ Υ 2 3 1 Υ = , 2 2 2 n−2 2 Υ = 4ψ (n − 1)r ψ + 1 (2λ + 1)r − 6βλ r ψ ψ Υ + 1 + β −6λ + × (3λ + 2)r ψ (−8n 12 3 2 2 3 2 ! !! n+1 2 2 n 7ψ Υ Υ Υ − r ψ ψ Υ Υ 3 4 3 2 3 4 2 2 2 1 1 1 × r ψ − 1 + 6nr ψ r ψ − 3 − − 8 , 2 2 (8 − 24λ)rψ (n − 1)r ψ + 1 3 2 2 2 n−2 Υ = 8β(3λ + 2) + 2λr r ψ ψ Υ + 1 + 2βrψ 13 2 3 2 3−n Υ + r ψ ψ Υ 2 3 1 3 2 2 2 n−2 2(9λ + 11)r ψ ψ (n − 1)r ψ + 1 2(3λ + 4)rψ Υ Υ 4n (−Υ ) (6λr + r) 2 2 3 3 15 1 −  + − , 3−n 2 2 2 2 Υ Υ Υ + r ψ ψ Υ 2 3 1 1 5 1 2 n−2 2 2 2 n Υ = r ψ ψ Υ + 1, Υ = (n − 1)r ψ + 1, Υ = Υ Υ − r ψ ψ Υ Υ , 14 2 3 15 2 16 4 3 2 3 1 1 1 n−2 4 3 2 3 24(λ + 1)ψ Υ Υ 4(3λ − 2)n (−Υ ) 12(3λ + 2)r ψ ψ Υ 3 16 1 1 3 2 15 Υ = − + − , 2 2 3−n 2 3 Υ Υ Υ + r ψ ψ Υ 5 1 2 3 1 n−3 n+2 2 2 2 2 2 2n 4 4(3λ + 4)nr ψ r ψ − 3 4(3λ + 2)r ψ ψ Υ r ψ ψ (Υ Υ ) Υ − ψ Υ Υ + (n − 2)Υ Υ 2 2 3 2 2 1 8 3 7 6 1 3 1 1 1 Υ = + , 3 3 Υ Υ 1 5 2 2 2 2 2 n−2 8(9λ + 8)r ψ ψ Υ 6(5λ + 6)nr ψ (−Υ ) 14r ψ ψ Υ Υ 2 2 1 3 2 16 3 15 Υ = + ψ Υ + (3λ + 2) + 4 , 19 3 15 3−n 2 2 2 2 Υ + r ψ ψ Υ Υ Υ 2 3 1 1 5 2rψ ψ Υ (r−2λr) 2 3 15 3 5 3 8βn r ψ 4λ + 4 7 4 3 3−n 2 2 8β(3λ + 1)n r ψ 2nr ψ (−Υ ) 24βrψ ψ Υ Υ +r ψ ψ Υ 2 3 1 2 1 2 3 15 2 1 Υ = + + λ + , 4 3 2 3−n Υ Υ Υ Υ + r ψ ψ Υ 1 1 1 2 3 1 2 2 2 n−2 4 2 3 3 2(5λ + 1)n r ψ (−Υ ) 4(7λ + 4)ψ Υ Υ 24λr ψ ψ Υ 2 1 3 16 2 3 15 Υ = − − + , 4 2 3−n 2 3 Υ Υ Υ + r ψ ψ Υ 1 5 2 3 1 4 2 2 2λ (11n+24)r ψ +(48−11n)r ψ +24 ( 2 ) 2 n−2 2 2 2 2λr ψ ψ Υ Υ 3 2 16 r ψ ψ Υ + 24 1 2 2 3 15 4n (−Υ ) 8λ + + 4 1 2 Υ = − , 3−n 2 2 2 Υ + r ψ ψ Υ Υ 2 3 1 1 2 4 2 2 2 Υ = 4(λ + 1)nΥ r ψ − 3 + λ 24nr ψ − 24nr ψ Υ + 3nΥ 23 1 2 2 1 2 1 n+1 2 2 2n n+2 4 2ψ Υ r ψ ψ (Υ Υ ) Υ − ψ Υ Υ + (n − 2)Υ Υ 3 2 1 8 3 7 6 3 1 1 1 1 + , 5 17 n+1 7λψ Υ Υ 3 16 2 1 4 2 2 2 2r ψ − + 4n r ψ − 1 + 6λnr ψ r ψ − 3 2 2 2 2 Υ = − 10λ, n−2 2 2 2 (λ − 2)ψ Υ Υ 2(2λ − 1)n (−Υ ) (7λ − 3)r ψ ψ Υ 3 16 1 2 3 15 2 1 Υ = β 4r ψ − − + 25 2 2 2 3−n 2 2 Υ Υ Υ + r ψ ψ Υ 5 1 2 3 1 4(λ + 3)ψ Υ 3 15 + − 8λ + 2λr Υ , 3−n Υ + r ψ ψ Υ 2 3 1 2 2 2 2 2 4 3 2 3 4(λ − 2)nr ψ r ψ − 3 8(λ − 3)r ψ ψ Υ 12λr ψ ψ Υ 2 2 2 3 15 3 2 15 Υ = + − , 3−n 3 3−n 2 2 2 3 Υ + r ψ ψ Υ Υ Υ + r ψ ψ Υ 2 3 1 1 2 3 1 1 1 n−3 n+2 2 2 2 2n 4 4λr ψ ψ Υ r ψ ψ (Υ Υ ) Υ − ψ Υ Υ + (n − 2)Υ Υ 3 2 2 1 8 3 7 6 1 3 1 1 1 Υ = , 7λψ Υ Υ 2 3 16 2r ψ + (13λ − 4)n (−Υ ) 2 2 1 Υ = 68λ + + 32, 2n−4 4 6 3 4 2 2 2 4β(λ + 1)n r ψ 44β(3λ + 2)nr ψ ψ (−Υ ) Υ Υ 2 3 2 15 1 Υ = − − Υ 29 2 4 n 2 2 4 2 Υ (r ψ (ψ Υ + 2) + r ψ + 1) 2 3 1 1 2 2 2 n 4 2 2(λ + 1)nr (−Υ ) r ψ (ψ Υ + 2) + r ψ + 1 1 2 3 1 2 − , 2(2λ+3)r ψ ψ Υ 3 4 2 3 2 15 4βn r ψ − 2 3−n 2 2 Υ +r ψ ψ Υ 2 3 1 n−2 2 n−2 1 Υ = 4βψ Υ 14λ + (λ + 1)r ψ ψ Υ + 8 − , 30 3 2 3 1 1 n+2 n−2 2λψ Υ Υ 2(3λ + 2)r ψ ψ Υ Υ 3 16 3 2 16 2 2 2 1 2 1 2 Υ = (7λ + 5)n r ψ (−Υ ) + + 2n (−Υ ) Υ 4λ + + 2 + 2r ψ 31 2 1 1 2 2 1 2 Υ Υ 5 5 2 4 2 2 2 × 2(4λ + 3)nΥ r ψ − 3 − (3λ + 2) 24nr ψ − 24nr ψ Υ + 3nΥ 1 2 2 1 2 1 !! n+1 n+2 2 2 2n 4 2ψ Υ r ψ ψ (Υ Υ ) Υ − ψ Υ Υ + (n − 2)Υ Υ 3 2 1 8 3 7 6 1 3 1 1 1 + , n+1 7ψ Υ Υ 3 16 2 1 2 2 4βr ψ 2(10λ + 7)n (−Υ ) Υ + (3λ + 2) − 6nr ψ r ψ − 3 2 1 1 2 2 2 Υ = 4β(9λ + 7) + − Υ 32 14 × 12β(λ + 1) + (2λ + 1)r , n−2 2 2 2 (7λ + 6)ψ Υ Υ 2(2λ + 3)n (−Υ ) (21λ + 17)r ψ ψ Υ 3 16 1 2 2 3 15 Υ = −8β(3λ + 2) + 4βr ψ − − + 33 2 2 2 3−n 2 2 Υ Υ Υ + r ψ ψ Υ 5 1 2 3 1 4(3λ + 4)ψ Υ 2(λ + 1)r ψ ψ Υ 3 15 3 2 15 2 2 − + 2(λ + 1)r Υ , Υ = − 2 12β(λ + 1) + r , 14 34 3−n 3−n 2 2 Υ + r ψ ψ Υ Υ + r ψ ψ Υ 2 3 1 2 3 1 1 1 2 n−3 2 2 2n n+2 4 Υ = 4(3λ + 2)r ψ ψ Υ r ψ ψ (Υ Υ ) Υ − ψ Υ Υ + (n − 2)Υ Υ , 35 3 2 2 1 8 3 7 6 3 1 1 1 1   (3λ+2)ψ Υ Υ 2 3 16 14r ψ (5λ + 4)(−n) (−Υ ) − 2 1 2   Υ = ψ Υ − 4(11λ + 6) . 36 3 15 1 18 Appendix (II) 8 16 n 7 14 n 6 12 n 5 10 n 4 8 n 3 6 n λ = −8ψ nψ r Υ − 56ψ nψ r Υ − 168ψ nψ r Υ − 280ψ nψ r Υ − 280ψ nψ r Υ − 168ψ nψ r Υ 1 3 3 3 3 3 3 2 1 2 1 2 1 2 1 2 1 2 1 2 4 n 2 n 16 8 n 14 7 n 12 6 n 10 5 n 8 4 n −56ψ nψ r Υ − 8ψ nψ r Υ + r ψ ψ Υ + 8r ψ ψ Υ + 28r ψ ψ Υ + 56r ψ ψ Υ + 70r ψ ψ Υ 3 3 2 3 3 3 3 3 2 1 1 2 1 2 1 2 1 2 1 2 1 6 3 n 4 2 n 2 n n 2 7 14 2n 2 6 12 2n 2 5 10 2n +56r ψ ψ Υ + 28r ψ ψ Υ + 8r ψ ψ Υ + ψ Υ − 12ψ nψ r Υ − 60ψ nψ r Υ − 120ψ nψ r Υ 3 3 2 3 3 2 1 2 1 1 1 3 2 1 3 2 1 3 2 1 2 4 8 2n 2 3 6 2n 2 2 4 2n 14 7 2 2n 12 6 2 2n 10 5 2 2n −120ψ nψ r Υ − 60ψ nψ r Υ − 12ψ nψ r Υ + 4r ψ ψ Υ + 24r ψ ψ Υ + 60r ψ ψ Υ 3 2 1 3 2 1 3 2 1 2 3 1 2 3 1 2 3 1 8 4 2 2n 6 3 2 2n 4 2 2 2n 2 2 2n 3 6 12 3n 3 5 10 3n +80r ψ ψ Υ + 60r ψ ψ Υ + 24r ψ ψ Υ + 4r ψ ψ Υ − 8ψ nψ r Υ − 24ψ nψ r Υ 2 3 1 2 3 1 2 3 1 3 1 3 2 1 3 2 1 3 4 8 3n 3 3 6 3n 12 6 3 3n 10 5 3 3n 8 4 3 3n 6 3 3 3n −24ψ nψ r Υ − 8ψ nψ r Υ + 6r ψ ψ Υ + 24r ψ ψ Υ + 36r ψ ψ Υ + 24r ψ ψ Υ 3 2 1 3 2 1 2 3 1 2 3 1 2 3 1 2 3 1 4 2 3 3n 6 3 4 4n 9 18 8 16 7 14 6 12 5 10 4 8 +6r ψ ψ Υ + 4r ψ ψ Υ − 2nψ r − 18nψ r − 72nψ r − 168nψ r − 252nψ r − 252nψ r 2 3 1 2 3 1 2 2 2 2 2 2 3 6 2 4 2 −168nψ r − 72nψ r − 18nψ r − 2n, 2 2 4 16 9 4 14 8 4 12 7 4 10 6 4 8 5 4 6 4 4 4 3 3 16 9 λ = 2βn r ψ + 12βn r ψ + 30βn r ψ + 40βn r ψ + 30βn r ψ + 12βn r ψ + 2βn r ψ + 4βn r ψ 2 2 2 2 2 2 2 2 3 14 8 3 12 7 3 10 6 3 8 5 3 6 4 3 4 3 3 4 3 n +20βn r ψ + 36βn r ψ + 20βn r ψ − 20βn r ψ − 36βn r ψ − 20βn r ψ + 76βn r ψ ψ Υ 2 2 2 2 2 2 2 1 3 2 2 2 2 16 9 2 14 8 2 12 7 2 10 6 2 8 5 2 6 −4βn r ψ − 22βn ψ + 18βn r ψ + 128βn r ψ + 368βn r ψ + 528βn r ψ + 340βn r ψ − 32βn r 2 2 2 2 2 2 4 2 4 3 2 4 3 n 2 2 2 2 2 2 n n ×ψ − 192βn r ψ + 660βn r ψ ψ Υ − 112βn r ψ + 130βn r ψ ψ Υ + 80βnψ − 160βψ ψ Υ + 144βnψ 3 3 2 3 2 2 2 2 2 1 2 2 1 1 n 16 9 14 8 12 7 10 6 10 4 5 4n 10 5 4 4n ×ψ Υ + 16βnr ψ + 192βnr ψ + 896βnr ψ + 2240βnr ψ − 2nr ψ ψ Υ + 4r ψ ψ Υ + 3360 1 2 2 2 2 3 2 1 2 3 1 8 5 8 4 4 4n 8 4 4 4n 8 4 5 5n 6 4 6 4 n 4 3 ×βnr ψ − 2nr ψ ψ Υ + 8r ψ ψ Υ + r ψ ψ Υ + 3136βnr ψ − 2720βr ψ ψ Υ + 1792βnr ψ 2 3 2 1 2 3 1 2 3 1 2 2 1 2 4 3 n 4 3 n 2 2 2 2 n 2 2 n −2208βr ψ ψ Υ + 1840βnr ψ ψ Υ + 576βnr ψ − 928βr ψ ψ Υ + 800βnr ψ ψ Υ 3 3 3 3 2 1 2 1 2 2 1 2 1 6 10 n 5 8 n 4 14 8 n 3 14 8 n 2 14 8 n 14 λ = −480ψ ψ r βΥ − 1760ψ ψ r βΥ + 18n r βψ ψ Υ + 36n r βψ ψ Υ + 18n r βψ ψ Υ + 32r β 3 3 3 3 3 3 2 1 2 1 2 1 2 1 2 1 8 n 4 12 7 n 3 12 7 n 2 12 7 n 12 7 n 12 7 n ×ψ ψ Υ + 72n r βψ ψ Υ + 220n r βψ ψ Υ + 212n r βψ ψ Υ + 80nr βψ ψ Υ + 32r βψ ψ Υ 3 3 3 3 3 3 2 1 2 1 2 1 2 1 2 1 2 1 4 10 6 n 3 10 6 n 2 10 6 n 10 6 n 4 8 5 n +108n r βψ ψ Υ + 520n r βψ ψ Υ + 798n r βψ ψ Υ + 544nr βψ ψ Υ + 72n r βψ ψ Υ + 600 3 3 3 3 3 2 1 2 1 2 1 2 1 2 1 3 8 5 n 2 8 5 n 8 5 n 4 6 4 n 3 6 4 n 2 6 ×n r βψ ψ Υ + 1432n r βψ ψ Υ + 1520nr βψ ψ Υ + 18n r βψ ψ Υ + 340n r βψ ψ Υ + 1358n r β 3 3 3 3 3 2 1 2 1 2 1 2 1 2 1 4 n 6 4 n 4 8 5 2 2n 3 8 5 2 2n 2 8 5 2 2n 8 5 ×ψ ψ Υ + 2240nr βψ ψ Υ + 16n r βψ ψ Υ + 192n r βψ ψ Υ + 440n r βψ ψ Υ + 240nr βψ 3 3 2 1 2 1 2 3 1 2 3 1 2 3 1 2 2 2n 3 6 4 2 2n 2 6 4 2 2n 6 4 2 2n 2 4 3 2 2n 4 3 2 2n ×ψ Υ + 80n r βψ ψ Υ + 424n r βψ ψ Υ + 464nr βψ ψ Υ + 144n r βψ ψ Υ + 392nr βψ ψ Υ 3 1 2 3 1 2 3 1 2 3 1 2 3 1 2 3 1 2 2 2 2n 2 6 10 2n 2 5 8 2n 2 4 6 2n 2 3 4 2n +120nr βψ ψ Υ − 284ψ ψ r βΥ − 1110ψ ψ r βΥ − 1800ψ ψ r βΥ − 1430ψ ψ r βΥ 2 3 1 3 2 1 3 2 1 3 2 1 3 2 1 2 2 2 2n 2 2n −540ψ ψ r βΥ − 74ψ ψ βΥ 3 2 1 3 1 4 12 7 2 2n 4 10 6 2 2n 3 12 7 2 2n 3 10 6 2 2n 2 12 7 2 2n λ = 16βn r ψ ψ Υ + 32βn r ψ ψ Υ + 32βn r ψ ψ Υ + 144βn r ψ ψ Υ + 24βn r ψ ψ Υ + 184 2 3 1 2 3 1 2 3 1 2 3 1 2 3 1 2 10 6 2 2n 2 10 6 3 3n 2 8 5 3 3n 2 6 3 4 3n 12 7 2 2n 12 ×βn r ψ ψ Υ + 24βn r ψ ψ Υ + 16βn r ψ ψ Υ − 8βn r ψ ψ Υ + 8βnr ψ ψ Υ − 10βr 2 3 1 2 3 1 2 3 1 3 2 1 2 3 1 2 7 2n 10 6 2 2n 10 6 3 3n 10 3 6 3n 8 5 3 3n 8 3 5 3n ×ψ ψ Υ + 56βnr ψ ψ Υ + 24βnr ψ ψ Υ − 62βr ψ ψ Υ + 56βnr ψ ψ Υ − 248βr ψ ψ Υ 3 2 1 2 3 1 2 3 1 3 2 1 2 3 1 3 2 1 8 4 5 4n 6 4 3 3n 6 3 4 3n 6 4 4 4n 6 5 4 5n 4 3 3 3n −22βr ψ ψ Υ + 8βnr ψ ψ Υ − 404βr ψ ψ Υ − 44βr ψ ψ Υ − 2βr ψ ψ Υ − 312βr ψ ψ Υ 3 2 1 2 3 1 3 2 1 3 2 1 3 2 1 3 2 1 4 3 3 3n 4 4 3 4n 2 3 2 3n −24βnr ψ ψ Υ − 22βr ψ ψ Υ − 94βr ψ ψ Υ 3 2 1 3 2 1 3 2 1 2 8 16 n 8 16 n 2 7 14 n 7 14 n 2 6 12 n 6 12 n λ = −3ψ n ψ r Υ − 2ψ ψ r Υ − 18ψ n ψ r Υ − 16ψ ψ r Υ − 45ψ n ψ r Υ − 56ψ ψ r Υ 5 3 3 3 3 3 3 2 1 2 1 2 1 2 1 2 1 2 1 2 5 10 n 5 10 n 2 4 8 n 4 8 n 2 3 6 n 3 6 n −60ψ n ψ r Υ − 112ψ ψ r Υ − 45ψ n ψ r Υ − 140ψ ψ r Υ − 18ψ n ψ r Υ − 112ψ ψ r Υ 3 3 3 3 3 3 2 1 2 1 2 1 2 1 2 1 2 1 2 2 4 n 2 4 n 2 n 16 8 n 14 7 n 12 6 n −3ψ n ψ r Υ − 56ψ ψ r Υ − 16ψ ψ r Υ + 11nr ψ ψ Υ + 71nr ψ ψ Υ + 195nr ψ ψ Υ + 295 3 3 3 2 3 3 3 2 1 2 1 1 2 1 2 1 2 1 10 5 n 8 4 n 6 3 n 4 2 n 2 n n 2 5 10 2n ×nr ψ ψ Υ + 265nr ψ ψ Υ + 141nr ψ ψ Υ + 41nr ψ ψ Υ + 5nr ψ ψ Υ − 2ψ Υ − 120ψ ψ r Υ 3 3 3 3 2 3 3 2 1 2 1 2 1 2 1 1 1 3 2 1 2 2 4 8 2n 2 4 8 2n 2 2 3 6 2n 2 3 6 2n 2 2 4 2n 2 2 2n −12ψ n ψ r Υ − 160ψ ψ r Υ − 3ψ n ψ r Υ − 120ψ ψ r Υ − 48ψ ψ r Υ − 8ψ ψ r Υ + 114 3 2 1 3 2 1 3 2 1 3 2 1 3 2 1 3 1 8 4 2 2n 6 3 2 2n 4 2 2 2n 2 9 18 2 8 16 2 7 14 2 6 12 2 5 ×nr ψ ψ Υ + 51nr ψ ψ Υ + 9nr ψ ψ Υ − n ψ r − 8n ψ r − 28n ψ r − 56n ψ r − 70n ψ 2 3 1 2 3 1 2 3 1 2 2 2 2 2 10 18 9 2 4 8 16 8 14 7 2 3 6 12 6 10 5 2 2 4 ×r + 3nr ψ − 56n ψ r + 25nr ψ + 92nr ψ − 28n ψ r + 196nr ψ + 266nr ψ − 8n ψ r 2 2 2 2 2 2 2 2 8 4 6 3 2 2 4 2 2 +238nr ψ + 140nr ψ − n ψ r + 52nr ψ + n + 11nr ψ 2 2 2 2 2 19 n n 2 2 7 14 2n 2 7 14 2n 2 2 6 12 2n 2 6 12 2n λ = −180ψ nψ βΥ + 200βψ ψ Υ − 3ψ n ψ r Υ − 8ψ ψ r Υ − 12ψ n ψ r Υ − 48ψ ψ r Υ 6 3 2 2 3 1 1 3 2 1 3 2 1 3 2 1 3 2 1 2 2 5 10 2n 14 7 2 2n 12 6 2 2n 10 5 2 2n 3 6 12 3n 3 2 6 12 3n −18ψ n ψ r Υ + 15nr ψ ψ Υ + 69nr ψ ψ Υ + 126nr ψ ψ Υ − 12ψ ψ r Υ − ψ n ψ r Υ 3 2 1 2 3 1 2 3 1 2 3 1 3 2 1 3 2 1 3 5 10 3n 3 2 5 10 3n 3 4 8 3n 3 2 4 8 3n 3 3 6 3n 3 2 4 3n 12 −48ψ ψ r Υ − 2ψ n ψ r Υ − 72ψ ψ r Υ − ψ n ψ r Υ − 48ψ ψ r Υ − 12ψ ψ r Υ + 9nr 3 2 1 3 2 1 3 2 1 3 2 1 3 2 1 3 2 1 6 3 3n 10 5 3 3n 8 4 3 3n 6 3 3 3n 4 5 10 4n 4 4 8 4n 4 3 6 4n ×ψ ψ Υ + 25nr ψ ψ Υ + 23nr ψ ψ Υ + 7nr ψ ψ Υ − 8ψ ψ r Υ − 16ψ ψ r Υ − 8ψ ψ r Υ 2 3 1 2 3 1 2 3 1 2 3 1 3 2 1 3 2 1 3 2 1 10 5 4 4n 8 4 4 4n 5 4 8 5n 2 10 6 2 8 5 2 6 4 2 4 3 +2nr ψ ψ Υ + 2nr ψ ψ Υ − 2ψ ψ r Υ + 84n r βψ + 560n r βψ + 700n r βψ + 420n r βψ 2 3 1 2 3 1 3 2 1 2 2 2 2 2 2 2 4 9 16 2 9 16 9 16 4 8 14 3 8 14 2 8 14 +124n r βψ − 6n ψ r β − 30n ψ r β − 36nψ r β − 36n ψ r β − 16n ψ r β − 140n ψ r β − 448n 2 2 2 2 2 2 2 8 14 4 7 12 3 7 12 2 7 12 7 12 4 6 10 3 6 10 ×ψ r β − 90n ψ r β − 96n ψ r β − 196n ψ r β − 2032nψ r β − 120n ψ r β − 240n ψ r β − 4800n 2 2 2 2 2 2 2 6 10 4 5 8 3 5 8 5 8 4 4 6 3 4 6 4 6 4 3 4 ×ψ r β − 90n ψ r β − 320n ψ r β − 6680nψ r β − 36n ψ r β − 240n ψ r β − 5696nψ r β − 6n ψ r β 2 2 2 2 2 2 2 2 3 3 4 3 4 3 2 2 2 2 2 −96n ψ r β − 2928nψ r β − 16n ψ r β − 832nψ r β − 100nψ β + 14n βψ 2 2 2 2 2 2 3 8 14 n 2 8 14 n 8 14 n 3 7 12 n 2 7 12 n λ = −104ψ n ψ r βΥ − 90ψ n ψ r βΥ − 88ψ ψ r βΥ − 376ψ n ψ r βΥ − 840ψ n ψ r βΥ 7 3 3 3 3 3 2 1 2 1 2 1 2 1 2 1 7 12 n 3 6 10 n 2 6 10 n 6 10 n 3 5 8 n −12ψ nψ r βΥ − 464ψ n ψ r βΥ − 2426ψ n ψ r βΥ − 1848ψ nψ r βΥ − 176ψ n ψ r βΥ − 3104 3 3 3 3 3 2 1 2 1 2 1 2 1 2 1 2 5 8 n 5 8 n 2 4 6 n 4 6 n 2 3 4 n ×ψ n ψ r βΥ − 5692ψ nψ r βΥ − 1806ψ n ψ r βΥ − 7608ψ nψ r βΥ − 344ψ n ψ r βΥ − 5124ψ n 3 3 3 3 3 3 2 1 2 1 2 1 2 1 2 1 3 4 n 2 2 n 4 14 8 n 14 8 n 4 12 7 n 12 7 n ×ψ r βΥ − 1640ψ nψ r βΥ + 2n r βψ ψ Υ + 88nr βψ ψ Υ + 8n r βψ ψ Υ + 56r βψ ψ Υ 3 3 3 3 3 2 1 2 1 2 1 2 1 2 1 2 1 4 10 6 n 10 6 n 4 8 5 n 8 5 n 4 6 4 n 3 6 4 +12n r βψ ψ Υ + 1800r βψ ψ Υ + 8n r βψ ψ Υ + 5080r βψ ψ Υ + 2n r βψ ψ Υ + 56n r βψ ψ 3 3 3 3 3 3 2 1 2 1 2 1 2 1 2 1 2 n 6 4 n 3 4 3 n 4 3 n 2 2 2 n 2 2 n 4 ×Υ + 6520r βψ ψ Υ + 40n r βψ ψ Υ + 4392r βψ ψ Υ + 34n r βψ ψ Υ + 1496r βψ ψ Υ + 3454r 3 3 3 3 3 1 2 1 2 1 2 1 2 1 2 1 3 2 2n 2 2 2 2n 2 2n 2 2 3 4 2n 2 3 4 2n 2 2 2 2n ×βψ ψ Υ + 1220r βψ ψ Υ + 106βψ ψ Υ − 468ψ n ψ r βΥ − 1408ψ nψ r βΥ − 520ψ nψ r βΥ 2 3 1 2 3 1 3 1 3 2 1 3 2 1 3 2 1 2 12 7 2 2n 12 7 2 2n 10 6 2 2n 8 5 2 2n 6 4 2 2n 2 4 λ = 60n r βψ ψ Υ + 114r βψ ψ Υ + 772r βψ ψ Υ + 2534r βψ ψ Υ + 4216r βψ ψ Υ − 48ψ n 2 3 1 2 3 1 2 3 1 2 3 1 2 3 1 3 7 12 2n 2 3 7 12 2n 2 7 12 2n 2 4 6 10 2n 2 3 6 10 2n 2 2 ×ψ r βΥ − 48ψ n ψ r βΥ − 216ψ nψ r βΥ − 96ψ n ψ r βΥ − 320ψ n ψ r βΥ − 96ψ n 2 1 3 2 1 3 2 1 3 2 1 3 2 1 3 6 10 2n 2 6 10 2n 2 4 5 8 2n 2 3 5 8 2n 2 2 5 8 2n 2 ×ψ r βΥ − 560ψ nψ r βΥ − 48ψ n ψ r βΥ − 496ψ n ψ r βΥ − 840ψ n ψ r βΥ − 840ψ n 2 1 3 2 1 3 2 1 3 2 1 3 2 1 3 5 8 2n 2 3 4 6 2n 2 2 4 6 2n 2 4 6 2n 10 6 3 3n 8 ×ψ r βΥ − 224ψ n ψ r βΥ − 1152ψ n ψ r βΥ − 1384ψ nψ r βΥ + 126r βψ ψ Υ + 528r β 2 1 3 2 1 3 2 1 3 2 1 2 3 1 5 3 3n 2 6 4 3 3n 6 4 3 3n 4 3 3 3n 4 3 3 3n 2 2 3 3n ×ψ ψ Υ + 16n r βψ ψ Υ + 924r βψ ψ Υ + 52nr βψ ψ Υ + 800r βψ ψ Υ + 278r βψ ψ Υ 2 3 1 2 3 1 2 3 1 2 3 1 2 3 1 2 3 1 3 2 6 10 3n 3 6 10 3n 3 2 5 8 3n 3 5 8 3n 3 4 6 3n 8 −48ψ n ψ r βΥ − 60ψ nψ r βΥ − 32ψ n ψ r βΥ − 164ψ nψ r βΥ − 52ψ nψ r βΥ + 42r β 3 2 1 3 2 1 3 2 1 3 2 1 3 2 1 5 4 4n 6 4 4 4n 4 3 4 4n 6 4 5 5n ×ψ ψ Υ + 108r βψ ψ Υ + 66r βψ ψ Υ + 6r βψ ψ Υ 2 3 1 2 3 1 2 3 1 2 3 1 References [1] Baade and Zwicky., PNAS 112, 1241 (2015). [2] K.R. Karmarkar, Proc. Indian Acad. Sci. A 27, 56 (1948). [3] L. Schlai, Ann. Mat. 5, 170 (1871). [4] J. Nash, Ann. of Math. 63, 20 (1956). [5] S.K. Maurya et al., Eur. Phys. J. C 75, 389 (2015). [6] S.K. Maurya et al., Eur. Phys. J. A 52 191 (2016). [7] P. Bhar et al., Eur. Phys. J. A 52, 312 (2016). [8] D. Deb, S.V. Ketov, S.K. Maurya, Mon. Not. R. Astron. Soc. 485, 5652 (2019). [9] S.K. Maurya, et al., Phys. Rev. D 100, (2019). [10] S.K. Maurya et al. , Ann. Physics 385, 532 (2017). [11] S.K. Maurya et al., Eur. Phys. J. C 76, 266 (2016). [12] S.K. Maurya et al., Eur. Phys. J. C 77, 1 (2016). [13] T. Harko et al., Phys. Rev. D 84, 024020 (2011). [14] M. Jamil, D. Momeni and R. Myrzakulov., Eur. Phys. J. C. 72, 1959 (2012). [15] M. Jamil, D. Momeni and R. Myrzakulov., Chin. Phys. Lett. 29, 109801 (2012). [16] H. Shabani and M. Farhoudi., Chin. Phys. Lett. 88, 044048 (2013). [17] H. Shabani and M. Farhoudi., Phys. Rev. D 90, 44031 (2014). [18] P.H.R.S. Moraes ., Eur. Phys. J. C 75, 168 (2015). [19] A. Alhamzawi and R. Alhamzawi., Int. J. Mod. Phys. D 25, 1650020 (2016). [20] P.H.R.S. Moraes, J.D. Arbanil and M. Malheiro., J. Cosmol. Astropart. Phys. 6, 5 (2016). [21] Das et al., Eur. Phys. J. C 76, 654 (2016). 20 [22] Moraes, al., arXiv:1806.04123v4. [23] Z. Yousaf, M.Z. Bhatti and M. Ilyas ., Eur. Phys. J. C 78, 307 (2018). [24] S.K.Mauryaa and F. T. Ortizb., Phys. Dark Universe 27, 100442 (2020). [25] S. Waheed., Symmetry 12, 962 (2020). [26] G. Mustafa et al., Eur. Phys. J. C 80, 26 (2020). [27] M. Ruderman., Annu Rev Astron Astr, 10, 427 (1972). [28] A. V. Astashenok, S. Capozziello and S. D. Odintsov., J. Cosmol. Astropart. Phys. 01, 001 (2015). [29] A. V. Astashenok, S. Capozziello and S. D. Odintsov., Phys. Lett. B 742, 160 (2015). [30] D. Momeni, P. H. R. S. Moraes and R. Myrzakulov., Astrophys. Space Sci 361, 228 (2018). [31] D. Momeni, M. Raza and R. Myrzakulov., Mod. Phys. Lett. A 31, 1650073 (2015). [32] S. Capozziello et al., Phys. Rev. D 93, 023501 (2016). [33] A. V. Astashenok, S. Capozziello and S. D. Odintsov., Phys. Rev. D 89, 103509 (2014). [34] L. Herrera, Phys. Lett. A 165, 206 (1992). [35] A. A. Starobinsky, Phys. Lett. B 91, 99 (1980). [36] M. Jamil et al., Eur. Phys. J. C 72, 1999 (2012). [37] H. Shabani and M. Farhoudi, Phys. Rev. D 88, 044048 (2013). [38] C.P. Singh and P. Kumar, Eur. Phys. J. C 74, 3070 (2014). [39] M. Sharif, Z. Yousaf, Astrophys. Space Sci. 354, 471 (2014). [40] I. Noureen and M. Zubair, Astrophys. Space Sci. 356, 103 (2015). [41] I. Noureen et al., Eur. Phys. J. C 75, 323 (2015). [42] S. Chandrasekhar, Astrophys. J. 140, 417 (1964). [43] H. Heintzmann and W. Hillebrandt, Astron. Astrophys. 38, 51 (1975). [44] W. Hillebrandt and K. O. Steinmetz, Astron. Astrophys. 53, 283 (1976). [45] D. Horvat, S. Ilijic and A. Marunovic, Class. Quantum Grav. 28, 025009 (2011). [46] D. D. Doneva and S. S. Yazadjiev, Phys. Rev. D 85, 124023 (2012). [47] H. O. Silva, C. F. B. Macedo, E. Berti and L. C. B. Crispino, Class. Quant. Grav. 32, 145008 (2015). [48] I. Bombaci, Astron. Astrophys. 305, 871 (1996). [49] M. Zubair, G. Abbas, and I. Noureen, Astrophys. Space Sci. 361,8 (2016). [50] A. K. Yadav, M. Mondal, and F. Rahaman, Pramana - J Phys 94, 90 (2020). http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png General Relativity and Quantum Cosmology arXiv (Cornell University)

Anisotropic spheres via embedding approach in $\mathcal{R}+\beta\mathcal{R}^{2}$ gravity with matter coupling

Download PDF
20 pages

Loading next page...
 
/lp/arxiv-cornell-university/anisotropic-spheres-via-embedding-approach-in-mathcal-r-beta-mathcal-r-AshfKqcruR

References

References for this paper are not available at this time. We will be adding them shortly, thank you for your patience.

ISSN
2212-6864
eISSN
ARCH-3332
DOI
10.1016/j.dark.2020.100747
Publisher site
See Article on Publisher Site

Abstract

Anisotropic spheres via embedding approach in R + βR gravity with matter coupling 1, ∗ 2, † 3, ‡ 4, § G. Mustafa, Xia Tie-Cheng, Mushtaq Ahmad, and M. Farasat Shamir Department of Mathematics, Shanghai University, Shanghai, 200444, Shanghai, People’s Republic of China Department of Mathematics, Shanghai University, Shanghai, 200444, Shanghai, People’s Republic of China. National University of Computer and Emerging Sciences, Chiniot-Faisalabad Campus, Pakistan. National University of Computer and Emerging Sciences, Lahore Campus, Pakistan. The manifesto of the current article is to investigate the compact anisotropic matter profiles in the context of one of the modified gravitational theories, known as f(R, T ) gravity, where R is a Ricci Scalar and T is the trace of the energy-momentum tensor. To achieve the desired goal, we capitalized on the spherical symmetric space–time and utilized the embedding class-1 solution via Karmarkar’s condition in modeling the matter profiles. To calculate the unidentified constraints, Schwarzschild exterior solution along with experimental statistics of three different stars LMC X-4, Cen X-3, and EXO 1785-248 are taken under consideration. For the evaluation of the dynamical equations, a unique model f(R, T ) = R + βR + λT has been considered, with β and λ being the real constants. Different physical aspects have been exploited with the help of modified dynamical equations. Conclusively, all the stars under observations are realistic, stable, and are free from all singularities. Keywords: Anisotropic spheres; f(R, T ) gravity; Compact stars; Embedding Class I. I. INTRODUCTION Late time evolution of stellar configurations, triggered by an immense gravitational pull has been anticipated to a great extent in the field of astrophysics and the modified gravitational theories. It expedites the examination of diverse attributes regarding the gravitating source by physical phenomena. Baade and Zwicky [1] forecast the inception of highly dense stellar objects inaugurating the debate that a supernova might be revolutionized into a highly dense star. This reality came into existence when exceptionally magnetized as well as rotating neutrons stars were detected. Therefore, a fundamental shift regarding normal stars to compact stars came into existence. By the newly discovered concept, the normal stars shifted into an extensive range, such as quark stars, neutron stars, gravastars, dark stars, and finally black holes. The actuality of the extensive range of these stars led the researchers to curiosity, regarding the formation of these stars. The stellar death of a normal star occurs, that is when the nuclear fusion reactions cease to act and burn all of their nuclear fuel results in the formation of new compact stars. The newly formed compact stars are primarily distinguished from the normal stars in two ways. Since all the fuel has been utilized by the star, hence the star cannot sustain against the gravitational collapse due to thermal pressure. Analogous to that, the white dwarf is stabilized due to strong degenerate electron pressure, while the neutron star is stabilized due to degenerate neutron pressure. Whereas black holes are entirely the collapsed remnants, therefore there is neither a thermal pressure nor a degenerate pressure sufficient enough to repress the centripetal pull of gravity; as a result, it leads towards the gravitational singularities and the event horizon. The formed compact stellar remnants consist of huge density and relatively small radii in contrast to the normal stars. The intention to investigate the physically stable models, leads us to an analytical approach regarding the Einstein field equations. One of the essential tools is to adopt the embedding class I space-time which transforms a four dimensional manifold into a Euclidean space of higher dimension. The conversion of curved embedding class space–time into higher dimensional space–time is substantial to develop exact new models in the field of astrophysics. The class I embedding condition leads towards a differential equation in the framework of spherically symmetric Electronic address: gmustafa3828@gmail.com Electronic address: xiatc@shu.edu.cn Electronic address: mushtaq.sial@nu.edu.pk Electronic address: farasat.shamir@nu.edu.pk; farasat.shamir@gmail.com arXiv:2101.00208v1 [gr-qc] 1 Jan 2021 2 space–time which connects the gravitational potentials i.e., g and g , the condition is also recognized as the rr tt Karmarkar condition [2]. The Karmarkar’s condition appears to be very influential in exploring new solutions for the astrophysical models. Schlia [3] was the pioneer in developing the Karmarkar condition for a spherically symmetric space–time. The embedding theorem based on the isometrics has been presented by Nash [4]. Maurya et al. [5]-[12] were the first explorers in the aspect of applying the embedding approach to the anisotropic matter configurations. After the new dawn of general relativity (GR), f(R) theory is considered to be quite a fascinating tool for the amplification of GR. Further, many researchers presented different versions of this theory, which were also very prosperous in diverse fields. The recent extension of this theory is regarded as f(R,T ) gravity, which was presented by Harko et al. [13]. The f(R,T ) theory has been the center of attention by many analysts and consequently many intriguing cosmological aspects have been unraveled [14]-[17]. The analysis regarding isotropic matter profile of the self-gravitating system and its stability has been done by Sharif et al. [18]. Alhamzawi and Alhamzawi [19] construed the occurrence of lensing of gravitation in the context of a modified f(R,T ) theory. Moraes et al. [20] numerically investigated the stability of the gravitational lensing by utilizing the Tolman–Oppenheimer–Volkov (TOV) equations in f(R,T ) gravity. Das et al. [21] formulated a family of solutions by characterizing the interior geometry of compact stars, permitting conformal motion under the influence of f(R,T ) gravity. Moraes et al. [22] investigated the configurations consisting of hydrostatic equilibrium along with fluids whose pressure was computed from equation of state (EoS) in the light of f(R,T ) gravitational theory. Yousaf et al. [23] investigated the formation of relativistic stellar profiles in the regime of f(R,T ) gravity by utilizing the Krori and Barura model. The study of dense anisotropic profiles consisting of charge has been investigated by Maurya and Aurtiz [24]. In this regard, they utilized the Durgapal–Fuloria model in the context of f(R,T ) gravity and applied gravitational decoupling utilizing geometric deformation. Waheed et al. [25] analyzed the existence of highly dense stellar configurations by utilizing Karmarkar along with the Pandey–Sharma condition. To do so, they used spherically symmetric space–time in the context of f(R,T ) gravity. Mustafa et al. [26] analyzed the Class 1 embedding condition in the presence of anisotropy matter profile and utilized the interior geometry of Schwarzschild along with Kohler–Chao solutions in modified gravity. The matter configuration consisting of nuclear density of 10 gm/cc exhibits the behavior of anisotropy i.e. p − p 6= 0 which exists due to certain factors involving magnetic flux, viscosity, phase transition, etc. t r In this regard, Ruderman [27] is the pioneer who argued about the anisotropy existing at the interior of the stars. Modified gravitational theories have provided an overwhelming approach in analyzing the anisotropic stellar con- figurations inheriting high matter profiles [28]-[33]. This work aims to analyze the modified f(R,T ) gravity to devise a realistic configuration which in nature is anisotropic. For this purpose, we take into account three different matter profiles i.e. LMC X-4, Cen X - 3, and EXO 1785–248, and apply a well-known embedding class 1 approach. Moreover, the structural aspect of anisotropic profiles has been examined by making use of spherically symmetric space–time along with the categorical f(R,T ) gravity model. The layout of this article is as follows: In Section 2, modified field equations have been formulated by utilizing the Karmarkar condition. Section 3 is to provide the matching condi- tions by considering Schwarzschild’s solution. The physical investigation has been done comprehensively in Section 4. Conclusive remarks have been provided in the last Section. II. f(R, T ) THEORY OF GRAVITY The modified form of Einstein-Hilbert action for extended f(R,T ) theory of gravity is defined as follows: S = [L + f(R,T )] −gd x, (1) where L and f(R,T ) denote the matter Lagrangian density and a modified function, respectively. Here, R and T are known as scalar curvature and trace of the energy-momentum tensor, respectively. Now, by varying Eq.(1), we get the following modified set of equations (1 − f (R,T ))T − f (R,T )Θ = − f(R,T )g + (R + (g ✷ − ∇ ∇ ))f (R,T ), (2) T νµ T νµ νµ νµ νµ µ ν R where, αβ 2 g δT ∂ L νµ m αβ Θ = = −2g − 2T + g L , νµ νµ νµ m νµ νµ αβ δg ∂g ∂g νµ ∂ ( −gg ∂ ) ∂f(R,T ) ∂f(R,T ) µ ν ✷ = √ , f (R,T ) = , f (R,T ) = , R T −g ∂R ∂T 3 with ∇, representing the covariant derivative. The energy-momentum tensor with the anisotropic matter source is defined as T = ρU U + p V V + p (U U − g − V V ), (3) νµ µ ν r µ ν t µ ν νµ µ ν where U represents the vector for 4-velocity and V is a vector in the direction of radial pressure. Further, the µ µ expressions, i.e., ρ, p and p are used to define define energy density, tangential and radial components of pressure, t r respectively. Using Eq.(3) in Eq.(2), we get the following set of equation: 1 1 G = (1 + f (R,T ))T + (∇ ∇ − g ✷)f (R,T ) + (f(R,T ) − Rf (R,T ))g T νµ µ ν νµ R R νµ νµ f (R,T ) 2 − ρg f (R,T )) . (4) νµ T We assume a static and spherically symmetric line element, which is defined as: 2 a(r) 2 b(r) 2 2 2 ds = −e dt + e dr + r dΩ , (5) 2 2 2 2 2 b(r) a(r) where the expression dΩ defines the g = r dθ and g = r sin θdφ components, e and e denote the θθ φφ gravitational components of stellar geometry. Further, we fix a quadratic model of f(R,T ), which is defined as: f(R,T ) = R + β × R + λ × T . (6) The considered model f(R,T ) = R + βR + λT involves a particular case of well-known Starobinsky model [35] with matter coupling. It is an interesting point that, in the Starobinsky model R + βR , a maximum value of M/M or beyond is reached when the value of the parameter β is selected to be negative. But, this leads to an issue; specifically, the Ricci scalar performs a damped oscillation. On the other hand, the Ricci scalar smoothly decreases to zero as we approach towards infinity for positive values of parameter β, for which the star can support a maximum mass lower than 2M/M . Now, we elaborate an eminent Karmarkar condition concisely which is the integral tool for current study. The infrastructure connecting the Karmarkar condition is established on the class I space of Riemannian geometry. A sufficient condition comprises of a second order symmetric tensor and the Riemann Christoffel tensor, given as Σ(Λ Λ − Λ Λ ) = R , ηµ υγ γµ νη υµ ηγ Λ ; n − Λ ; ν = 0. νµ νη Here ; stands for covariant derivative whereas Σ = ±1. These values signify a space-like or time-like manifold relying on the sign considered as − or +. Now, the Karmarkar condition is defined as R R = R R + R R . (7) 1414 2323 1224 1334 1212 3434 These Riemann tensor components are given below as follows. a(r) ′′ ′ 2 ′ ′ 2 2 b(r) e (2a (r) + a (r) − a (r)b (r)) r sin θ(e − 1) R = , R = , 1414 2323 b(r) 4 e ′ 2 ′ a(r) rb (r) rsin θb (r)e − b(r) R = , R = , 1212 3434 2 2 R = R sin θ, R = 0, 1334 1224 1224 where, R 6= 0. A differential equation can be achieved by utilizing the Karmarkar condition using Eq. (7) as ′ ′ a (r)b (r) ′ ′ ′ 2 ′′ ′ 2 b(r) − a (r)b (r) + a (r) − 2 a (r) + a (r) = 0, e 6= 1. (8) b(r) 1 − e Integration of Eq. (8) provides a connection between two main gravitational components of the space-time as follows b(r) a(r) ′ 2 e = e × a (r) + 1 + K, (9) 4 where K is a constant of integration. We choose a specific model for a g component which is expressed as tt a(r) 2 e = ψ r ψ + 1 , (10) 1 2 where ψ , ψ are assumed as constants, n is an integer. By plugging Eq. (10) in Eq. (9), we get the g component, 1 2 rr which is calculated as b(r) 2 2 n−2 e = r ψ ψ r ψ + 1 + 1, (11) 2 3 2 where ψ = 4n × ψ × ψ × K. It is mentioned here that we get realistic results for n > 2. Now, we are able to 3 2 1 calculate the following set of modified field equations for the anisotropic stellar configuration. 2 2 2 1 ψ 12β(λ + 1)n r ψ (−Υ ) 2 2 1 ρ = × − + Υ 2 3−n 4 n−2 2 3 Υ Υ + r ψ ψ Υ 2 3 1 2(λ + 1)(2λ + 1) r ψ ψ Υ + r 1 1 2 3 4 6 3 n+2 4β(3λ + 1)n r ψ 8β(3λ + 2) ψ Υ Υ 3 16 2 2 4 2 2 1 2 + − × − + ψ r 4nΥ r ψ − 3 − 3n r ψ − 6r ψ + 1 2 1 2 2 4 4 2 Υ Υ Υ 1 1 5 n+1 2 2ψ Υ 2nr (−Υ ) Υ 3 1 10 n+2 1 2 2 2n 4 − × r ψ ψ (Υ Υ ) Υ − ψ Υ Υ + (n − 2)Υ Υ − − Υ − Υ + Υ 2 1 8 3 7 6 9 11 12 3 1 1 1 3 2 Υ Υ 5 1 2 2 4 n r ψ Υ n 2r ψ ψ Υ 2 13 3 2 15 n−2 2 + + 2ψ Υ 2β(9λ + 5) + Υ (2λ + 1)r − 2(3βλ + β) − × λΥ 3 14 14 2 1 3−n Υ Υ 1 Υ + r ψ ψ Υ 1 2 3 1 2 2 − 4 6β + r + β 24λ + 2r ψ + Υ − Υ , (12) 2 17 18 3−n Υ + r ψ ψ Υ 2 3 1 −1 2ψ Υ Υ 4r ψ Υ 3 15 24 2 23 p = − 4βrψ + + Υ + Υ − 2Υ Υ + 4rψ r 2 21 22 14 20 3 2 3−n 4 4(λ + 1)(2λ + 1)r Υ Υ Υ + r ψ ψ Υ 14 2 3 1 1 2 3 2 2n r ψ Υ 2nrψ 25 2 n−2 2 2 2 × ψ Υ + + Υ (2λ + 1)r − 2(3βλ + β) − 2β(15λ + 7) β 8(5λ + 4) + 2r ψ 14 2 1 2 Υ Υ n−2 8(λ + 2)ψ Υ Υ 4(5λ + 6)n (−Υ ) ψ Υ Υ 3 16 1 3 15 28 × + + + Υ − Υ − Υ 26 27 14 2 2 3−n Υ Υ Υ + r ψ ψ Υ 5 1 2 3 1 2λr ψ ψ Υ 2 3 15 × + 4 (λ + 1)r − 6βλ , (13) 3−n Υ + r ψ ψ Υ 2 3 1 2 2 ψ 4βΥ n r ψ Υ 2ψ Υ Υ 1 2 31 2 33 3 15 32 p = + + − × n (β (−8λ 2 4 2 3−n n−2 3 Υ Υ Υ Υ + r ψ ψ Υ 2 3 1 1 2(λ + 1)(2λ + 1) r ψ ψ Υ + r 1 1 2 3 1 n−2 2 2 2 2 2 4(8λ + 7)ψ Υ Υ 4(3λ + 5)n (−Υ ) 4(7λ + 6)nr ψ r ψ − 3 4(24λ + 19)r ψ ψ Υ 2 2 3 16 1 2 2 1 3 15 + 2r ψ + − − 2 2 3 3−n 2 2 Υ Υ Υ Υ + r ψ ψ Υ 5 1 1 2 3 1 ! ! ! 4 2 3 3 Υ 12(3λ + 2)r ψ ψ Υ Υ 36 35 2 3 15 + + + + 8 + Υ Υ + Υ − Υ , (14) 14 34 29 30 3−n 3−n 3 2 3 Υ + r ψ ψ Υ Υ + r ψ ψ Υ Υ 2 3 1 2 3 1 5 1 1 where Υ , {i = 1, ..., 36}, are given in the Appendix (I). III. COMPARISON OF EXTERIOR AND INTERIOR SOLUTION By considering the Jebsen-Birkhoff’s theorem, the spherically symmetric vacuum solution of GR field equations must be asymptotically flat. In particular, the spacetime is of the form 2 2 −1 2 2 2 2 2 ds = εdt − ε dr − r (dθ + sin θdφ ), (15) 2M where ε = 1 − . Here, M denotes the stellar mass of the star. Now, considering the constraint p (r = R ) = 0 at r ǫ the boundary r = R and using metric coefficients g and g from Eq. (5) and Eq. (15), we calculate the following ǫ tt rr 5 TABLE I: Predicted values of ψ , ψ , ψ , and λ with β = 2 and and (Radii=9.711 km, & Mass =1.29M/M ). 1 2 3 ⊙ LMC X-4 n ψ ψ ψ λ 1 2 3 3 0.4329 0.001273 4.7820 0.066676 5 0.4363 0.000729 7.6640 0.040511 10 0.4387 0.000352 14.900 0.021340 20 0.4399 0.000173 29.392 0.011919 50 0.4406 0.000068 72.882 0.006322 100 0.4408 0.000034 145.373 0.004466 −6 500 0.4410 6.827742×10 725.285 0.002985 TABLE II: Predicted values of ψ , ψ , ψ , and λ with β = 2 and (Radii= 10.136 km, & Mass =1.49M/M ). 1 2 3 ⊙ Cen X - 3 n ψ ψ ψ λ 1 2 3 3 0.3765 0.001421 4.5683 0.055841 5 0.3807 0.000805 7.2751 0.032774 10 0.3838 0.000386 14.082 0.015565 20 0.3852 0.000186 28.267 0.007027 50 0.3861 0.000074 68.656 0.007027 100 0.3864 0.000036 138.251 0.000239 −6 500 0.3866 7.447003×10 682.743 0.000113 expressions 2M ψ ψ R + 1 = 1 − , (16) 1 2 −1 n−2 2M 2 2 ψ ψ R ψ R + 1 + 1 = 1 − , (17) 3 2 2 ǫ ǫ n−1 M nψ ψ ψ R + 1 = , (18) 1 2 2 p (r = R ) = 0. (19) r ǫ Utilizing these boundary conditions from Eqs. (16-19), we get the following relations −n (R − 2M) 1 − 2nM−nR +M ψ = , (20) ψ = , (21) R (nR − (2n + 1)M) 1−n ψ = 2n 1 − , (22) 2nM − nR + M λ + λ + λ + λ 1 2 3 4 λ = , (23) λ + λ + λ + λ 5 6 7 8 where λ , {i = 1, ..., 8} are given in the Appendix (II). The estimated values of the above parameters, i.e., ψ , ψ , ψ , & λ are given in Table I, Table II and Table III. 1 2 3 6 TABLE III: Predicted values of ψ , ψ , ψ , and λ with β = 2 and (Radii=8.849 km, & Mass =1.30M/M ). 1 2 3 ⊙ EXO 1785-248 n ψ ψ ψ λ 1 2 3 3 0.3768 0.001862 4.5698 0.076448 5 0.3811 0.001055 7.2777 0.044225 10 0.3842 0.000496 14.361 0.020799 20 0.3856 0.000248 27.733 0.009361 50 0.3864 0.000098 68.685 0.002592 100 0.3867 0.000048 136.944 0.000352 −6 500 0.3869 9.759977×10 683.029 0.000134 IV. PHYSICAL ANALYSIS In this section, we briefly present the results by analyzing the physical attributes along with different aspects of the stellar configurations under the acquired f(R,T ) gravity model. In order to fulfill the purpose, experimental data of distinct stars i.e. LMC X-4, Cen X – 3 and EXO 1785–248 is used. All the attributes of the stellar configurations are depicted in tabular form as well as graphically. A. Gravitational Metric Potential The existence of anomalies within the sphere such as geometric singularities are contemplated to be an essential peculiarity in the investigation of stellar spheres. In order to unravel the existence of singularities, we examine the a b nature of gravitational potential g = e and g = e at the core r = 0 of the sphere. Physical essence and endurance tt rr of the models rely upon the gravitational metric potentials and it should be decreasing on regular intervals within the spherical structures. It can be observed from Fig. 1 that the metric potential with in the interior of the sphere b(r=0) a(r=0) exhibits the behavior of e = 1 and e 6= 0 which is consistent and physically valid. It can also be observed that both of the metric potentials exhibit minimum values at the center and show non-linear increasing behavior towards the boundary. FIG. 1: Visual representation of gravitational potentials with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) B. Energy Density and Pressure Evolutions Prior to analysis of the anisotropy, we investigate the evolutional change of the matter profiles in connection to the energy density ρ along with anisotropic stresses such as p and p . The energy density along with pr and pt exhibits r t the exceptional behavior of high density of matter configuration. The phenomenal high density is due to the strong forces of attraction which are regarded as dipole interactions and intermolecular forces. Numerical values of density and the components of the pressure for the three compact spheres are provided in Tables IV-VI. All of the physical attributes remain positive and appear to be finite at the core. It confirms that the current system is independent 7 of all singularities. From the Figs. 2-4, it is evident that the matter configuration under consideration attains the maximum mass at the core and tends to zero at the boundary of the star, which depicts the high compactness of the stellar spheres. These graphical plots establish the presence of anisotropy of the compact sphere under the influence of our f(R,T ) model. FIG. 2: Visual representation of ρ with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) FIG. 3: Visual representation of p with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) FIG. 4: Visual representation of p with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) t 8 C. Anisotropy and Gradients In order to model the interior geometry of the relativistic stellar configuration under current circumstances, the role of anisotropy is crucial for the compact sphere modeling and it is represented as Δ = p − p . (24) t r It depicts the information regarding the anisotropic nature of the stellar configuration. If p > p then the anisotropy t r is considered to be non negative and is drawn outwards and depicted as Δ < 0. Whereas if p > p then the anisotropy r t turns out to be negative and this shows that anisotropy is drawn inwards. From 5 it is observed that for our ongoing study anisotropy remains positive, hence, directed outwards. The deviation of radial derivatives of the energy density dρ dp dp r t and pressure components , i.e., , and , are shown in Figs. 6-8 such that dr dr dr FIG. 5: Visual representation of Δ with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) dρ dp dp r t < 0, < 0, < 0. (25) dr dr dr It can be observed from the second order derivatives that the pressure components and energy density show the maximum value at the core r = 0, i. e., 2 2 2 d ρ d p d p r t > 0, > 0, > 0. (26) 2 2 2 dr dr dr FIG. 6: Visual representation of gradient of ρ with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) 9 FIG. 7: Visual representation of gradient of p with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) FIG. 8: Visual representation of gradient of p with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) D. Energy Conditions Energy conditions appear to be quite helpful in analyzing the realistic distribution of matter. These attributes play a decisive role to classify the exotic and normal mater distribution within the stellar model. The energy conditions have been crucially important in debating the issues related to cosmology and astrophysics. The energy condition are classified as NEC : ρ > 0, ρ + p ≥ 0, ρ + p ≥ 0, DEC : ρ ≥ |p |, ρ ≥ |p |, r t r t WEC : ρ − p ≥ 0, ρ − p ≥ 0, r t SEC : ρ − p ≥ 0, ρ − p ≥ 0, ρ − p − 2p ≥ 0. (27) r t r t TABLE IV: Predicted values of physical parameters at center and boundary. LMC X-4 a(r=0) b(r=0) 3 2 3 n e e ρ (g/cm ) p = p (dyne/cm ) ρ (g/cm ) p /ρ = p /ρ R r t 0 r 0 t 0 0 0 0 0 14 34 14 3 0.431 1.0 0.7546 × 10 13.5220 × 10 1.8378 × 10 0.074 14 34 14 5 0.432 1.0 0.8479 × 10 14.1351 × 10 1.7467 × 10 0.081 14 34 14 10 0.433 1.0 0.9735 × 10 15.9780 × 10 1.6629 × 10 0.094 14 34 14 20 0.434 1.0 1.0004 × 10 16.1932 × 10 1.5492 × 10 0.110 14 34 14 50 0.440 1.0 1.0875 × 10 17.1260 × 10 1.4207 × 10 0.124 14 34 14 100 0.441 1.0 1.1787 × 10 18.1910 × 10 1.3328 × 10 0.135 14 34 14 500 0.442 1.0 1.2433 × 10 18.9280 × 10 1.2558 × 10 0.136 10 TABLE V: Predicted values of physical parameters at center and boundary. Cen X - 3 (mass =1.49M/M & radii= 10.136 km) a(r=0) b(r=0) 3 2 3 n e e ρ (g/cm ) p = p (dyne/cm ) ρ (g/cm ) p /ρ = p /ρ R r t 0 r 0 t 0 0 0 0 0 14 34 14 3 0.376 1.0 1.7236 × 10 16.4320 × 10 1.9578 × 10 0.095 14 34 14 5 0.377 1.0 1.8429 × 10 17.1451 × 10 1.8769 × 10 0.103 14 34 14 10 0.378 1.0 1.9325 × 10 18.9380 × 10 1.7629 × 10 0.125 14 34 14 20 0.379 1.0 2.0137 × 10 19.1632 × 10 1.5092 × 10 0.140 14 34 14 50 0.380 1.0 2.0945 × 10 20.8460 × 10 1.4007 × 10 0.164 14 34 14 100 0.381 1.0 2.3527 × 10 21.1220 × 10 1.4428 × 10 0.169 14 34 14 500 0.382 1.0 2.5723 × 10 21.9180 × 10 1.4958 × 10 0.172 TABLE VI: Predicted values of physical parameters at center and boundary. EXO 1785-248 (Mass =1.30M/M & Radii=8.849 km) a(r=0) b(r=0) 3 2 3 n e e ρ (g/cm ) p = p (dyne/cm ) ρ (g/cm ) p /ρ = p /ρ R r t 0 r 0 t 0 0 0 0 0 14 34 14 3 0.376 1.0 0.8876 × 10 15.4634 × 10 1.9684 × 10 0.082 14 34 14 5 0.381 1.0 0.9574 × 10 16.7635 × 10 1.8073 × 10 0.099 14 34 14 10 0.384 1.0 1.0165 × 10 17.3650 × 10 1.7871 × 10 0.114 14 34 14 20 0.385 1.0 1.0969 × 10 18.7532 × 10 1.6063 × 10 0.131 14 34 14 50 0.386 1.0 1.1607 × 10 19.8591 × 10 1.4534 × 10 0.146 14 34 14 100 0.387 1.0 1.2485 × 10 20.1813 × 10 1.3115 × 10 0.153 14 34 14 500 0.388 1.0 1.4524 × 10 20.9994 × 10 1.3089 × 10 0.156 Here, NEC stands for null energy condition, SEC for strong energy condition , DEC for dominant energy condition and WEC for week energy condition. It can be seen from Fig. 9 that all the energy bounds exhibit decreasing behavior with the increase in radii of the compact stellar sphere. FIG. 9: Evolution of energy bounds with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) E. Analysis of Stability The stability of the stellar configuration plays a decisive role in analyzing the consistency of the acquired model. Many analytical discussions have been done in order to find the stability of the matter configuration but Herrera’s cracking conception emerged to be very effective [34]. The radial and tangential speed of sound are defined as dp dp r t 2 2 v = , and v = . (28) sr st dρ dρ For the conservation of causality condition, components of the speed of sound must be with the bounds of the interval 2 2 2 2 [0, 1] i.e 0 ≤ v and v ≤ 1. It can be observed from the Figs. 10-11 that the condition i.e 0 ≤ v and v ≤ 1 is r t r t 2 2 satisfied by both velocity components. Apart from that Abreu condition i.e −1 ≤ |v − v | ≤ 0 has also been satisfied t r 11 and can be observed from the Fig. 11. The validity of both of the aspects affirm the viability and the effectiveness of 2 2 our model. Further, it can also be observed that inverse Abreu condition i.e 0 ≤ |v − v | ≤ 1 is also satisfied. t r 2 2 FIG. 10: Visual representation of v and v with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and r t n = 500(⋆) 2 2 2 2 FIG. 11: Visual representation of |v − v | and |v − v | with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), t r r t n = 100(⋆), and n = 500(⋆) F. Equilibrium Analysis for Modified f(R, T ) Gravity In this section, we will analyze the equilibrium condition by considering the stability of the acquired solution of the three different stellar configuration. For the purpose, we make use of the TOV equation [36]-[41] dp dp dρ r t ′ λ − − 2 + 3 2Δ dp a dr dr dr − − (ρ + p ) + = 0. (29) r dr 2 3(2λ + 1) The above equation characterizes the necessary and sufficient condition for the hydrostatic-equilibrium. It comprises of four different forces dp dp dρ r t ′ λ − − 2 + 3 dr dr dr a dp 2Δ F = (ρ + p ), F = , F = , F = . (30) g r h a e 2 dr r 3(2λ + 1) • F represents the anisotropy force. • F represents the hydrostatic force. • F represents the gravitational force. g 12 • F represents the extra force. Consequently, the TOV equation can also be written as F + F + F + F = 0. From the attained graph as shown g h a e in Fig. 12, it is deduced that all of the forces sum up to neutralize the total effect, and this confirms the existence of the stable stellar structures. FIG. 12: Visual representation of F , F , F , and F with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), g h a e n = 100(⋆), and n = 500(⋆) G. Evolution of Adiabatic Index For the energy density, EoS stiffness can be better described by the adiabatic index. The stability of relativistic as well as the non-relativistic stellar structures can be explained through the adiabatic index. The concept of the dynamical stability via the radial adiabatic index was presented by Chandrasekhar [42]. This was further utilized by many authors [43]-[48]. For the system to be dynamically stable, the adiabatic index must go beyond 4/3. Adiabatic index corresponding to the radial stress is given as ρ + p dp ρ + p r r r Γ = ( ) = v . (31) p dρ p r r One of the quite fascinating fact of the above equation is that the stability of the Newtonian matter configuration is achieved when Γ > 4/3. While if Γ = 4/3 then a neutral equilibrium is achieved whereas if Γ < 4/3 then an r r r unstable matter configuration consisting of anisotropy is achieved. Anisotropic matter profile via adiabatic index can be elaborated as 4 ρ p 4 p − p i ri ti ri Γ = + ( r + ). (32) ′ ′ 3 2 | p | 3 | p | r ri ri From the Fig. 13, the graphical behavior of Γ with respect to increasing radii can be observed. It is noted that Γ r r shows the monotonically increasing conduct for all the stellar spheres and Γ is always greater than 4/3. Hence, Γ r r is consistent for the stability of our model in f(R,T ) gravity. H. Equation of state The evolution of the emergence of compact stars can be determined by EoS of matter. Moreover, EoS has a strong impact on the conditions of nucleosynthesis. Therefore, EoS is a vital tool in many astrophysical simulations. The EoS is considered to be a ratio of the pressure terms p and p with density. The components of EoS for the study r t of stellar configuration are i.e., ω and ω and are mathematically connected as r t p p r t ω = , ω = . (33) r t ρ ρ From Fig. 14, it can be observed that with the increase in radii, the components of the EoS shows monotonically decreasing behavior and are always less than 1. Moreover, the positive nature is observed for both of the components of EoS i.e. ω and ω with in the matter configuration. The accomplishment of the condition i.e. 0 ≤ ω and ω < 1 r t r t unveils that the our obtained solutions are valid and legitimate. 13 FIG. 13: Evolution of Γ with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) FIG. 14: Evolution of energy EoS with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) I. Compactness factor and Surface Redshift For the existence of any matter configuration mass function, compactness factor along with surface redshift function are contemplated to be an essential constituent. The fundamental relation for the mass function is given as m(r) = 4πρr dr. (34) FIG. 15: Evolution of m(r) with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) 14 FIG. 16: Evolution of μ(r) with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) FIG. 17: Evolution of z with n = 3(⋆), n = 5(⋆), n = 10(⋆), n = 20(⋆), n = 50(⋆), n = 100(⋆), and n = 500(⋆) The generalize compactness factor i.e. μ(r) is represented as m 1 μ(r) = = 4πρr dr. (35) r r The strong intermolecular interaction forces with in the stellar matter configuration and its corresponding EoS can be characterized by the term surface redshift i.e., Z . The generalized relation is Z = − 1. (36) 2 n 1/2 (ψ (ψ r + 1) ) 1 2 Figs. 15-17 represent the evolution of redshift function along with compactness factor corresponding to the increasing radii. It can be observed that the surface redshift is always less than 5 i.e., 5 i.e. Z ≤ 5 and the compactness factor remains less than 0.30 i.e., μ(r) ≤ 0.30. All of the mentioned functions are positive throughout the configuration. Hence, our models are stable. V. CONCLUSION The manifesto of the current study is to identify the realistic and stable configuration for the stellar sphere in the modified f(R,T ) gravitational theory. For the analysis of the stellar matter configuration, a viable model f(R,T ) = R + βR + λT is considered along with spherically symmetric space–time. In order to achieve the current objective, observational data of three compact spheres i.e. LMC X-4, Cen X-3 and EXO 1785–248 inheriting an anisotropic matter distribution has been utilized. This realistic range for masses of stars under this study is 1.29 to 1.5 solar mass. In this context, the considered models LMC X-4 (mass 1.29 M/M ), Cen X-3 (mass 1.49 M/M ) and EXO 1785-248 ⊙ ⊙ (mass 1.29 M/M ) are with in suitable range under this study. Moreover, the radii are also with in prescribed range 8.849 km to 10.136 km. In general, the study is valid for other models of stars and those within the given ranges of 15 mass and radius under this study. Embedding class 1 condition is used to find the potential i.e. g by considering rr a(r) 2 the primary potential as i.e g = e = χ r χ + 1 . To find the unknown constraints, Schwarzschild’s exterior tt 1 2 solution has been utilized. All the obtained results can be summarized as: • Metric potentials: The presence of singularities within the stellar configuration is an essential topic, worthy of debate. The stability of the stellar matter configuration depends on it. Therefore, the aspects of metric potential b(r=0) play a decisive role. The graphical behavior of Fig. 1 depicts that the fundamental condition i.e. e = 1 a(r=0) and e 6= 0 has been encompassed by the gravitational potential. Increasing attribute of the potentials is also observed throughout the configuration. Therefore, the potentials are free from any singularity and so is our model. • Energy density and stress constraints: From the Figs. 2-4 the evolution of density and components of stresses i.e., p and p can be observed. The density along with stress components show decreasing behavior and are non- r t negative throughout the configuration. The peak value is accomplished at the core while decreasing evolution is observed with the increase in radii and it tends to 0 towards the boundary. • Anisotropy and gradients: Fig. 5 depicts anisotropy behavior for our current configuration. It is noted that p 6= p and p > p , therefore, △ > 0 so the anisotropy is positive and is directed outwards. From Fig. 6-8 t r t r gradient of density and stress components are reviewed and it is noticed that all the gradients are negative and dρ dp dp r t exhibit decreasing behavior i.e., < 0, < 0, < 0. Since the non-positive behavior of all gradients along dr dr dr with their vanishing attribute at r = 0 is observed, therefore, our stellar configuration is stable. • Energy bounds: From the Figs. 2-4 and Fig. 9 the characteristic behavior of functions ρ, p , p , ρ − p , ρ − p , r t r t ρ−p − 2p are depicted. Decreasing characteristic behavior is observed with the increase in radii which further r t fulfills the bounds of energy such as (NEC), (SEC), (DEC) and (WEC). Hence the matter profile is realistic and viable. • Causality analysis The behavior of the constraints of sound speed is depicted in Figs. 10 and Fig. 11. The decreasing attribute is observed for both the components of the speed and it is shown that they are always 2 2 within the limits i.e., 0 ≤ v ≤ 1 and 0 ≤ v ≤ 1. The fulfillment of Aberu condition is also observed i.e., r t 2 2 2 2 −1 ≤ |v − v | ≤ 0. For the current model v > v . Therefore, fulfillment of all the conditions confirms the t r t r viability of the stellar sphere. • Equilibrium and EoS analysis: The balancing nature of all the forces, i.e., F , F and F is depicted in Fig. a h g 12. As all these forces add up to 0 and balance the effect of each other, therefore, the equilibrium condition is satisfied. From Fig. 14, the attributes of the parameters of EoS are observed. It is concluded that constraints i.e., ω and ω of EoS are positive in the interior of stellar profile and are in the stability bounds of 0 ≤ ω and r t r ω < 1. • Adiabatic index stability analysis: The behavior of adiabatic index Γ can be seen in Fig. 13 showing that Γ > 4/3. It also depicts the positive and decreasing nature, justifying the effectiveness of our system in the framework of f(R,T ) theory. • Redshift, mass function and compactness factor: Fig. 15-17 exhibit the compactness factor, mass function along with gravitational redshift. It can be seen that both μ(r) ≤ 0.30 and m(r) show increasing behavior. Apart from this, Z shows the decreasing attribute and Z ≤ 5 which is in alliance with the stability of the configuration. s s It is worth mentioning here that our obtained solutions in current study represent more dense stellar structures as compared to past related works on compact objects in f(R,T ) gravity [23, 25, 26, 49, 50]. 16 Appendix (I) 4 3 2 2 3 48β(3λ + 2)r ψ ψ (n − 1)r ψ + 1 3 2 2 2 Υ = r ψ + 1, Υ = , Υ = Υ (n − 1)r ψ − n + 4 , 1 2 2 3 1 2 3−n 2 3 Υ + r ψ ψ Υ 2 3 1 2 2 n 2 2 2 Υ = Υ (n − 1)(2n − 3)r ψ + 5n − 8 , Υ = r ψ ψ Υ + Υ , Υ = r ψ (n − 1)(2n − 3)r ψ + 9n − 15 + 6, 4 1 2 5 2 3 6 2 2 1 1 2 4 2 2 Υ = r ψ (n − 1)(n + 1)(2n − 3)r ψ + 3(n(n + 4) − 9)r ψ + 6n − 3 + 3, 7 2 2 22(3λ+2)nr ψ (−Υ ) 2 2 2 2 2 1 2βr ψ ψ (n − 1)r ψ + 1 48λ − + 24 2 2 2 2 2 1 Υ = r ψ (n − 1)r ψ − 3n + 9 + 3, Υ =  , 8 2 2 9 3−n 2 2 Υ + r ψ ψ Υ 2 3 1 n−2 2 2 n 8β(3λ + 2)ψ ψ Υ Υ Υ − r ψ ψ Υ Υ 3 2 4 3 2 3 1 1 1 2 n−2 Υ = λ + + λr ψ ψ Υ , 10 2 3 2rψ ψ (6λr+r) (n−1)r ψ +1 3 2 ( 2 ) 3 4 2 4βn r ψ − 4(3λ + 1) 3−n 2 2 Υ +r ψ ψ Υ 2 3 1 Υ = , 2 2 2 n−2 2 Υ = 4ψ (n − 1)r ψ + 1 (2λ + 1)r − 6βλ r ψ ψ Υ + 1 + β −6λ + × (3λ + 2)r ψ (−8n 12 3 2 2 3 2 ! !! n+1 2 2 n 7ψ Υ Υ Υ − r ψ ψ Υ Υ 3 4 3 2 3 4 2 2 2 1 1 1 × r ψ − 1 + 6nr ψ r ψ − 3 − − 8 , 2 2 (8 − 24λ)rψ (n − 1)r ψ + 1 3 2 2 2 n−2 Υ = 8β(3λ + 2) + 2λr r ψ ψ Υ + 1 + 2βrψ 13 2 3 2 3−n Υ + r ψ ψ Υ 2 3 1 3 2 2 2 n−2 2(9λ + 11)r ψ ψ (n − 1)r ψ + 1 2(3λ + 4)rψ Υ Υ 4n (−Υ ) (6λr + r) 2 2 3 3 15 1 −  + − , 3−n 2 2 2 2 Υ Υ Υ + r ψ ψ Υ 2 3 1 1 5 1 2 n−2 2 2 2 n Υ = r ψ ψ Υ + 1, Υ = (n − 1)r ψ + 1, Υ = Υ Υ − r ψ ψ Υ Υ , 14 2 3 15 2 16 4 3 2 3 1 1 1 n−2 4 3 2 3 24(λ + 1)ψ Υ Υ 4(3λ − 2)n (−Υ ) 12(3λ + 2)r ψ ψ Υ 3 16 1 1 3 2 15 Υ = − + − , 2 2 3−n 2 3 Υ Υ Υ + r ψ ψ Υ 5 1 2 3 1 n−3 n+2 2 2 2 2 2 2n 4 4(3λ + 4)nr ψ r ψ − 3 4(3λ + 2)r ψ ψ Υ r ψ ψ (Υ Υ ) Υ − ψ Υ Υ + (n − 2)Υ Υ 2 2 3 2 2 1 8 3 7 6 1 3 1 1 1 Υ = + , 3 3 Υ Υ 1 5 2 2 2 2 2 n−2 8(9λ + 8)r ψ ψ Υ 6(5λ + 6)nr ψ (−Υ ) 14r ψ ψ Υ Υ 2 2 1 3 2 16 3 15 Υ = + ψ Υ + (3λ + 2) + 4 , 19 3 15 3−n 2 2 2 2 Υ + r ψ ψ Υ Υ Υ 2 3 1 1 5 2rψ ψ Υ (r−2λr) 2 3 15 3 5 3 8βn r ψ 4λ + 4 7 4 3 3−n 2 2 8β(3λ + 1)n r ψ 2nr ψ (−Υ ) 24βrψ ψ Υ Υ +r ψ ψ Υ 2 3 1 2 1 2 3 15 2 1 Υ = + + λ + , 4 3 2 3−n Υ Υ Υ Υ + r ψ ψ Υ 1 1 1 2 3 1 2 2 2 n−2 4 2 3 3 2(5λ + 1)n r ψ (−Υ ) 4(7λ + 4)ψ Υ Υ 24λr ψ ψ Υ 2 1 3 16 2 3 15 Υ = − − + , 4 2 3−n 2 3 Υ Υ Υ + r ψ ψ Υ 1 5 2 3 1 4 2 2 2λ (11n+24)r ψ +(48−11n)r ψ +24 ( 2 ) 2 n−2 2 2 2 2λr ψ ψ Υ Υ 3 2 16 r ψ ψ Υ + 24 1 2 2 3 15 4n (−Υ ) 8λ + + 4 1 2 Υ = − , 3−n 2 2 2 Υ + r ψ ψ Υ Υ 2 3 1 1 2 4 2 2 2 Υ = 4(λ + 1)nΥ r ψ − 3 + λ 24nr ψ − 24nr ψ Υ + 3nΥ 23 1 2 2 1 2 1 n+1 2 2 2n n+2 4 2ψ Υ r ψ ψ (Υ Υ ) Υ − ψ Υ Υ + (n − 2)Υ Υ 3 2 1 8 3 7 6 3 1 1 1 1 + , 5 17 n+1 7λψ Υ Υ 3 16 2 1 4 2 2 2 2r ψ − + 4n r ψ − 1 + 6λnr ψ r ψ − 3 2 2 2 2 Υ = − 10λ, n−2 2 2 2 (λ − 2)ψ Υ Υ 2(2λ − 1)n (−Υ ) (7λ − 3)r ψ ψ Υ 3 16 1 2 3 15 2 1 Υ = β 4r ψ − − + 25 2 2 2 3−n 2 2 Υ Υ Υ + r ψ ψ Υ 5 1 2 3 1 4(λ + 3)ψ Υ 3 15 + − 8λ + 2λr Υ , 3−n Υ + r ψ ψ Υ 2 3 1 2 2 2 2 2 4 3 2 3 4(λ − 2)nr ψ r ψ − 3 8(λ − 3)r ψ ψ Υ 12λr ψ ψ Υ 2 2 2 3 15 3 2 15 Υ = + − , 3−n 3 3−n 2 2 2 3 Υ + r ψ ψ Υ Υ Υ + r ψ ψ Υ 2 3 1 1 2 3 1 1 1 n−3 n+2 2 2 2 2n 4 4λr ψ ψ Υ r ψ ψ (Υ Υ ) Υ − ψ Υ Υ + (n − 2)Υ Υ 3 2 2 1 8 3 7 6 1 3 1 1 1 Υ = , 7λψ Υ Υ 2 3 16 2r ψ + (13λ − 4)n (−Υ ) 2 2 1 Υ = 68λ + + 32, 2n−4 4 6 3 4 2 2 2 4β(λ + 1)n r ψ 44β(3λ + 2)nr ψ ψ (−Υ ) Υ Υ 2 3 2 15 1 Υ = − − Υ 29 2 4 n 2 2 4 2 Υ (r ψ (ψ Υ + 2) + r ψ + 1) 2 3 1 1 2 2 2 n 4 2 2(λ + 1)nr (−Υ ) r ψ (ψ Υ + 2) + r ψ + 1 1 2 3 1 2 − , 2(2λ+3)r ψ ψ Υ 3 4 2 3 2 15 4βn r ψ − 2 3−n 2 2 Υ +r ψ ψ Υ 2 3 1 n−2 2 n−2 1 Υ = 4βψ Υ 14λ + (λ + 1)r ψ ψ Υ + 8 − , 30 3 2 3 1 1 n+2 n−2 2λψ Υ Υ 2(3λ + 2)r ψ ψ Υ Υ 3 16 3 2 16 2 2 2 1 2 1 2 Υ = (7λ + 5)n r ψ (−Υ ) + + 2n (−Υ ) Υ 4λ + + 2 + 2r ψ 31 2 1 1 2 2 1 2 Υ Υ 5 5 2 4 2 2 2 × 2(4λ + 3)nΥ r ψ − 3 − (3λ + 2) 24nr ψ − 24nr ψ Υ + 3nΥ 1 2 2 1 2 1 !! n+1 n+2 2 2 2n 4 2ψ Υ r ψ ψ (Υ Υ ) Υ − ψ Υ Υ + (n − 2)Υ Υ 3 2 1 8 3 7 6 1 3 1 1 1 + , n+1 7ψ Υ Υ 3 16 2 1 2 2 4βr ψ 2(10λ + 7)n (−Υ ) Υ + (3λ + 2) − 6nr ψ r ψ − 3 2 1 1 2 2 2 Υ = 4β(9λ + 7) + − Υ 32 14 × 12β(λ + 1) + (2λ + 1)r , n−2 2 2 2 (7λ + 6)ψ Υ Υ 2(2λ + 3)n (−Υ ) (21λ + 17)r ψ ψ Υ 3 16 1 2 2 3 15 Υ = −8β(3λ + 2) + 4βr ψ − − + 33 2 2 2 3−n 2 2 Υ Υ Υ + r ψ ψ Υ 5 1 2 3 1 4(3λ + 4)ψ Υ 2(λ + 1)r ψ ψ Υ 3 15 3 2 15 2 2 − + 2(λ + 1)r Υ , Υ = − 2 12β(λ + 1) + r , 14 34 3−n 3−n 2 2 Υ + r ψ ψ Υ Υ + r ψ ψ Υ 2 3 1 2 3 1 1 1 2 n−3 2 2 2n n+2 4 Υ = 4(3λ + 2)r ψ ψ Υ r ψ ψ (Υ Υ ) Υ − ψ Υ Υ + (n − 2)Υ Υ , 35 3 2 2 1 8 3 7 6 3 1 1 1 1   (3λ+2)ψ Υ Υ 2 3 16 14r ψ (5λ + 4)(−n) (−Υ ) − 2 1 2   Υ = ψ Υ − 4(11λ + 6) . 36 3 15 1 18 Appendix (II) 8 16 n 7 14 n 6 12 n 5 10 n 4 8 n 3 6 n λ = −8ψ nψ r Υ − 56ψ nψ r Υ − 168ψ nψ r Υ − 280ψ nψ r Υ − 280ψ nψ r Υ − 168ψ nψ r Υ 1 3 3 3 3 3 3 2 1 2 1 2 1 2 1 2 1 2 1 2 4 n 2 n 16 8 n 14 7 n 12 6 n 10 5 n 8 4 n −56ψ nψ r Υ − 8ψ nψ r Υ + r ψ ψ Υ + 8r ψ ψ Υ + 28r ψ ψ Υ + 56r ψ ψ Υ + 70r ψ ψ Υ 3 3 2 3 3 3 3 3 2 1 1 2 1 2 1 2 1 2 1 2 1 6 3 n 4 2 n 2 n n 2 7 14 2n 2 6 12 2n 2 5 10 2n +56r ψ ψ Υ + 28r ψ ψ Υ + 8r ψ ψ Υ + ψ Υ − 12ψ nψ r Υ − 60ψ nψ r Υ − 120ψ nψ r Υ 3 3 2 3 3 2 1 2 1 1 1 3 2 1 3 2 1 3 2 1 2 4 8 2n 2 3 6 2n 2 2 4 2n 14 7 2 2n 12 6 2 2n 10 5 2 2n −120ψ nψ r Υ − 60ψ nψ r Υ − 12ψ nψ r Υ + 4r ψ ψ Υ + 24r ψ ψ Υ + 60r ψ ψ Υ 3 2 1 3 2 1 3 2 1 2 3 1 2 3 1 2 3 1 8 4 2 2n 6 3 2 2n 4 2 2 2n 2 2 2n 3 6 12 3n 3 5 10 3n +80r ψ ψ Υ + 60r ψ ψ Υ + 24r ψ ψ Υ + 4r ψ ψ Υ − 8ψ nψ r Υ − 24ψ nψ r Υ 2 3 1 2 3 1 2 3 1 3 1 3 2 1 3 2 1 3 4 8 3n 3 3 6 3n 12 6 3 3n 10 5 3 3n 8 4 3 3n 6 3 3 3n −24ψ nψ r Υ − 8ψ nψ r Υ + 6r ψ ψ Υ + 24r ψ ψ Υ + 36r ψ ψ Υ + 24r ψ ψ Υ 3 2 1 3 2 1 2 3 1 2 3 1 2 3 1 2 3 1 4 2 3 3n 6 3 4 4n 9 18 8 16 7 14 6 12 5 10 4 8 +6r ψ ψ Υ + 4r ψ ψ Υ − 2nψ r − 18nψ r − 72nψ r − 168nψ r − 252nψ r − 252nψ r 2 3 1 2 3 1 2 2 2 2 2 2 3 6 2 4 2 −168nψ r − 72nψ r − 18nψ r − 2n, 2 2 4 16 9 4 14 8 4 12 7 4 10 6 4 8 5 4 6 4 4 4 3 3 16 9 λ = 2βn r ψ + 12βn r ψ + 30βn r ψ + 40βn r ψ + 30βn r ψ + 12βn r ψ + 2βn r ψ + 4βn r ψ 2 2 2 2 2 2 2 2 3 14 8 3 12 7 3 10 6 3 8 5 3 6 4 3 4 3 3 4 3 n +20βn r ψ + 36βn r ψ + 20βn r ψ − 20βn r ψ − 36βn r ψ − 20βn r ψ + 76βn r ψ ψ Υ 2 2 2 2 2 2 2 1 3 2 2 2 2 16 9 2 14 8 2 12 7 2 10 6 2 8 5 2 6 −4βn r ψ − 22βn ψ + 18βn r ψ + 128βn r ψ + 368βn r ψ + 528βn r ψ + 340βn r ψ − 32βn r 2 2 2 2 2 2 4 2 4 3 2 4 3 n 2 2 2 2 2 2 n n ×ψ − 192βn r ψ + 660βn r ψ ψ Υ − 112βn r ψ + 130βn r ψ ψ Υ + 80βnψ − 160βψ ψ Υ + 144βnψ 3 3 2 3 2 2 2 2 2 1 2 2 1 1 n 16 9 14 8 12 7 10 6 10 4 5 4n 10 5 4 4n ×ψ Υ + 16βnr ψ + 192βnr ψ + 896βnr ψ + 2240βnr ψ − 2nr ψ ψ Υ + 4r ψ ψ Υ + 3360 1 2 2 2 2 3 2 1 2 3 1 8 5 8 4 4 4n 8 4 4 4n 8 4 5 5n 6 4 6 4 n 4 3 ×βnr ψ − 2nr ψ ψ Υ + 8r ψ ψ Υ + r ψ ψ Υ + 3136βnr ψ − 2720βr ψ ψ Υ + 1792βnr ψ 2 3 2 1 2 3 1 2 3 1 2 2 1 2 4 3 n 4 3 n 2 2 2 2 n 2 2 n −2208βr ψ ψ Υ + 1840βnr ψ ψ Υ + 576βnr ψ − 928βr ψ ψ Υ + 800βnr ψ ψ Υ 3 3 3 3 2 1 2 1 2 2 1 2 1 6 10 n 5 8 n 4 14 8 n 3 14 8 n 2 14 8 n 14 λ = −480ψ ψ r βΥ − 1760ψ ψ r βΥ + 18n r βψ ψ Υ + 36n r βψ ψ Υ + 18n r βψ ψ Υ + 32r β 3 3 3 3 3 3 2 1 2 1 2 1 2 1 2 1 8 n 4 12 7 n 3 12 7 n 2 12 7 n 12 7 n 12 7 n ×ψ ψ Υ + 72n r βψ ψ Υ + 220n r βψ ψ Υ + 212n r βψ ψ Υ + 80nr βψ ψ Υ + 32r βψ ψ Υ 3 3 3 3 3 3 2 1 2 1 2 1 2 1 2 1 2 1 4 10 6 n 3 10 6 n 2 10 6 n 10 6 n 4 8 5 n +108n r βψ ψ Υ + 520n r βψ ψ Υ + 798n r βψ ψ Υ + 544nr βψ ψ Υ + 72n r βψ ψ Υ + 600 3 3 3 3 3 2 1 2 1 2 1 2 1 2 1 3 8 5 n 2 8 5 n 8 5 n 4 6 4 n 3 6 4 n 2 6 ×n r βψ ψ Υ + 1432n r βψ ψ Υ + 1520nr βψ ψ Υ + 18n r βψ ψ Υ + 340n r βψ ψ Υ + 1358n r β 3 3 3 3 3 2 1 2 1 2 1 2 1 2 1 4 n 6 4 n 4 8 5 2 2n 3 8 5 2 2n 2 8 5 2 2n 8 5 ×ψ ψ Υ + 2240nr βψ ψ Υ + 16n r βψ ψ Υ + 192n r βψ ψ Υ + 440n r βψ ψ Υ + 240nr βψ 3 3 2 1 2 1 2 3 1 2 3 1 2 3 1 2 2 2n 3 6 4 2 2n 2 6 4 2 2n 6 4 2 2n 2 4 3 2 2n 4 3 2 2n ×ψ Υ + 80n r βψ ψ Υ + 424n r βψ ψ Υ + 464nr βψ ψ Υ + 144n r βψ ψ Υ + 392nr βψ ψ Υ 3 1 2 3 1 2 3 1 2 3 1 2 3 1 2 3 1 2 2 2 2n 2 6 10 2n 2 5 8 2n 2 4 6 2n 2 3 4 2n +120nr βψ ψ Υ − 284ψ ψ r βΥ − 1110ψ ψ r βΥ − 1800ψ ψ r βΥ − 1430ψ ψ r βΥ 2 3 1 3 2 1 3 2 1 3 2 1 3 2 1 2 2 2 2n 2 2n −540ψ ψ r βΥ − 74ψ ψ βΥ 3 2 1 3 1 4 12 7 2 2n 4 10 6 2 2n 3 12 7 2 2n 3 10 6 2 2n 2 12 7 2 2n λ = 16βn r ψ ψ Υ + 32βn r ψ ψ Υ + 32βn r ψ ψ Υ + 144βn r ψ ψ Υ + 24βn r ψ ψ Υ + 184 2 3 1 2 3 1 2 3 1 2 3 1 2 3 1 2 10 6 2 2n 2 10 6 3 3n 2 8 5 3 3n 2 6 3 4 3n 12 7 2 2n 12 ×βn r ψ ψ Υ + 24βn r ψ ψ Υ + 16βn r ψ ψ Υ − 8βn r ψ ψ Υ + 8βnr ψ ψ Υ − 10βr 2 3 1 2 3 1 2 3 1 3 2 1 2 3 1 2 7 2n 10 6 2 2n 10 6 3 3n 10 3 6 3n 8 5 3 3n 8 3 5 3n ×ψ ψ Υ + 56βnr ψ ψ Υ + 24βnr ψ ψ Υ − 62βr ψ ψ Υ + 56βnr ψ ψ Υ − 248βr ψ ψ Υ 3 2 1 2 3 1 2 3 1 3 2 1 2 3 1 3 2 1 8 4 5 4n 6 4 3 3n 6 3 4 3n 6 4 4 4n 6 5 4 5n 4 3 3 3n −22βr ψ ψ Υ + 8βnr ψ ψ Υ − 404βr ψ ψ Υ − 44βr ψ ψ Υ − 2βr ψ ψ Υ − 312βr ψ ψ Υ 3 2 1 2 3 1 3 2 1 3 2 1 3 2 1 3 2 1 4 3 3 3n 4 4 3 4n 2 3 2 3n −24βnr ψ ψ Υ − 22βr ψ ψ Υ − 94βr ψ ψ Υ 3 2 1 3 2 1 3 2 1 2 8 16 n 8 16 n 2 7 14 n 7 14 n 2 6 12 n 6 12 n λ = −3ψ n ψ r Υ − 2ψ ψ r Υ − 18ψ n ψ r Υ − 16ψ ψ r Υ − 45ψ n ψ r Υ − 56ψ ψ r Υ 5 3 3 3 3 3 3 2 1 2 1 2 1 2 1 2 1 2 1 2 5 10 n 5 10 n 2 4 8 n 4 8 n 2 3 6 n 3 6 n −60ψ n ψ r Υ − 112ψ ψ r Υ − 45ψ n ψ r Υ − 140ψ ψ r Υ − 18ψ n ψ r Υ − 112ψ ψ r Υ 3 3 3 3 3 3 2 1 2 1 2 1 2 1 2 1 2 1 2 2 4 n 2 4 n 2 n 16 8 n 14 7 n 12 6 n −3ψ n ψ r Υ − 56ψ ψ r Υ − 16ψ ψ r Υ + 11nr ψ ψ Υ + 71nr ψ ψ Υ + 195nr ψ ψ Υ + 295 3 3 3 2 3 3 3 2 1 2 1 1 2 1 2 1 2 1 10 5 n 8 4 n 6 3 n 4 2 n 2 n n 2 5 10 2n ×nr ψ ψ Υ + 265nr ψ ψ Υ + 141nr ψ ψ Υ + 41nr ψ ψ Υ + 5nr ψ ψ Υ − 2ψ Υ − 120ψ ψ r Υ 3 3 3 3 2 3 3 2 1 2 1 2 1 2 1 1 1 3 2 1 2 2 4 8 2n 2 4 8 2n 2 2 3 6 2n 2 3 6 2n 2 2 4 2n 2 2 2n −12ψ n ψ r Υ − 160ψ ψ r Υ − 3ψ n ψ r Υ − 120ψ ψ r Υ − 48ψ ψ r Υ − 8ψ ψ r Υ + 114 3 2 1 3 2 1 3 2 1 3 2 1 3 2 1 3 1 8 4 2 2n 6 3 2 2n 4 2 2 2n 2 9 18 2 8 16 2 7 14 2 6 12 2 5 ×nr ψ ψ Υ + 51nr ψ ψ Υ + 9nr ψ ψ Υ − n ψ r − 8n ψ r − 28n ψ r − 56n ψ r − 70n ψ 2 3 1 2 3 1 2 3 1 2 2 2 2 2 10 18 9 2 4 8 16 8 14 7 2 3 6 12 6 10 5 2 2 4 ×r + 3nr ψ − 56n ψ r + 25nr ψ + 92nr ψ − 28n ψ r + 196nr ψ + 266nr ψ − 8n ψ r 2 2 2 2 2 2 2 2 8 4 6 3 2 2 4 2 2 +238nr ψ + 140nr ψ − n ψ r + 52nr ψ + n + 11nr ψ 2 2 2 2 2 19 n n 2 2 7 14 2n 2 7 14 2n 2 2 6 12 2n 2 6 12 2n λ = −180ψ nψ βΥ + 200βψ ψ Υ − 3ψ n ψ r Υ − 8ψ ψ r Υ − 12ψ n ψ r Υ − 48ψ ψ r Υ 6 3 2 2 3 1 1 3 2 1 3 2 1 3 2 1 3 2 1 2 2 5 10 2n 14 7 2 2n 12 6 2 2n 10 5 2 2n 3 6 12 3n 3 2 6 12 3n −18ψ n ψ r Υ + 15nr ψ ψ Υ + 69nr ψ ψ Υ + 126nr ψ ψ Υ − 12ψ ψ r Υ − ψ n ψ r Υ 3 2 1 2 3 1 2 3 1 2 3 1 3 2 1 3 2 1 3 5 10 3n 3 2 5 10 3n 3 4 8 3n 3 2 4 8 3n 3 3 6 3n 3 2 4 3n 12 −48ψ ψ r Υ − 2ψ n ψ r Υ − 72ψ ψ r Υ − ψ n ψ r Υ − 48ψ ψ r Υ − 12ψ ψ r Υ + 9nr 3 2 1 3 2 1 3 2 1 3 2 1 3 2 1 3 2 1 6 3 3n 10 5 3 3n 8 4 3 3n 6 3 3 3n 4 5 10 4n 4 4 8 4n 4 3 6 4n ×ψ ψ Υ + 25nr ψ ψ Υ + 23nr ψ ψ Υ + 7nr ψ ψ Υ − 8ψ ψ r Υ − 16ψ ψ r Υ − 8ψ ψ r Υ 2 3 1 2 3 1 2 3 1 2 3 1 3 2 1 3 2 1 3 2 1 10 5 4 4n 8 4 4 4n 5 4 8 5n 2 10 6 2 8 5 2 6 4 2 4 3 +2nr ψ ψ Υ + 2nr ψ ψ Υ − 2ψ ψ r Υ + 84n r βψ + 560n r βψ + 700n r βψ + 420n r βψ 2 3 1 2 3 1 3 2 1 2 2 2 2 2 2 2 4 9 16 2 9 16 9 16 4 8 14 3 8 14 2 8 14 +124n r βψ − 6n ψ r β − 30n ψ r β − 36nψ r β − 36n ψ r β − 16n ψ r β − 140n ψ r β − 448n 2 2 2 2 2 2 2 8 14 4 7 12 3 7 12 2 7 12 7 12 4 6 10 3 6 10 ×ψ r β − 90n ψ r β − 96n ψ r β − 196n ψ r β − 2032nψ r β − 120n ψ r β − 240n ψ r β − 4800n 2 2 2 2 2 2 2 6 10 4 5 8 3 5 8 5 8 4 4 6 3 4 6 4 6 4 3 4 ×ψ r β − 90n ψ r β − 320n ψ r β − 6680nψ r β − 36n ψ r β − 240n ψ r β − 5696nψ r β − 6n ψ r β 2 2 2 2 2 2 2 2 3 3 4 3 4 3 2 2 2 2 2 −96n ψ r β − 2928nψ r β − 16n ψ r β − 832nψ r β − 100nψ β + 14n βψ 2 2 2 2 2 2 3 8 14 n 2 8 14 n 8 14 n 3 7 12 n 2 7 12 n λ = −104ψ n ψ r βΥ − 90ψ n ψ r βΥ − 88ψ ψ r βΥ − 376ψ n ψ r βΥ − 840ψ n ψ r βΥ 7 3 3 3 3 3 2 1 2 1 2 1 2 1 2 1 7 12 n 3 6 10 n 2 6 10 n 6 10 n 3 5 8 n −12ψ nψ r βΥ − 464ψ n ψ r βΥ − 2426ψ n ψ r βΥ − 1848ψ nψ r βΥ − 176ψ n ψ r βΥ − 3104 3 3 3 3 3 2 1 2 1 2 1 2 1 2 1 2 5 8 n 5 8 n 2 4 6 n 4 6 n 2 3 4 n ×ψ n ψ r βΥ − 5692ψ nψ r βΥ − 1806ψ n ψ r βΥ − 7608ψ nψ r βΥ − 344ψ n ψ r βΥ − 5124ψ n 3 3 3 3 3 3 2 1 2 1 2 1 2 1 2 1 3 4 n 2 2 n 4 14 8 n 14 8 n 4 12 7 n 12 7 n ×ψ r βΥ − 1640ψ nψ r βΥ + 2n r βψ ψ Υ + 88nr βψ ψ Υ + 8n r βψ ψ Υ + 56r βψ ψ Υ 3 3 3 3 3 2 1 2 1 2 1 2 1 2 1 2 1 4 10 6 n 10 6 n 4 8 5 n 8 5 n 4 6 4 n 3 6 4 +12n r βψ ψ Υ + 1800r βψ ψ Υ + 8n r βψ ψ Υ + 5080r βψ ψ Υ + 2n r βψ ψ Υ + 56n r βψ ψ 3 3 3 3 3 3 2 1 2 1 2 1 2 1 2 1 2 n 6 4 n 3 4 3 n 4 3 n 2 2 2 n 2 2 n 4 ×Υ + 6520r βψ ψ Υ + 40n r βψ ψ Υ + 4392r βψ ψ Υ + 34n r βψ ψ Υ + 1496r βψ ψ Υ + 3454r 3 3 3 3 3 1 2 1 2 1 2 1 2 1 2 1 3 2 2n 2 2 2 2n 2 2n 2 2 3 4 2n 2 3 4 2n 2 2 2 2n ×βψ ψ Υ + 1220r βψ ψ Υ + 106βψ ψ Υ − 468ψ n ψ r βΥ − 1408ψ nψ r βΥ − 520ψ nψ r βΥ 2 3 1 2 3 1 3 1 3 2 1 3 2 1 3 2 1 2 12 7 2 2n 12 7 2 2n 10 6 2 2n 8 5 2 2n 6 4 2 2n 2 4 λ = 60n r βψ ψ Υ + 114r βψ ψ Υ + 772r βψ ψ Υ + 2534r βψ ψ Υ + 4216r βψ ψ Υ − 48ψ n 2 3 1 2 3 1 2 3 1 2 3 1 2 3 1 3 7 12 2n 2 3 7 12 2n 2 7 12 2n 2 4 6 10 2n 2 3 6 10 2n 2 2 ×ψ r βΥ − 48ψ n ψ r βΥ − 216ψ nψ r βΥ − 96ψ n ψ r βΥ − 320ψ n ψ r βΥ − 96ψ n 2 1 3 2 1 3 2 1 3 2 1 3 2 1 3 6 10 2n 2 6 10 2n 2 4 5 8 2n 2 3 5 8 2n 2 2 5 8 2n 2 ×ψ r βΥ − 560ψ nψ r βΥ − 48ψ n ψ r βΥ − 496ψ n ψ r βΥ − 840ψ n ψ r βΥ − 840ψ n 2 1 3 2 1 3 2 1 3 2 1 3 2 1 3 5 8 2n 2 3 4 6 2n 2 2 4 6 2n 2 4 6 2n 10 6 3 3n 8 ×ψ r βΥ − 224ψ n ψ r βΥ − 1152ψ n ψ r βΥ − 1384ψ nψ r βΥ + 126r βψ ψ Υ + 528r β 2 1 3 2 1 3 2 1 3 2 1 2 3 1 5 3 3n 2 6 4 3 3n 6 4 3 3n 4 3 3 3n 4 3 3 3n 2 2 3 3n ×ψ ψ Υ + 16n r βψ ψ Υ + 924r βψ ψ Υ + 52nr βψ ψ Υ + 800r βψ ψ Υ + 278r βψ ψ Υ 2 3 1 2 3 1 2 3 1 2 3 1 2 3 1 2 3 1 3 2 6 10 3n 3 6 10 3n 3 2 5 8 3n 3 5 8 3n 3 4 6 3n 8 −48ψ n ψ r βΥ − 60ψ nψ r βΥ − 32ψ n ψ r βΥ − 164ψ nψ r βΥ − 52ψ nψ r βΥ + 42r β 3 2 1 3 2 1 3 2 1 3 2 1 3 2 1 5 4 4n 6 4 4 4n 4 3 4 4n 6 4 5 5n ×ψ ψ Υ + 108r βψ ψ Υ + 66r βψ ψ Υ + 6r βψ ψ Υ 2 3 1 2 3 1 2 3 1 2 3 1 References [1] Baade and Zwicky., PNAS 112, 1241 (2015). [2] K.R. Karmarkar, Proc. Indian Acad. Sci. A 27, 56 (1948). [3] L. Schlai, Ann. Mat. 5, 170 (1871). [4] J. Nash, Ann. of Math. 63, 20 (1956). [5] S.K. Maurya et al., Eur. Phys. J. C 75, 389 (2015). [6] S.K. Maurya et al., Eur. Phys. J. A 52 191 (2016). [7] P. Bhar et al., Eur. Phys. J. A 52, 312 (2016). [8] D. Deb, S.V. Ketov, S.K. Maurya, Mon. Not. R. Astron. Soc. 485, 5652 (2019). [9] S.K. Maurya, et al., Phys. Rev. D 100, (2019). [10] S.K. Maurya et al. , Ann. Physics 385, 532 (2017). [11] S.K. Maurya et al., Eur. Phys. J. C 76, 266 (2016). [12] S.K. Maurya et al., Eur. Phys. J. C 77, 1 (2016). [13] T. Harko et al., Phys. Rev. D 84, 024020 (2011). [14] M. Jamil, D. Momeni and R. Myrzakulov., Eur. Phys. J. C. 72, 1959 (2012). [15] M. Jamil, D. Momeni and R. Myrzakulov., Chin. Phys. Lett. 29, 109801 (2012). [16] H. Shabani and M. Farhoudi., Chin. Phys. Lett. 88, 044048 (2013). [17] H. Shabani and M. Farhoudi., Phys. Rev. D 90, 44031 (2014). [18] P.H.R.S. Moraes ., Eur. Phys. J. C 75, 168 (2015). [19] A. Alhamzawi and R. Alhamzawi., Int. J. Mod. Phys. D 25, 1650020 (2016). [20] P.H.R.S. Moraes, J.D. Arbanil and M. Malheiro., J. Cosmol. Astropart. Phys. 6, 5 (2016). [21] Das et al., Eur. Phys. J. C 76, 654 (2016). 20 [22] Moraes, al., arXiv:1806.04123v4. [23] Z. Yousaf, M.Z. Bhatti and M. Ilyas ., Eur. Phys. J. C 78, 307 (2018). [24] S.K.Mauryaa and F. T. Ortizb., Phys. Dark Universe 27, 100442 (2020). [25] S. Waheed., Symmetry 12, 962 (2020). [26] G. Mustafa et al., Eur. Phys. J. C 80, 26 (2020). [27] M. Ruderman., Annu Rev Astron Astr, 10, 427 (1972). [28] A. V. Astashenok, S. Capozziello and S. D. Odintsov., J. Cosmol. Astropart. Phys. 01, 001 (2015). [29] A. V. Astashenok, S. Capozziello and S. D. Odintsov., Phys. Lett. B 742, 160 (2015). [30] D. Momeni, P. H. R. S. Moraes and R. Myrzakulov., Astrophys. Space Sci 361, 228 (2018). [31] D. Momeni, M. Raza and R. Myrzakulov., Mod. Phys. Lett. A 31, 1650073 (2015). [32] S. Capozziello et al., Phys. Rev. D 93, 023501 (2016). [33] A. V. Astashenok, S. Capozziello and S. D. Odintsov., Phys. Rev. D 89, 103509 (2014). [34] L. Herrera, Phys. Lett. A 165, 206 (1992). [35] A. A. Starobinsky, Phys. Lett. B 91, 99 (1980). [36] M. Jamil et al., Eur. Phys. J. C 72, 1999 (2012). [37] H. Shabani and M. Farhoudi, Phys. Rev. D 88, 044048 (2013). [38] C.P. Singh and P. Kumar, Eur. Phys. J. C 74, 3070 (2014). [39] M. Sharif, Z. Yousaf, Astrophys. Space Sci. 354, 471 (2014). [40] I. Noureen and M. Zubair, Astrophys. Space Sci. 356, 103 (2015). [41] I. Noureen et al., Eur. Phys. J. C 75, 323 (2015). [42] S. Chandrasekhar, Astrophys. J. 140, 417 (1964). [43] H. Heintzmann and W. Hillebrandt, Astron. Astrophys. 38, 51 (1975). [44] W. Hillebrandt and K. O. Steinmetz, Astron. Astrophys. 53, 283 (1976). [45] D. Horvat, S. Ilijic and A. Marunovic, Class. Quantum Grav. 28, 025009 (2011). [46] D. D. Doneva and S. S. Yazadjiev, Phys. Rev. D 85, 124023 (2012). [47] H. O. Silva, C. F. B. Macedo, E. Berti and L. C. B. Crispino, Class. Quant. Grav. 32, 145008 (2015). [48] I. Bombaci, Astron. Astrophys. 305, 871 (1996). [49] M. Zubair, G. Abbas, and I. Noureen, Astrophys. Space Sci. 361,8 (2016). [50] A. K. Yadav, M. Mondal, and F. Rahaman, Pramana - J Phys 94, 90 (2020).

Journal

General Relativity and Quantum Cosmology – arXiv (Cornell University)

Published: Jan 1, 2021

There are no references for this article.