Strange quark stars: the role of excluded volume effects

We study cold strange quark stars employing an enhanced version of the quark-mass density-dependent model which incorporates excluded volume effects to address non-perturbative QCD repulsive interactions. We provide a comparative analysis of our mass formula parametrization with previous models from the literature. We identify the regions within the parameter space where three-flavor quark matter is more stable than the most tightly bound atomic nucleus (stability window). Specifically, we show that excluded volume effects do not change the Gibbs free energy per baryon at zero pressure, rendering the stability window unaffected. The curves of pressure versus energy density exhibit various shapes -- convex upward, concave downward, or nearly linear -- depending on the mass parametrization. This behavior results in different patterns of increase, decrease, or constancy in the speed of sound as a function of baryon number density. We analyze the mass-radius relationship of strange quark stars, revealing a significant increase in maximum gravitational mass and a shift in the curves towards larger radii as the excluded volume effect intensifies. Excluded volume effects render our models compatible with all modern astrophysical constraints, including the properties of the recently observed low-mass compact object HESSJ1731.


Introduction
Strange quark matter (SQM) is a type of quark matter composed of a roughly equal mixture of up, down, and strange quarks.This matter form is theoretically unique due to its absolute stability, rendering it self-bound at the quantum level.SQM's stability arises from its lower energy per baryon compared to nuclear matter, suggesting that it could be the true ground state of hadronic matter under certain conditions.
The concept of SQM can be traced back to seminal papers by Itoh [1], Bodmer [2], and Witten [3] who proposed that quark matter containing strange quarks could be more stable than ordinary nuclear matter.Their hypothesis, often referred to as the 'strange matter hypothesis', laid the groundwork for subsequent research in this field.These pioneering works suggested that, under extreme conditions like those in the core of neutron stars, quarks may deconfine and form a stable mixture.
To describe the properties of SQM, various effective models have been employed.The most prominent among these is the MIT Bag model, which treats quarks as free particles inside a hypothetical 'bag', confining them with a constant energy density (the bag constant).This model effectively balances the quarks' freedom with confinement, allowing for a simplified yet insightful exploration of SQM characteristics [4].Other popular models include the Nambu-Jona-Lasinio (NJL) model, which incorporates chiral symmetry breaking but does not allow the existence of self-bound matter for the most accepted parametrizations [5,6].In contrast, the quark-mass density-dependent model (QMDDM) offers an alternative approach, where the masses of quarks depend on the baryon number density [7][8][9][10][11][12][13][14][15][16][17][18][19][20].As first shown in [11], this model stands out for achieving the confinement property through quark-mass density-dependence, without the need for an explicit bag constant.

arXiv:2405.18249v1 [nucl-th] 28 May 2024
Strange quark stars, if they exist, are compact objects entirely composed of SQM.These stars significantly deviate from conventional neutron stars in terms of many of their physical properties [21].Theoretical calculations show them to be smaller and denser, characterized by mass-radius relationships that set them apart from their hadronic counterparts.Another key feature distinguishing strange quark stars from neutron stars is their predicted cooling behaviors.Unlike neutron stars, which cool primarily through neutrino emission and photon radiation, strange quark stars could exhibit significantly more rapid cooling rates, a consequence of the direct Urca process being unimpeded in quark matter [22,23].Furthermore, strange quark stars are theorized to exhibit unique gravitational wave signatures.Their different internal structures could result in oscillation modes that, when perturbed, would emit gravitational waves with characteristic frequencies [24][25][26].These could potentially be detected by advanced gravitational wave observatories, offering a novel method to probe the existence of these hypothetical stars.
Before the discovery of pulsars with masses around 2M ⊙ [27][28][29][30][31], simple EOSs for SQM, such as the standard MIT bag model and the standard QMDDM, were consistent with the astrophysical data available at the time, specifically the observed masses and radii of compact stars.However, with the observation of pulsars with ∼ 2M ⊙ the EOSs provided by the simplest versions of those effective models were found to be too soft, meaning they could not account for the high masses of these newly discovered pulsars.To reconcile these models with observations, the introduction of new elements is necessary.In the case of the MIT bag model, the incorporation of color superconductivity [32,33] and repulsive vector interactions [34,35] allowed the construction of strange quark star models with masses above 2M ⊙ .Similarly, changes were made to the QMDDM to account for the observational data of the more massive pulsars.These modifications often involved changes to the mass parameterization of the quasiparticles.For instance, introducing a mass term that increases with the baryon number density allowed the model to achieve M max > 2M ⊙ [20,[36][37][38].
In recent work, we implemented modifications to the QMDDM, with a specific focus on the role of repulsive interactions among quarks [39].This is achieved through the excluded volume effect, which prevents quark quasiparticles from being arbitrarily close to each other.In a dense quark matter environment, as is expected in the core of neutron stars, this effect renders the EOS substantially stiffer.In this study, we will employ the model described in Ref. [39] to explore the properties of strange quark stars.
The paper is organized as follows.Section 2 introduces the mass formula incorporating excluded volume effects, followed by a comparative examination of this parameterization against others reported in the literature.In Section 3, we present an overview of the general aspects of the QMDDM, both with and without the consideration of excluded volume, and derive the EOS for SQM that will be used in the rest of this work.Section 4 is dedicated to presenting our results.In Section 5, we succinctly summarize our key findings and present the conclusions derived from our study.

Mass formula
The QMDDM represents quark matter as a system of quasiparticles in such a way that the energy density and particle number densities maintain the same functional form as in a non-interacting system.This is achieved by using effective masses that incorporate interaction effects.Accordingly, the equivalent mass is modeled as comprising two components: a flavor-dependent (constant) current mass m 0i and an interacting term m I that depends on the medium properties: To accurately represent strongly interacting quark matter as a quasiparticle system, m I must be consistent with some key properties of QCD.This formula must reflect quark confinement and asymptotic freedom, requiring m I to fulfill the following criteria: lim where n B reprresents the baryon number density.An explicit mass form, derived from these principles, has been proposed.This involves expanding the mass expression into a Laurent series with respect to the Fermi momentum p F (centered around p F = 0), retaining leading terms of both positive and negative exponents [40]: Since the Fermi momentum scales with n 1/3 B , the following mass formula is derived [20]: As discussed in [20] and related references, the term that scales with n 1/3 B is critical for the maximum gravitational mass of stellar configurations to reach a value above 2M ⊙ , as required by present pulsar observations.Simply using the first term of Eq. ( 5) proves insufficient for stellar configurations to fulfill pulsar's constraints.However, there is an issue with Eq. ( 5) as the term increasing with n 1/3 B contradicts the asymptotic freedom condition because it diverges when n B → ∞.This is not surprising since Eq. ( 5), being a series expansion, makes sense only for low n B .Despite this, Eq. ( 5) can still be employed because the undesirable behavior associated with this term is not apparent at the typical densities present in the cores of compact stars but becomes significant at much higher densities.This is due to the very small numerical value adopted for the parameter C in Eq. (5).
Another modification to the mass formula was adopted in Ref. [41] to take into account that the quark-quark effective interaction in quark matter should be isospin dependent.According to this, the equivalent quark mass should vary with isospin, an aspect not taken into account in the previous formulas.The form of this isospin dependence remains uncertain and ideally requires determination through nonperturbative QCD calculations.Ref. [41] proposed the following phenomenological parameterization: Here, D I , α, and β are free parameters and τ q represents the isospin quantum number of quarks-with τ q = 1 for u quarks, τ q = −1 for d quarks, and τ q = 0 for s quarks.The isospin asymmetry is parametrized by δ = 3(n d − n u )/(n d + n u ), being n i the particle number density of quarks of the flavor i.Including isospin dependence in the mass significantly affects the symmetry energy of quark matter and the characteristics of strange quark stars, thereby enabling these objects to satisfy the 2M ⊙ constraint [41].
Our approach to the mass formula takes a different path.For simplicity, we neglect the contribution of the isospin contribution.Then, we begin with an expression that explicitly satisfies Eqs.(2) and (3): The latter expression is 'flavor blind' because its dependence on the medium is uniform across all flavors [15].As previously mentioned in references [42] and [43], the parameter a correlates with the form of the potential v(r) governing the quark-quark interaction, ex-pressed as v(r) ∝ r a .While many research papers adopt a = 1, in this study, we opt to treat it as a free parameter and investigate the ramifications of its variation.We will revisit this issue in the discussion of Figure 3.In this formula, we will incorporate all non-perturbative and higher-order perturbative effects of QCD exclusively through the 'excluded volume' effect, rather than introducing additional terms into the mass formula.As discussed in Ref. [39], the underlying concept is that particles, when in close proximity, effectively 'exclude' a certain region around them.This exclusion is due to repulsive interactions among the particles, stemming from a combination of non-perturbative phenomena associated with confinement, higher-order perturbative QCD interactions, and the dynamic behavior of the strong coupling constant.
In practice, the excluded volume effect is incorporated by replacing the system's volume V with the available volume Ṽ, accounting for the space b(n B ) occupied by each quasibaryon: being N B the total baryon number of the system.When implementing this prescription, various thermodynamic quantities are altered, such as the Fermi momenta [39].Additionally, since the quasiparticle mass depends on the baryon number density, Eq. ( 7) is reformulated as: Introducing the factor q, defined as: the mass formula with excluded volume effects then becomes: This mass formula was introduced in Ref. [39] and will be used throughout the rest of this work.In Eq. ( 10), the ratio q must necessarily lie between 0 and 1, as the available volume can never be larger than the total system's volume.This means that the parametrization of the excluded volume per particle, b(n B ), cannot be chosen arbitrarily.Instead, it must be carefully chosen to ensure that the restriction 0 < q < 1 is always satisfied.The mass parametrization we have adopted includes non-perturbative QCD effects while maintaining agreement with the asymptotic freedom constraint detailed in Eq. (3).

General overview of the QMDDM
In this section, we present a synopsis of the QMDDM EOS as developed in Refs.[42] and [39].In Sec.3.1 we summarize the EOS for point-like particles and in Sec.3.2 we show how to incorporate excluded volume effects in that EOS.In Sec.3.3 we use the results of Secs.3.1 and 3.2 to construct the EOS for strange quark matter.

The EOS for point-like particles
We will describe the system as a mixture of non-interacting quarks species with effective masses M i and free electrons.The total Helmholtz free energy depends on both the system's volume and the number of particles {N i }, and it is the sum of the contributions from each species: where, e χ(x e ) (electrons), (13) with g = 6 and g e = 2.The function χ(x) is defined by: where x is the dimensionless Fermi momentum given by: x e = 1 m e 6π 2 n e g e 1/3 (electrons), and n i = N i /V is the particle number density of the i-species.
At T = 0, the energy density coincides with the Helmholtz free energy per unit volume.Therefore: where: e χ(x e ) (electrons).
The total pressure is obtained from p = − ∂F/∂V| {Nj} , and is the sum of the partial pressures of each species: being e ϕ(x e ) (electrons), where ϕ(x) is defined by: The 'bag constant' B i is given by: with β(x) defined by: In the flavor blind model the mass has the form M i = m 0i + m I (n B ), i.e. the interaction term m I of the mass depends only on n B .Therefore: Thus, for flavor blind mass formulas, Eq. ( 22) reads: In our previous publication [42], an error was identified in the definition of the bag constant which was further corrected in [39].Eq. ( 22) represents the accurate expression for any mass formula, while Eq. ( 25) applies specifically to the flavor blind scenario.By substituting Eq. ( 22) into (20), and subsequently replacing Eq. ( 20) into (19), we arrive at the same expression for the pressure as given in Refs.[19,20].Finally, the chemical potential of the i-species is derived from For any mass formula, the chemical potential of u, d, and s quasiparticles is given by the following expression (cf.Appendix B of [39]): which coincides with the expressions given in Refs.[19,20].In the flavor blind scenario, the latter equation simplifies to: For the electron component, the chemical potential is given by: It should be noted that the second term of Eq. ( 27) differs from the expression in Eq. (88) of Ref. [42].The version presented above is indeed the accurate one.

Excluded volume effects: general formulas
As shown in Ref. [39], the incorporation of excluded volume effects into any zerotemperature EOS originally formulated for point-like particles in the Helmholtz representation is straightforward.The procedure unfolds as follows.Initially, we select an ansatz for the excluded volume per particle, b(n B ), from which we derive the following quantities: It's worth noting that the functional form of b(n B ) must be selected to ensure that the constraint 0 < q < 1 is always verified.Subsequently, we adapt the expressions for the energy density ϵ pl , pressure p pl , and chemical potentials µ pl,i of point-like particles.This involves rewriting ϵ pl , p pl , and µ pl,i in terms of the modified variable set {n j /q}, and then multiplying them by the factors q, δ, and λ as outlined below: In this manner, we obtain expressions that incorporate the effect of excluded volume, derived from the known expressions for point-like particles.

Strange quark matter EOS
Now, we adopt an ansatz for b consistent with the asymptotic freedom behavior of QCD, i.e., we adopt a formula that allows b to approach zero at asymptotically large densities.Based on the prescription of Ref. [39], the excluded volume is defined as being κ a positive constant.Using Eqs. ( 29), ( 30) and ( 31), the functions q, δ and λ are: Replacing these parameters in Eqs. ( 32), ( 33) and ( 34) we obtain: where As our focus is on describing quark matter within the interiors of cold compact stars, we will enforce the conditions of local electric charge neutrality and chemical equilibrium under weak interactions.Under these conditions, neutrinos are free to exit the system (µ ν e = 0), leading to the following chemical equilibrium conditions: Meanwhile, charge neutrality can be expressed as:

Numerical results for the EOS under compact star conditions
In this section, we will use the EOS introduced in Section 3.3 to explore the properties of strange quark matter and strange quark stars.We will explore various values for the parameters a, C, and κ of the model.The parameter κ is dimensionless.The units of a and C are not fixed and are chosen such that n B is in fm −3 and M i is in MeV in Eq. (11).
Figure 1 illustrates the Gibbs free energy per baryon, G/n B = (ϵ + p)/n B , plotted against pressure.Depending on the chosen parameters a, C, and κ, the value of G/n B at zero pressure may either be above or below the energy per nucleon in the most tightly bound atomic nucleus, 62 Ni, which is approximately 930 MeV.For parameter sets where G/n B is less than 930 MeV at zero pressure and temperature, we are in the case of selfbound quark matter.In these cases, bulk quark matter is stable in vacuum and does not convert into hadronic matter.If this self-bound matter includes all three flavors (u, d, and s) and leptons at p = T = 0, it is known as SQM.In such circumstances, it is conceivable for compact stars composed entirely of quark matter, termed self-bound quark stars, to exist.On the other hand, if G/n B exceeds 930 MeV, the matter is classified as hybrid matter, which is characterized by a transition from a hadronic state at lower pressures to a deconfined state at higher pressures.In such cases, stars containing quark matter would be hybrid stars, featuring a core of quark matter surrounded by hadronic matter.
In Figure 1, we selected parameter sets to ensure that the G/n B curves for uds quark matter are positioned below 930 MeV at zero temperature and pressure.It is important to note that these G/n B curves are sensitive to variations in the parameter κ.Despite this, the value of G/n B at zero pressure-a critical factor in determining the self-bound nature of quark matter-is unaffected by changes in κ.This is evident in the three panels of Figure 1, where the curves converge at the same point within each panel.
Figure 2 presents the stability window of quark matter as determined by the parameters a and C. As previously explained, the value of the Gibbs free energy per baryon at zero pressure does not depend on the value of κ.Therefore, the stability window shown in the figure holds for any value of the excluded volume.The plot includes two significant curves.The solid black curve labeled with e 0 ud = 930MeV represents parameter values at which G/n B at zero pressure is exactly 930 MeV for two-flavor quark matter.The dashed red curve with e 0 uds = 930MeV corresponds to parameter values where G/n B at p = 0 is exactly 930 MeV for three-flavor quark matter.The light green shaded area in between indicates the parameter values for a and C where three-flavor quark matter is more stable than the most tightly bound atomic nucleus, and two-flavor quark matter is less stable than this nucleus (strange matter hypothesis).The EOS models represented by these parameters support the possibility of self-bound quark matter without contradicting the existence of atomic nuclei.Within this region, we have selected three sets of parameters that will be used for the calculations in the following figures, which are identified by red dots in the stability window.The dark green region corresponds to parameter values where two-flavor quark matter (in bulk) is more stable than the most tightly bound atomic nucleus.This area represents a scenario where atomic nuclei could decay into a single bag of deconfined u and d quarks, which is not observed in Nature.Additionally, it is important to note that throughout the stability window, the energy per baryon for matter with three flavors Parameter space for quark matter in bulk.The shaded light green region denotes the parameter combinations for which SQM is the true ground state of matter.The region above the light green area, labeled as 'hybrid matter', corresponds to a scenario where hadronic matter is energetically favored at low pressures, whereas quark matter becomes preferred at higher pressures.The dark green region below, marked as 'nonexistence of nuclei', represents parameter choices where nuclei would be unstable.The same parameter space is valid for any choice of the excluded volume parameter κ.Note that the stability window presented here has some quantitative differences from that shown in [42].The reason for this discrepancy is that Ref. [42] contained an error in the expression of the chemical potential, as previously mentioned at the end of Section 3.1.always remains lower than for two-flavor matter, meaning that the present model does not support the existence of absolutely stable two-flavor matter.This means that in the current model, the dark green region represents a scenario where nuclei would ultimately decay into small bags of quark matter consisting of u, d, and s quarks.One might speculate that the parameters of the dark green region are physically viable because the existence of such small hypothetical configurations might be practically suppressed by large surface and curvature tensions.In such a scenario, self-bound quark matter would be possible in bulk, but suppressed within nuclei.Yet, this paper opts for a more restrictive approach, considering the dark green area as unphysical.In the white region, the parameters a and C render both two-and three-flavor quark matter less stable than the most tightly bound atomic nucleus.Here, quark matter could only exist within the core of compact stars, leading to the formation of what are known as hybrid stars.The range of a values depicted ensures that the speed of sound is always lower than the speed of light throughout the figure.
In Figure 3, we present the total pressure p as a function of the energy density ϵ using the same parameter choices as in Figure 1.In each case, the pressure turns negative at a finite energy density due to the effect of the bag constant B. The EOS is highly sensitive to the parameter κ, with an increase in κ consistently leading to a stiffer EOS.As κ increases, the bag constant diminishes, causing the point at which the pressure of matter vanishes to move to lower energy densities.Although not apparent in the figure due to the chosen scale, the excluded volume effect disappears at high densities, which results in the convergence of the curves for different κ values.
An interesting feature of the parameterizations used in each panel is the change in the concavity of the curves.In panel (a), the curves are concave downwards, in panel (b) they are almost linear, and in panel (c) they are concave upwards.As discussed in Refs.[42] and [43], the parameter a is related to the form of the potential v(r) of the quark-quark interaction, such that v(r) ∝ r a (this parameter a appears in the mass as m I ∝ n −a/3

B
). Accordingly, panel (a) corresponds to a cubic potential, panel (b) to a quadratic potential, and panel (c) to a linear potential.It is important to note that the change in concavity of the  EOS occurs for quadratic interactions.Most research papers use a = 1, corresponding to a linear interaction potential.The speed of sound c s is shown in Figure 4 for the same parameter choices of previous figures.In every scenario, c s approaches the conformal limit, where c s = 1/ √ 3 asymptotically.Since the pressure approaches zero at a finite density, we truncate the curves at the density where p = 0. Notice that for a = 3, c s decreases with increasing n B ; for a = 2, it remains nearly constant; and for a = 1, it increases with n B .We have numerically confirmed that this change in behavior is directly associated with the parameter a (linked to the quark-quark interaction potential), while the parameter C plays a minor role in this feature.This transition in behavior is observed very close to a = 2.When c s is a decreasing function of density (panel a), we observe that as the value of the parameter κ grows, c s is smaller for a given density.Conversely, when c s rises with density, an increase in the value of κ leads to a higher speed of sound at a given density.These differences disappear at asymptotically high densities, where all curves ultimately converge.
In Figure 5, we present the gravitational mass M G versus radius R relationship for strange quark stars, using the same parameter sets as in the previous figures.The following  characteristics are consistent across all panels of Figure 5.For point-like quarks (κ = 0), the maximum M G does not satisfy the 2M ⊙ constraint, and the mass-radius curve does not match the astrophysical constraints established by NICER [30,31,44,45] observations.When taking into account the excluded volume of quasiparticles, the EOS becomes stiffer, as previously mentioned.This leads to an increased maximum gravitational mass and larger stellar radii at a given mass, and the resulting curves satisfy the main astrophysical constraints.Notably, the maximum gravitational mass can reach nearly 2.5M ⊙ when κ = 0.5, indicating a significant rise due to the effects of the excluded volume.Finally, the curves with the excluded volume effect shown in Figure 5 also satisfy the constraints imposed by the observation of HESSJ1731 [46].This object has both a very small mass and radius, making it difficult to describe using hadronic matter models.
In Figure 6 we show the behavior of the tidal deformability Λ as a function of strange star gravitational mass M G across different sets of model parameters.Neutron stars experience significant tidal deformation during the inspiral phase of a neutron star-neutron star merger leaving a detectable imprint on the observed gravitational waveform of the merger.The extent of this effect is expressed in terms of the star's dimensionless tidal deformability, which characterizes the induced mass quadrupole moment in response to an external tidal field.Specifically, the induced mass-quadrupole moment Q ij and the applied tidal field ϵ ij are related to this parameter by The dimensionless tidal deformability parameter can be expressed as: where C ≡ M G /R is the compactness of the star and k 2 is the tidal Love number [47].In the three panels of Figure 6, a general feature observed is the steep decrease in Λ with increasing M G .Alongside this trend, as the parameter κ increases, there is an intensification of repulsive forces within strange quark matter, leading to greater tidal deformability at a fixed M G .The error bar represents the estimated tidal deformability of a 1.4 M ⊙ neutron star from GW170817, which is Λ 1.4 = 190 +390 −120 at the 90% confidence level.Notably, the model with a = 1 aligns more closely with the GW170817 error bar, particularly when compared to the model with a = 3, where compatibility is only observed for κ = 0.5.This difference can be understood through the impact of the parameter C on tidal deformability.As seen in Figure 5, decreasing a shifts the curves toward regions of larger radii, reducing compactness and enhancing tidal deformability.
Finally, Figure 7 displays the gravitational mass of the star as a function of its baryon mass (M B ).The baryon number for a canonical neutron star (M G = 1.4M ⊙ ) can be estimated as A B ≃ 1.4M ⊙ /m n = 1.7 × 10 57 .It is useful to define the equivalent mass of non-interacting baryons, termed the baryon mass (or rest mass) of the star, denoted as M B ≡ A B m B , where m B ≈ m n represents the mass of one baryon.A common choice for m B is to set it equal to the neutron mass, m n .This assumption leads to the equation: where A B57 ≡ A B /10 57 .The values of A B and M B remain constant throughout the evolutionary path of an isolated star.This constancy in baryon mass is essential for understanding the potential transformation of a metastable hadronic star into a strange quark star.In this proposed scenario, it is hypothesized that the strange star does not form directly during a supernova explosion but rather emerges through a phase conversion from another neutron star.During such a transformation, the baryon mass must remain unchanged between the initial and final states.To illustrate this, we have included in Figure 7 the curve for the widely used SLY4 hadronic EOS [48] for comparison.As previously mentioned, the phase conversion of a hadronic star occurs along a vertical trajectory in Figure 7.A critical aspect of this figure is that the SLY4 curve is consistently above those of the strange matter objects, with few exceptions.This positioning indicates that the conversion of a hadronic star into a strange star is an exothermic process.This finding qualitatively aligns with the results obtained using different EOS for strange quark matter [49][50][51].In all panels of Figure 7, the curves demonstrate that as κ increases, the slope of the curve becomes steeper.This suggests that a higher κ value, indicative of stronger repulsive interactions within the neutron star, leads to a higher gravitational mass for the same baryon mass.Generally, the curves in panel (a) exhibit a steeper slope than those in panel (b), which in turn are steeper than those in panel (c).This pattern is associated with the fact that stiffer EOSs (see Figure 3) result in a greater gravitational mass for a fixed baryon mass.Consequently, as the parameter a decreases, the energy released decreases, as evidenced by the narrowing gap between the SLY4 and the SQM curves.

Summary and Conclusions
In this work, we examined the properties of strange quark matter using the recently developed QMDDM with excluded volume effects, as detailed in [39].Our study focuses on the properties of self-bound compact stars formed by this matter.
High-density quark matter at zero temperature is represented as a system of quasiparticles whose effective masses depend on the baryon number density.In the initial version of the model [19,20,42], the mass formula for these quasiparticles phenomenologically incorporates key quark properties such as color confinement, asymptotic freedom, and chiral symmetry breaking/restoration.This is achieved through an effective mass that diminishes with increasing n B .A recent improvement in this model accounts for potential repulsive interactions between quasiparticles by considering them as surrounded by an 'excluded volume'.In our current study, we contrasted our mass formula with those proposed in earlier studies, particularly those that incorporate an additional term that grows with n B [20].It is argued that this additional term represents leading-order perturbative interactions in QCD.However, such an additional term contradicts the property of asymptotic freedom which requires that the mass of quasiparticles tend to the current quark mass at high n B .In practice, this inconsistency does not pose significant issues for astrophysical applications, as the EOS is not applied in regimes of extremely high density.However, in a more consistent approach it is crucial to guarantee that the effects of any interaction vanish at asymptotically high densities.Thus, this article focused on a mass formula incorporating an excluded volume that decreases with n B and vanishes as n B → ∞.
Our results show that the Gibbs free energy per baryon at zero pressure remains unaffected by the κ parameter associated with the excluded volume per quasiparticle.Consequently, the stability window shown in Figure 2 is the same across different κ values.Within this stability window, we focused on the light green shaded area, associated with parameter values for a and C where three-flavor quark matter is more stable than the most tightly bound atomic nucleus, while two-flavor quark matter falls short of this stability (commonly referred to as the SQM hypothesis).The EOS models represented by these parameters support the possibility of strange quark stars.
From the SQM region, we selected three parameter sets for our calculations, indicated by red dots in Figure 2. The pressure curves in Figure 3 exhibit a strong dependence on the parameter a, which is manifested in variations of concavity.Specifically, our numerical findings indicate that the curves are convex upward for a ≲ 2, linear for a ≈ 2, and concave downward for a ≳ 2. Additionally, an increase in the κ parameter leads to an upward shift of the entire curve, i.e., at any given energy density, a higher κ value corresponds to increased pressure.Such a behavior indicates a rise in the stiffness of the EOS for larger κ values.Correspondingly, the shape of the pressure-energy density curves directly influences the behavior of the speed of sound as a function of baryon number density, as apparent from Figure 4.For a ≲ 2, c s increases with n B ; for a ≈ 2, the curves are almost horizontal; and for a ≳ 2, c s decreases with n B .In all cases, we observe that the curves approach the conformal limit at asymptotically large densities.The speed of sound curves for a = 1 exhibit the same qualitative behavior as those presented in [20].It is noteworthy that the constant speed of sound model, widely used in the literature (e.g., [52][53][54]), corresponds to a = 2, which is associated with a quadratic interaction among quarks.
We examined the mass-radius relationship of strange quark stars.A key feature of the curves is the significant increase in maximum gravitational mass as the effect of excluded volume is enhanced.Also, for fixed values of the parameter a, the excluded volume effect substantially shifts the curves towards larger radii.The curve without excluded volume effects is significantly below the observational limit of 2 M ⊙ and does not meet the observational constraints imposed by recent NICER [30,31,44,45] and LIGO/Virgo [55] observations.However, by incorporating a significant fraction of excluded volume into the model, the EOS becomes stiff enough to attain the 2 M ⊙ limit.Additionally, with the general increase in the radii of all objects, the other observational limits are also fulfilled.Another interesting aspect is that models incorporating the excluded volume effect naturally satisfy observations of the recently observed object by HESS [46].This object is uniquely characterized by its very small mass and radius, which challenges theoretical models of objects containing hadronic matter [56].As shown in Figure 5, models of strange quark matter agree more naturally with the observation of HESSJ1731.Another noteworthy observation is GW190814, which was detected by LIGO/Virgo [55].This event involved the merger of a 23 M ⊙ black hole with a compact object of 2.6 M ⊙ .Should the mass of the latter be confirmed, it could either be the smallest known black hole or the most massive known neutron star.An alternative interpretation arises from our curves, suggesting that this object could be a strange quark star.We have also examined the tidal deformability Λ of strange quark stars in relation to the GW170817 event (Figure 6).The model with a = 1 more closely matches the GW170817 constraints compared to the a = 3 model, where compatibility is limited to κ = 0.5.This is because decreasing a shifts the curves toward larger radii, reducing compactness and increasing tidal deformability, as illustrated in Figure 5. Finally, we analyzed the relationship between the baryon and gravitational masses of strange quark stars (Figure 7) and compared these findings with those derived from the commonly used SLY4 hadronic EOS.The consistent positioning of the SLY4 curve above that of the strange matter curves, suggests that the phase conversion of a metastable isolated hadronic star into a strange quark star is exothermic, in qualitative agreement with findings from previous works [49][50][51]57].
To conclude, it is worth comparing our results with those obtained in previous studies that used the QMDDM but employed different parameterizations of the quasiparticle mass formula.Specifically, prior studies introduced an additional term that increases with density (see e.g.Refs.[20,[36][37][38] and references therein).This inclusion enables the maximum gravitational mass of stellar configurations to fulfill the 2 M ⊙ constraint.The reason for this is explained by the stability window depicted in Figure 4 of Xia et al. [20].In this figure, it becomes apparent that increasing the C constant in Eq. ( 5) requires a shift towards lower values of the D constant of the same equation to remain within the stability window for SQM.By doing this, the effective mass of the quasiparticles is reduced in the density ranges typical of compact stars, thereby enhancing the EOS stiffness.However, this mass reduction effect vanishes at higher densities, as the mass formula of Eq. ( 5) diverges for large n B .By taking a different strategy, the effect of the mass formula used in our study ultimately yields similar outcomes.When we consider the excluded volume effect, the interaction mass m I for quasiparticles takes the specific form m I = Cq a/3 /n a/3 B .Since 0 < q < 1, it follows that the coefficient within the numerator of m I decreases relative to the model without excluded volume effects.Thus, a similar effect to that obtained by [20] is ultimately achieved, but with the advantage that in our case, we do not find any issues with the behavior of the mass formula at asymptotically high densities.

Figure 1 .
Figure 1.Gibbs free energy per baryon G/n B as a function of pressure p for different parameter sets labeled (a), (b), and (c).Within each panel, three values of the excluded volume parameter κ are considered: solid lines represent point-like particles (κ = 0), dashed lines correspond to κ = 0.3, and dotted lines to κ = 0.5.The horizontal line at 930 MeV indicates the energy per baryon in 62 Ni, serving as a benchmark for stability against decay into normal nuclear matter.Notice that within each panel all curves coincide at p = 0 (see discussion in the text for further details on this issue).

Figure 2 .
Figure 2.Parameter space for quark matter in bulk.The shaded light green region denotes the parameter combinations for which SQM is the true ground state of matter.The region above the light green area, labeled as 'hybrid matter', corresponds to a scenario where hadronic matter is energetically favored at low pressures, whereas quark matter becomes preferred at higher pressures.The dark green region below, marked as 'nonexistence of nuclei', represents parameter choices where nuclei would be unstable.The same parameter space is valid for any choice of the excluded volume parameter κ.Note that the stability window presented here has some quantitative differences from that shown in[42].The reason for this discrepancy is that Ref.[42] contained an error in the expression of the chemical potential, as previously mentioned at the end of Section 3.1.

Figure 3 .
Figure 3. Pressure p as a function of the energy density ϵ for the same sets of model parameters as shown in Figure 1.In panel (a), the curves exhibit concave downward behavior; in panel (b), they are straight lines; and in panel (c), they are concave upward.As the value of κ increases, the EOS becomes stiffer.

Figure 4 .
Figure 4. Speed of sound c s , normalized to the speed of light, as a function of the baryon number density normalized to nuclear saturation density n B /n 0 .In panel (a) the curves show a decreasing trend which corresponds to the concave curvature seen in the pressure versus energy density plot in the previous figure.Panel (b) presents nearly horizontal curves, consistent with the near-linear relation between p and ϵ observed previously.Finally, in panel (c) the curves are ascending reflecting the convex curvature of the curves in the previous figure.The gray horizontal line indicates the conformal limit of c s = 1/ √ 3.It is observed that c s approaches this limit at high densities.

Figure 5 .
Figure 5. Mass-radius relationship for strange quark stars calculated under the same parameter sets of previous figures.The solid, dashed, and dotted lines within each panel have the same meaning as in previous figures.The different shaded regions correspond to observational data from various sources, providing constraints on the mass and radius of observed NSs.

Figure 6 .
Figure 6.Variation of tidal deformability (Λ) with the star gravitational mass (M G ) for the same model parameters as in previous figures.The error bar is the tidal deformability of a 1.4 M ⊙ neutron star estimated from GW170817 to be Λ 1.4 = 190 +390 −120 at the 90% level.

Figure 7 .
Figure 7. Relationship between gravitational mass (M G ) and baryon mass (M B ) for the same model parameters as in previous figures.For comparative purposes, we have included the curve representing the SLY4 hadronic EOS.Each panel illustrates how M G increases roughly linearly with M B across different model parameters.