---
title: 'Strange Quark Stars: Hypotheses, EOS, and Phenomenology'
url: https://www.emergentmind.com/topics/strange-quark-stars-sqss
type: topic
---

# Strange Quark Stars: Hypotheses, EOS, and Phenomenology

Strange quark stars (SQSs) are hypothetical compact stars entirely composed of deconfined quark matter, usually taken to be three-flavor strange quark matter (SQM) containing \(u\), \(d\), and \(s\) quarks, together with electrons when needed for charge neutrality. In the strange-matter hypothesis, bulk three-flavor SQM can have an energy per baryon lower than that of the most stable atomic nucleus, so the star is self-bound by the strong interaction rather than only by gravity. SQSs have been theoretically proposed since the 1970s, can have masses and compactness similar to neutron stars (NSs), and remain observationally unconfirmed despite decades of multiwavelength data [2507.22033][1210.1910][2404.00363].

## 1. Strange-matter hypothesis and defining properties

The defining microscopic assumption is absolute stability of bulk three-flavor SQM. A commonly used criterion is that the zero-pressure energy per baryon,
\[
e^0_{uds}\equiv \frac{\epsilon}{n_B}\bigg|_{p=0},
\]
satisfies \(e^0_{uds}<930\ \mathrm{MeV}\), while two-flavor matter remains unstable; this is the standard way to realize the Bodmer–Witten picture without destabilizing ordinary nuclei [2503.11515][1210.1910]. In this regime, the pressure vanishes at a finite, nonzero density rather than at zero density, so SQSs are self-bound objects with a sharp surface [2503.11515][2102.07721].

Self-boundness leads to several structural consequences repeatedly emphasized across the literature. Bare quark stars exhibit the approximate scaling \(M\propto R^3\), may have no strict minimum mass in the slowly rotating case, and typically occupy radii of order \(R\sim 10\)–\(12\) km with maximum masses around \(\sim 2\,M_\odot\) in representative models [1210.1910]. This differs qualitatively from hadronic NS sequences, which are gravitationally bound and organized essentially by a single parameter, the central density [1210.1910].

The same self-bound character underlies modern scaling analyses. In both the quark-mass density-dependent model with excluded-volume corrections and the vector MIT bag model, the maximum mass is controlled by the zero-pressure binding energy and by repulsive interactions. One explicit fit is
\[
M_{\mathrm{max}} = \frac{1}{\sqrt{q}\, \Bigl[5.1098 - 3.5154 \Bigl(\frac{e^0_{uds}}{930\,\mathrm{MeV}}\Bigr)\Bigr]}\, M_\odot,
\]
accurate to better than \(0.2\%\) in the quoted model class; a related effective form is used in the vector MIT bag realization [2503.11515]. This suggests that observational upper limits on compact-star masses can directly constrain the depth of self-binding and the strength of repulsive quark interactions.

## 2. Equation-of-state frameworks

Because first-principles QCD at compact-star densities remains incomplete, SQS research is organized around phenomenological and QCD-inspired equations of state (EOSs). The literature represented here uses several distinct frameworks, all aiming to realize self-bound three-flavor matter while preserving thermodynamic consistency and stability constraints.

| Framework | Key ingredients | Notable consequence |
|---|---|---|
| MIT bag model | bag pressure \(B\), sometimes density dependent | simple self-bound EOS; constant-\(B\) models are often too restrictive for the heaviest compact objects |
| NJL / MNJL | chiral dynamics; Fierz-weighted interaction via \(\alpha\) | increasing \(\alpha\) stiffens the EOS and raises \(M_{TOV}\) and \(\Lambda\) |
| QMDD / DDQM / CIDDM | density-dependent quark masses; sometimes excluded volume or isospin dependence | finite-density self-binding, simple scaling laws, and strong sensitivity to repulsion parameters |
| Dyson–Schwinger model | dressed propagators and quark–gluon vertex truncations | stable SQSs exist only in a narrow allowed parameter region |
| Perturbative QCD + bag term | running \(\alpha_s\), running \(m_s\), effective bag constant | density-dependent \(B\) can satisfy both \(M_{TOV}>2\,M_\odot\) and tidal constraints |

The MIT bag model remains the canonical baseline. In its simplest form the EOS can be written linearly, but recent work increasingly replaces a fixed bag constant by a density-dependent one. A representative prescription is
\[
{\mathcal{B}_{\rm bag}(\rho)}=\mathcal{B}_{\infty}+(\mathcal{B}_{0}-\mathcal{B}_{\infty})e^{-\alpha(\rho/\rho_0)^2},
\]
with \(\rho_0 = 0.17\,{\rm fm^{-3}}\), \(\alpha = 0.17\), \(\mathcal{B}_{0} = 400\,{\rm MeV/fm^3}\), and \(\mathcal{B}_{\infty} = 8.99\,{\rm MeV/fm^3}\) in one rotating-star study [2601.09532]. In perturbative-QCD-based models, a different density-dependent ansatz,
\[
B = B_0 \exp\!\left[-a\left(\frac{n_B}{n_0}-1\right)^2\right],
\]
is used to emulate nonperturbative vacuum physics, and the constant-\(B\) limit is recovered at \(a=0\) [2411.12048].

NJL-like models incorporate chiral dynamics more explicitly. In the modified NJL construction,
\[
\mathcal{L}=(1-\alpha)\mathcal{L}_{NJL}+\alpha \mathcal{L}_F,
\]
increasing \(\alpha\) makes the EOS stiffer, increases \(M_{TOV}\), and also increases the tidal deformability \(\Lambda\); the SQM stability requirement restricts \(\alpha\lesssim 0.94\) in the quoted parameter set [2402.14040]. In QCD-motivated APT and BPT descriptions, the infrared running coupling remains finite, producing maximum rotating SQS masses of \(3.02\)–\(3.94\,M_\odot\) depending on the onset density and coupling prescription [2104.00544].

Density-dependent-mass models provide another major branch. The CIDDM model uses
\[
m_q = m_{q_0} + \frac{D}{n_B^z} - \tau_q\,\delta\, D_I\, n_B^\alpha e^{-\beta n_B},
\]
which enforces confinement as \(n_B\to 0\) and asymptotic freedom as \(n_B\to\infty\) for \(\beta>0\) [1406.5610]. The older quark mass-density-dependent model uses
\[
m_q = m_{q0} + \frac{C}{n_b^x},
\]
and finds maximum masses in the range \(1.5\text{--}1.8\,M_\odot\) across \(x=1/10,\ 1/5,\ 1/3,\ 1,\ 2,\ 3\), with radii decreasing as \(x\) increases [1007.4737]. More recent DDQM treatments add full thermodynamic self-consistency at finite temperature and lepton fraction, enabling proto-SQS evolution calculations [2311.12511][2601.11106].

## 3. Surfaces, crusts, and equilibrium structure

A central distinction from NSs is the surface structure. A bare quark star has quark matter exposed at the surface, and the quark surface itself is very thin, of order \(\sim 1~\mathrm{fm}\). Electrons extend above the positively charged quark surface in a layer of thickness of several hundred fermis, producing a dipole layer and ultra-high electric fields of order \(10^{18}\ \mathrm{V/cm}\), or \(E\gtrsim 10^{19}\ \mathrm{V/cm}\) if the matter is in a color-superconducting state [1210.1910]. The field energy density can be comparable to the energy density of strange matter itself, changing the mass by about \(15\%\) and the radius by about \(5\%\) [1210.1910].

SQSs may also be crusted rather than bare. In a dressed star, the nuclear crust is electrostatically suspended above the quark surface and is not in direct contact with the quark core. The crust consists of heavy ions immersed in an electron gas, and the density at its base cannot exceed neutron drip,
\[
\rho_{\rm base} \approx 4.3\times 10^{11}\ \mathrm{g/cm^3},
\]
so only the outer crust exists, not the inner crust familiar from NSs [1210.1910]. This gives quark stars with crusts a genuine two-parameter equilibrium sequence, determined by the central density and the crust-base density [1210.1910].

These surface properties feed directly into observables. Because the surface constituents of a bare SQS are not gravitationally bound in the usual way, the standard Eddington luminosity limit does not apply. Photon luminosities can exceed \(10^{38}\ \mathrm{erg/s}\), and pair-plasma luminosities can reach \(\sim 3\times 10^{51}\ \mathrm{erg/s}\) for surface temperatures around \(10^{11}\) K [1210.1910]. By contrast, crusted SQSs are Eddington-limited at the outer nuclear surface [1210.1910].

The same sharp surface complicates numerical relativity. Fully general-relativistic merger simulations therefore introduced a very thin crust via a polytrope,
\[
p = k \rho^\Gamma,\qquad k=8.12,\ \Gamma=1.90,
\]
together with an enthalpy-continuity condition and a baryon-mass rescaling \(\chi=940/850\simeq 1.11\). In the MIT2cfl EOS, the physical zero-pressure baryon mass is \(m_B\simeq 850\ \mathrm{MeV}\), and the crust added for numerical regularization is only about two grid cells wide, with spatial width \(\sim 240\) m, mass \(\sim 5\times 10^{-3}M_\odot\), and tidal deformability changed by only about \(1.3\%\) [2102.07721].

## 4. Formation channels and thermal evolution

A specific late-stage binary-supernova channel for SQS formation has recently been explored in detail. In this scenario, an evolved carbon-oxygen or Wolf-Rayet star collapses in a compact binary with an already existing NS companion. The pre-supernova orbital period is only a few minutes. Core collapse produces a newborn NS and a supernova, while fallback onto the newborn object and ejecta capture by the companion drive hypercritical, highly super-Eddington accretion onto both compact stars [2507.22033].

The accreted material circularizes into disks around both stars within about one orbital period. The compact-star evolution is then followed by imposing
\[
\dot{J}=\tau_{\rm acc},
\]
with
\[
\tau_{\rm acc}= \begin{cases} \chi\,l\,\dot{m}_b, & \text{onto the newborn NS},\\ \chi\,l\,\dot{M}_b, & \text{onto the NS companion}, \end{cases}
\]
where \(\chi\le 1\) is the angular-momentum transfer efficiency and \(l\) is the specific angular momentum at the disk inner radius [2507.22033]. In the two-family scenario adopted there, quark deconfinement is tied to the strangeness fraction rather than an arbitrary mass limit; the qualitative nucleation threshold is
\[
Y_S\equiv n_s/n_B \sim 0.2\text{--}0.3.
\]
Once a seed appears, the hadronic star rapidly converts into an SQS if absolutely stable SQM exists in bulk [2507.22033].

The outcome depends strongly on progenitor mass, explosion energy, and spin. For the hadronic EOSs used in that study, the static maximum masses are only about \(1.58\,M_\odot\), rotation raises the supportable mass to \(\sim 1.85\text{--}1.88\,M_\odot\), and the critical mass for deconfinement is around \(1.52\text{--}1.57\,M_\odot\). A \(15\,M_\odot\) progenitor does not push the companion to \(n_{\rm crit}\), but \(25\,M_\odot\) and \(30\,M_\odot\) progenitors can do so for the weakest explosion energies. For the newborn NS, rapid initial rotation with \(j_{\rm ns,0}\gtrsim 1.2\) tends to lead to the Keplerian mass-shedding limit before deconfinement, whereas slower spins favor compression and conversion [2507.22033].

When conversion occurs in that channel, the star is re-evaluated at the same baryon number and angular momentum with a quark EOS. The resulting SQS has a smaller gravitational mass and a larger equatorial radius than the progenitor hadronic star. The gravitational-mass change is about \(0.02\text{--}0.025\,M_\odot\) for the unpaired bag model and about \(0.22\text{--}0.25\,M_\odot\) for the CFL bag model, corresponding to conversion energies of roughly \(4\times 10^{52}\ \mathrm{erg}\) and \(4\times 10^{53}\ \mathrm{erg}\), respectively [2507.22033]. Bound remnants can emerge as NS-NS, NS-SQS, or SQS-SQS binaries and may merge on timescales of order \(\sim 10\ \text{kyr}\) [2507.22033].

Separate proto-SQS studies track stars assumed to be quark stars from birth. In one DDQM sequence, the star evolves through neutrino-trapped and neutrino-transparent stages into a cold \(T=0\) object, with maximum masses decreasing from \(2.20\,M_\odot\) to \(2.09\,M_\odot\); higher neutrino concentration makes the star slightly more massive because trapped neutrinos suppress strange-quark production and stiffen the EOS [2311.12511]. A closely related self-consistent thermodynamic treatment gives a similar decrease, from \(2.21\,M_\odot\) and \(R=14.13\ \mathrm{km}\) in the earliest neutrino-trapped stage to \(2.07\,M_\odot\) and \(R=13.22\ \mathrm{km}\) in the final cold star, along a fixed baryon-mass evolutionary picture [2601.11106].

## 5. Rotation, magnetic fields, and oscillation spectra

Rapid rotation is particularly natural for self-bound stars. For a typical pulsar-mass strange star (\(\sim 1.45\,M_\odot\)), the Kepler period may be as short as \(P_K \sim 0.55\text{--}0.8\ \mathrm{ms}\), compared with about \(\sim 1\) ms for an NS of the same mass [1210.1910]. More recent fully relativistic rotating models with a density-dependent bag constant find that the static maximum mass is \(2.35\,M_\odot\), while rotating sequences in the range \(1100\)–\(1300\) Hz reach \(M^{\rm max}_{\rm g}\approx 2.55\text{--}2.87\,M_\odot\); in the summary of that study the explored interval is quoted as \(2.55\)–\(2.78\,M_\odot\) [2601.09532].

Rotation also introduces a minimum mass at sufficiently high frequency. While slowly rotating self-bound stars have no minimum mass, the same study finds that above about \(1150\) Hz the sequence develops a minimum-mass configuration, with \(M^{\rm min}_{\rm g}\approx 0.54\text{--}2.13\,M_\odot\) over \(1175\)–\(1300\) Hz. Both the minimum-mass sequence and the mass-shedding sequence obey an approximately linear relation \(f=\tilde a M_g+\tilde b\), with slopes \(\tilde a=0.078\) and \(0.077\ {\rm kHz}/M_\odot\) and \(\tilde b=1.13\ {\rm kHz}\), motivating the statement that the Keplerian frequency depends almost linearly on mass with slope about \(0.08~{\rm kHz}/M_\odot\) [2601.09532].

Strong magnetic fields make the pressure anisotropic. In the CIDDM model,
\[
P_\parallel = \sum_i \mu_i n_i - \mathcal{E}_{\rm tot},\qquad
P_\perp = \sum_i \mu_i n_i - \mathcal{E}_{\rm tot} + B^2 - M B,
\]
so \(P_\parallel\) decreases and \(P_\perp\) increases with field strength [1406.5610]. The orientation of the internal field then becomes decisive: radial orientation reduces the maximum mass, while transverse orientation increases it. For the fast \(B\)-profile, the largest masses are \(M_{\max}\approx 2.13\,M_\odot\) in the transverse case and \(1.57\,M_\odot\) in the radial case; for the slow profile they become \(2.18\,M_\odot\) and \(1.29\,M_\odot\), respectively [1406.5610].

Axisymmetric magnetized rotating equilibria using the density-dependent MIT bag model and fields up to \(5\times10^{17}\,\mathrm{G}\) reach even larger masses when rapid spin is included. In the strongest-field, \(f=1200\,\mathrm{Hz}\) configuration, the maximum mass is \(2.80\,M_\odot\), the circumferential radius is \(14.6\,\mathrm{km}\), and the deformation parameter \(a=R_{\rm eq}/R_{\rm pol}\) reaches \(1.55\). The binding energy per baryon lies between \(171\) and \(184\ \mathrm{MeV}\), and the compactness is quoted as roughly \(\beta\approx 0.25\,M_\odot/{\rm km}\) [2601.09529].

Oscillation spectra provide a complementary diagnostic. For newly born SQSs, the quadrupolar \(g\)-mode eigenfrequencies are about one order of magnitude lower than in newborn NSs, independent of whether the quark matter is described by the MIT bag model or the NJL model. In the MIT calculation, the \(l=2,\ n=1\) \(g\)-mode is \(82.3,\ 78.0,\ 63.1\,\mathrm{Hz}\) at \(t=100,\ 200,\ 300\,\mathrm{ms}\), whereas the corresponding NS values are \(717.6,\ 774.6,\ 780.3\,\mathrm{Hz}\). By contrast, SQS \(f\)- and \(p\)-modes are much higher, with \(f\)-mode frequencies around \(2980\)–\(3016\,\mathrm{Hz}\) and \(p_1\)-mode frequencies around \(16792\)–\(18282\,\mathrm{Hz}\) in the same study [1701.00418].

## 6. Binary dynamics, gravitational waves, and electromagnetic transients

Fully general-relativistic simulations of equal-mass \(1.35+1.35\,M_\odot\) strange-star binaries show both similarities to and differences from hadronic mergers. For the MIT2cfl quark EOS, the stellar radius is \(R\simeq 11.81\) km and the tidal deformability is \(\Lambda\simeq 789\), while the comparison DD2 hadronic stars have \(R\simeq 13.21\) km and \(\Lambda\simeq 858\) [2102.07721]. Despite similar inspiral phasing when expressed through \(\Lambda\), dynamical mass loss is reduced: \(M_{\rm ej}=2.68\times10^{-3}M_\odot\) for quark stars versus \(2.94\times10^{-3}M_\odot\) for hadronic stars, and the quark-star ejecta exhibit much smaller high-velocity and very-high-entropy tails [2102.07721].

The same simulations show that merger and post-merger frequencies obey the same quasi-universal relations derived from hadronic binaries when written in terms of tidal deformability, but not when written in terms of average compactness. This is important because it implies that gravitational-wave frequency information alone may not cleanly distinguish SQSs from NSs unless independent radius information is available [2102.07721]. Low-frequency newborn-star \(g\)-modes add another potential discriminator because they are tied directly to the relativistic composition of quark matter [1701.00418].

Several recent compact-binary observations have been interpreted in this context. Using infrared-finite QCD couplings in APT and BPT, one study infers rotating SQS maximum masses of \(3.02\)–\(3.94\,M_\odot\) and compares them with a GW190425 remnant estimate of \(3.11\)–\(3.54\,M_\odot\), concluding that the remnant might be a strange quark star [2104.00544]. In the two-families picture of the binary-supernova formation channel, mixed NS-SQS systems are explicitly mentioned as relevant to events like GW170817 [2507.22033].

Electromagnetic transients have long been linked to SQS physics. A recent review emphasizes that strong gravitational-wave emission may arise from mergers, mode excitation, continuous-wave instabilities, or even close-in strange quark planets, while fierce electromagnetic bursts may be powered by neutron-star to strange-star conversion, binary strange-star mergers, or crust collapse [2404.00363]. The same review notes that such mechanisms have been discussed in connection with short gamma-ray bursts and fast radio bursts [2404.00363]. In the supernova-binary conversion channel, the quoted conversion energies up to \(4\times10^{53}\ \mathrm{erg}\) further reinforce this possibility, although detailed radiation-hydrodynamic modeling is left open [2507.22033].

## 7. Constraints, extensions, and unresolved issues

No unambiguous SQS identification has yet emerged, and much of current work focuses on discriminants that differ from hadronic-star systematics. One robust line of attack is the existence of SQS-specific universal relations. Across both the QMDD and vector MIT bag models, the moment of inertia, tidal deformability, and compactness obey relations distinct from hadronic stars. Examples include
\[
\frac{I}{M_G^3} = b_4 \mathcal{C}_G^4 + b_3 \mathcal{C}_G^3 + b_2 \mathcal{C}_G^2 + b_1 \mathcal{C}_G + b_0,
\]
and
\[
\log_{10}\Lambda = a_3 \mathcal{C}_G^3 + a_2 \mathcal{C}_G^2 + a_1 \mathcal{C}_G + a_0,
\]
with quoted coefficients that differ from the usual hadronic fits [2503.11515]. This suggests that combined measurements of \(I\), \(\Lambda\), \(M\), and \(R\) could separate self-bound from gravity-bound stars without detailed EOS reconstruction.

Another active issue is how to reconcile heavy compact objects with tidal constraints. In perturbative QCD plus an effective bag term, constant-\(B\) models cannot comfortably exceed about \(2.03\,M_\odot\) while satisfying \(70<\Lambda_{1.4M_\odot}<580\), whereas a density-dependent \(B\) can raise \(M_{TOV}\) to \(2.16\), \(2.21\), \(2.30\), and \(2.43\,M_\odot\) for \(a=0.1,\ 0.2,\ 0.4,\ 0.8\), respectively [2411.12048]. The same framework is used to argue that PSR J0952-0607, PSR J2215+5135, PSR J0740+6620, and even the GW190814 secondary object can be interpreted as SQSs [2411.12048].

More speculative extensions push SQSs into the mass-gap regime. In dRGT-like massive gravity with an MNJL EOS, increasing the bag constant softens the EOS and lowers \(\Lambda\), while modified hydrostatic balance allows self-bound SQSs in the interval \(2.5M_\odot \lesssim M \lesssim 5M_\odot\) that still satisfy \(\Lambda_{1.4M_\odot}\lesssim 580\). The analysis further requires \(R>R_{Sch}\), sub-horizon compactness, and nonzero \(\Lambda_{M_{TOV}}\) to exclude black-hole behavior [2402.14040]. A different extension, involving fermionic dark matter admixed through a vector dark boson, finds that dark matter lowers mass and radius relative to the no-DM case, but can still yield configurations consistent with GW190814 [2209.09021].

SQSs have also been proposed as dark-matter detectors. By scanning 1403 solitary pulsar-like compact stars, one study identifies PSR J1801-0857D as the strongest source of limits on scalar-DM scattering, with age \(9.71\) Gyr and Galactic-center distance \(3.06\) kpc. If that object is assumed to be an SQS rather than an NS, the inferred DM–proton cross-section limits become much weaker and can be comparable to direct-detection bounds; this leads to the stated implication that if scalar dark matter were observed in future terrestrial experiments, old pulsars would be favored to be SQSs rather than NSs [1603.07518].

The central controversy remains unchanged: SQSs are theoretically consistent in many EOS frameworks and now embedded in detailed calculations of formation, evolution, oscillations, magnetized rotation, merger dynamics, and multimessenger phenomenology, yet no single observational signature has isolated them beyond ambiguity. The literature therefore increasingly treats SQS identification as a problem of joint inference across mass, radius, spin, \(\Lambda\), mode spectroscopy, ejecta properties, and transient energetics rather than as a one-observable classification problem [2507.22033][2404.00363].

Source: https://www.emergentmind.com/topics/strange-quark-stars-sqss