Quark Star

One squeeze too far

Start with a neutron star — already the densest object that holds a recognisable surface. A sugar-cube of its core would weigh as much as a mountain range, the whole 1.4-solar-mass object packed into a sphere about 12 kilometres across. Now imagine pushing it harder: dropping in more mass, or simply asking what happens at the very centre where the pressure is highest. The neutrons are crammed so tightly that they begin to overlap. A neutron is not a point; it is itself a bag of three quarks bound by the strong force, roughly a femtometre wide. When neutrons are squeezed until they touch and merge, the natural question becomes: do the quarks still belong to individual neutrons at all?

Quantum chromodynamics — the theory of the strong interaction — answers that above a critical density the quarks deconfine. They stop being prisoners of separate baryons and instead form one continuous Fermi sea filling the whole star. That object is a quark star. In the most striking version of the idea, the matter is not just deconfined but contains a third flavour, the strange quark, in roughly equal numbers to up and down quarks. Such three-flavour strange quark matter may be more stable than ordinary nuclear matter, in which case the star is bound by the strong force itself and would not fall apart even with gravity switched off. It is the only star we can imagine that is, in a literal sense, one gigantic hadron.

The QCD phase transition

The defining physics is deconfinement. At everyday densities, quarks are confined: you cannot isolate a single quark because the strong force between them grows with separation, and the energy poured into pulling them apart instead creates new quark-antiquark pairs. But QCD also predicts asymptotic freedom — at very short distances, or equivalently very high densities and temperatures, the coupling weakens and quarks behave almost as free particles. Squeeze nuclear matter past a few times nuclear saturation density,

ρ₀ ≈ 2.7 × 10¹⁴ g/cm³   (nuclear saturation density)
ρ_transition ≈ 5 – 10 ρ₀  (deconfinement, model-dependent)

and the energetically favoured state flips from a gas of nucleons to a soup of deconfined quarks. The transition is the high-density, low-temperature corner of the QCD phase diagram — the same diagram whose high-temperature corner (the quark-gluon plasma) is probed at RHIC and the LHC by smashing heavy nuclei together. Heavy-ion collisions explore hot, nearly baryon-free matter; quark stars probe the opposite limit, cold and extremely baryon-rich. Lattice QCD, the workhorse for the hot regime, fails at high baryon density because of the notorious "sign problem," so the dense corner relevant to quark stars is the most uncertain region of the entire phase diagram.

Why "strange," and why it might be stable

Naively, adding a third, heavier flavour ought to cost energy: the strange quark has a current mass of about 95 MeV/c², versus a few MeV for up and down. But there is a compensating gain. In two-flavour quark matter the Pauli principle forces up and down quarks into ever-higher momentum states as you add more of them. Converting some down quarks into strange quarks opens a new Fermi sea to fill, lowering the average energy per particle. Arnold Bodmer noted in 1971, and Edward Witten argued forcefully in 1984, that the net effect could push the energy per baryon of three-flavour strange quark matter below that of the most tightly bound nucleus, iron-56:

E/A (iron-56)        ≈ 930.4 MeV   (most bound ordinary nucleus)
E/A (strange matter) < 930.4 MeV   (Bodmer-Witten hypothesis)

If that inequality holds, strange quark matter is the true ground state of baryonic matter, and ordinary nuclei — including everything you are made of — are only metastable, kept from decaying because converting a nucleus into strange matter would require simultaneously creating many strange quarks via the weak interaction, an astronomically slow process. The hypothesis is unproven and remains a genuine "maybe," but it is what makes strange stars conceptually clean: it lets quark matter be self-bound at zero external pressure, so the star is held together by the strong force, with gravity merely confining the dilute outer skin.

The bag model and the equation of state

The simplest workable description is the MIT bag model. It treats the quark matter as a gas of nearly free, relativistic quarks confined inside a region of perturbed QCD vacuum, with a constant "bag pressure" B pushing inward to represent the energy cost of that vacuum. The equation of state is remarkably simple:

P = (ε − 4B) / 3        (massless-quark MIT bag EoS)
B^(1/4) ≈ 145 – 165 MeV  (bag constant, for stable strange matter)

Two features follow immediately. First, the pressure vanishes not at zero energy density but at a finite surface density ε_s = 4B — the matter holds together on its own and ends in a sharp surface, with the density dropping from nuclear values to zero over roughly a femtometre. Second, the sound speed approaches the conformal limit c_s² → c²/3, a stiff but causal equation of state. The window of bag constants that makes strange matter stable (E/A below iron) yet two-flavour matter unstable (so ordinary nuclei survive) is narrow, which is part of why the hypothesis is testable in principle.

A self-bound star built on this equation of state behaves, at low mass, like an incompressible drop of constant density ρ ≈ 4B/c². Its radius then scales with mass like a planet rather than a neutron star:

M = (4/3) π R³ ρ    →    R ∝ M^(1/3)   (low-mass strange-star branch)

This is the opposite of the neutron-star mass-radius relation, where heavier means smaller. Only as the mass climbs toward the general-relativistic maximum does the strange-star curve bend back and turn over.

How you would tell one apart

Because no quark star is confirmed, the interesting question is observational: what would betray one? Several signatures have been proposed.

• A small radius at a given mass. A strange star is more compact. NICER and gravitational-wave radius measurements at the 1-km level can, in principle, separate a soft self-bound branch from ordinary nuclear models — though current data sit in the overlap region.
• The bare quark surface. With no nuclear crust, a strange star's surface can sustain enormous electric fields, of order 10¹⁷–10¹⁸ V/cm, in a thin electron layer a few hundred femtometres thick. This is unlike anything on a neutron star and would alter pair production, polar-cap physics, and the thermal spectrum.
• Fast, clean cooling. Unpaired quark matter allows the direct-Urca-like "quark Urca" process, dumping energy into neutrinos extremely fast, so a young quark star could look anomalously cold for its age. Conversely, the colour-flavour-locked phase gaps every quark and slows cooling dramatically — so the thermal history is a sensitive, if model-dependent, probe.
• Rapid, sustained spin. Self-bound matter resists the r-mode oscillation instability differently from fluid neutron stars, so a strange star might survive at sub-millisecond spin periods that would shake an ordinary neutron star apart.
• Distinctive seismology. A bare surface and a different shear modulus change the spectrum of magnetar quasi-periodic oscillations seen after giant flares.

Quark star versus neutron star, by the numbers

Property	Neutron star	Strange quark star
Constituents	Neutrons, protons, electrons	Deconfined u, d, s quarks + electrons
What binds it	Gravity + nuclear repulsion	Strong force (self-bound)
Typical mass	1.1 – 2.3 M☉	up to ~2 M☉ (model-dependent)
Typical radius (1.4 M☉)	~11 – 13 km	~8 – 11 km
Central density	~5 – 10 ρ₀	~4 – 8 ρ₀
Surface	~1 km crust, smooth density drop	Bare quark surface, drop over ~1 fm
Mass–radius trend	R decreases with M	R ∝ M¹ᐟ³ at low mass
Minimum mass	~0.1 M☉ (set by crust)	essentially none (strange dwarfs allowed)
Surface electric field	negligible	10¹⁷–10¹⁸ V/cm

The overlap in mass and radius is exactly the problem: across most of the observed population, a quark star and a neutron star look almost the same from the outside. The clean discriminators — a sub-8-km radius at 1.4 M☉, a bare-surface spectrum, or a sub-millisecond spin — all live at the edges of what current instruments can measure.

Real numbers and real candidates

Several observed compact objects have, at one time or another, been floated as quark-star candidates:

• RX J1856.5-3754. A nearby (~120 pc) isolated, thermally emitting neutron star. Early single-temperature blackbody fits implied a radiation radius of only 5–6 km, which would be hard for a normal neutron star and prompted strange-star papers around 2002. Two-temperature atmosphere models later revised the radius up to ~14 km (or ~17 km for the canonical mass), defusing the claim — a textbook lesson in how atmosphere modelling, not exotic physics, drove the small radius.
• 3C 58. The compact object in this young supernova remnant (often linked to SN 1181) appeared too cold for standard neutron-star cooling, suggesting enhanced neutrino emission of the kind quark matter provides. The interpretation is not unique — exotic nuclear processes can cool a neutron star similarly fast.
• Massive pulsars. PSR J0740+6620 weighs about 2.08 M☉ and PSR J0348+0432 about 2.01 M☉. A heavy star is sometimes invoked for quark matter, but heavy stars actually constrain exotic matter: any equation of state, quark or nuclear, must be stiff enough to support 2 M☉, which rules out the softest quark models. NICER's radius for J0740+6620 is about 12.4 km.
• GW170817. The 2017 binary neutron-star merger measured the tidal deformability Λ of ~1.4 M☉ stars to be Λ(1.4) ≲ 580–800, favouring relatively compact, soft stars and tightening the allowed equation-of-state band — a constraint every quark-star model must now respect.

To anchor the scale: a strange star of 1.4 M☉ with R ≈ 10 km has a mean density near 7 × 10¹⁴ g/cm³, a surface gravity around 1.9 × 10¹⁴ cm/s² (≈ 2 × 10¹¹ g), and a compactness GM/Rc² ≈ 0.21 — high enough that its own surface, seen from far away, is gravitationally redshifted by roughly 30 percent. A teaspoon (≈ 5 cm³) of its matter would weigh about 3.5 billion tonnes.

Inside: colour superconductivity

The interior is not a featureless quark gas. At the densities and modest temperatures of a settled compact star, quarks near the Fermi surface attract through the strong force and form Cooper pairs, exactly as electrons do in a metallic superconductor. Because quarks carry colour charge, the paired state is a colour superconductor, and several distinct pairing patterns compete:

• Colour-flavour-locked (CFL) phase. The most symmetric and likely densest phase. All three colours and all three flavours pair symmetrically; every quark at the Fermi surface is gapped (gaps of order 10–100 MeV), the matter is electrically neutral with no electrons needed, and it is rigid. CFL matter may be the densest stable matter anywhere in the cosmos.
• 2SC phase. Two-flavour colour superconductivity, where only up and down quarks of two colours pair; arises at lower density where strange quarks are scarcer.
• Crystalline / LOFF phases. When Fermi surfaces are mismatched, pairing can occur in a spatially modulated, crystalline pattern with potentially large shear strength — a candidate for storing the stresses that power pulsar glitches or magnetar flares.

These phases differ in how they emit neutrinos, conduct heat, and respond to rotation, so the thermal and timing history of a real quark star would be a fingerprint of which colour-superconducting condensate fills its core.

Where they might come from

• Conversion of a neutron star. A neutron star that accretes enough mass, or spins down enough to raise its central density past the deconfinement threshold, could undergo a sudden phase transition. The conversion releases ~10⁵²–10⁵³ erg of binding energy — a "quark nova" has been proposed as a second explosion days after a supernova, and as a possible engine for some superluminous supernovae and gamma-ray bursts.
• Direct collapse. A sufficiently massive stellar core could in principle collapse straight to a quark star if the deconfinement transition happens during the bounce.
• Seeded by a strangelet. If a pre-existing strangelet — perhaps a relic from the early universe — is captured by a neutron star, and the Bodmer-Witten hypothesis holds, it can convert the whole star to strange matter from the inside out.
• Mergers. The hot, hyper-dense remnant of a binary neutron-star merger briefly reaches densities well above any isolated star, a natural place to cross into the quark phase; its gravitational-wave ringdown could carry the imprint of a transition.

Common misconceptions and edge cases

• "A quark star is just a very heavy neutron star." No — the distinction is the phase of matter, not the mass. A quark star can be light; the point is that its interior is deconfined and (for a strange star) self-bound, regardless of where it sits on the mass scale.
• "Quark stars are confirmed to exist." They are not. Every candidate to date has an ordinary-neutron-star explanation. The field is healthy precisely because the question is open.
• "All quark stars are strange stars." Not necessarily. A hybrid star has a normal nuclear envelope and crust over a deconfined-quark core, with a phase boundary inside. A pure strange star is deconfined all the way to a bare surface. Hybrids are bound by gravity; only the self-bound strange variety needs Bodmer-Witten.
• "The surface is a hard shell like a planet's crust." The surface of a bare strange star is a quark-matter boundary a femtometre thick, possibly capped by a thin (~hundreds of metres) crust of ordinary nuclei held up off the quark surface by an intense electric field — a "strange dwarf crust." It is nothing like rock.
• "A strangelet would eat the Earth." The ice-nine fear was taken seriously and dismissed: stable strangelets are most likely positively charged and Coulomb-repelled from nuclei, and cosmic rays far more energetic than any collider have struck the Moon for billions of years without converting it.
• "Quark matter violates the speed of light because it is so stiff." The bag-model sound speed tops out at c_s² = c²/3, comfortably causal. Some interacting models push c_s higher to support 2 M☉, but never past c.