QCD Magnetars: When Mass, Not Spin, Unleashes a Superstrong Magnetic Field ✛

Original: "Beyond Spin: QCD Magnetars"
· Rachid Ouyed
arXiv:2607.17556v1 · 2026-07-20 · CC BY 4.0 · 6 min
A unified model explains anomalous X-ray pulsars, superluminous supernovae, and fast radio bursts as stages in the life of a hybrid neutron star with a ferromagnetic quark core.

Context

In June 2024, the CHIME radio telescope caught the hyperactive repeater FRB 20240114A firing off over 700 bursts in a single hour—an event that unleashed nearly all the magnetic energy of an ordinary magnetar. These fast radio bursts and dozens of other enigmatic transients—anomalous X-ray pulsars, superluminous supernovae, and luminous fast blue transients—have puzzled astrophysicists for years. Where do the fields come from? Traditional models tie the field to rapid spin of a newborn neutron star, but the birth rate of such “spun‑up” objects is a hundred times lower than the observed number of magnetars. Massive star evolution (stellar evolution) produces neutron stars, but only a negligible fraction have millisecond periods. Ouyed (2026) offers a different view: the field is set not by spin, but by mass. A neutron star heavier than ~2,1 M☉ carries a dormant ferromagnetic volcano from birth—a quark core ready to awaken as soon as spin‑down blows the fuse. The magnetic energy in a single cubic centimeter of such a core (4×$10^{34}$ erg) equals the energy the Sun radiates in a century. Fritz Zwicky already suspected that neutron stars are born in supernovae, and now we see that their true power may slumber for ages.

Methods

To test the hypothesis, we employed direct Monte‑Carlo population synthesis. From log‑normal distributions of mass, spin period, and magnetic field—consistent with the observed radio pulsar population (median values: mass 1,35 M☉, period 200 ms, field $10^{12}$,5 G)—we randomly generated $10^7$ neutron stars. For each, we computed the spin‑down rate and the time until a quark nova would trigger. The outcome was sorted by observable channel: whether the supernova ejecta remained optically thick at the time of the quark nova; whether the period exceeded the critical value $P_{fast}$ (~5,5 ms); whether a solid crust had time to form, and so forth. Matching the measured relative rates of five event types—the fraction of magnetars among supernova remnants, the fraction of Type I superluminous supernovae, the fraction of double‑humped superluminous supernovae, luminous fast blue transients, and the association of magnetars with supernova remnants—allowed a Bayesian sampler (MCMC) to reconstruct the posterior distributions of the key parameters $M_{dec}$ and $P_{fast}$. Direct counting of outcomes in a large sample sidesteps the need for approximations.

Results

The analysis yields $M_{dec}$ ≈ $2,05^{+0}$,$30_{-0}$,32 M☉ and $P_{fast}$ ≈ $5,45^{+4}$,$86_{-3}$,11 ms (68% confidence intervals). The widths of the natal distributions—0,15 dex in mass, 0,72 dex in period, and 0,61 dex in field—lie within the range of independent estimates. The critical mass is about 2,1 M☉. The model reproduces all five observational constraints simultaneously: the fraction of magnetars among all neutron stars is ~11%, close to the observed 5–15%; the fraction of Type I superluminous supernovae is about 3,2×10^−4 of the core‑collapse rate; the double‑humped fraction is less than 10% of those; the fraction of luminous fast blue transients is ~1,1×10^−3; and roughly 30% of magnetars are associated with supernova remnants, as required by the statistics of the McGill catalog (8 out of 30). Much like the Chandrasekhar limit for white dwarfs, a critical mass here separates the fates of compact objects.

Implications

If mass, not spin, decides the field’s destiny, then every magnetar is a former massive neutron star, not a lucky one born spinning exceptionally fast. This resolves a half‑century‑old contradiction. Moreover, the delay between core collapse and quark nova, set by the spin‑down rate, stretches from days to millions of years, naturally explaining the appearance of fast radio bursts in the most unexpected places—from the outskirts of elliptical galaxies to globular clusters. The quark nova itself becomes a new source of r‑process nucleosynthesis: its neutron‑rich ejecta synthesize lanthanides and actinides, producing a kilonova without a neutron‑star merger. The crust‑free phase of a QCD magnetar, lasting centuries until a crust builds up, explains X‑ray‑quiet repeating radio bursts. In an extended picture, a quark nova can also launch gamma‑ray bursts (a relativistic jet), and a gravitational‑wave burst makes it a target for multi‑messenger astronomy. Decades after Jocelyn Bell Burnell discovered the first pulsars, we are beginning to see that neutron stars are far more diverse.

Future development

The model makes a series of testable predictions. First, for all magnetars, the field should correlate with mass, not with spin period. Second, luminous fast blue transients must show an infrared excess—the afterglow of a quark‑nova kilonova. Third, repeating fast radio bursts in their early stages should lack X‑ray counterparts, in contrast to the standard crust‑fracture model. Finally, old superluminous supernovae discovered in the 2000s (SN 2005ap, SN 2006gy) are now entering their transparency window for gigahertz radiation, and targeted monitoring in the next few years could catch the awakening of a hidden FRB repeater. On the flip side, finding even one magnetar with a millisecond birth period and a hadronic mass above 2,1 M☉ would put the model in peril.

Impact

The findings touch multiple fields at once—from the physics of dense quark matter to the chemical evolution of galaxies. The scenario imposes concrete requirements on the equation of state of the hadron‑quark transition and on ferromagnetic ordering of quarks. In astrophysics, it ties together magnetar and supernova populations into a single knot, predicting a new class of isolated kilonovae and a legacy Faraday rotation in radio bursts.

Next steps

The immediate next steps include measuring magnetar masses through pulsar‑pair timing or X‑ray pulse profiling, and searching for infrared kilonovae in samples of luminous fast blue transients. In parallel, multi‑dimensional MHD simulations of current‑driven instabilities at the boundary between quark core and hadronic envelope are needed.

Key open problems

The work directly ties into several unsolved problems in physics: the nature of the hadron‑quark phase transition inside neutron stars, the origin of relativistic jets and magnetic fields in collapsars, and the question of astrophysical sites for heavy‑element synthesis. The conclusions about delayed deconfinement can be tested by future gravitational‑wave detectors capturing non‑equilibrium oscillations of a hybrid star. Someday, a new hyper‑active repeater caught by CHIME will deliver not only radio bursts but also the signature of a kilonova—and then the dormant volcano of the quark star will finally awaken in our imagination.

The knowledge machine's view

This section is addressed to the author of the paper and their colleagues. Neighbouring works were found by meaning-based search across our archive, not by word overlap. We do not review: we show the field around the work.

What the work does

The paper proposes the quark-nova hypothesis as a common mechanism for five types of transients (FRBs, magnetars, superluminous supernovae, etc.), where the critical mass of a neutron star determines the birth of an ultra-strong magnetic field. Using Monte Carlo population synthesis with Bayesian sampling, they obtain estimates of M_dec ≈ 2.05 M☉ and P_fast ≈ 5.45 ms, which reproduce the observational constraints.

Where its strength is

The strength lies in unifying several classes of transients under a single physical mechanism—a phase transition to quark matter—that naturally links stellar mass to magnetic field generation, bypassing the problem of a shortage of rapidly rotating progenitors.

Where it could go next

Worth trying

  • Test how well the model's predicted magnetic field distributions agree with estimates from turbulent dynamo in a related work.
  • Apply the method of analyzing FRB time structure (the Pincus–Lyapunov diagram) to synthetic bursts generated within the quark-nova scenario.
  • Consider the possibility of joint constraints on the triggering mass of quark-novae from population synthesis and independent neutron star mass data.

Why this may matter

If the quark-nova indeed accounts for several classes of transients, it provides a unified observational signature to test the existence of quark matter in neutron stars, which is important for the physics of strong interactions and the astrophysics of compact objects.

Original abstract

as written by the authors arXiv:2607.17556v1

We present a unified framework in which anomalous X-ray pulsars (AXPs), soft gamma-ray repeaters (SGRs), superluminous supernovae (SLSNe-I), luminous fast blue optical transients (LFBOTs), and fast radio bursts (FRBs) originate from quark deconfinement in the core of a massive neutron star (NS). Spontaneous ferromagnetism in the deconfined phase generates ~10^18 G core fields, producing a highly magnetized hybrid star (HS) - the "QCD magnetar" - whose surface field is set by NS mass, not birth spin. The Quark-Nova (QN) that forms the HS ejects ~0.01M_sun of neutron-rich outer layers, powering a kilonova and leaving the HS crustless for centuries. During this phase, magnetic instabilities at the hadron-quark interface release rising flux ropes that power X-ray-quiet FRBs; as the crust reforms, the source evolves into an X-ray-loud AXP&SGR. Two parameters govern the model: a critical mass M_dec triggering deconfinement, and a critical period P_fast separating fast and slow rotators. Fast rotators inject spin-down energy into the QN ejecta, producing an LFBOT - directly observable once the SN ejecta is optically thin, or via binary accretion with no preceding SN; otherwise the LFBOT is reprocessed by the SN ejecta into an SLSN-I. A Bayesian Monte Carlo population synthesis reproduces the observed rates of AXPs&SGRs, SLSNe-I, and LFBOTs with M_dec ~2.1M_sun and P_fast ~5.5 ms, and predicts non-merger r-process signatures and kilonovae from isolated NSs, with or without an LFBOT. The QN ejecta also carries its own DM and RM, independent of environment, predicted to appear as excess dispersion and rotation measure in FRBs once Galactic, host, and intergalactic contributions are removed. These provide direct observational tests of the hadron-quark phase transition and ferromagnetic ordering in dense quark matter.

The magnetic field of $10^{18}$ G generated in the quark core is so immense that if it filled the volume of Earth, its total energy would equal all the sunlight the Sun has emitted over a billion years.

t_{\rm NS,SpD} \simeq 6.3\times 10^{4}\, {\rm yr} \times \left(\frac{P_0}{0.2\, {\rm s}}\right)^2 \left(\frac{B_0}{10^{12.5}\, {\rm G}}\right)^{-2}
The time it takes for the spin period to double under magnetic dipole radiation.
t'_{\rm HS,SpD} \sim 12.7\, {\rm yr} \times \left(\frac{P_{\rm HS}}{0.28\, {\rm s}}\right)^2 \left(\frac{B_{\rm HS}}{10^{15}\, {\rm G}}\right)^{-2}
The time for a QCD magnetar to exhaust its rotational energy after the quark nova.

Key numbers

  • Critical mass: 2,1 M☉
  • Critical period: 5,5 ms
  • Magnetic field in core: ~1018 G
  • Surface field: ~1015 G
  • Crust-free phase duration: centuries
Scientists
Albert EinsteinBarry BarishKip ThorneRichard FeynmanRussell HulseSubrahmanyan Chandrasekhar
Concepts
hyperon magnetorotational supernova quark gluon plasma kilonova supernova nucleosynthesis neutron star light emission quark star quantum chromodynamics gamma ray burst merger magnetar pulsations in grb baryon radio pulsar pulsar stellar nucleosynthesis X-ray radiation nuclear spin crab nebula strong interaction hadron relic neutrino