---
title: 'Strangeon-Star Model: Compact Star Physics'
url: https://www.emergentmind.com/topics/strangeon-star-model
type: topic
---

# Strangeon-Star Model: Compact Star Physics

Searching arXiv for recent and foundational papers on strangeon stars to ground the article in the literature.
Strangeon stars are a proposed class of pulsar-like compact stars in which the bulk matter is not ordinary nucleonic matter and not deconfined strange quark matter, but a self-bound condensed phase composed of localized multi-quark clusters called **strangeons**. In this framework, strangeons are “strange nucleons,” i.e. quark clusters with approximate three-light-flavor symmetry involving \(u\), \(d\), and \(s\) quarks, and a strangeon star is treated as a **3-flavored gigantic nucleus** rather than a gravity-bound neutron star [1701.08463; 1706.04707]. The model is motivated by the possibility that, at densities of a few times nuclear saturation density and energy scales of order \(\sim 0.5\ \mathrm{GeV}\), QCD remains strongly coupled, so quarks may cluster rather than deconfine into a weakly interacting Fermi liquid [1706.04707]. A central consequence is that strangeon stars are **self-bound by the strong interaction**, can have a **sharp finite-density surface**, may become **globally solid** at low temperature, and can therefore exhibit phenomenology that differs qualitatively from both conventional neutron stars and MIT-bag strange stars [1701.08463; 1706.04707; 2509.13732].

## 1. Conceptual definition and microphysical basis

The strangeon-star model generalizes Witten’s strange-matter idea by allowing the stable bulk phase of dense matter to consist of **localized quark clusters** rather than free deconfined quarks [1701.08463]. In this terminology, a strangeon is a cluster containing roughly equal numbers of \(u\), \(d\), and \(s\) quarks; specific cluster sizes mentioned in the literature include \(6, 9, 12, 18,\) or more quarks, although the precise number is not fixed [1706.04707; 1701.08463]. The intended distinction is threefold: nucleons are ordinary three-quark baryons without built-in strangeness, quark stars assume deconfined quarks, and strangeon stars assume **confined but clustered** quarks as the relevant effective degrees of freedom [1706.04707].

The microphysical motivation rests on the claim that matter in compact stars occupies a nonperturbative QCD regime. The model emphasizes that at compact-star densities the strong coupling can be large, making a weak-coupling quark description unreliable and suggesting that localization into clusters is plausible [1701.08463]. In parallel, the model assumes approximate **3-flavor symmetry restoration**: dense matter becomes energetically favorable to include strange quarks, but the energy scale is still insufficient for full quark deconfinement [1706.04707]. This is the defining contrast with strange quark stars, where bulk stability is attributed to deconfined quark matter.

A complementary phenomenological construction appears in the liquid-drop treatment of strangeon matter, where strangeon matter is modeled as a 3-flavor analog of ordinary nuclei with energy per baryon
\[
E/A = M + b_\mathrm{vol,\,s} + b_\mathrm{surf,\,s} A^{-1/3}.
\]
In that formulation, stability improves with baryon number, and for \(M\sim \text{GeV}\) and \(\epsilon\sim 100\ \text{MeV}\), the critical baryon number can be as low as \(A_{\rm c}\approx 300\) [1708.03908]. This does not constitute a first-principles derivation, but it provides a phenomenological basis for the claim that bulk strangeon matter could be absolutely stable.

## 2. Self-bound structure, equation of state, and stellar properties

A defining structural assumption is that strangeon matter is **self-bound**, so the stellar surface has a finite density rather than the smooth density decline characteristic of neutron stars [1706.04707; 2109.13234]. In recent applications to accretion physics, this sharp interface is made explicit: the bulk strangeon matter at the surface is described as having density \(\gtrsim 10^{14}\,\mathrm{g\,cm^{-3}}\), while accreted normal matter above it in a polar mound can remain at only \(\sim 10^{11}\,\mathrm{g\,cm^{-3}}\) [2509.13732]. This sharp surface density discontinuity is one of the key physical differences from a neutron star.

Several phenomenological EOS constructions are used in the literature. One widely used family adopts a **Lennard-Jones model** for the strangeon-strangeon interaction,
\[
u(r)=4\epsilon\left[\left(\frac{\sigma}{r}\right)^{12}-\left(\frac{\sigma}{r}\right)^6\right],
\]
or equivalently
\[
u(r)=4U_0\left[\left(\frac{r_0}{r}\right)^{12}-\left(\frac{r_0}{r}\right)^6\right].
\]
With a simple-cubic lattice, this yields energy density and pressure of the form
\[
\rho c^2=2U_0\left(6.2\,r_0^{12}n^5-8.4\,r_0^6n^3\right)+nmc^2,
\]
\[
P=4U_0\left(12.4\,r_0^{12}n^5-8.4\,r_0^6n^3\right),
\]
or equivalent expressions parameterized by \(\epsilon\), \(\sigma\), and \(N_q\) [2401.11754; 2109.13234]. The physical rationale is that the repulsive core makes the EOS very stiff.

A second phenomenological route is the **linked bag model**, which extends bag-model ideas to condensed matter in both 2-flavor and 3-flavor sectors [2008.06932]. In this model the maximum mass of strangeon stars can be as large as \(\sim 2.5M_\odot\), while the tidal deformability of a \(1.4M_\odot\) strangeon star lies in the range of \(180\lesssim \Lambda_{1.4} \lesssim 340\) [2008.06932]. The same work argues that increasing the number of valence quarks per strangeon lowers the energy per baryon and stiffens the EOS, thereby increasing the maximum mass.

The macroscopic consequences of self-binding and stiffness recur across the literature. Strangeon stars are predicted to support high masses, with early overview papers quoting
\[
M_{\max}\sim 2\text{--}3\,M_\odot
\]
[1706.04707], while specific Lennard-Jones calculations report maximum stable baryonic masses as large as
\[
M_{b,\max}^{\rm stable}\big|_{\rm SS}\simeq 4.4\,M_\odot
\]
for the chosen EOS in the nonrotating case [2401.11754]. Because the stars are self-bound, the mass-radius behavior differs from neutron stars: the radius can increase with mass over much of the sequence, and the minimum mass can be very small, even down to **planet-like masses** \(\sim 10^{-2}M_\odot\) in the overview discussion [1706.04707]. Low-mass realizations are used explicitly in the proposed interpretation of Calvera, where a strangeon-star atmosphere fit yields \(R\lesssim 4\,\mathrm{km}\), \(B\lesssim 10^{11}\,\mathrm{G}\), and \(M\lesssim 0.1\,M_\odot\) [1709.07679].

## 3. Rotation, tidal deformability, and merger remnants

The rotating and tidally deformed strangeon-star problem has been studied in perturbative general relativity with particular attention to the role of the **finite surface density** [2109.13234]. In the Hartle-Thorne slow-rotation treatment, the finite surface density introduces crucial surface corrections into the perturbation equations and matching conditions. This distinguishes self-bound strangeon stars from ordinary neutron stars and affects the inferred moment of inertia, quadrupole moment, and tidal Love number [2109.13234].

The same study concludes that strangeon stars satisfy quasi-universal I-Love-Q relations in terms of
\[
\bar I\equiv \frac{I}{M^3},\qquad \bar Q \equiv -\frac{Q_{\rm r}}{M^3\chi^2},\qquad \Lambda \equiv \frac{\lambda}{M^5},
\]
with relative deviations from standard fits generally below \(1\%\) [2109.13234]. It also notes that the conservative constraint \(\Lambda(1.4M_\odot)\lesssim 800\) implies \(M_{\rm max}\lesssim 4.2\,M_\odot\) for the strangeon-star model within the scanned parameter space [2109.13234].

Merger-related calculations exploit the unusually large maximum masses implied by strangeon-matter EOSs. One merger study reports that the tidal deformability for a \(1.4M_\odot\) strangeon star can be
\[
\Lambda(1.4)\approx 381.9,
\]
consistent with GW170817 bounds [1710.04964]. In that picture, two \(\sim 1.4M_\odot\) strangeon stars can merge to produce a **hyper-massive strangeon star** of mass around \(\sim 2.6M_\odot\), with multimessenger signatures distinct from neutron-star mergers [1710.04964]. The same work proposes a “strangeon kilonova” in which early blue emission is powered by decay of unstable strangeon nuggets and later red emission by remnant spin-down [1710.04964].

A complementary post-merger scenario focuses on slowly rotating massive strangeon stars as engines of short-GRB X-ray plateaus. Using the Lennard-Jones EOS parameter set LX3630 with
\[
n_{b,s}=0.36\ {\rm fm}^{-3},\qquad U_0=30\ {\rm MeV},\qquad N_q=18,
\]
the calculation finds that rigid rotation increases the maximum supported gravitational mass along constant-baryon-number sequences by about **9.3%** over \(M_{\rm TOV}\), compared with about **5.6%** for the neutron-star comparison model [2401.11754]. In that framework, spin-down-induced contraction releases gravitational energy,
\[
L_{\rm grav}=\dot E_{\rm grav}=\frac{dE_{\rm grav}}{d\Omega}\frac{d\Omega}{dt},
\]
which is proposed as an alternative power source for short-GRB plateau emission [2401.11754]. The fit results imply magnetic dipole fields in the range
\[
B_p \sim 0.88\text{ to }8.73 \times 10^{15}\ {\rm G}
\]
with efficiencies
\[
\eta_g \sim 0.11\text{ to }0.69,
\]
and initial periods from about \(2.05\) to \(10.94\) ms [2401.11754].

## 4. Surface physics, atmospheres, and the strangeness barrier

The strangeon-star model assigns central importance to the **surface**, because the transition from normal matter to strangeon matter is regulated not only by electromagnetic effects but also by flavor conversion. The relevant concept is the **strangeness barrier**: ordinary 2-flavor matter cannot simply become part of the 3-flavor strangeon phase without weak-interaction conversion [1701.08463; 1706.04707]. This barrier is repeatedly invoked to explain why a strangeon star can be bare, or instead be covered by a corona, atmosphere, or crust of normal matter [1706.04707].

In accretion settings, the barrier is more important than the Coulomb barrier. In the ULXP accretion-column calculation, the Coulomb penetration probability is written as
\[
\tau = \exp\!\left[-\left(\frac{E_G}{E}\right)^{1/2}\right],
\]
with \(E_G = 0.49\,\mathrm{MeV}\) and \(E = k_B T_C\), and the conclusion is that the Coulomb barrier alone is not a strong obstacle [2509.13732]. Instead, because weak conversion is slow, most inflowing material is “bounced back” and accumulates above the surface, producing a thermal mound [2509.13732].

The same surface physics underlies atmosphere models for X-ray dim isolated neutron stars (XDINSs). In the first bremsstrahlung-atmosphere model, a strangeon star is surrounded by a thin, two-temperature plasma atmosphere formed and maintained by ISM-accreted matter because the strangeness barrier inhibits direct incorporation into the star [1603.08288]. The emergent flux is written as
\[
F_{\nu}^{\infty}\simeq \pi\left(\frac{R_{\mathrm{opt}^{\infty}}}{d}\right)^2 B_{\nu}\left(1-e^{-\tau_{\infty}(\nu)}\right),
\]
with optical depth
\[
\tau_{\infty}(\nu)= 3.92\times10^{-45} \frac{n_{\mathrm{i0}^{2}}(kT_{\mathrm{i}})_{\mathrm{keV}}R_{\mathrm{km}}}{(h\nu)_{\mathrm{keV}^{3.5}\left(1-e^{-h\nu/(kT_{\mathrm{e}})}\right)}}.
\]
The fitted radiation radii of XDINSs are from 7 to 13 km, while the modelled electron temperatures are between 50 and 250 eV, except RX J0806.4-4123 with a radiation radius \(\sim 3.5\) km [1603.08288].

The later nonuniform-atmosphere extension attributes the optical/UV excess and spectral deviation to **bremsstrahlung emission from a nonuniform plasma atmosphere** supplied by accretion that is funneled to the poles and then spreads over the surface [1705.03763]. The model uses
\[
n_e = n_i = n_{i0}(\theta)\exp\!\left(-\frac{m_i g z}{kT_i}\right),\qquad
n_{i0}(\theta)=\frac{n_0}{1+\xi\,\theta^\gamma},
\]
and reports electron temperatures \(\sim100-200\) eV with radiation radii \(R_{\rm opt}^{\infty}\sim5-14\) km [1705.03763]. For five XDINSs—RX J0720.4−3125, RX J0806.4−4123, RX J1308.6+2127, RX J1605.3+3249, and RX J1856.5−3754—the spectra from optical/UV to X-ray bands could be well fitted and exhibit gaussian absorption lines at \(\sim 100-500\) eV [1705.03763].

The same atmosphere framework is applied to Calvera. A joint fit to **XMM-Newton** and **Chandra** data with the strangeon-star atmosphere model yields \(kT_e \approx 0.61\text{--}0.67\,\mathrm{keV}\), reduced chi-square values around \(\chi_\nu^2 \simeq 1.00\text{--}1.09\), and supports the interpretation of Calvera as a **low-mass strangeon star** with inactive magnetosphere, small radius, and likely fallback-disk braking [1709.07679].

## 5. Solid-state dynamics: glitches, oscillations, and free energy

A major claim of the strangeon-star framework is that the star can be **globally solid**, not merely crusted, and that this opens a unified class of starquake phenomena [1701.08463; 1706.04707]. Glitch modeling is a central application. In the original glitch formulation, the starquake is decomposed into **plastic flow** in a fractured outer layer and **elastic motion** in the inner region [1707.07471]. The moment of inertia is parameterized as
\[
I=I_0(1+\epsilon)(1+\eta),
\]
and the recovery coefficient is defined as
\[
Q=\frac{\Delta\Omega_d}{\Delta\Omega_g}.
\]
The model proposes an empirical relation
\[
Q=\exp\left(-a\frac{\Delta \Omega_g}{\Omega}\right),
\]
with observationally reasonable range \(a \sim 10^6 - 10^7\) [1707.07471]. It also derives the inter-glitch waiting time
\[
t_q\simeq \frac{|\Delta \sigma|}{\dot{\sigma}}
\]
and argues that Crab- and Vela-like glitches can be explained in a unified bulk-invariable starquake picture with small radiative output [1707.07471].

The activity analysis extends this framework statistically. Using the observed relation
\[
\dot{\nu}_g \simeq 0.01\,|\dot{\nu}|,
\]
the strangeon-star glitch-activity model infers that the shear modulus must be
\[
\mu \simeq 3\times 10^{34}\ {\rm erg\,cm^{-3}}
\]
for a canonical star of mass \(1.4M_\odot\) and radius \(10\,{\rm km}\) [2011.01496]. It further argues that the observed glitch activity requires **about ten times the oblateness shift accumulated during the glitch interval**, i.e. \(\Delta \epsilon_m\sim 10\,\Delta\epsilon\), and estimates energy releases of
\[
\Delta E \sim 2.4\times 10^{40}\,{\rm erg}
\]
for a Vela-like glitch with \(\Delta \nu/\nu \sim 10^{-6}\), and
\[
\Delta E \sim 4.2\times 10^{41}\,{\rm erg}
\]
for a Crab-like glitch with \(\Delta \nu/\nu \sim 10^{-8}\) [2011.01496].

The recovery model further interprets post-glitch relaxation as a pressure-restoring viscous inflow of fragments into pressure-deficient equatorial cracks [2301.09088]. The outflow velocity is assumed to decay as
\[
\upsilon(t)=\upsilon_0 e^{-t/\tau},
\]
leading to the usual exponential frequency recovery
\[
\nu(t)\simeq \nu_{0-}+\Delta\nu_p+\Delta\nu_d\,e^{-t/\tau}.
\]
Fits to five glitches give representative cracking depths and viscous timescales such as \(h/R=0.005,\ \tau=120\) d for B1838-04 and \(h/R=0.1,\ \tau=50\) d for B1800-21 [2301.09088]. The reported trend is that \(h/R\) increases with glitch size \(\Delta\nu/\nu\) [2301.09088].

The solid-state interpretation is also applied to magnetar giant-flare quasi-periodic oscillations. In the strangeon-star torsional-mode study, the large shear modulus
\[
\mu \sim 10^{32}\,{\rm erg\,cm^{-3}}
\]
pushes the fundamental \(\ell=2\) torsional mode to roughly
\[
{}_{2}f_0 \simeq 145\text{--}277~\mathrm{Hz},
\]
while first overtones are roughly
\[
{}_{\ell}f_1 \sim 300\text{--}1700~\mathrm{Hz}
\]
[2309.09847]. This enables identifications such as \(150\) Hz with \({}_{2}f_0\), \(625\) Hz with \({}_{2}f_1\), and \(1837\) Hz with \({}_{3}f_6\) [2309.09847]. Low-frequency QPOs are then attributed to ocean-crust interface modes with ocean densities \(\rho \sim 10^6\)–\(10^9\,\mathrm{g\,cm^{-3}}\), temperatures \(T \sim 10^8\)–\(10^9\) K, and thickness \(\sim 10\)–\(50\) m [2309.09847].

A related line of work studies stored elastic/gravitational free energy through pressure anisotropy in general relativity. Using the anisotropy ratio
\[
\frac{p_t-p_r}{p_r},
\]
the calculation finds that values of order
\[
\frac{p_t-p_r}{p_r} \sim 10^{-4}
\]
are already sufficient for
\[
\Delta E_b \gtrsim 10^{46}\ \mathrm{erg},
\]
with the broader range \(\Delta E_b \sim 10^{44} - 10^{47}\ \mathrm{erg}\) for \(\eta_{1,2}\sim 10^{-4}-10^{-3}\) [2305.19687]. This is presented as an alternative energy reservoir for SGR giant flares, FRBs, and related transients without requiring extremely strong magnetic fields [2305.19687].

## 6. Thermal evolution, neutrinos, and accretion-powered high-energy phenomena

The strangeon-star model has also been developed as a thermal-evolution framework. In the supernova context, a newborn strangeon star is assumed to be born as a hot liquid with internal energy of order
\[
E_{\rm bind} \simeq 1.5\times10^{53}\left(\frac{M}{M_\odot}\right)\ {\rm erg},
\]
and to cool through neutrino and photon emission until it reaches a melting temperature
\[
T_m \sim 1\!-\!2\ {\rm MeV}
\]
[1705.08188]. The total internal energy is split into
\[
U = U_{\rm s}+U_e+U_{\rm pion}+U_\nu+U_\gamma,
\]
with pion excitation emphasized as a major high-temperature contribution [1705.08188]. Because the neutrino mean free path is estimated as
\[
l \sim 10^3\ {\rm cm},
\]
the star is opaque and cooling is diffusive [1705.08188]. The model identifies the sharp drop in SN1987A neutrino emission with the liquid-solid transition of strangeon matter [1705.08188]. This suggests that a future core-collapse supernova neutrino burst showing a cut-off or sharp drop would be a characteristic test of the model.

A much more recent extension applies the strangeon-star surface physics to **ultraluminous X-ray pulsars**. There the accretion column is solved in a one-dimensional steady model, and the weak-conversion probability is taken as \(\eta \sim 10^{-13}\) with \(t^*_{\rm weak}=10^{-10}\,\mathrm{s}\), yielding a mound-base density \(n_{\rm base}\sim 10^{35}\,\mathrm{cm^{-3}}\) and \(\rho_{\rm base}\sim 10^{11}\,\mathrm{g\,cm^{-3}}\) [2509.13732]. The resulting thermal mound can reach
\[
H \simeq 0.7\text{--}0.95\,\mathrm{km},
\]
with base temperature
\[
T \gtrsim 10^9\,\mathrm{K},
\]
and efficient neutrino cooling via **electron–positron pair annihilation** [2509.13732].

The heating rate entering the neutrino luminosity is written as
\[
L_\nu \approx S\int_0^H Q_{\rm acc}(H)\,dH, \qquad Q_{\rm acc}(H)=\frac{\dot M}{2S}\left[\frac{GM}{(R+H)^2}\right].
\]
The luminosity decomposition is
\[
L_{\rm tot}=\frac{GM\dot M}{R}=L_\nu+L_{\rm ph}+L_{\rm heat},
\]
with a surface heat-transport contribution
\[
L_{\rm heat} \simeq 1.2\times 10^{36}\,\mathrm{erg\,s^{-1}},
\]
emphasized to be essentially independent of the accretion rate [2509.13732]. Two regimes are distinguished:
\[
\dot M < 10^{20}\,\mathrm{g\,s^{-1}}
\]
where photons dominate, and
\[
\dot M > 10^{21}\,\mathrm{g\,s^{-1}}
\]
where photon trapping makes neutrinos the main cooling channel and the photon luminosity saturates near
\[
L_{\rm ph} \sim 10^{41}\,\mathrm{erg\,s^{-1}}.
\]
For detectability, only Swift J0243.6\(+\)6124 is estimated to reach a marginally interesting Earth flux,
\[
F_\nu \sim 2.3\times 10^{-2}\,\mathrm{MeV\,cm^{-2}\,s^{-1}},
\]
while extragalactic ULXPs remain far below current backgrounds [2509.13732].

## 7. Observational status, comparisons, and open issues

The strangeon-star model is presented throughout the literature as a **third alternative** to conventional neutron stars and strange quark stars [1706.04707]. Compared with neutron stars, it posits self-bound matter, a finite surface density, stronger surface binding, a potentially globally solid interior, and a typically stiffer EOS [1706.04707; 1701.08463]. Compared with quark stars, it replaces deconfined quarks with localized quark clusters and often emphasizes solid-state behavior [1706.04707; 2509.13732].

This distinction matters because several observational signatures invoked in the model depend specifically on clustering and solidity rather than on generic self-binding. Examples include: starquake-driven glitches and their recovery [1707.07471; 2301.09088], high-frequency magnetar QPOs from torsional oscillations [2309.09847], phase-transition cutoffs in supernova neutrino bursts [1705.08188], and accretion thermodynamics shaped by the **strangeness barrier** [2509.13732]. By contrast, properties such as high compactness, large surface redshift, or self-bound mass-radius behavior are shared more broadly with MIT-bag strange-star models [1406.4118; 1509.00401]. This suggests that some observables constrain self-binding in general, whereas others are intended to probe specifically strangeon clustering or solid-state physics.

Several studies explicitly state that the model remains **phenomenological**. The liquid-drop and linked-bag papers stress that the detailed microscopic dynamics are uncertain and not derived from first-principles QCD [1708.03908; 2008.06932]. The rotational and tidal calculations likewise note that perfect-fluid perturbation theory neglects the true shear modulus of a solid strangeon star, which may matter quantitatively for deformations and oscillations [2109.13234]. The QPO analysis adopts a phenomenological Lennard-Jones EOS rather than deriving the interaction from QCD [2309.09847]. The accretion-column work in ULXPs makes specific choices for weak-conversion efficiency and timescale, which control the mound structure and neutrino luminosity [2509.13732]. These are not contradictions within the framework, but they delimit the current status of the model.

The observational program proposed in the literature is correspondingly broad. FAST and SKA are cited as probes of drifting subpulses, glitches, precession, and other timing phenomena [1701.08463; 1706.04707]. eXTP, NICER-like pulse-profile modeling, and future X-ray polarimetry are invoked as tests of surface composition, atmosphere geometry, and compactness [1701.08463; 2109.13234; 1705.03763]. Gravitational-wave measurements of \(\Lambda\), \(\tilde\Lambda\), and post-merger dynamics are expected to constrain or falsify regions of the strangeon-matter parameter space [1710.04964; 2109.13234; 2401.11754]. Future supernova neutrino observations could test the predicted liquid-solid transition cutoff [1705.08188], and ULXP neutrino searches, although presently unfavorable for all but the nearest Galactic case, provide a further possible diagnostic of strangeon-star surface physics [2509.13732].

Taken together, the strangeon-star model constitutes a coherent but still conjectural program: dense matter is assumed to form a self-bound, clustered 3-flavor phase; compact stars built from this matter acquire a sharp surface, very stiff EOS, and often a globally solid interior; and these properties are then used to reinterpret a wide range of compact-star phenomena—from glitches and magnetar QPOs to XDINS spectra, short-GRB plateaus, merger remnants, supernova neutrino bursts, and ULXP accretion columns [1701.08463; 1706.04707; 2509.13732]. A plausible implication is that the most discriminating tests will be those that directly probe the conjunction of **self-binding**, **surface conversion barriers**, and **solid-state response**, rather than compactness alone.

Source: https://www.emergentmind.com/topics/strangeon-star-model