Anisotropic strange stars in Tolman–Kuchowicz spacetime

We attempt to study a singularity-free model for the spherically symmetric anisotropic strange stars under Einstein’s general theory of relativity by exploiting the Tolman–Kuchowicz (Tolman in Phys Rev 55:364, 1939; Kuchowicz in Acta Phys Pol 33:541, 1968) metric. Further, we have assumed that the cosmological constant Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varLambda $$\end{document} is a scalar variable dependent on the spatial coordinate r. To describe the strange star candidates we have considered that they are made of strange quark matter distribution, which is assumed to be governed by the MIT bag equation of state. To obtain unknown constants of the stellar system we match the interior Tolman–Kuchowicz metric to the exterior modified Schwarzschild metric with the cosmological constant, at the surface of the system. Following Deb et al. (Ann Phys 387:239, 2017) we have predicted the exact values of the radii for different strange star candidates based on the observed values of the masses of the stellar objects and the chosen parametric values of the Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varLambda $$\end{document} as well as the bag constant B\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathcal {B}}$$\end{document}. The set of solutions satisfies all the physical requirements to represent strange stars. Interestingly, our study reveals that as the values of the Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varLambda $$\end{document} and B\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathcal {B}}$$\end{document} increase the anisotropic system become gradually smaller in size turning the whole system into a more compact ultra-dense stellar object.


Introduction
Einstein's general theory of relativity (GR) [4] represents a grand tool of gravitation to understand uniquely the febric of the space-time and therefore the behavioral pattern of the cosmic bodies as well as phenomena. The first successful a e-mail: mahmoodkhalid@unizwa.edu.om b e-mail: ddeb.rs2016@physics.iiests.ac.in c e-mail: saibal@associates.iucaa.in d e-mail: kumar001947@gmail.com e e-mail: sourav.rs2016@physics.iiests.ac.in field theoretical application of GR to black hole solution by Schwarzschild [5] and to the structure of the universe by Einstein himself [6] made pavement for a new avenue in the research fields of astrophysics and cosmology.
In the aspect of Cosmology, the cosmological constant Λ as introduced by Einstein in GR to match the Mach principle and to have a non-expanding static solution of the universe, becomes very significant. But thirteen years later the pioneering observations by Edwin Hubble [7] divulge that actually universe is expanding, which forced Einstein to abandon the concept of the cosmological constant as a 'blunder'. Later, in 1998 scientists [8,9] came up with some ground breaking scientific evidences by observing the high redshift Type Ia supernova led to the decision that expansion of the universe is actually accelerating in manner, which reinstated the idea of Einstein's cosmological constant. Although the erstwhile cosmological constant was conjectured as a constant quantity but gradually it appears that Λ is actually a scalar variable. This changes with time and the obtained decreasing value of Λ (≤ 10 −56 cm −2 ) [8,9] as obtained from the observational evidences agree with the concept of its variable nature. Now, as Λ has been used in the field equations as a scalar variable, hence it should be a function of the spatial coordinate too. Importantly, likewise the framework of cosmology where variation of Λ with the time has shown a significant effect, many researchers [10][11][12][13][14][15] have successfully studied the effect of the variation of Λ in the realm of astrophysics with respect to the spatial coordinate.
As a consequence of gravitational collapse we get a panorama of stellar formations which include white dwarfs, neutron stars and black holes (in the observational level) and also quark stars (in the theoretical level). In the present work, we are interested to study the strange stars, which are made of strange quark matter (SQM). The possibility of the existence of the hypothetical ultra-dense strange stars was first conjectured by several researchers [16][17][18][19]. The SQM, which is made of the equal number of up (u), down (d) and strange (s) quarks, is assumed to be the true ground state for the confined hadrons [17,18,20,21], as predicted by the strange quark matter hypothesis [22][23][24][25][26][27]. Though it is quite difficult to distinguish neutron stars from the strange stars in the basis of mass and radius but the strange stars, which are ultra-dense and massive stellar objects, have formed a different hypothetical branch for the compact stars. However, likewise the neutron stars and white dwarfs [28][29][30], strange stars do not exhibit continuum in the equilibrium configurations. Interestingly, the neutron star equation of state (EOS) failed to explain the compactness of the recently observed compact stellar objects like 4U 1820 − 30, S AX J 1808.4 − 3658, 4U 1728 − 34, Her X − 1, R X J 185635 − 3754 and P S R 0943 + 10, etc., whereas SQM EOS [18,19,[31][32][33][34] has satisfactorily explained the compactness of the stellar candidates.
The compact stars are basically considered as a spherically symmetric and isotropic ultra-dense stellar objects. However, isotropy may be a favoured feature but need not be a general characteristics of the stellar objects. Anisotropic factor (Δ = p t − p r ) which deals with the inhomogeneity in pressures into two components -radial pressure and tangential pressure, is actually a function to consider the internal situation of the star compared to the idealized isotropic case. The concept of anisotropy was originally proposed by Ruderman [35] and later on by several other scientists [36][37][38]. As the reasons behind the pressure anisotropy different factors are thought to be responsible, such as very high density region in the core region, various condensate states (like pion condensates, meson condensates etc.), superfluid 3A, mixture of fluids of different types, rotational motion, presence of magnetic field, phase transition etc. The basic ideas involved in the anisotropy and their applications in the diverse fields are available in the recent literature [3,[39][40][41][42][43][44][45][46][47][48][49][50][51][52][53][54][55][56][57].
It is argued by Maurya [58] that the chances of having anisotropy is much higher in the compact stars because the interaction among the particles is highly relativistic and they become very random to maintain any uniform distribution throughout the region. This relativistic nature of particles in the compact stars could be one of the possible reasons for giving birth of significant anisotropy in the compact stars. It is therefore understood that the anisotropic force inside the stellar systems make the compact objects more compact than the isotropic condition which eventually makes the possible transition of a neutron star to strange star.
Under the above background, our motivation in the present paper is to study a spherically symmetric anisotropic stellar model and the effects of the cosmological constant, which is assumed to vary with the spatial coordinate, on the stellar system. The outline of this study is as follows: We provide the basic stellar equations to describe the assumed spherically symmetric anisotropic system and the Einstein field equa-tions for the anisotropic stellar source in Sect. 2. In Sect. 3 we have explained the necessity of the choice of appropriate EOS for SQM, whereas the solutions of the Einstein field equations and the expression for the different physical parameters has been featured in Sect. 4. In Sect. 5 we have discussed junction conditions to derive the unknown constants. Several tests, viz., energy conditions, stability (via conservation equation, Herrera condition, and adiabatic index), compactification factor and redshift, have been conducted to study the physical validity of the stellar model in Sect. 6. The last Sect. 7 is devoted as a platform for providing some salient features and concluding remarks along with a brief discussion.

Basic stellar equations
To describe the spacetime of a compact stellar objects, we consider the spherically symmetric line element as Literature survey [59] shows that some interesting static spherically-symmetric perfect fluid solutions in canonical coordinates for the above metric were proposed by Tolman [1], Patwardhan and Vaidya [60] and Mehra [61] where λ(r ) = ln(1 + ar 2 + br 4 ), and by Kuchowicz [2] and Leibovitz [62] where ν(r ) = Br 2 + 2 ln C with a, b, B and C as positive constants. The metric, therefore, with the mentioned metric potentials will be referred later on simply as the Tolman-Kuchowicz (TK) metric.
The Einstein field equations for the present anisotropic stellar system is given as where T where p r , p t and ρ represent the radial pressure, tangential pressure and matter-energy density, respectively. Here, u i is the 4-velocity and v i is the radial 4-vector of the fluid element. We define the mass function of the present spherically symmetric stellar system as follows We have considered the erstwhile Einsteinian cosmological constant as radially dependent, i.e. Λ r = Λ(r ) and hence the Einstein field equations for the metric (1) with energymomentum tensor T i j are obtained as where ' ' denotes the differentiation with respect to the radial coordinate r . In above field equations we have considered geometrized units G = c = 1.
Now, the equation of continuity, T i j;i = 0 of the present anisotropic system reads

Equation of state
Now to describe the interior structure of the stellar system we assume the radial pressure is directly related to the energy density by a relation, known as the equation of state (EOS). Once we define EOS, immediately the Einstein field equations (6)-(8) can easily be tackled. In the present investigation, we have assumed the MIT bag EOS to describe anisotropic SQM, made of up, down and strange quark only [23]. Many researchers have successfully used the MIT bag EOS in their recent studies [63][64][65][66][67][68][69][70] for SQM. The MIT bag EOS is given as follows where α is a constant and it has value 1/3 for the massless strange quarks whereas it is equal to 0.28 for the massive strange quarks having mass 250 MeV [71]. Here B represents the bag constant. It is worth mentioning that to include all the corrections due to the energy and pressure functions of SQM, in the simplified and linear MIT bag model EOS an ad hoc bag function has been introduced. In the present study, the mathematical analysis has been carried out by considering α = 0.28 and B = 60 MeV/fm 3 .

Solution of the Einstein field equations
By using the metric (1) and Eqs. (6)- (8) and (3), we can write the expressions for radial pressure p r , tangential pressure p t and energy density ρ in the following form as where β = 32 Bπ α. We have shown behaviour of the density function in Fig. 1. Again, the variation of the pressures viz., p r and p t against the radial coordinate have been featured in Fig. 2. From Figs. 1 and 2 we find ρ, p r and p t decrease monotonically from the maximum value at the centre to the minimum value at the surface and confirm regularity of the achieved solutions. We have presented variation of the varying cosmological constant against the radial coordinate in Fig. 3 Hence, the anisotropic stress Δ(r ) = p t − p r for the present system is as follows We have shown the variation of the anisotropic stress Δ with respect to the fractional radial coordinate r/R in Fig. 4.

Junction conditions
The first fundamental form of the boundary surface involved in the metric (1) should be the same whether obtained from the interior or exterior space-time and this guarantees that for some coordinate system the metric components g i j will be continuous across the surface. Therefore, we consider that the exterior system is equivalent to the metric which is the modified Schwarzschild metric with the cosmological constant and hence turns into the usual Schwarzschild metric for vanishing Λ.
The continuity of the metric functions, involved in the interior and the exterior metrics, at the outer boundary of the fluid sphere (r = R), as well as the requirement of matching condition for the radial pressure, i.e. p r = 0 at the surface, immediately provides the constants a, b, B and C in the following forms In the later part of the article it can be observed that tuning of these constants are very important to determine different physical properties of compact stars. The variations of the metric potentials viz., e ν and e λ based on the obtained values of the constants, against the radial coordinate is shown in Fig. 5. Now, many scientists [72,73] have shown that anisotropy is minimum at the surface, on the other hand, many others [15,47,74] have also shown that anisotropy is maximum at the surface. To resolve the confusion that what should be the real feature of the anisotropic stress at the surface Deb et al. [3] have maximized anisotropy at the surface without pre-assuming it as maximum or minimum at r = R and conclusively showed that for a stable and spherically symmetric anisotropic stellar system the anisotropic stress must be max-imum at the surface as an inherent property of the compact stars. In the present study, following Deb et al. [3] we have predicted different physical parameters for the strange star candidates, including the exact values of radius by considering the observed values of mass of the strange star candidates for B = 60 MeV/fm 3 (see Table. 1).

The physical properties of the stellar model
In this section we will discuss physical validity of the achieved solutions based on the few physical tests viz., energy conditions, equilibrium of forces, Herrera cracking concept, etc.

Energy conditions
The energy conditions, viz., Null Energy Condition (NEC), Weak Energy Condition (WEC) and Strong Energy Condition (SEC) will be valid only when the following inequalities hold simultaneously for our system given as S EC : ρ + p r ≥ 0, ρ + p r + 2 p t ≥ 0, The energy conditions are plotted in Fig. 6 which features that as our system is consistent with all the inequalities simultaneously, it confirms that the achieved solution is physically viable.
6.2 Stability of the stellar model

TOV equation
The Tolaman-Oppenhiemer-Volkoff (TOV) ( [1,81]) equation for the anisotropic matter with variable cosmological constant (Λ r ) in radial direction are given as It is Obviously that the above generalized TOV equation describes the equilibrium condition for the strange star subject to the gravitational (F g ) and hydrostatic (F h ) plus another anisotropic force (F a ) due to the anisotropic nature of the stellar object. Now, the above equation can be written via different forces as with (28) Figure 7 shows the counterbalancing pattern of different forces to attain equilibrium of the spherical configuration. Figure 7 also features that in each cases the repulsive anisotropic force F a acts along the outward direction to counter balance the combined effect of the attractive forces viz., gravitational force F g and the hydrodynamic force F h , which act along the inward direct.

Herrera condition
To verify stability of our anisotropic model we have plotted the radial (v 2 r ) and transverse (v 2 t ) sound speeds in Fig. 7. It is observed that the inequalities 0 ≤ v 2 r ≤ 1 and 0 ≤ v 2 t ≤ 1 are satisfied everywhere within the stellar system and obeys the causality condition. Now the radial and transverse sound speeds are given by where On the other hand, following the cracking concept [82,83] we checked stability of the local anisotropic configuration by finding out the potentially stable region where the radial speed of sound is greater than the transverse speed of sound, i.e. v 2 t − v 2 r ≤ 1. From Fig. 8 one can note that 0 ≤ v 2 r ≤ 1 and 0 ≤ v 2 t ≤ 1 so that | v 2 t − v 2 r |≤ 1. Thus the present compact stellar model provides a stable configuration of the spherical system by satisfying both the causality condition and the Herrera cracking concept.

Adiabatic index
To check the dynamical stability of the stellar system against an infinitesimal radial adiabatic perturbation we follow the pioneering work by Chandrasekhar [84] where he predicted that for relativistic stellar system to be stable the adiabatic index Γ should be greater than 4/3. Note that later on many researchers [85][86][87][88] have successfully tested the prediction by Chandrasekhar for both the isotropic and anisotropic stellar objects. The adiabatic index is defined as In Fig. 9 we have successfully shown the behaviour of the adiabatic index Γ which features that in the present case Γ is well above the critical limit 4/3 and confirms stability of the system against an infinitesimal radial adiabatic perturbation.   The mass function of the present spherical stellar system is given as Fig. 9 Behavior of the adiabatic index Γ as a function of the fractional radial distance for the strange star candidate L MC X − 4 is shown where In terms of the mass function, let us now define the compactification factor as follows Hence the surface redshift is given by We have plotted variation of the compactification factor and redshift function in the left and right panel of Fig. 10.

Discussion and Conclusion
In the present study, we have considered an anisotropic fluid sphere in the framework of Einstein's general theory of relativity. The spherical spacetime is assumed to be of the Tolman-Kochowicz type and solutions of the Einstein field equations are found by using this metric. One can observe that all the parameters involved in the solutions set are viable within the specified physical conditions. Further, following the work of Stergioulas [71] we have analysed the present system based on the phenomenological MIT bag model EOS which is suitable for the massive strange quarks. To derive the unknown constants and radius of the strange star candidates in the present study we have used the observed values of the mass of the strange stars candidates as presented in Table 1 and also assumed parametric values of B and Λ given as B = 50, 60 and 70 MeV/fm 3 , and Λ = 0.004, 0.005 and 0.006. For further analysis of the achieved solutions, we have considered L MC X − 4 as the representative of the strange star candidates.
However, we would like to highlight some of the salient features of the present models which are as follows: 1. Density and pressure: We have shown the profile of density in Fig. 1, which features that for B = 60 MeV/fm 3 and the chosen values of Λ, the nuclear density inside the compact stellar model is much higher than the normal nuclear density ρ normal = 2.3 × 10 14 g/cm 3 . Hence, it confirms that the stellar system corresponds to the ultra-dense strange stars [35,89,90]. In Fig. 2 we have shown the variation of p r and p t in the left and right panel, respectively. On the other hand, we have shown variation of the anisotropic stress in Fig. 4, which indicates that the anisotropy for our system is minimum at the centre and maximum at the surface as predicted by Deb et al. [3]. Interestingly, Fig. 4 also features that as the values of Λ increase anisotropic stress of the systems increase consequently. We have presented the variation of the cosmological constant which is a scalar variable dependent on the spatial coordinate in Fig. 3. The variations of the metric potentials, i.e., e ν and e λ are depicted in the left and right panel of Fig. 5, respectively. Further, Figs. 1,2 and 5 confirm that our system is free from all sorts of singularities viz., physical as well as geometrical singularities.

Energy conditions:
We find that for both the chosen cases the energy conditions, viz., NEC, WEC, SEC and DEC are consistent with our system. We have shown variation of the energy conditions with respect to the radial coordinate in Fig. 6, which confirms physical acceptability of the obtained solutions.
3. Equilibrium of forces: In Fig. 7 we observe that for all the chosen values of Λ the equilibrium of the forces is achieved. In every cases the anisotropic force F a which acts along the outward direction (repulsive nature) counterbalances the combined effect of the gravitational force F g and the hydrodynamic force F h which act along the inward direction (attractive nature).

Stability:
To discuss stability of the system we have examined both the causality condition and the Herrera cracking condition. Fig. 8 exhibits that for all the cases the inequalities 0 ≤ v 2 r ≤ 1, 0 ≤ v 2 t ≤ 1 and | v 2 t − v 2 r |≤ 1 are valid simultaneously and hence confirms stability of the stellar system in terms of the sound velocity of the system. Again, Fig. 9 reveals that for our system the adiabatic index Γ is less than the critical limit 4/3 which confirms that our system is  completely stable against the infinitesimal radial adiabatic perturbation.

Compactification factor and redshift:
We have featured variation of the compactification factor and redshift function in the left and right panel in Fig. 10, respectively. For the chosen values of B and Λ we have derived the factor 2M/R and surface redshift Z s for different strange star candidates and presented their values in Table 1. We find from these Tables 2 and 3 that the values of 2M/R and Z s increase gradually with the increasing values of Λ and B. For all the cases as 2M/R < 8/9, so our system is consistent with the Buchdahl condition [75]. Present investigation reveals that the surface redshift of the different stellar candidates are within the range 0.58 − 0.41 and such high redshift values indicate that the stellar candidates under study are actually strange stars [91]. The obtained values of Z s due to different strange star candidates are well within the provided limits given as 0 < Z s ≤ 1 [46,47,[92][93][94] and Z s ≤ 2 [75,95,96] which provides physical consistency of the proposed stellar system.
6. Comparative discussion: With the motivation to analysis the effects of the variable cosmological constant Λ(r ) and bag constant B on the anisotropic stellar system, we have presented Tables 2 and 3. For the fixed chosen value of B as B = 60 MeV/fm 3 we find from Table 2 that as the values of Λ increases the different physical parameters viz., ρ c , p c , 2M/R and Z s increase gradually, whereas the radii of the stellar candidates decrease consequently. This observation leads to the conclusion that as Λ increases the stellar systems become smaller in size turning the strange star candidates into a more compact stellar object. The same effects can be observed for the increasing values of B as presented in Table 3. Hence, with the increasing values of Λ and B the systems become gradually diminutive along with the rising values of density and surface redshift, which are the perfect situations to describe the strange star candidates.
The pioneering detection of gravitational waves (GW170817) [97] from the coalescence of a neutron star binary system, has given a realistic opportunity to study the properties of the matter inside the stars in such extreme condition which cannot be recreated in the terrestrial laboratories. In their recent study, Abbott et al. [98] have presented the radii of the two neutron stars. The authors have shown that based on only the LIGO and Virgo data the two neutron star radii are R 1 = 10.8 +2.0 −1.7 km and R 2 = 10.7 +2.1 −1.5 km for the heavier and lighter stars, respectively at the 90% confidence level. However, for the neutron stars having masses more than 1.97 M which includes the electromagnetic observa-tions Abbott et al. [98] predicted that radii of the two neutron stars are exactly the same, i.e. R 1 = 11.9 +1.4 −1.4 km and R 2 = 11.9 +1.4 −1.4 km at the 90% confidence level. They also predicted that at twice the nuclear saturation density the pressure should be 3.5 +2.7 −1.7 ×10 34 dyne/cm 2 . In the present study, our model reveals that for a hypothetical strange quark star having mass 2.2 M (say) the radius would be 9.981 km and central pressure would be 4.836 × 10 35 dyne/cm 2 , where the value of B is chosen as 60 MeV/fm 3 . It is to note that our model predicts the stars having more dense and smaller in size compared to the neutron stars, however, this situation can easily be well accommodated within the range of data as predicted by Abbott et al. [98].
For the last few years, the scalar-tensor theories of gravity become popular as an extended theory of General Relativity (GR) to address several puzzles, such as the dark energy problem in the realm of Cosmology and the dark matter problem in the realm of astrophysics. In this line, f (R) gravity theory has appeared as one of the promising alternative theories to address those unsolved issues which are studied in detail in the following literature [99][100][101][102][103][104][105][106][107][108]. Now, our present study deals with the cosmological constant Λ which is assumed to be a scalar variable dependent on the spatial coordinate r . Interestingly, we find that the violation of standard GR predictions has much in common with the results achieved by Resco et al. [108] in their work where they have explored the main features and also studied the existence of neutron stars in the framework of f (R) paradigmatic models. Hence, the achieved results from our model have actually encouraged the extended gravity theories as a competent and alternative route to the cosmological constant or dark energy. However, following Dombriz and Dobado [99] we are looking forward to further detailed study on strange stars based on the assumed Tolman-Kuchowicz metric in the framework of extended gravity theory without considering the cosmological constant.
In this context, it is important to mention that though in the present study under the framework of GR we have used the standard asymptotic (or stellar) mass definition as in Eq. (5) [m (r ) = 4π r 0 ρ (r ) r 2 dr] but for modified f (R) gravity theory this standard definition of mass should be avoided. Astashenok et al. [109] in their study have shown that in the case of modified f (R) gravity theory as at the surface of the stars the scalar curvature R does not vanish hence it invites not only the existence of exterior non-Schwarzschild solutions for realistic equations of state but also a new definition for the gravitational mass of the stellar system. As a result, the exterior spacetime region (gravitational sphere) of the stars also contribute to the total gravitational mass which can be perceived by a distant observer.
As a final comment, in this article employing the Tolman-Kuchowicz ( [1,2]) metric and using parametric values of Λ(r ) and B we have successfully presented singularity free and completely stable stellar system which is suitable to describe the strange star candidates exhibiting anisotropic nature.