---
title: 'Strange Quark Stars: Self-Bound Compact Stars'
url: https://www.emergentmind.com/topics/strange-quark-stars
type: topic
---

# Strange Quark Stars: Self-Bound Compact Stars

Searching arXiv for recent and foundational papers on strange quark stars to ground the article.
Strange quark stars are hypothetical compact stars composed of absolutely stable three-flavor strange quark matter, i.e. deconfined up, down, and strange quarks, possibly with electrons, and in some models with either no crust, a thin nuclear crust, or a strangelet crust. Their existence is tied to the strange quark matter hypothesis, according to which strange quark matter may be the true ground state of hadronic matter [1210.1910]. In contrast to neutron stars, they are self-bound by the strong interaction, a property that leads to distinctive structural, dynamical, and observational consequences, including the scaling \(M \propto R^3\) for low-mass self-bound configurations, the possible existence of strange dwarfs and strange quark planets, and, for bare stars, ultra-high surface electric fields on the order of \(10^{18}\) to \(10^{19}\,\mathrm{V/cm}\) [1210.1910; 1501.02122]. Contemporary work studies these objects across multiple theoretical settings, including the MIT bag model, Dyson-Schwinger and contact-interaction approaches to dense QCD, modified-gravity frameworks, dark-matter admixture, stellar evolution, and binary-merger simulations [2006.00479; 1603.05755; 2604.20159; 2404.00363].

## 1. Definition, self-binding, and stellar taxonomy

Strange quark stars are proposed to be made of absolutely stable strange quark matter. In the overview of quark-star structure, such objects can be either bare or enveloped in thin nuclear crusts, and their structure differs qualitatively from that of neutron stars because self-binding by the strong force, rather than gravity alone, determines the equilibrium over a wide mass range [1210.1910]. A central consequence is the low-mass scaling \(M \propto R^3\), which also underlies the prediction of a continuous hydrostatically stable sequence extending from \(1\)–\(2\,M_\odot\) strange stars down to strange dwarfs and strange planets [1210.1910; 1501.02122].

Bare quark stars are described as having deconfined quark matter up to the surface, with the surface itself extremely thin, of order \(1\) fm, and an electron layer several hundred fermis thick above it [1210.1910]. Dressed, or crusted, quark stars instead possess a thin nuclear crust made of heavy ions in an electron gas, held out of contact with the quark core by an intense electric dipole layer. The maximum base density of such a crust is set by the neutron drip point, approximately \(4.3\times10^{11}\,\mathrm{g/cm^3}\), because free neutrons would otherwise be absorbed into the quark matter core [1210.1910].

A distinctive structural feature is that quark stars with crusts form two-parameter sequences. Their structure depends on the central star density and the density at the base of the crust, whereas neutron stars are primarily described by one-parameter sequences governed by central density [1210.1910]. The same self-bound character also implies that there is no minimum mass for strange stars, in contrast to neutron stars, and that compact SQM objects can extend into the planetary regime [1210.1910; 1501.02122].

For bare stars, the electron cloud produces ultra-high electric fields \(E \sim 10^{18}-10^{19}\,\mathrm{V/cm}\), with the associated energy density \(\mathcal{E}=E^2/(8\pi)\) comparable to the quark matter energy density; the review of quark-star structure notes that this can affect the stellar mass and radius by up to \(15\%\) and \(5\%\), respectively [1210.1910]. Bare stars are also not subject to the Eddington luminosity limit, which has been invoked in discussions of strong pair-creation-driven surface emission [1210.1910]. These properties collectively motivate treating strange quark stars not merely as neutron-star analogues, but as a broader class of self-bound compact objects with a richer taxonomy.

## 2. Equations of state and microscopic descriptions of strange quark matter

A large fraction of the literature models strange quark stars using the MIT bag model or closely related variants. In the 4D Einstein-Gauss-Bonnet study, two equations of state are used: a massless-quark approximation with
\[
\epsilon = 3P + 4B
\]
and a cold-star approximation,
\[
\epsilon = 3.05\,P + 368,
\]
quoted for \(B = 70\, \mathrm{MeV}/\mathrm{fm}^3\) [2006.00479]. In other MIT-based treatments, quark matter is written as
\[
P_s = \frac{1}{3}(\epsilon_s-4B),
\]
or equivalently
\[
p_r = \frac{1}{3}(\rho-4B),
\]
with the bag constant fixed according to the model under study [1801.05031; 2509.10583; 2105.13899]. The short review on recent progress presents the MIT form schematically as
\[
P = k(\epsilon c^2 - 4B),
\]
with \(k=1/3\) for massless quarks [2404.00363].

More microscopic approaches attempt to encode nonperturbative QCD effects directly. In the Dyson-Schwinger quark model, the strange quark matter pressure is written as
\[
P(\mu_u,\mu_d,\mu_s) = - DS + \sum_{q=u,d,s}\int_{\mu_q^0}^{\mu_q} d\mu \, n_q(\mu),
\]
with model constraints imposed by the stability of ordinary nuclear matter and by the requirement that strange quark matter be more bound than iron if the strange-matter hypothesis is to hold [1603.05755]. Respecting these constraints, the maximum mass is about \(1.9\) solar masses and typical radii are \(9\)–\(11\) km [1603.05755]. In a later Poincaré-covariant study using a symmetry-preserving vector\(\otimes\)vector contact interaction, the equation of state is constructed from a momentum-independent quark propagator at finite chemical potential, with
\[
P = -B + \sum_{i=u,d,s,e} P_i(\mu_i), \qquad
\epsilon = -P + \sum_{i=u,d,s,e}\mu_i n_i,
\]
subject to beta equilibrium and charge neutrality [2604.20159]. That analysis finds that reducing the coupling constant stiffens the EOS, whereas increasing the ultraviolet cutoff softens it, and identifies parameter sets \(\alpha_{ir}=0.735\pi,\, \Lambda_{uv}=0.905\,\mathrm{GeV}\) and \(\alpha_{ir}=0.588\pi,\, \Lambda_{uv}=0.9955\,\mathrm{GeV}\), together with \(B \approx (0.106\,\mathrm{GeV})^4\), as compatible with multi-messenger constraints [2604.20159].

Other effective descriptions include the quasi-particle model and the extended MIT bag model. In the dark-matter-admixed study, the quasi-particle model uses effective masses and a running coupling, with pressure
\[
P(\mu) = -B + \sum_{i=u,d,s,e} P_i(\mu_i),
\]
while the extended MIT bag model incorporates medium effects through density-dependent effective masses and writes
\[
P = -\sum_i (\Omega_i + B_i(\mu_i)) - B, \qquad
\epsilon = \sum_i \mu_i \rho_i - P
\]
[2409.09275]. The review on recent progress further lists the NJL model and quasi-particle model as two of the three popular phenomenological descriptions besides MIT, emphasizing that the resulting mass-radius relations depend sensitively on model assumptions [2404.00363].

A recent scaling analysis abstracts away from specific microphysics and links macroscopic observables to the energy per baryon at zero pressure \(e_{uds}^0\) and repulsive interactions. In that work, the maximum mass obeys
\[
M_{\mathrm{max}} = \frac{1}{\sqrt{q^*}} \left[5.1098 - 3.5154 \left( \frac{e_{uds}^0}{930\,\mathrm{MeV}} \right) \right]\, M_{\odot},
\]
with \(q^*=q\) in the QMDD model and \(1/\sqrt{q^*}=1+1.3\,G_V^{0.8}\) in the vector MIT model [2503.11515]. This suggests that self-binding and repulsion can be parameterized in a way that exposes robust macroscopic trends across substantially different microscopic models.

## 3. Structure equations, stability, and compactness

Given an equation of state, stellar structure is obtained by solving the Tolman-Oppenheimer-Volkoff equations or modified variants. In the review of recent progress, the standard GR equations are written as
\[
\frac{dp(r)}{dr} = - \frac{G \left[ \epsilon + p \right]\left[ m(r) + 4\pi r^3 p \right]}{r(r-2G m(r))}, \qquad
\frac{dm(r)}{dr} = 4\pi r^2 \epsilon
\]
[2404.00363]. These equations also appear in the Dyson-Schwinger study as the basis for deriving mass-radius curves [1603.05755].

The self-bound character of strange quark matter changes both boundary conditions and stability diagnostics. In the 4D Einstein-Gauss-Bonnet treatment, the pressure gradient and mass equation are modified by the Gauss-Bonnet coupling \(\alpha\),
\[
\frac{dP}{dr} = -\frac{G \epsilon(r) m(r)}{c^2 r^2}
\frac{
\left[1 + \frac{P(r)}{\epsilon(r)}\right]\left[1 + \frac{4\pi r^3 P(r)}{c^2 m(r)} - \frac{2G\alpha m(r)}{c^2 r^3}\right]
}{
\left[1 + \frac{4G\alpha m(r)}{c^2 r^3} \right]\left[1- \frac{2G m(r)}{c^2 r}\right]
},
\]
\[
\frac{dm(r)}{dr} = \frac{6\alpha G m(r)^2 + 4\pi r^6 \epsilon(r)}{4\alpha G r m(r) + c^2 r^4},
\]
with the GR limit recovered as \(\alpha \to 0\) [2006.00479]. Positive \(\alpha\) increases the maximum mass and radius, while negative \(\alpha\) reduces them [2006.00479]. For the massless quark EOS with \(B=70\,\mathrm{MeV}/\mathrm{fm}^3\), the quoted maximum mass changes from \(1.82\,M_\odot\) and \(R=9.93\,\mathrm{km}\) in GR to \(2.04\,M_\odot\) and \(R=10.12\,\mathrm{km}\) for \(\alpha=5\,\mathrm{km}^2\) [2006.00479].

The same study defines compactness as
\[
C = r_g/R = \frac{2GM}{c^2R},
\]
and reports values in the range \(0.5 < C < 0.6\), with compactness slightly higher for negative \(\alpha\) and lower for positive \(\alpha\) at fixed mass [2006.00479]. Stability is tested by the condition \(dM/d\epsilon_c>0\) and by the adiabatic index
\[
\gamma = \left( 1 + \frac{\epsilon}{P} \right)\left(\frac{dP}{d\epsilon}\right),
\]
with all studied configurations satisfying \(\gamma>\gamma_{\rm crit}\) [2006.00479].

In mimetic gravity, the equilibrium equation acquires an additional force \(F_E\),
\[
-\frac{dp_r}{dr} - \frac{j'(r)}{2} (\rho + p_r) + \frac{2}{r}(p_t - p_r) + F_E = 0,
\]
and equilibrium is reported not to be achieved unless the mimetic extra force is included [2509.10583]. The same paper checks NEC, WEC, DEC, and SEC, quotes \(\Gamma > 4/3\) as the adiabatic stability criterion, and finds \(0 \leq V_s^2 \leq 1\) together with surface redshift below \(2\) in all configurations studied [2509.10583].

Universal and scaling relations specific to strange quark stars have also been proposed. The 2025 scaling study gives an EOS-independent relation for the normalized moment of inertia,
\[
\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,
\]
with \(b_0 = 172.89\), \(b_1=-2306.9\), \(b_2=12816\), \(b_3=-33089\), \(b_4=32699\), as well as
\[
\log_{10}\Lambda = a_3 \mathcal{C}_G^3 + a_2 \mathcal{C}_G^2 + a_1 \mathcal{C}_G + a_0
\]
with \(a_0=5.9886\), \(a_1=-38.861\), \(a_2=123.93\), \(a_3=-184.19\), and
\[
\mathcal{C}_G = d_2\mathcal{C}_B^2 + d_1\mathcal{C}_B
\]
with \(d_2=-0.682\), \(d_1=0.954\) [2503.11515]. That work emphasizes that these relations differ significantly from those previously established for hadronic stars.

## 4. Oscillations, tidal deformability, and response to perturbations

Tidal deformability is a central observable in strange-star phenomenology. In the recent review, it is defined through
\[
Q_{ij} = -\lambda E_{ij}, \qquad
\Lambda = \frac{\lambda}{m^5}, \qquad
k_2 = \frac{3}{2}\lambda R^{-5},
\]
with \(k_2\) obtained from the internal perturbation problem coupled to the background stellar structure [2404.00363]. The Poincaré-covariant contact-interaction study uses
\[
\Lambda = \frac{2}{3}k_2 \left(\frac{R}{M}\right)^5
\]
and explicitly accounts for nonzero surface energy density in bare self-bound stars [2604.20159]. For \(1.4\,M_\odot\) stars it quotes \(\Lambda_{1.4\,M_\odot}\approx699\) for \(\alpha_{ir}=0.735\pi\), \(B=(0.106\,\mathrm{GeV})^4\), and \(\Lambda_{1.4\,M_\odot}\approx598\) for \(\alpha_{ir}=0.588\pi\), \(\Lambda_{uv}=0.9955\,\mathrm{GeV}\), \(B=(0.106\,\mathrm{GeV})^4\), describing both as fully compatible with GW170817 constraints [2604.20159].

Low-mass strange objects exhibit especially distinctive tidal properties. For bare strange quark planets and dwarfs, one study finds \(k_2 \approx 0.75\) and the analytic scalings
\[
\lambda \propto m^{5/3}, \qquad \Lambda \propto m^{-10/3},
\]
adding that bare strange quark planets are effectively rigid to tides and have extremely small tidal deformabilities [2105.13899]. The same work states that, for a typical \(0.6\,M_\odot\) compact star, the tidal deformability of a strange dwarf is about \(1.4\) times less than that of a normal white dwarf, and that the distinction between strange planets and normal planets is even more pronounced [2105.13899].

Oscillation spectra provide a complementary diagnostic. In the Newtonian Cowling treatment of strange stars with a strangelet crust, the \(l=2\) spheroidal \(f\)-mode and \(p\)-mode frequencies are calculated for a two-component core-plus-crust model [1702.05691]. For homogeneous strange quark stars, \(f\)-modes lie in the \(4\)–\(5\) kHz range and \(p_1\)-modes in the \(25\)–\(35\) kHz range, higher than for typical neutron stars, while the addition of a realistic strangelet crust changes the spheroidal frequencies only very slightly, by a few Hz for the \(f\)-mode and up to about \(2\%\) for the \(p\)-modes [1702.05691]. The paper stresses that Newtonian gravity overestimates the frequencies compared with relativistic calculations [1702.05691].

Fully general-relativistic simulations have extended the analysis from isolated stars to mergers. A novel smooth-crust prescription was introduced to regularize the enthalpy discontinuity at the surface, using a very thin crust with mass \(\sim5\times10^{-3}M_\odot\), thickness \(\sim240\) m, and only \(\sim1\%\) effect on tidal deformability [2102.07721]. The oscillation frequencies of isolated simulated stars agree with perturbative predictions, with the \(F_0\) mode within \(1\)–\(3\%\) of linear perturbation theory [2102.07721]. In equal-mass \(1.35\,M_\odot\) binary mergers, the inspiral and post-merger frequencies follow the same quasi-universal relations derived from hadronic stars when expressed in terms of the tidal deformability \(\Lambda\), including
\[
f_2 \approx 5.832 - 0.800\,\Lambda^{1/5}\,\text{kHz},
\]
but not when expressed in terms of the average compactness [2102.07721]. This is a direct statement that universal relations framed in \(\Lambda\) can mask the self-bound nature of strange stars.

Radial oscillations have also been proposed as a dissipation channel specific to strange quark matter objects. One recent study assumes a critical surface density
\[
\rho_0 = 4.7\times10^{14}\,\mathrm{g/cm^3},
\]
below which strange quark matter becomes unstable and decays into photons, hadrons, and leptons [2507.10617]. For small oscillation amplitudes \(x=\Delta R/R\), it gives the total energy to be radiated for small constant-density objects as
\[
T = \frac{9}{2}\beta_s M c^2 x^2 \approx 7.4\times10^{48}\,\mathrm{erg}\,\frac{M}{M_\odot}x^2,
\]
and quotes \(T\approx 6.6\times10^{45}\,\mathrm{erg}\) for a \(1.4\,M_\odot\) star with \(x=0.001\), with energy emitted on timescales of order \(1\) ms [2507.10617]. The same work argues that larger amplitudes may lead to fragmentation or dissolution of the surface layers [2507.10617].

## 5. Formation, evolution, and astrophysical environments

Several channels for strange-star formation have been proposed. In low-mass X-ray binaries, one study assumes that quark deconfinement occurs when the neutron-star core density reaches \(\sim5\rho_0\), where \(\rho_0\sim2.7\times10^{14}\,\mathrm{g\,cm^{-3}}\), corresponding to a critical density of \(\sim1.35\times10^{15}\,\mathrm{g\,cm^{-3}}\) [1303.2458]. With a standard equation of state, the corresponding critical gravitational mass is \(\sim1.8\,M_\odot\), implying that a typical \(1.4\,M_\odot\) neutron star must accrete \(\sim0.4\,M_\odot\) before conversion [1303.2458]. Population synthesis in that work gives a conversion fraction of about \(1\permil\)–\(10\%\) of LMXBs, depending strongly on accretion efficiency and initial neutron-star mass, and an estimated Galactic birthrate of isolated strange stars or submillisecond pulsars of \(\sim5\)–\(70\) per Myr if the conversion disrupts the binary [1303.2458].

The two-families scenario embeds strange-star formation in binary evolution more broadly. There, below a critical gravitational mass \(M_\mathrm{max}^H-\Delta M\) only hadronic stars exist; in the range \((M_\mathrm{max}^H-\Delta M)\leq M\leq M_\mathrm{max}^H\) hadronic stars and strange quark stars coexist; and above \(M_\mathrm{max}^H\) all compact objects are strange quark stars, with \(M_\mathrm{max}^H\sim1.5\)–\(1.6\,M_\odot\) and \(\Delta M\sim0.1\)–\(0.15\,M_\odot\) [1707.01586]. Large-scale population synthesis identifies accretion from a secondary onto a neutron star as the main channel for strange-star formation in binaries, accounting for \(72\)–\(96\%\) of cases, while strange quark stars make up \(1\)–\(4\%\) of all compact objects in binaries and \(3\)–\(18\%\) in LMXBs [1707.01586]. Double strange-quark-star systems are rare, with mergers at a rate of \(\sim12\) Gyr\(^{-1}\) in the Milky Way [1707.01586].

A more recent binary-supernova scenario considers a compact system containing a neutron-star companion and an evolved carbon-oxygen or Wolf-Rayet star. Three-dimensional SPH simulations show that fallback onto the newborn neutron star and hypercritical accretion onto the companion can raise the central density enough to trigger deconfinement [2507.22033]. That work states that the accretion rates reach \(10^{-3}\)–\(10^{-1}\,M_\odot\,\mathrm{s}^{-1}\), and uses the strangeness threshold \(Y_S=n_s/n_B\gtrsim0.2\)–\(0.3\) as the deconfinement criterion [2507.22033]. It further argues that the conversion can release \(\sim10^{52}\)–\(10^{53}\) erg and may leave NS–SQS or SQS–SQS binaries [2507.22033].

Thermal evolution from proto-strange stars to cold stable strange stars has also been modeled. A self-consistent thermodynamic treatment based on the baryon density-dependent quark mass model follows the star through four stages: neutrino-trapped birth, early deleptonization, maximum heating/neutrino transparency, and final cooling to \(T=0\) after about \(100\) years [2601.11106]. The maximum mass and radius decrease monotonically along this sequence, from \(2.21\,M_\odot\), \(14.13\,\mathrm{km}\) in the first stage to \(2.07\,M_\odot\), \(13.22\,\mathrm{km}\) for the cold stable strange star [2601.11106]. The same study reports core temperatures below \(15\) MeV at birth, below \(23\) MeV in the second stage, below \(34\) MeV during maximum heating, and \(0\) MeV in the final state, while noting that cold strange-star configurations are consistent with HESS J1731-347, PSR J1231-1411, PSR J0030+0451, PSR J0348+0432, and PSR J0740+6620 [2601.11106].

## 6. Observational probes, multimessenger signatures, and current constraints

The chief observational difficulty is that strange stars and neutron stars can have similar masses and compactness. Several papers explicitly state that distinguishing them is extremely difficult using traditional mass-radius measurements, spin, or cooling alone [2102.07721; 2109.15161]. The general-relativistic merger simulations sharpen this point: quark-star and hadronic-star binaries with comparable tidal deformabilities produce very similar inspiral gravitational-wave signals, even when their radii differ by over \(1\) km, and become distinguishable only if independent radius or compactness information is available [2102.07721].

One proposed avenue is the search for strange quark planets. Because SQM objects are self-bound and extraordinarily dense, the tidal disruption radius is much smaller than for normal planets. In the strange-planet merger study, the disruption radius is written as
\[
r_{\rm td} \approx 5.1 \times 10^{10}
\left(\frac{M}{1.4 M_\odot}\right)^{1/3}
\left(\frac{\rho_0}{10\, {\rm g\,cm}^{-3}}\right)^{-1/3} \text{ cm},
\]
which for a strange planet with \(\rho_0=4\times10^{14}\,\mathrm{g\,cm^{-3}}\) becomes \(r_{\rm td}\approx1.5\times10^{6}\) cm [1501.02122]. The same work gives the GW power
\[
P = \frac{32 G^4 M^2 m^2 (M + m)}{5c^5 a^5},
\]
and the strain amplitude for a circularized binary,
\[
h = 5.1 \times 10^{-23}
\left(\frac{\mathcal{M}}{1~M_\odot}\right)^{5/3}
\left(\frac{P_{\rm orb}}{1~{\rm hr}}\right)^{-2/3}
\left(\frac{d}{10~{\rm kpc}}\right)^{-1},
\]
reporting \(h\approx1.7\times10^{-22}\) at \(10\) kpc for a strange planet of mass \(10^{-4}M_\odot\) near disruption [1501.02122]. A related observational study proposes that pulsar planets with orbital periods \(P_{\rm orb}<6100\) s are strong SQM candidates because normal planets could not survive so close to the host [2109.15161].

Electromagnetic manifestations have likewise been discussed. Bare stars can exceed the Eddington limit and may have pair-creation-driven surface emission [1210.1910]. Accretion-induced or spontaneous conversion of a neutron star into a strange star has been associated with large energy releases: the Dyson-Schwinger study gives
\[
E_r = [M_{\mathrm{NS}}(M_B)-M_{\mathrm{SQS}}(M_B)]c^2
\]
and quotes \(E_r \simeq (0.8\,–\,2.5)\times10^{53}\,\text{erg}\) for a typical \(1.4\,M_\odot\) neutron star and \(E_r^{\max}=3.6\times10^{53}\,\text{erg}\) for the heaviest stars considered [1603.05755]. The two-families population-synthesis work links such transitions to GRB-like phenomena not associated with supernovae or double-neutron-star mergers [1707.01586]. The recent review also notes proposals connecting strange-star crust collapse and binary-merger or conversion events with short gamma-ray bursts and fast radio bursts [2404.00363].

Current observational constraints do not yet compel a unique interpretation. The Poincaré-covariant contact-interaction EOSs are reported to match pulsar mass measurements, GW170817, and NICER-related constraints for suitable parameter choices [2604.20159]. The evolutionary study argues that cold strange-star sequences can reproduce observed masses and radii of several compact objects that are difficult to explain in standard neutron-star models [2601.11106]. A Bayesian ranking analysis based on NICER mass-radius measurements goes further, stating that for PSR J0614-3329 the measured \(M = 1.44^{+0.06}_{-0.07} M_{\odot}\) and \(R_{eq} = 10.29^{+1.01}_{-0.86}\) km provide a strong case for strange quark stars over physically motivated neutron-star models compatible with such a low radius [2508.02652]. That study reports that many nucleonic EOSs are decisively disfavored relative to a broad class of strange-star EOSs, although some softer or crossover models remain comparatively competitive [2508.02652].

Controversies remain. The 4D Einstein-Gauss-Bonnet paper itself notes that the claim that the resulting 4D theory is of pure graviton has been cast in doubt on several grounds [2006.00479]. In merger phenomenology, the fate of ejected strange quark matter remains unsettled: fully relativistic simulations find that dynamical ejecta are \(\sim20\%\) less for strange-star binaries than for hadronic ones, with much smaller high-velocity and high-entropy tails, but also conclude that predicting kilonova or nucleosynthetic signatures requires an accurate model for evaporation of ejected quarks into nucleons [2102.07721]. This suggests that the principal challenge is no longer just constructing viable strange-star models, but identifying observables that are simultaneously robust against EOS degeneracy and sensitive to self-bound quark matter.

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