Abstract

The possibility of strange stars mixed with dark energy to be one of the candidates for dark energy stars is the main issue of the present study. Our investigation shows that quark matter atcs as dark energy after a certain yet unknown critical condition inside the quark stars. Our proposed model reveals that strange stars mixed with dark energy feature a physically acceptable stable model and mimic characteristics of dark energy stars. The plausible connections are shown through the mass-radius relation as well as the entropy and temperature. We particularly note that a two-fluid distribution is a major reason for the anisotropic nature of the spherical stellar system.

1. Introduction

The study of the different aspects of quark matter has drawn attention among astrophysicists and particle physicists in the last two decades. Bhattacharyya and his coworkers [13] proposed that after the few microseconds of the big bang, the universe undergoes a quark-gluon phase transition which may be a process of origin and survival of the quark matter. The nature of the confining force which triggered this phase transition is almost unknown. Witten [4] first considered that at the critical temperature  MeV, a small portion of the colored objects like quark and gluon leaves hadronization through a phase transition to form the colored particles called quark nuggets (QNs). These QNs are made of , , and quarks and have a density that is a few times higher than the normal nuclear density. This was further investigated by [57]. For both cases of the compact stars and the early universe, Ghosh [8] investigated the role of quark matter in the phase transition. It is believed that quark matter exists in the core of neutron stars [9], in strange stars [10], and as small pieces of strange matter [11]. Investigations by Rahaman et al. [12] and Brilenkov et al. [13] led to an interesting and important result according to which quark matter plays the same role as dark energy on the global level. It is worth mentioning that for the last several years, LHC at CERN has been trying to recreate the situation encountered before and early in the hadronization period, performing collisions of the relativistic nuclei [14].

Chapline [15] proposed that a gravitationally collapsing compact star with a mass greater than a few solar masses, a quantum critical surface for spacetime, and an interior region consisting of a large amount of dark energy compared to the ordinary spacetime can be defined as a dark energy star. He also predicted that this surface of a compact dark energy star is a quantum critical shell [16]. When ordinary matter which has energy beyond the critical energy enters the quantum critical region, it decays into constituent products and corresponding radiation is emitted in the outward direction perpendicular to that quantum critical surface. For the matter having energy less than , these constituent products and the radiation can pass through that critical surface and follow a diverging geodesic inside the star. For compact objects and compact stars at the center of galaxies, the energy of the quarks and gluons inside the nucleons is higher than the critical energy [17]. According to the Georgi-Glashow grand unified model, nucleons decay in a process in which a quark decays into a positron and two antiquarks. So the observation of the excess positron in the center of the galaxy may validate the presence of dark energy stars.

In the present article, we have tried to investigate the possible connection between the proposed quark star model mixed with dark energy and dark energy stars. Following the works of Rahaman et al. [12] and Brilenkov et al. [13], we propose that quark matter is one of the possible candidates for dark energy. We are considering an anisotropic quark star model where we assume that the dark energy density is linearly proportional to the quark matter density. The proposed stellar configuration consists of two kinds of fluid: (i) quark nuggets (QNs) and (ii) dark energy having a repulsive nature. We have avoided any interaction between the fluids for the sake of the simplicity of the model. To describe the equation of state (EOS) for the effective fluid of the stellar model, we have used the MIT bag EOS. At this point, a short comment on the mechanism that allows for the accumulated dark energy inside the star is needed. A possibility would be that a mechanism similar to the one responsible for dark matter accretion might be active. In such a case, a possible candidate would be weakly interactive massive particles (WIMPs) that can accrue. At the same time, they have their antiparticles so that they can annihilate to create a heat source. Since dark energy density is expected to be very high, this heating would control the internal structure of the stars. However, one should keep in mind that WIMPs interact weakly with baryons, whereas a dark energy particle most likely exhibits only gravitational interactions with baryons. Therefore, eventually, the accretion rate, as mentioned earlier, would be much lower than expected. In any case, all this is a theoretical conjecture. In reality, the mechanism that allows for the accumulated dark energy inside the stars remains unknown, which is also beyond the scope of the present work to discuss and needs further investigation.

It is worth mentioning that the anisotropy in compact stars may arise due to phase transition, the mixture of two fluids, the existence of type 3A superfluid, bosonic composition, rotation, pion condensation, etc., at the microscopic level. In the present study of the proposed anisotropic two-fluid model, we consider the compact stars PSRJ 1614-2230, VelaX-1, PSRJ 1903+327, Cen X-3, and SMCX-4 as testing candidates.

2. Field Equations

To describe the interior of a relativistic compact star mixed with dark energy, we are considering the following spacetime line element where and depend only on the radial coordinate .

The energy-momentum tensor components of the proposed two-fluid model are given by where

, , and represent the quark matter density, the radial pressure, and the tangential pressure, respectively, whereas and represent, respectively, the dark energy density and the radial pressure within the star. On the other hand, by , , and , we represent, respectively, the effective energy density, the effective radial pressure, and the effective tangential pressure of the matter distribution of the stellar system.

Using Equations (1) and (2), the Einstein field equations for the present spherically symmetric anisotropic compact star read

Now, to solve the Einstein equations for our spherical distribution, we consider the following ansatz: (i) the strange quark matter (SQM) distribution obeys the equation of state (EOS) of the phenomenological MIT bag model, i.e., where is the bag constant, and (ii) the dark energy radial pressure is related to the dark energy density as

Here, the second ansatz represents the EOS of the matter distribution, which is called “degenerate vacuum” or “false vacuum” [1821].

To close the system of equations following Mak and Harko [22], we consider a functional form for and as follows: where both the functional forms for and are considered in such a way that they have a maximum value at the stellar core and decrease gradually to attain their minimum values at the surface. Note that the density for strange quark matter (SQM) at the surface has a nonzero finite value, i.e., , which has been taken care of in the assumed functional form of (see Equation (11)). denotes the central density of SQM, and controls the amount of DE matter distribution corresponding to SQM, which also plays an important role in in determining of phase transition from the quark matter to the dark energy. denotes the total radius of the stellar object.

Now, from Equation (6), we have where is the mass function of the star, which is defined as follows:

The effective gravitational mass of the star can be found using Equations (11), (12), and (14) as

Taking a vanishing radial pressure at the surface, we find from the ansatzes (i) and (ii) and Equations (4), (11), and (12)

Using Equations (6)–(8), the complete set of structure equations is given by

Now, substituting Equations (7) and (9)–(13) into (8), one can easily get an expression for which yields the expected form of anisotropy () as follows:

3. Physical Features of the Model

From the above expression for the physical parameter (see Figure 1), following the method of [23] and after using Equations (15) and (16), we obtain the following equation in order to get the maximum anisotropy at the surface:

Solving Equation (20), using different observational mass values for the various stars considered and with the choice of the parametric values of the bag constant as 83 MeV/(fm)3 [24] and , we get different values for the radius of the star. We choose only that value of which is physically valid and consistent with the Buchdahl condition [25] and find that the anisotropy is maximum at the surface of the star. It is found that the central pressure () is  dyne/cm2 for SMC X-1 due to  MeV/fm3 and .

According to Buchdahl [25], the maximum allowed mass-radius ratio for a static spherically symmetric compact star is . Later Mak and Harko [26] came up with a more generalized expression for the same mass-radius ratio. Now in our model, the effective gravitational mass which is defined by Equation (15) is given as . The variation of the effective mass with the radius of the star is shown in Figure 2.

Now the maximum value of the effective mass of the star corresponding to can be derived using the equation . Similarly, using the equation , we can derive the maximum radius for . Due to  MeV/fm3 and , the maximum effective mass and the corresponding radius are 3.582 M and  km, respectively.

To derive the entropy and the temperature of the stellar model, we are turning to the Gibbs relation, , where is the local entropy density, is the local temperature, is the chemical potential, and is the number density of the matter distribution inside the ultradense star. Let us assume for the sake of simplicity that the matter distribution inside the stellar configuration is isotropic in nature while the value of is negligible. Hence, the Gibbs relation becomes

Now using the first and the second laws of thermodynamics along with the first ansatz, we have the following relation: where is the volume of the stellar configuration. Since and is a total differential, one may find from Equation (22) where is the integration constant and . Here, represents the famous Stefan-Boltzmann constant.

Hence, after some algebra, one can find the entropy density as which provides a basic idea about the total entropy of the compact stellar system.

The variation of temperature in the interior region of the different strange star candidates is shown in Figure 3. From this figure, we see that the temperature in the central region is maximum while it decreases with the radial coordinate and becomes minimum at the surface, which is physically acceptable. For the stellar configurations, we find that the temperature is higher than the Fermi melting point ( K) for quarks. Hence, all the quark matter in the ultradense compact star remains in the form of quark-gluon plasma.

For the anisotropic static stellar configuration, though radial pressure vanishes at the surface, the tangential pressure does not. However, as the radial pressure is continuous at the boundary, we already satisfy Synge’s junction condition [27] in the case of static spherical symmetry. In the boundary, the interior solution and the exterior Schwarzschild solution should match in order to satisfy the fundamental junction condition. The metric coefficients are continuous at the surface where . The second fundamental form is also continuous on the boundary surface. Now the intrinsic stress-energy tensor at the boundary surface (where ) can be defined as the surface stresses, i.e., the surface energy and the surface tangential pressures which in the present situation are given as and . Thus, the complete spacetime is given by our interior metric and the exterior Schwarzschild metric, which are matched smoothly on the boundary surface .

4. Discussion and Conclusions

Under the proposed, model we have presented a data set for the physical parameters of some strange star candidates in Table 1.

Let us highlight the major results of the proposed model: (i) the quark matter is converting into dark energy under a certain critical condition; (ii) with the presence of dark energy inside a strange star, the latter behaves like a dark energy star; (iii) the high temperature distribution (>Fermi melting point for quarks) inside the star confirms the presence of quark matter in the form of quark-gluon plasma; (iv) all the physical and structural features of the proposed ultradense strange star model match well with a dark energy star as suggested by Chapline [15]; (v) some of the physical tests, such as the energy conditions, the TOV equations, and the sound speed constraint, are found to be satisfied in the presented model, and thus, the model has a stable configuration in all respects. According to Herrera [28] and Andréasson [29], in order to form a physically acceptable matter distribution, the quark matter also has to respect the condition , where represents the radial sound speed of the quark matter. This leads to the result that the acceptable value of lies in the range .

An obvious query regarding the present investigation may be as follows: are there any current or future missions that could provide us with some data that would allow us to probe the inner structure of relativistic compact stars? The recent observation of two gravitational waves (GW) events, GW170817 [30] and GW190425 [31], put a well-defined constraint on the EOS of the neutron stars as , where is the effective tidal deformability of the binary neutron star system, which leads to ruling out of the EOS that supports radii above 13 km for the  M neutron stars. In fact, observational data from GW favoured soft EOS and ruled out stiff EOS, such as ms1b, h4, etc. On the other hand, the observation data from the Neutron Star Interior Composition Explorer (NICER) [32] sets stringent constraints on the radius of the neutron stars, which supports the presence of stiffer EOS. Interestingly, this dichotomy of the results from GW observations and NICER data becomes helpful to understand the true nature of EOS of the compact stars and their interior structure. Recent observational constraints suggest that primordial black holes (PBHs) with mass scales~10-12 M can account for the vast majority of dark matter (DM) in the universe. We need an increase in primordial scalar curvature perturbations to the order of at the scale to produce these PBHs [33]. Based on the analysis of years of data from the Kepler satellite, searching for short-duration bumps induced by gravitational microlensing, unique limits on the acceptable masses of a DM halo composed of PBHs and/or any other massive compact halo object have been imposed [34]. The masses range from  M to 10-7 M. Prospective analyses of the whole Kepler data set should find PBH DM or rule out some of this spectrum, and space missions like WFIRST Wide-Field Infrared Survey Telescope (WFIRST) have the power to reach an order of magnitude more.

Let us comment on some shortcomings of the proposed model, which do exist as follows: (i) we proposed that quark matter is converting into dark energy, but we cannot predict under which condition this is actually happening, and (ii) in the range , we consider as a constant parametric term. Therefore, it is difficult to predict through this model whether the dark energy in the stellar configuration remains constant or it varies with time within the provided range of . In connection to this comment, one can specifically note that the metric coefficients are considered independent of time in the present model. However, all these issues can be considered in a future investigation.

Data Availability

No new data has been generated for this manuscript.

Disclosure

The presentation of the manuscript in “Can strange stars mimic dark energy stars?” with arXiv: 1611.02253 is a part of the present article with revised authorships and version [35].

Conflicts of Interest

There is no conflict of interest in connection to the present work.

Acknowledgments

We all are grateful to Professor A. DeBenedictis, Simon Fraser University, for several valuable suggestions. FR and SR are thankful to the Inter-University Centre for Astronomy and Astrophysics (IUCAA), Pune, India, for providing Visiting Associateship under which a part of this work was carried out. Research of DD is funded by the C.V. Raman Postdoctoral Fellowship (Reg. No. R(IA)CVR-PDF/2020/222) from the Department of Physics, Indian Institute of Science. SRC would like to thank the Southern Federal University (SFedU) for financial support (Grant No. P-VnGr/21-05-IF).