Next Article in Journal
Magic Three-Qubit Veldkamp Line and Veldkamp Space of the Doily
Previous Article in Journal
Effect of Nano-Graphene Oxide and n-Butanol Fuel Additives Blended with Diesel—Nigella sativa Biodiesel Fuel Emulsion on Diesel Engine Characteristics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Physically Acceptable Embedded Class-I Compact Stars in Modified Gravity with Karmarkar Condition

1
Faculty of CORE Humanities, Prince Mohammad Bin Fahd University, Al Khobar 31952, Saudi Arabia
2
Department of Mathematics, Shanghai University, Shanghai 200444, China
3
Department of Mathematics, COMSATS University Islamabad, Lahore Campus, Punjab 5459, Pakistan
4
Jiangsu Key Labortary for NSLSCS, School of Mathematical Sciences, Nanjing Normal University, Nanjing 210023, China
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(6), 962; https://doi.org/10.3390/sym12060962
Submission received: 27 March 2020 / Revised: 30 April 2020 / Accepted: 7 May 2020 / Published: 5 June 2020

Abstract

:
The present study is devoted to explore the existence of a new family of compact star solutions by adopting the Karmarkar as well as Pandey–Sharma condition in the background of f ( R , T ) modified gravitational framework. For this purpose, we consider static spherically symmetric spacetime with anisotropic fluid distribution in absence of electric charge. In respect of Karmarkar condition, we assume a specific model of g r r metric potential representing a new family of solutions which is also compatible with the Pandey–Sharma condition. This assumed model permits us to calculate the g t t component of metric tensor by making the use of Karmarkar condition. Further, we investigate the interior solutions for V e l a X 1 model of compact star by utilizing this new family of solutions for different values of parameter λ . We have tuned the solution for V e l a X 1 so that the solutions matches the observed mass and radius. For the same star we have extensively discussed the behavior of the solutions. It is found that these solutions fulfill all the necessary conditions under the observational radii and mass attribute data for small values of parameter λ and hence physically well-behaved and promising. Through graphical analysis, it is observed that our obtained analytical solutions are physically acceptable with a best degree of accuracy for n [ 1.8 , 7 ) { 2 , 4 , 6 } , where parameter n is involved in the discussed model. It is also noticed the causality condition is violated for all n 7 and the tangential sound velocity v t is observed as complex valued for all 0 < n < 1.8 . Likewise, we explore these properties by considering large parameter λ values. It is seen that the presented model violates all the physical conditions for n { 2 , 4 , 6 } , while some of these for large values of λ . Consequently, it can be concluded that the parameters n and λ have a strong impact on the obtained solutions.

1. Introduction

The study of relativistic compact stars in the background of general relativity (GR) and extended theories of gravity is considered as one of the most promising subjects among the researchers during last few decades. The exact solution of Einstein’s field equations (EFEs) for the first time was calculated by Schwarzschild [1] in 1916 for the interior of a compact object. Later various researchers made efforts to obtain the new exact solutions of EFEs in reference to different astrophysical and cosmological aspects. Tolman [2] and Oppenheimer [3] both studied some more realistic and convoluted models of non-traversable relativistic objects in proximity with data-based actuality. In 1939, Oppenheimer and Snyder [4] examined the phenomenon of gravitational collapse with homogeneity based dust sphere and they presented an interesting and simple collapse mode in this respect. In 1922, Jeans [5] introduced the anisotropic pressure in self gravitational bodies in the context of Newtonian gravity. Further, in 1933, Lemaitre [6] used the anisotropic model for the first time in the background of GR by taking constant energy density with a radial pressure source. Later on, Bowers and Liang [7], in 1974, employed anisotropic spheres in GR. Ruderman [8] argued that the existence of nonzero anisotropy is regarded as a significant component of a stellar object without an electric charge. Gleiser and Dev [9] argued that pressure anisotropy impacts the structure, mass and other physical properties of extremely compact spheres. Herrera and Santos [10] described the effect of anisotropy on the characteristics of the stellar stars. The MIT bag model arising from observations in fundamental particle physics relates the pressure to the density of the star via a linear relation of the form p = α ρ B where B is the Bag constant. This equation of state has been successfully used to model compact objects in general relativity ranging from neutron stars through to strange star candidates.
The most interesting fact of modern cosmology is the accelerated expanding nature of our cosmos which is also supported by the suggestions coming from different recent observational experiments like LSS, CMB, BAO and Supernovae-Ia [11,12,13,14]. It is argued that the standard GR theory fails to describe the hidden exotic nature of its responsible factor labeled as dark energy (DE). It is realized that exotic cosmological astral structures may be determined by reviewing the relativistic astral objects in the framework of modified theory of gravity. It is also argued that recent data from solar system tests of gravity might also be compatible with some other theories of gravity [15,16]. In the literature, many alternative theories of gravity are available including f ( R ) theory, f ( T ) gravity, f ( G ) theory, f ( T , T G ) modified framework, braneworld scenario, Gauss–Bonnet gravity, Kalb–Ramond background and scalar-tensor theories especially Brans–Dicke gravity [16,17,18,19,20,21,22,23]. It is argued that GR theory faced stringent constraints for describing the galactic scale, the strong lensing regime and CMBR tests. In this respect, its recently proposed extended version namely f ( R ) theory is considered as quite successful and interesting because of its correspondence to scalar-tensor gravity, i.e., the GR with an extra scalar-field candidate, whose scalar-potential is in double correspondence to the generic function f [24].
Recently, another interesting modification of f ( R ) theory has been proposed in the literature [25] by taking an arbitrary function into account namely f ( R , T ) , where T represents the trace of energy-momentum matter of the ordinary matter lagrangian part of the action. This modification is acknowledged as f ( R , T ) modified theory. The f ( R , T ) theory has drawn much attention of the investigators and many interesting results have been found in the cosmological aspects [26,27,28,29,30]. Sharif et al. [31] explored f ( R , T ) theory of gravity to investigate some features related to the stability of a locally isotropic self-gravitating system. Zubair et al. [32] discussed the possible formation of compact objects by employing the Krori and Barua metric in the background of f ( R , T ) theory. Alhamzawi and Alhamzawi [33] described the phenomenon of gravitational lensing in the framework of f ( R , T ) theory. The topic of stability has been investigated numerically by Moraes et al. [34] using TOV equation in f ( R , T ) theory. Das et al. [35] calculated an exact solution of compact objects using modified field equations by taking the conformal symmetry motion into account in the background of f ( R , T ) theory.
In the literature, numerous researchers have recommended an embedding spacetime metric for the basic structure of compact stars which is consistent with Karmarkar condition defined by R 1414 = ( R 1212 R 3434 + R 1224 R 1334 ) / R 2323 [36]. The Karmarkar condition provides a connection between the two main parts of metric potentials of a spherically symmetric spacetime. By adopting Karmarkar condition, one of the metric components can be calculated which is further capable of generating the second metric component [37,38,39]. In this respect, various families of solutions are generated by Bhar et al. [40,41], Maurya et al. [42,43,44] and Newton Singh et al. [45,46,47,48] which are physically acceptable for embedding spacetimes and can be used to describe the interior solution of stellar objects. These new families of solutions satisfying the EFEs along with Karmarkar condition are termed as embedding class-I solutions. Further, many different stellar models have been assumed to investigate the new well-behaved classes of solutions by considering the connection of metric potentials in embedding class-I solutions [49,50,51]. Recently, Mustafa et al. [51,52] have explored the realistic stellar anisotropic model with well behaved nature solutions satisfying Karmarker condition in f ( R , T ) gravity.
The present study investigates the existence of spherical symmetric solutions for V e l a X 1 model of compact stars in the framework of f ( R , T ) theory of gravity. For this purpose, we shall adopt the Karmarkar as well as Pandey Sharma condition for simplification purposes. The paper is designed in the following order of sections. In the next section, we shall present the mathematical formulations of modified f ( R , T ) theory of gravity along with its field equations for a spherically symmetric geometry with anisotropic matter source. Here we shall discuss a new configuration of embedding class-I solutions in f ( R , T ) theory and the physical boundary conditions in the same section. Section 3 will provide a detailed graphical analysis of some physical properties and the stability of V e l a X 1 model of compact stars for different values of λ . Final section concludes the whole work by focusing on major findings.

2. Basic Structure of f ( R , T ) Gravitational Framework

In this part, we shall explain the essential structure of modified f ( R , T ) gravitational framework with all the necessary assumptions taken for this work. The action of f ( R , T ) modified theory is defined by [24]:
S = 1 16 π f ( R , T ) g d 4 x + L m g d 4 x ,
where L m symbolizes the ordinary matter Lagrangian density. By taking variation of the above action in Equation (1), we get the following set of dynamical equations
8 π T μ ν f T ( R , T ) T μ ν f T ( R , T ) Θ μ ν = 1 2 f ( R , T ) g μ ν + f R ( R , T ) R μ ν + ( g μ ν μ ν ) f R ( R , T ) .
In this set of field equations, the notations ∇ and □ are, respectively introduced for the covariant derivative and the d’Alembert operator defined by = μ ( g g μ ν ν ) / g . Further, the function f with subscripts refer to the derivatives given by f R ( R , T ) = f ( R , T ) / R and f T ( R , T ) = f ( R , T ) / T . In addition, the term Θ μ ν is defined by
Θ μ ν = g α β δ T μ ν δ g μ ν = 2 T μ ν + g μ ν L m 2 g α β 2 L m g μ ν g α β .
For this work, we assume the matter contents as anisotropic fluid distribution defined by the following energy-momentum tensor
T μ ν = ( ρ + p t ) υ μ υ ν p t g μ ν + ( p r p t ) ξ μ ξ ν ,
where υ μ and ξ μ represent the 4-velocity vectors defined by υ μ = e μ δ 0 μ and ξ μ = e ν δ 1 μ which satisfy the relations υ μ υ μ = ξ μ ξ μ = 1 . Here ρ , p r and p t are the notations for ordinary matter density, radial and tangential pressures, respectively.
The set of field Equation (2) can be further arranged to the following form:
G μ ν = 1 f R ( R , T ) ( 8 π + f T ( R , T ) ) T μ ν + ( μ ν g μ ν ) f R ( R , T ) 1 2 ( f ( R , T ) R f R ( R , T ) ) g μ ν ρ g μ ν f T ( R , T ) .
The geometry of an object exhibiting static and spherically symmetry properties is defined by the line element given by
d s 2 = e μ ( r ) d t 2 + e ν ( r ) d r 2 + r 2 d Ω 2 ,
where d Ω 2 = d θ 2 + s i n 2 θ d ϕ 2 . Further, we consider an interesting and simple linear model of f ( R , T ) modified theory which is defined by
f ( R , T ) = R + λ T .
Utilizing the metric components of Equation (5) in Equation (4), the expressions of energy density function and pressure components take the following form
8 π ρ = e ν ( r ) 4 ( λ + 1 ) ( 2 λ + 1 ) r 2 4 λ + 4 λ r 2 μ ( r ) + ( 2 λ + 1 ) r 2 μ ( r ) 2 ( λ + 1 ) 2 r 2 μ ( r ) r 2 μ ( r ) ν ( r ) + r 2 μ ( r ) 2 + 4 r μ ( r ) 4 r ν ( r ) 4 e ν ( r ) + 4 + 2 r 2 μ ( r ) ( 2 λ + 1 ) r μ ( r ) r ν ( r ) 4 + 4 λ r ν ( r ) + 4 λ e ν ( r ) ,
8 π p r = e ν ( r ) 4 2 λ 2 + 3 λ + 1 r 2 2 ( 2 λ + 1 ) r 2 μ ( r ) + 2 λ e ν ( r ) 1 + ( 2 λ + 1 ) r 2 μ ( r ) ν ( r ) ( 2 λ + 1 ) r 2 μ ( r ) 2 + ( λ + 1 ) 2 r 2 μ ( r ) r 2 μ ( r ) ν ( r ) + r 2 μ ( r ) 2 + 4 r μ ( r ) 4 r ν ( r ) 4 e ν ( r ) + 4 + 4 ( λ + 1 ) r ν ( r ) ,
8 π p t = e ν ( r ) 4 2 λ 2 + 3 λ + 1 r 2 ( λ + 1 ) 2 r 2 μ ( r ) r 2 μ ( r ) ν ( r ) + r 2 μ ( r ) 2 + 4 r μ ( r ) 4 r ν ( r ) 4 e ν ( r ) + 4 + 2 ( 2 λ r + r ) μ ( r ) + 2 ( λ + 1 ) e ν ( r ) 1 + r ν ( r ) ,
8 π ( p t p r ) = e ν ( r ) 4 ( λ + 1 ) r 2 r 2 r μ ( r ) + μ ( r ) r μ ( r ) r ν ( r ) 2 2 ν ( r ) + 4 e ν ( r ) 1 .
Here, the expression 8 π ( p t p r ) defines the measurement of anisotropy, and is represented by ▵, i.e., = 8 π ( p t p r ) . The isotopy of the configuration can be achieved in the limit of zero anisotropy factor, i.e., = 0 .

2.1. The Karmarker Condition

Here we shall define the well-known Karmarker condition and the resulting expression of unknown functions. The Karmarkar condition is defined by the relation
R 1414 = R 1212 R 3434 + R 1224 R 1334 R 2323 ,
where R 2323 0 . In this respect, the Pandey–Sharma criterion, i.e., R 2323 0 defines a new spacetime with an embedding class-I. Here by correlating the space-time and Karmarker as well as Pande–Sharma condition, we get the following differential equation in terms of metric potentials
ν ( r ) μ ( r ) + μ 2 ( r ) 2 ( μ ( r ) + μ 2 ( r ) ) = μ ( r ) ν ( r ) 1 e ν ,
where we have assumed that e μ ( r ) 1 . Consequently, the solution to the above equation can be written as
e μ ( r ) = B e ν ( r ) 1 d r + A 2
with B and A as arbitrary parameters. Further, we shall find the solution of field Equations (7)–(9) by using the relationship of metric coefficients g t t and g r r given by Equation (13). Let us consider an interesting and significant ansatz for the metric coefficient g r r [53] which is given by
e ν ( r ) = 1 + ( a r 2 + 1 ) n c r 2 ( b r 2 + 1 ) 2 ,
where a , b , c and n are all arbitrary parameters. Utilizing Equation (14) in Equation (13), we obtain the metric coefficient g t t as follows
e μ ( r ) = B b r 2 + 1 a b r 2 + b a b r 2 + a n 2 c r 2 a r 2 + 1 n b r 2 + 1 2 2 F 1 n 2 , n 2 ; 1 n 2 ; a b a b r 2 + a b n r + A 2 .
Further, we substitute the values of e μ ( r ) and e ν ( r ) from Equations (14) and (15) in the field Equations (7) and (8) and consequently, we get the following expressions for the energy density function, the radial and tangential pressure components and measure of anisotropy parameter:
8 π ρ = b r 2 + 1 ( λ + 1 ) ( 2 λ + 1 ) r a r 2 + 1 c r 2 a r 2 + 1 n + b r 2 + 1 2 2 × c r 2 a r 2 + 1 n b r 2 + 1 2 A b n r a b r 2 + b a b r 2 + a n / 2 + B b r 2 + 1 c r 2 a r 2 + 1 n b r 2 + 1 2 2 F 1 n 2 , n 2 ; 1 n 2 ; a b a b r 2 + a × b n b r 2 + 1 a b r 2 + b a b r 2 + a n / 2 A ( 2 λ + 1 ) b r 2 + 1 r 2 a b ( 2 n 1 ) r 2 + 2 n + 3 b + 3 + c r 2 a r 2 + 1 n + 1 c r 2 a r 2 + 1 n b r 2 + 1 2 + B λ r b r 2 + 1 r 2 a b ( n + 1 ) r 2 + n + 3 + b + 3 + 2 c r 2 a r 2 + 1 n + 1 + B c ( 2 λ + 1 ) r a r 2 + 1 n b r 2 + 1 r 2 a b × ( 2 n 1 ) r 2 + 2 n + 3 b + 3 + c r 2 a r 2 + 1 n + 1 2 F 1 n 2 , n 2 ; 1 n 2 ; a b a b r 2 + a ,
8 π p r = ( b r 2 + 1 ) ( λ + 1 ) ( 2 λ + 1 ) r a r 2 + 1 c r 2 a r 2 + 1 n + b r 2 + 1 2 2 × c r 2 a r 2 + 1 n b r 2 + 1 2 A b n r a b r 2 + b a b r 2 + a n / 2 + B b r 2 + 1 c r 2 a r 2 + 1 n b r 2 + 1 2 2 F 1 n 2 , n 2 ; 1 n 2 ; a b a b r 2 + a × b n b r 2 + 1 a b r 2 + b a b r 2 + a n / 2 A ( 2 λ + 1 ) a r 2 + 1 c r 2 a r 2 + 1 n + b r 2 + 1 2 × c r 2 a r 2 + 1 n b r 2 + 1 2 + B r b r 2 + 1 r 2 λ a b ( n 3 ) r 2 + n 1 3 b 2 a b r 2 + a + b λ 2 2 B c ( λ + 1 ) r 3 a r 2 + 1 n + 1 + B c ( 2 λ + 1 ) r a r 2 + 1 n + 1 c r 2 a r 2 + 1 n + b r 2 + 1 2 2 F 1 n 2 , n 2 ; 1 n 2 ; a b a b r 2 + a ,
8 π p t = c r ( λ + 1 ) ( 2 λ + 1 ) c r 2 a r 2 + 1 n b r 2 + 1 2 c r 2 a r 2 + 1 n + b r 2 + 1 2 2 × a r 2 + 1 n 1 A b n r a b r 2 + b a b r 2 + a n / 2 + B b r 2 + 1 c r 2 a r 2 + 1 n b r 2 + 1 2 2 F 1 n 2 , n 2 ; 1 n 2 ; a b a b r 2 + a × b n a b r 2 + b a b r 2 + a n / 2 A ( 2 λ + 1 ) b r 2 + 1 a b ( n 1 ) r 4 + r 2 ( a n + a b ) + 1 c r 2 a r 2 + 1 n b r 2 + 1 2 + B r b r 2 + 1 r 2 a λ b ( n 1 ) r 2 + n + 1 + b n r 2 + n + 2 b λ + λ + 2 + B c r 3 a r 2 + 1 n + 1 B c ( 2 λ + 1 ) r a r 2 + 1 n a b ( n 1 ) r 4 + r 2 ( a n + a b ) + 1 2 F 1 n 2 , n 2 ; 1 n 2 ; a b a b r 2 + a ,
= r b r 2 + 1 c r 2 a r 2 + 1 n b r 2 + 1 2 8 π ( λ + 1 ) a r 2 + 1 c r 2 a r 2 + 1 n + b r 2 + 1 2 2 × c a r 2 + 1 n + 1 b r 2 + 1 b a ( n 2 ) r 2 2 + a n A b n r a b r 2 + b a b r 2 + a n / 2 + B b r 2 + 1 c r 2 a r 2 + 1 n b r 2 + 1 2 2 F 1 n 2 , n 2 ; 1 n 2 ; a b a b r 2 + a × b n b r 2 + 1 a b r 2 + b a b r 2 + a n / 2 A c r 2 a r 2 + 1 n b r 2 + 1 2 B r + B c r a r 2 + 1 n × 2 F 1 n 2 , n 2 ; 1 n 2 ; a b a b r 2 + a .

2.2. Boundary Conditions

In this subsection, we shall define the boundary conditions required by taking the obtained solutions as interior geometry while the exterior geometry is described by the Schwarzschild line-element which is defined as
d s 2 = Ψ d t 2 + Ψ 1 d r 2 + r 2 d Ω 2 ,
where Ψ = ( 1 2 M / r ) and consequently, Ψ 1 = ( 1 2 M / r ) 1 . In addition, d Ω 2 = d θ 2 + s i n 2 θ d ϕ 2 . Here we shall apply a constraint on the radial coordinate given by r > r s , where r s is the Schwarzschild radius. It is worthwhile to specify here that without applying this condition, black hole solutions can be obtained.
Using the metric components g t t and g r r from Equations (14) and (15) over the boundary r = R , we obtain the following relations
1 2 M R = B b R 2 + 1 a b R 2 + b a b R 2 + a n 2 c R 2 a R 2 + 1 n b R 2 + 1 2 2 F 1 n 2 , n 2 ; 1 n 2 ; a b a b R 2 + a b n R + A 2 ,
1 2 M R 1 = 1 + ( a R 2 + 1 ) n c R 2 ( a R 2 + 1 ) 2 ,
p r ( r = R ) = 0 .
Further, by taking boundary conditions in Equations (21) and (22) into account, expressions for the unknowns can be obtained as follows
c = 2 M b r 2 + 1 2 a r 2 + 1 n r 2 ( 2 M r ) ,
A = r 2 M r B F a b r 2 + b a b r 2 + a n 2 c a r 2 + 1 n b n ,
B = b B 1 ( 2 λ + 1 ) n r r 2 M r b B 2 n r b r 2 + 1 B 3 B 4 F ( 2 λ + 1 ) a r 2 + 1 ,
where
F = 2 F 1 n 2 , n 2 ; 1 n 2 ; a b a b R 2 + a , B 1 = c a r 2 + 1 n 2 + 1 a b r 2 + b a b r 2 + a n / 2 c r 2 a r 2 + 1 n + b r 2 + 1 2 , B 2 = a b r 2 + b a b r 2 + a n / 2 b r 2 + 1 r 2 λ a b ( n 3 ) r 2 + n 1 3 b 2 a b r 2 + a + b λ 2 2 c ( λ + 1 ) r 2 a r 2 + 1 n + 1 , B 3 = b r 2 + 1 c a r 2 + 1 n c r 2 a r 2 + 1 n b r 2 + 1 2 c r a r 2 + 1 n , B 4 = c r 2 ( 1 + a r 2 ) n + ( 1 + b r 2 ) 2 .
Other remaining parameters a , b and n can be considered as free and dimensionless whereas notations M and R refer to the mass and radius of an object, respectively. In this study, for further calculations and discussions, we consider the V e l a X 1 ( 1.77 M , 9.56 km) compact star model along with a = 0.001 km 2 and b = 0.0015 km 2 .

3. Graphical Analysis of Some Physical Properties of Presented Model in f ( R , T ) Theory

In the present section, we shall investigate graphically our presented stellar object model in f ( R , T ) modified gravity by focusing on some necessary physical properties and conditions satisfied by an astronomical compact star object. For this purpose, we shall concentrate on some physical features of the stellar interior like metric potential, energy density distribution, pressure profile, Zeldovich’s condition and gradients of energy density and pressure profile with respect to radial coordinate. In this respect, we shall discuss an astronomical object like V e l a X 1 ( 1.77 M , 9.56 km) under the present model. Since the parameter λ , which is present in the linear form of f ( R , T ) function, is very significant as it indicates the deviation from GR, therefore we shall investigate the physical properties of this model by taking two cases of parameter λ into account: its small values like λ = 0.35 and large values like λ = 2 .

3.1. Analysis of Some Features of the Model for λ = 0.35

3.1.1. Metric Potential

Herein this segment, we shall describe the graphical behavior of metric potentials for an astronomical object, i.e., g t t = e μ , and g r r = e ν . From Equations (31) and (33), it can be observed that the metric functions satisfy e ν ( r = 0 ) = 1 and e μ ( r = 0 ) 0 , which indicated that our obtained model is physically interesting and acceptable. Further, it can be seen from the left panel of Figure 1 that both metric components exhibit regular, increasing and hence monotonic behavior in the inner side of astronomical object. Therefore, it is concluded that the metric potential components satisfy the required necessary condition which ensures physical viability of the obtained model.

3.1.2. Zeldovich’s Condition

Here, we explore the energy density distribution and radial pressure profile at central values r = 0 , i.e., ρ ( r = 0 ) and p r ( r = 0 ) , which are defined as
ρ c = 3 c 8 π ( λ + 1 ) ( 2 λ + 1 ) ,
p r c = c b n b a n / 2 B ( λ + 2 ) A c ( 2 λ + 1 ) B c ( 2 λ + 1 ) 2 F 1 n 2 , n 2 ; 1 n 2 ; 1 b a 8 π ( λ + 1 ) ( 2 λ + 1 ) A b n b a n / 2 + B c 2 F 1 n 2 , n 2 ; 1 n 2 ; 1 b a .
The Zeldovich’s condition is defined as a ratio of density and pressure given by Equations (27) and (28), i.e., p r c / ρ c , which turns out be
b n b a n / 2 B ( λ + 2 ) A c ( 2 λ + 1 ) B c ( 2 λ + 1 ) 2 F 1 n 2 , n 2 ; 1 n 2 ; 1 b a 3 c A b n b a n / 2 + B c 2 F 1 n 2 , n 2 ; 1 n 2 ; 1 b a 1 .
It is evident from the Equation (29) and Table 1 that the Zeldovich’s condition is satisfied for all suitable values of involved parameters in f ( R , T ) gravity under the present model. From the consistency of Zeldovich’s condition, it is confirmed that our presented model is physically acceptable for a stellar object.

3.1.3. Energy Density

In this segment, we shall explore the graphical behavior of energy density for our obtained model under the Karmarker condition for a compact star V e l a X 1 in f ( R , T ) modified theory. The graphical illustration of energy density has been provided in the right plot of Figure 1. Here, it is noted that the density function exhibits positive decreasing behavior versus radial coordinate for r 9.56 along with 1.8 n < 7 excluding the values n = 2 , 4 and n = 6 . In this plot, the solid graphs illustrate the current development of the energy density function for 1.8 n < 7 and the dashed graphs refer to its behavior for n = 7 .

3.1.4. Pressure

Here, we shall provide the graphical behavior of radial and tangential pressure profiles for a V e l a X 1 compact star model representing a new family of solutions under the Karmarker condition in f ( R , T ) theory. The graphical development of radial and pressure components have been presented in Figure 2. It can be seen from these plots of Figure 2 that both radial and tangential pressures admit positive decreasing behavior for the radial coordinate satisfying r 9.56 with 1.8 n < 7 (excluding even numbers between 1.8 n < 7 ).

3.1.5. Anisotropy Distribution

In this part, we shall describe the graphical behavior of anisotropy measurement parameter denoted by ▵ versus radial coordinate which is provided in the left panel of Figure 3. It indicates that the measure of anisotropy distribution, i.e., = 8 π ( p t p r ) is positive throughout the region of the compact star V e l a X 1 in f ( R , T ) modified theory of gravity. It is observed that for > 0 , i.e., p t > p r , the anisotropic force is seen as outward-directed while for < 0 , this force is inward-directed. Since, in the present study, as we have > 0 and consequently, it indicates that the anisotropic pressure is repulsive and hence supports to the structure of compact object.

3.1.6. Gradients

Here, we shall present the graphical illustration of radial derivatives of energy density and radial pressure profile, i.e., d ρ / d r , d p r / d r . The corresponding graphical illustration is provided in the right panel of Figure 3. Here seperate visible curves refer to the derivatives of the energy density function, while the overlapping curves correspond to the derivative of pressure profiles with respect to radial coordinate for a V e l a X 1 compact object. From this graph, it can be easily observed that all curves representing derivatives show negative decreasing response and consequently, satisfy the primary conditions, i.e., d ρ / d r < 0 , d p r / d r < 0 and d p t / d r < 0 . It is also seen that these derivatives are zero at the central radius, i.e., d ρ ( r = 0 ) / d r = 0 , d p r ( r = 0 ) / d r = 0 , and d p t ( r = 0 ) / d r = 0 . Hence it can be concluded that these gradients meet the required constraints in this modified theory.

3.1.7. Equilibrium Condition

In this section, we shall explore the stability of the obtained new family of solutions by taking the equilibrium equation into account. For this purpose, we shall employ a well-known and fundamental equation, known as Tolman–Oppenheimer–Volkov (TOV) equation for checking the stability of the present configuration. This equation is given as follows
2 r ( p t p r ) d p r d r μ ( r ) 2 ( ρ + p r ) = 0 ,
The above equation describes the equilibrium condition by considering three forces namely anisotropic, hydrostatic and gravitational forces into account. The forces are defined by the following mathematical relations:
F a = 2 r ( p t p r ) , F h = d p r d r , F g = μ ( r ) ( ρ + p r ) 2
and consequently, Equation (30) takes the form
F a + F h + F g = 0 .
Further, we discuss the graphical response of these three compulsory forces in reference to the equilibrium equation. In this respect, the behavior of anisotropic force is provided in the left panel of Figure 4, while the plot of gravitational force is given in its right panel. Likewise, the graphical response of hydrostatic force is presented in the left panel of Figure 5, while the stability of the configuration can be observed by the graphical representation of equilibrium condition provided in the right panel of Figure 5. It can be easily observed that these forces balance each other’s effect for V e l a X 1 in f ( R , T ) modified theory of gravity and hence the stability of the model has been achieved under the mass and radii which are described in Table 1. Thus, another necessary required condition is satisfied in our case.

3.1.8. Energy Conditions

The study of energy constraints for checking the viability of a constructed model in GR as well as modified theories is regarded as one of the most interesting topic. These constraints also help to fix the values of arbitrary constants. These constraints are defined in terms four different inequalities which are labeled as D E C , S E C , W E C and N E C and are defined as
N E C : ρ 0 , W E C : ρ p t 0 , ρ p r 0 , S E C : ρ p r 2 p t 0 , D E C : ρ > | p r | , ρ > | p t | .
Here, we discuss the graphical response of these different energy conditions in the connection of N E C , W E C , S E C and D E C . The N E C and D E C expressions have been already discussed in the previous section of energy density and pressure profile which indicated that these constraints are consistent in our case. Further, the graphical illustration of remaining W E C and S E C constraints are provided in the right and left panels of Figure 6 and Figure 7. It can be easily observed that these all energy bounds remain satisfied in the current scenario which is a necessary condition for the physical acceptance of compact stars under the observational data (as given in Table 1) in f ( R , T ) theory.

3.1.9. Mass Function, Compactness Parameter and Surface Red-Shift Function

In this segment, we shall discuss the relation of mass and radius of a compact star graphically. In this connection, we shall discuss mass function, mass-radius ratio and surface red-shift function, which are listed below in detail.
The expression of mass function for a compact star with radius r is defined as
m ( r ) = 4 π 0 r ( r 2 ρ ) d r .
In the current scenario, we use the radius under our presented model and the corresponding graph is presented in the left part of Figure 8. It can be easily observed that the mass function exhibits positive increasing behavior throughout the interior of stellar object, i.e., m ( r ) > 0 , for r 9.56 . Further, it can be noticed that m ( r ) 0 as r 0 and consequently, it can be concluded that this condition leads to the regularity of the mass function at the center.
Next we shall discuss graphically the mass-radius function, i.e., compactness parameter, which is defined by the following mathematical expression
u ( r ) = 2 m ( r ) r .
From the right part of Figure 8, the graphical development of compactness parameter can be noticed which indicates that parameter u ( r ) remains positive and increasing. It is also easy to verify from this graph that the mass-radius function satisfy the Bhuchdahl [54] condition, i.e., u ( r ) < 8 / 9 .
Now, with the help of mass-radius function, we investigate the graphical response of surface red-shift function for the considered stellar object, which is defined as
Z s = 1 1 2 u ( r ) 1 .
The graphical response of surface red-shift function, i.e., Z s against the radial coordinate r is presented in the left panel of Figure 9. From this graph, it is found that Z s vanishes at r = 0 and gradually increases with the increasing values of r. It is also worthwhile to mention here that in our case, Z s satisfy the Bohmer and Harko [55] condition under the anisotropic configuration, i.e., Z s < 5 .

3.1.10. Equation of State

In this segment, we calculate the equations of state for both radial and tangential components of pressure, which are denoted by w r and w t , respectively. These parameters are defined in the normal and transverse directions as follows
w r = p r ρ ,
w t = p t ρ .
These two ratios define the equations of state. Here, we explore the graphical response of these two ratios, i.e., w r and w t , in the reference of the equation of state parameter. It is noticed from the right part of Figure 9, that both these ratios are less than 1, which is a mandatory condition for a stellar object. It is also shown that both ratios exhibit decreasing behaviors under the current scenario. Hence it can be concluded that the equations of state satisfy the necessary requirement in f ( R , T ) theory for V e l a X 1 model.

3.1.11. Causality Stability Analysis

In this segment, we shall describe the causality stability analysis for the obtained solutions. In this reference, we shall discuss two sorts of the speed of sound parameter (radial and tangential) for V e l a X 1 compact star model. These speed of sound parameters are represented by v t 2 and v r 2 for radial and tangential components, respectively and can be calculated as
v t 2 = d p t d ρ , v r 2 = d p r d ρ .
It is worthwhile to mention here that the casuality stability of a configuration can be achieved if both the tangential and radial speed of sound parameters take values less than 1. In addition, the difference between these velocities denoted by v t 2 v r 2 must satisfy Abrea condition [56,57], which is given by 1 v t 2 v r 2 0 .
For the present case, the graphs of these parameters are provided in Figure 10 and Figure 11. It is observed form the left part of Figure 10 that the expression v r 2 attains values less than 1. It is also noted from the right part of Figure 10 that the expression v t 2 takes values less than 1. Further, it is seen that for 1.8 n < 7 , both the radial and tangential velocities satisfy the necessary condition, but for n = 7 , both radial and tangential velocities take values greater than 1 and hence violate the necessary condition in this case. Furthermore, it can be noticed from Figure 10 that our model is not compatible with the causality stability criterion for n = 7 , which is represented by dashed graphs for both velocities. From the left panel of Figure 11, it can be observed that the difference v t 2 v r 2 is consistent with the condition 1 v t 2 v r 2 0 . Hence it can be concluded that our obtained model is potentially stable.

3.1.12. The Adiabatic Index Stability Analysis

In this portion, we shall explore an important and necessary ratio of two specific heats given by
Γ r = v r 2 ρ + p r p r .
It is significant to mention here that the stability of Newtonian sphere can be obtained under the restriction of Γ r > 4 / 3 . In literature, it is pointed out that this expression is used to define a neutral equilibrium when the condition Γ r = 4 / 3 holds. In the present case, the graphical development of Γ r parameter is provide in the right panel of Figure 11. It is easy to check that the relation Γ r shows monotonic increasing behavior for n < 7 and also is compatible with the condition Γ r > 4 / 3 for V e l a X 1 compact star, while it exhibits decreasing trend for n = 7 . Hence, the adiabatic index is satisfied with the stability condition in the reference of f ( R , T ) theory for n < 7 . In Table 2, we have presented a summary of our obtained solutions. All the physical properties can be seen from Table 2.

3.2. Graphical Analysis of Some Physical Properties for λ = 2

Here we shall investigate the behavior of presented model by exploring different physical properties graphically (as discussed in previous subsections) for large values of parameter λ . For this purpose, we fix the parameter as λ = 2 and the corresponding calculated values of unknowns c , A and B under V e l a X 1 ( 1.77 M , 9.56 km) are provided in Table 3. Further the graphical analysis of different properties of our obtained solution is given in Figure 12, Figure 13, Figure 14, Figure 15, Figure 16, Figure 17, Figure 18, Figure 19, Figure 20, Figure 21 and Figure 22. From Figure 12, it is easy to verify that the metric potentials satisfy the necessary condition required for a realistic compact star model. Further, the density function shows positive decreasing behavior versus radial coordinate which is also physically acceptable. In Figure 13, it can be seen that radial pressure exhibits physically valid behavior while the tangential pressure shows negative behavior near the boundary and hence refers to non-realistic compact star model. The measure of anisotropy parameter takes negative values as shown in Figure 14 which corresponds to p t < p r and hence the force on the compact star will be inward directed. In addition, the radial derivatives of density and radial pressure exhibit negative decreasing behavior which also favor to a realistic model. The graphical response of gravitational, hydrostatic and anisotropic forces is provided in Figure 15 and Figure 16 which indicate the stability of compact star configuration.
The graphical behavior of energy constraints is provided in Figure 17 and Figure 18 which showed that these constraints are valid for the present compact star model. The curves for mass function, compactness parameter and surface red-shift function are provided in Figure 19 and Figure 20, which indicate positive increasing behavior and hence favors to the physically acceptable model of compact star. Figure 21 and Figure 22 provide the graphical response of speed of sound parameters (radial and tangential) as well as adiabatic index. It can be easily observed that for large values of λ , the conditions 1 < v t 2 v r 2 < 0 and Γ r > 4 3 fail to hold and hence refer to a non-realistic compact star model.

4. Conclusions

In the present paper, we have investigated a new family of solutions of embedding class-I model to find the physically acceptable solutions for V e l a X 1 compact star in modified f ( R , T ) theory of gravity. For this purpose, we have adopted the well-known Karmarkar condition in the reference of Pandey–Sharma condition. Herein this study, we have explored some interesting physical properties related to the compact star configuration graphically and analytically. For the present study, we have taken spherically symmetric geometry with the matter contents involving anisotropic pressure profile without any net electric charge. Here we have assumed e ν ( r ) = 1 + ( a r 2 + 1 ) n c r 2 ( b r 2 + 1 ) 2 model for g r r metric component. It is seen that the parameter n has an important role in the discussion as the the compact star solutions for the choices n = 2 , n = 4 , n = 6 can not be calculated. It is also observed that some important physical properties are also not satisfied for n = 7 . Likewise, it is found that the solutions for n > 7 can not be calculated. In addition, we have used a simple linear model given by f ( R , T ) = R + λ T for this study. In this model, the parameter λ has also an important role for the present discussion. In order to check its impact on the obtained solution, we have considered its two values given by λ = 0.35 and λ = 2 . It is seen that the obtained solutions exhibited stable and realistic behavior for n [ 1.8 , 7 ) 2 , 4 , 6 . For λ = 0.35 , the graphical illustration of different physical properties have been presented in Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11. In the following, we shall summarize our major findings.
  • From Figure 1, it is seen that the metric components, i.e., g t t = e μ , and g r r = e ν , are calculated as e ν ( r = 0 ) = 1 and e μ ( r = 0 ) 0 , which shows that this model is physically realistic and acceptable for V e l a X 1 compact star;
  • The Zeldovich’s condition, i.e., p r c / ρ c 1 , has been observed to be satisfied which can be confirmed From the Table 1;
  • The energy density function has remained positive throughout in the stellar interior V e l a X 1 compact object for this current model which can be seen from Figure 1;
  • Both the radial and tangential pressure components has remained positive throughout in the stellar interior for V e l a X 1 compact star. This behavior of both components can be verified from Figure 2;
  • It is found that the measure of anisotropic pressure ( r ) is positive throughout the region of the compact star and consequently, it supports to the structure of compact star which has been provided in Figure 3;
  • The derivative of energy density function and radial pressure with respect to radial coordinate are calculated as negative and it is also confirmed form their graphical behavior given in Figure 3;
  • From the Table 2 and Figure 4 and Figure 5, it is found that the forces F a , F h and F g are consistent with the equilibrium Condition for V e l a X 1 model. It is also observed that these forces balanced each other’s effect and hence left the whole configuration stable;
  • The energy conditions, i.e., ρ , p r , p t , ρ p r , ρ p t , ρ p r 2 p t , are observed to be satisfied for this model. Their graphically illustration has been provided in Figure 6 and Figure 7;
  • The mass-radii function m ( r ) has remained positive, increasing and regular. Its graphical behavior can be observed from Figure 8;
  • The compactness parameter u ( r ) has remained positive and it also satisfied the Buchdahl limit, i.e., u ( r ) 8 / 9 which can be seen from Figure 8;
  • From the graph, it is observed that the surface red-shift function, i.e., Z s turned to be zero at r = 0 and gradually increased with the increasing radial coordinate. In addition, it satisfied the Bohmer and Harko condition under the anisotropic configuration, i.e., Z s < 5 ;
  • In radial and tangential directions, we have defined two equation of state (EoS), namely w r and w t . It is noticed that the values of these EoS parameters remained positive inside the stellar object and also both assumed values less than 1 which can be seen from Figure 9;
  • It is expressed From the Table 2 and Figure 10 and Figure 11 that the radial and tangential speeds of sound for compact stars, which are denoted by v r 2 and v t 2 satisfied the necessary condition for 1.8 n < 7 , while for n = 7 , these radial and tangential velocities are observed greater than 1, and consequently, violated the necessary condition. Further, the causality stability condition for n = 7 is also violated for the present model;
  • In reference of adiabatic index, it is found From the Table 2 and Figure 11 that the adiabatic index Γ r satisfied the inequality Γ r > 4 / 3 and showed increasing behavior.
For the second choice λ = 2.00 , the following properties are failed to meet the basic criteria for the existence of realistic compact stars.
  • Tangential pressure profile has turned out as negative near the boundary, i.e., r = R of astral objects. The negative tangential pressure components showed that our calculated solutions are non-realistic.
  • The anisotropy function exhibited negative behavior throughout the configuration. The non-zero anisotropy function is a necessary requirement for the existence of realistic objects.
  • The expression v t 2 v r 2 turned out to be inconsistent with the Abrea condition, i.e., 1 v t 2 v r 2 0 .
  • Likewise, the Adiabatic index, i.e., Γ r has also been incompatible with the condition Γ r > 4 / 3 .
It is found that these solutions fulfill all the necessary conditions under the observational radii and mass attribute data for small values of parameter λ = 0.35 and hence physically well-behaved and promising. The graphical analysis of different properties of our obtained solution is given in Figure 12, Figure 13, Figure 14, Figure 15, Figure 16, Figure 17, Figure 18, Figure 19, Figure 20, Figure 21 and Figure 22 for λ = 2 . In fact, the obtained solutions are not confirmed satisfied entirely for λ > 0.35 . It is not possible to define a limit on parameter λ with observational data in this present study. The smaller values of λ have provided us a well-behaved results with good accuracy. Hence, to sum-up, it can be concluded that our obtained model exhibited well-behaved nature and physically considerable for 1.8 n < 7 . For other choices of parameter n and λ , the solutions do not correspond to a realistic compact star model. If we put λ = 0 , we can get the results for GR case, which is already discussed by Singh et al. [53].

Author Contributions

Conceptualization, S.W. and M.Z.; methodology, G.M., S.W., and A.A.; software, G.M. and A.A.; validation, S.W., G.M., M.Z., and A.A.; formal analysis, S.W. and M.Z.; investigation, S.W., G.M., M.Z., and A.A.; data curation, S.W. and G.M.; writing—original draft preparation, S.W., A.A., and G.M.; writing—review and editing, M.Z.; visualization, M.Z.; supervision, S.W. and M.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research has received APC fee from Prince Mohammad Bin Fahd University, Al Khobar 31952, Saudi Arabia.

Acknowledgments

We are grateful to the anonymous referees for their valuable comments, which improved the presentation of this paper. Saira Waheed thanks to Prince Mohammad Bin Fahd University, Al Khobar 31952, Saudi Arabia for financial support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Schwarzschild, K. Über das Gravitationsfeld einer Kugel aus inkompressibler Flüssigkeit nach der Einsteinschen Theorie. Sitz. Dtsch. Akad. Wiss. Math. Phys. Berl. 1916, 24, 424. [Google Scholar]
  2. Tolman, R.C. Static solutions of Einstein’s field equations for spheres of fluid. Phys. Rev. 1939, 55, 364. [Google Scholar] [CrossRef] [Green Version]
  3. Oppenheimer, J.R.; Volkoff, G.M. On massive neutron cores. Phys. Rev. 1939, 55, 374. [Google Scholar] [CrossRef]
  4. Oppenheimer, J.R.; Snyder, H. On continued gravitational contraction. Phys. Rev. 1939, 56, 455. [Google Scholar] [CrossRef] [Green Version]
  5. Jeans, J. The motions of stars in a Kapteyn universe. Mon. Not. R. Astron. Soc. 1922, 82, 122–132. [Google Scholar] [CrossRef]
  6. Lemaitre, G. L’univers en expansion. Ann. Soc. Sci. Brux. A 1937, 10, 1–19. [Google Scholar]
  7. Bowers, R.L.; Liang, E.P.T. Anisotropic spheres in general relativity. Astrophys. J. 1974, 188, 657. [Google Scholar] [CrossRef]
  8. Ruderman, R. Pulsars: Structure and dynamics. Annu. Rev. Astron. Astrophys. 1972, 10, 427–476. [Google Scholar] [CrossRef]
  9. Dev, K.; Gleiser, M. Anisotropic Stars: Exact Solutions. Gen. Relat. Gravit. 2002, 34, 1793–1818. [Google Scholar] [CrossRef]
  10. Herrera, L.; Santos, N.O. Local anisotropy in self-gravitating systems. Phys. Rep. 1997, 286, 53–130. [Google Scholar] [CrossRef]
  11. Riess, A.G.; Filippenko, A.V.; Challis, P.; Clocchiatti, A.; Diercks, A.; Garnavich, P.M.; Gilliland, R.L.; Hogan, C.J.; Jha, S.; Kirshner, R.P.; et al. Observational evidence from supernovae for an accelerating universe and a cosmological constant. Astron. J. 1998, 116, 1009. [Google Scholar] [CrossRef] [Green Version]
  12. Perlmutter, S.; Aldering, G.; Goldhaber, G.; Knop, R.A.; Nugent, P.; Castro, P.G.; Deustua, S.; Fabbro, S.; Goobar, A.; Groom, D.E.; et al. Measurements of Ω and Λ from 42 high-redshift supernovae. Astrophys. J. 1999, 517, 565. [Google Scholar] [CrossRef]
  13. de Bernardis, P.; Ade, P.A.; Bock, J.J.; Bond, J.R.; Borrill, J.; Boscaleri, A.; Coble, K.; Crill, B.P.; De Gasperis, G.; Farese, P.C.; et al. A flat Universe from high-resolution maps of the cosmic microwave background radiation. Nature 2000, 6781, 955–959. [Google Scholar] [CrossRef] [PubMed]
  14. Knop, R.A.; Aldering, G.; Amanullah, R.; Astier, P.; Blanc, G.; Burns, M.S.; Conley, A.; Deustua, S.E.; Doi, M.; Ellis, R.; et al. New constraints on ΩM, ΩΛ, and w from an independent set of 11 high-redshift supernovae observed with the Hubble Space Telescope. Astrophys. J. 2003, 598, 102. [Google Scholar] [CrossRef] [Green Version]
  15. Capozziello, S.; De Laurentis, M. Extended theories of gravity. Phys. Rep. 2011, 509, 167–321. [Google Scholar] [CrossRef] [Green Version]
  16. Capozziello, S.; De Laurentis, M. The dark matter problem from f(, 𝒯) gravity viewpoint. Ann. Phys. 2012, 524, 545–578. [Google Scholar] [CrossRef]
  17. Starobinksy, A.A. A New Type of Isotropic Cosmological Models Without Singularity. Phys. Lett. B 1980, 91, 99–102. [Google Scholar] [CrossRef]
  18. Bhawal, B.; Kar, S. Lorentzian wormholes in Einstein-Gauss-Bonnet theory. Phys. Rev. D 1992, 46, 2464. [Google Scholar] [CrossRef]
  19. Rahaman, F.; Kalam, M.; Ghosh, A. Existence of wormholes in Einstein-Kalb-Ramond space time. arXiv 2006, arXiv:gr-qc/0605095. [Google Scholar]
  20. Chakraborty, S.; Bandyopadhyay, T. Wormhole and its analogue in brane world. arXiv 2008, arXiv:0707.1181. [Google Scholar]
  21. Ferraro, R.; Fiorini, F. Modified teleparallel gravity: Inflation without an inflaton. Phys. Rev. D 2007, 75, 084031. [Google Scholar] [CrossRef] [Green Version]
  22. Sharif, M.; Waheed, S. Energy conditions in a generalized second-order scalar-tensor gravity. Adv. High Energy Phys. 2013, 2013. [Google Scholar] [CrossRef]
  23. Kofinas, G.; Saridakis, N.E. Teleparallel equivalent of Gauss-Bonnet gravity and its modifications. Phys. Rev. D 2014, 90, 084044. [Google Scholar] [CrossRef] [Green Version]
  24. Zubair, M.; Kousar, F. Cosmological reconstruction and energy bounds in f(R, RαβRαβ, ϕ) gravity. Eur. Phys. J. C 2016, 76, 254. [Google Scholar] [CrossRef] [Green Version]
  25. Ketov, S.V. Modified supergravity and early universe: The meeting point of cosmology and high-energy physics. Int. J. Mod. Phys. A 2013, 28, 1330021. [Google Scholar] [CrossRef] [Green Version]
  26. Harko, T.; Lobo, F.S.N.; Nojiri, S.; Odintsov, S.D. f(, 𝒯) gravity. Phys. Rev. D. 2011, 84, 024020. [Google Scholar] [CrossRef] [Green Version]
  27. Jamil, M.; Momeni, D.; Myrzakulov, R. Attractor solutions in f(, 𝒯) cosmology. Eur. Phys. J. C. 2012, 72, 1959. [Google Scholar] [CrossRef]
  28. Jamil, M.; Momeni, D.; Myrzakulov, R. Violation of the First Law of Thermodynamics in f(, 𝒯) Gravity. Chin. Phys. Lett. 2012, 29, 109801. [Google Scholar] [CrossRef] [Green Version]
  29. Shabani, H.; Farhoudi, M. f(, 𝒯) cosmological models in phase space. Phys. Rev. D 2013, 88, 044048. [Google Scholar] [CrossRef] [Green Version]
  30. Shabani, H.; Farhoudi, M. Cosmological and solar system consequences of f(, 𝒯) gravity models. Phys. Rev. D 2014, 90, 44031. [Google Scholar] [CrossRef] [Green Version]
  31. Moraes, P.H.R.S. Cosmological solutions from induced matter model applied to 5D gravity and the shrinking of the extra coordinate. Eur. Phys. J. C 2015, 75, 168. [Google Scholar] [CrossRef] [Green Version]
  32. Sharif, M.; Yousaf, Z. Dynamical analysis of self-gravitating stars in f(, 𝒯) gravity. Astrophys. Space Sci. 2014, 354, 471–479. [Google Scholar] [CrossRef]
  33. Zubair, M.; Abbas, G.; Noureen, I. Possible formation of compact stars in gravity. Astrophys. Space Sci. 2016, 361, 8. [Google Scholar] [CrossRef]
  34. Alhamzawi, A.; Alhamzawi, R. Gravitational lensing by f(, 𝒯) gravity Int. J. Mod. Phys. D 2016, 25, 1650020. [Google Scholar] [CrossRef]
  35. Moraes, P.H.R.S.; Arbañil, J.D.; Malheiro, M. Stellar equilibrium configurations of compact stars in f(R, T) theory of gravity. J. Cosmol. Astropart. Phys. 2016, 6, 5. [Google Scholar] [CrossRef] [Green Version]
  36. Das, A.; Rahaman, F.; Guha, B.K.; Ray, S. Compact stars f(, 𝒯) gravity. Eur. Phys. J. C 2016, 76, 654. [Google Scholar] [CrossRef] [Green Version]
  37. Karmarkar, K.R. Gravitational metrics of spherical symmetry and class one. Proc. Indian Acad. Sci. A 1948, 27, 56. [Google Scholar] [CrossRef]
  38. Komathiraj, K.; Maharaj, S.D. Classes of exact Einstein–Maxwell solutions. Gen. Relat. Gravit. 2007, 39, 2079–2093. [Google Scholar] [CrossRef] [Green Version]
  39. Ivanov, B.V. Collapsing shear-free perfect fluid spheres with heat flow. Gen. Relat. Gravit. 2012, 44, 1835–1855. [Google Scholar] [CrossRef]
  40. Ivanov, B.V. A conformally flat realistic anisotropic model for a compact star. Eur. Phys. J. C 2018, 78, 332. [Google Scholar] [CrossRef] [Green Version]
  41. Bhar, P.; Singh, K.N.; Pant, N. Compact star modeling with quadratic equation of state in Tolman VII space–time. Indian J. Phys. 2017, 91, 701–709. [Google Scholar] [CrossRef]
  42. Bhar, P.; Singh, K.N.; Rahaman, F.; Pant, N.; Banerjee, S. A charged anisotropic well-behaved Adler–Finch–Skea solution satisfying Karmarkar condition. Int. J. Mod. Phys. D 2017, 26, 1750078. [Google Scholar] [CrossRef] [Green Version]
  43. Maurya, S.K.; Gupta, Y.K.; Ray, S.; Deb, D. Generalised model for anisotropic compact stars. Eur. Phys. J. C 2016, 76, 693. [Google Scholar] [CrossRef] [Green Version]
  44. Maurya, S.K.; Gupta, Y.K.; Smith, T.T.; Rahaman, F. A new exact anisotropic solution of embedding class one. Eur. Phys. J. A 2016, 52, 191. [Google Scholar] [CrossRef]
  45. Maurya, S.K.; Gupta, Y.K.; Ray, S.; Deb, D. A new model for spherically symmetric charged compact stars of embedding class 1. Eur. Phys. J. C 2017, 77, 45. [Google Scholar] [CrossRef] [Green Version]
  46. Singh, K.N.; Bhar, P.; Pant, N. A new solution of embedding class I representing anisotropic fluid sphere in general relativity. Int. J. Mod. Phys. D 2016, 25, 1650099. [Google Scholar] [CrossRef] [Green Version]
  47. Singh, K.N.; Pant, N.; Govender, A.M. Anisotropic compact stars in Karmarkar spacetime. Chin. Phys. C 2017, 41, 015103. [Google Scholar] [CrossRef] [Green Version]
  48. Singh, K.N.; Pant, N. A family of well-behaved Karmarkar spacetime describing interior of relativistic stars. Astrophys. Space Sci. 2016, 361, 173. [Google Scholar] [CrossRef]
  49. Singh, K.N.; Bhar, P.; Pant, N. Solutions of the Einstein’s field equations with anisotropic pressure compatible with cold star model. Astrophys. Space Sci. 2016, 361, 339. [Google Scholar] [CrossRef]
  50. Rahaman, F.; Chakraborty, K.; Kuhfittig, P.K.F.; Shit, G.C.; Rahman, M. A new deterministic model of strange stars. Eur. Phys. J. C 2014, 74, 3126. [Google Scholar] [CrossRef] [Green Version]
  51. Malaver, M. Strange Quark Star Model with Quadratic Equation of State. Math. Appl. 2014, 1, 9–15. [Google Scholar]
  52. Herrera, L.; Barreto, W. General relativistic polytropes for anisotropic matter: The general formalism and applications. Phys. Rev. D 2013, 88, 084022. [Google Scholar] [CrossRef] [Green Version]
  53. Mustafa, G.; Zubair, M.; Waheed, S.; Xia, T. Realistic stellar anisotropic model satisfying Karmarker condition in f(R, T) gravity. Phys. J. C 2020, 80, 26. [Google Scholar]
  54. Singh, K.N.; Pant, N.; Tewari, N.; Aria, A.K. Embedded class solutions compatible for physical compact stars in general relativity. Eur. Phys. J. A 2018, 54, 77. [Google Scholar] [CrossRef]
  55. Buchdahl, H.A. General relativistic fluid spheres. Phys. Rev. D 1959, 116, 1027. [Google Scholar] [CrossRef]
  56. Bohmer, C.G.; Harko, T. Bounds on the basic physical parameters for anisotropic compact general relativistic objects. Class. Quantum Gravit. 2006, 23, 6479. [Google Scholar] [CrossRef] [Green Version]
  57. Abreu, H.; Hernandez, H.; Nunez, L.A. Sound speeds, cracking and the stability of self-gravitating anisotropic compact objects. Class. Quantum Gravit. 2007, 24, 4631. [Google Scholar] [CrossRef]
Figure 1. It shows the evolution of metric potentials (left) and energy density function (right) versus r for V e l a X 1 by using the parameters values from the Table 1.
Figure 1. It shows the evolution of metric potentials (left) and energy density function (right) versus r for V e l a X 1 by using the parameters values from the Table 1.
Symmetry 12 00962 g001
Figure 2. Provides the graphical response of radial (left) and tangential pressures (right) versus r for V e l a X 1 , using values of parameters From the Table 1.
Figure 2. Provides the graphical response of radial (left) and tangential pressures (right) versus r for V e l a X 1 , using values of parameters From the Table 1.
Symmetry 12 00962 g002
Figure 3. Describes the graphical behavior of anisotropy function (left) and radial derivatives of ρ and p r (right) versus for V e l a X 1 , by using the parameters values from the Table 1.
Figure 3. Describes the graphical behavior of anisotropy function (left) and radial derivatives of ρ and p r (right) versus for V e l a X 1 , by using the parameters values from the Table 1.
Symmetry 12 00962 g003
Figure 4. Represents the behavior of forces F a (left) and F g (right) versus r for V e l a X 1 by using the parameters values from the Table 1.
Figure 4. Represents the behavior of forces F a (left) and F g (right) versus r for V e l a X 1 by using the parameters values from the Table 1.
Symmetry 12 00962 g004
Figure 5. Illustrates the behavior of F h (left) and the balanced combine effect of F a , F h and F g (right) versus r for V e l a X 1 by using the parameters values from the Table 1.
Figure 5. Illustrates the behavior of F h (left) and the balanced combine effect of F a , F h and F g (right) versus r for V e l a X 1 by using the parameters values from the Table 1.
Symmetry 12 00962 g005
Figure 6. Shows the evolution of ρ p r (left) and ρ p t (right) versus r by using the parameters values from the Table 1.
Figure 6. Shows the evolution of ρ p r (left) and ρ p t (right) versus r by using the parameters values from the Table 1.
Symmetry 12 00962 g006
Figure 7. Corresponds to the behavior of ρ p r 2 p t (left) and all energy constraints together (right) against radial coordinate for V e l a X 1.
Figure 7. Corresponds to the behavior of ρ p r 2 p t (left) and all energy constraints together (right) against radial coordinate for V e l a X 1.
Symmetry 12 00962 g007
Figure 8. Shows the graphical behavior of mass-radius relation function (left) and compactness parameter (right) versus r for V e l a X 1.
Figure 8. Shows the graphical behavior of mass-radius relation function (left) and compactness parameter (right) versus r for V e l a X 1.
Symmetry 12 00962 g008
Figure 9. Provides the graph of surface red-shift function (left) and equations of state (right) versus r for V e l a X 1.
Figure 9. Provides the graph of surface red-shift function (left) and equations of state (right) versus r for V e l a X 1.
Symmetry 12 00962 g009
Figure 10. Represents the behavior of v r 2 (left) and v t 2 (right) against r for V e l a X 1.
Figure 10. Represents the behavior of v r 2 (left) and v t 2 (right) against r for V e l a X 1.
Symmetry 12 00962 g010
Figure 11. Shows the plots of v t 2 v r 2 (left) and Γ r (right) versus r by using the parameters values from the Table 1.
Figure 11. Shows the plots of v t 2 v r 2 (left) and Γ r (right) versus r by using the parameters values from the Table 1.
Symmetry 12 00962 g011
Figure 12. Shows the evolution of metric potential (left) and energy density function (right) against r for V e l a X 1 , by using the parameters values from the Table 3.
Figure 12. Shows the evolution of metric potential (left) and energy density function (right) against r for V e l a X 1 , by using the parameters values from the Table 3.
Symmetry 12 00962 g012
Figure 13. Provides the behavior of radial pressure (left) and tangential pressure (right) versus r for V e l a X 1 by using the parameters values from the Table 3.
Figure 13. Provides the behavior of radial pressure (left) and tangential pressure (right) versus r for V e l a X 1 by using the parameters values from the Table 3.
Symmetry 12 00962 g013
Figure 14. Describes the graphical response of anisotropy function (left) and radial derivative of ρ and p r (right) versus r for V e l a X 1, by using the parameters values from the Table 3.
Figure 14. Describes the graphical response of anisotropy function (left) and radial derivative of ρ and p r (right) versus r for V e l a X 1, by using the parameters values from the Table 3.
Symmetry 12 00962 g014
Figure 15. Refers to the plots of forces F a (left) and F g (right) versus r for V e l a X 1, by using the parameters values from the Table 3.
Figure 15. Refers to the plots of forces F a (left) and F g (right) versus r for V e l a X 1, by using the parameters values from the Table 3.
Symmetry 12 00962 g015
Figure 16. The behavior of F h (left) and the combine effect of F a , F h and F g (right) versus r by using the parameters values from the Table 3.
Figure 16. The behavior of F h (left) and the combine effect of F a , F h and F g (right) versus r by using the parameters values from the Table 3.
Symmetry 12 00962 g016
Figure 17. Shows the evolution of ρ p r (left) and ρ p t (right) against r.
Figure 17. Shows the evolution of ρ p r (left) and ρ p t (right) against r.
Symmetry 12 00962 g017
Figure 18. Describes the graphs of ρ p r 2 p t (left) and all energy conditions (right) against r.
Figure 18. Describes the graphs of ρ p r 2 p t (left) and all energy conditions (right) against r.
Symmetry 12 00962 g018
Figure 19. Shows the graphical behavior of mass-radius relation (left) and compactness parameter (right) versus r.
Figure 19. Shows the graphical behavior of mass-radius relation (left) and compactness parameter (right) versus r.
Symmetry 12 00962 g019
Figure 20. Plots correspond to the behavior of surface red-shift function (left) and equation of state (EoS) parameters (right) versus r.
Figure 20. Plots correspond to the behavior of surface red-shift function (left) and equation of state (EoS) parameters (right) versus r.
Symmetry 12 00962 g020
Figure 21. Shows the graphical response of v r 2 (left) and , v t 2 (right) against r.
Figure 21. Shows the graphical response of v r 2 (left) and , v t 2 (right) against r.
Symmetry 12 00962 g021
Figure 22. Shows the behavior of v t 2 v r 2 (left) and Γ r (right) against r.
Figure 22. Shows the behavior of v t 2 v r 2 (left) and Γ r (right) against r.
Symmetry 12 00962 g022
Table 1. Calculated values of c , A and B under V e l a X 1 ( 1.77 M , 9.56 km) with a = 0.001 km 2 b = 0.0015 km 2 and λ = 0.35 for different suitable values of n
Table 1. Calculated values of c , A and B under V e l a X 1 ( 1.77 M , 9.56 km) with a = 0.001 km 2 b = 0.0015 km 2 and λ = 0.35 for different suitable values of n
ncAB p rc / ρ c
1.80.007113.485300.043290.48013
2.20.00686−4.276400.043420.50382
2.60.00663−1.368600.043550.52824
3.00.00639−0.623930.043680.55319
3.40.00618−0.121250.043810.55319
3.80.005960.924850.043940.60458
4.20.00576−1.416300.044070.63089
4.60.00556−0.461140.044210.65790
5.00.00537−0.185170.044330.68546
5.40.005180.009220.044470.71351
5.80.005000.366620.044660.74225
6.20.00483−0.330380.044740.77137
6.60.00467−0.008880.044870.80142
7.00.004510.098960.045020.83202
Table 2. Calculated results for V e l a X 1 compact star under Karmarkar condition in f ( R , T ) theory for λ = 0.35 .
Table 2. Calculated results for V e l a X 1 compact star under Karmarkar condition in f ( R , T ) theory for λ = 0.35 .
Expressions 1.8 n < 7 n = 7
ρ ρ > 0 , satisfied ρ > 0 , satisfied
p r p r > 0 , satisfied p r > 0 , satisfied
p t p t > 0 , satisfied p t > 0 , satisfied
p r c / ρ c p r c / ρ c 1 p r c / ρ c 1
( r ) ( r ) > 0 , satisfied ( r ) > 0 , satisfied
d ρ / d r d ρ / d r < 0 , satisfied d ρ / d r < 0 , satisfied
d p r / d r d p r / d r < 0 , satisfied d p r / d r < 0 , satisfied
F a , F h and F g forces are Balanceforces are Balance
ρ p r ρ p r > 0 , satisfied ρ p r > 0 , satisfied
ρ p t ρ p t > 0 , satisfied ρ p t > 0 , satisfied
ρ p r 2 p t ρ p r 2 p t > 0 , satisfied ρ p r 2 p t > 0 , satisfied
m ( r ) m ( r ) > 0 , satisfied m ( r ) > 0 , satisfied
u ( r ) 0 < u ( r ) < 8 / 9 satisfied 0 < u ( r ) < 8 / 9 satisfied
Z s 0 < Z s < 5 , satisfied 0 < Z s < 5 , satisfied
w r 0 < w r < 1 , satisfied 0 < w r < 1 , satisfied
w t 0 < w t < 1 , satisfied 0 < w t < 1 , satisfied
v r 2 v r 2 < 1 , satisfied v r 2 > 1 , not-satisfied
v t 2 v t 2 < 1 , satisfied v t 2 > 1 , not-satisfied
v t 2 v r 2 1 v t 2 v r 2 0 , satisfied 1 v t 2 v r 2 0 , satisfied
Γ r 0 < Γ r < 4 / 3 , satisfied 0 < Γ r < 4 / 3 , satisfied
Table 3. Calculated values of c , A and B under V e l a X 1 ( 1.77 M , 9.56 km) with a = 0.001 km 2 and b = 0.0015 km 2 for different suitable values of n with λ = 2 .
Table 3. Calculated values of c , A and B under V e l a X 1 ( 1.77 M , 9.56 km) with a = 0.001 km 2 and b = 0.0015 km 2 for different suitable values of n with λ = 2 .
ncAB p rc / ρ c
1.80.007104.8700260.0655619.2046
2.20.00686−6.9292100.0661359.5856
2.60.00662−2.5192580.0667209.9767
3.00.00639−1.3911690.06731610.3777
3.40.00617−0.6248520.06792210.7924
3.80.005960.9983920.06853911.2200
4.20.00576−2.6746530.06916711.6628
4.60.00556;−1.1878360.06980712.1187
5.00.00537−0.7617500.07045912.5898
5.40.00518−0.4609700.07112313.0764
5.80.005000.0718000.10633413.5785
6.20.00483−1.0275550.07249014.0983
6.60.00467−0.5153080.07319314.6350
7.00.00451−0.3467810.07391115.1897

Share and Cite

MDPI and ACS Style

Waheed, S.; Mustafa, G.; Zubair, M.; Ashraf, A. Physically Acceptable Embedded Class-I Compact Stars in Modified Gravity with Karmarkar Condition. Symmetry 2020, 12, 962. https://doi.org/10.3390/sym12060962

AMA Style

Waheed S, Mustafa G, Zubair M, Ashraf A. Physically Acceptable Embedded Class-I Compact Stars in Modified Gravity with Karmarkar Condition. Symmetry. 2020; 12(6):962. https://doi.org/10.3390/sym12060962

Chicago/Turabian Style

Waheed, Saira, Ghulam Mustafa, Muhammad Zubair, and Asifa Ashraf. 2020. "Physically Acceptable Embedded Class-I Compact Stars in Modified Gravity with Karmarkar Condition" Symmetry 12, no. 6: 962. https://doi.org/10.3390/sym12060962

APA Style

Waheed, S., Mustafa, G., Zubair, M., & Ashraf, A. (2020). Physically Acceptable Embedded Class-I Compact Stars in Modified Gravity with Karmarkar Condition. Symmetry, 12(6), 962. https://doi.org/10.3390/sym12060962

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop