Magnetars are neutron stars endowed with magnetic fields exceeding roughly 10^14 gauss, a quadrillion times stronger than Earth’s field and hundreds to thousands of times more intense than the fields of ordinary radio pulsars. Their magnetic energy reservoirs power some of the most spectacular transients known: superluminous supernovae, whose light curves are widely modeled as magnetar-driven; long gamma-ray bursts, for which millisecond protomagnetars have been proposed as central engines; and at least some fast radio bursts, the millisecond-duration radio flashes whose repeaters have been firmly linked to a Galactic magnetar. Despite this central role in high-energy astrophysics, the rate at which magnetars actually form relative to ordinary neutron stars has remained stubbornly uncertain, preventing direct, quantitative comparisons between magnetar birth rates and the observed rates of these cosmic beacons.

The new study, led by Celsa Pardo-Araujo of the Institute of Space Sciences (ICE, CSIC) and the Institute of Space Studies of Catalonia, with Nanda Rea, Michele Ronchi and Vanessa Graber, attacks the problem from an unexpected direction: by taking a complete census of the very youngest neutron stars in the Milky Way. The team focused on isolated neutron stars younger than two thousand years, a sample in which observational incompleteness is minimized and each object retains the magnetic and thermal fingerprints of its birth. In this young Galactic population, only nine of the twenty-three detected sources, or about forty percent, are classical rotation-powered pulsars, the lighthouse-like radio beacons that have historically dominated neutron star catalogs. The remaining sixty percent are magnetars or so-called central compact objects, dim X-ray sources embedded in young supernova remnants whose surfaces show puzzling thermal properties.

To translate this snapshot of the youngest population into a birth fraction, the researchers constructed a population synthesis of isolated neutron star classes in the Galaxy using the ML-Poppyns code, a framework previously refined with simulation-based inference techniques for radio pulsar populations. The synthesis tracks each simulated star from birth through its full dynamical trajectory through the Galactic gravitational potential, incorporating natal kicks imparted by the asymmetric supernova explosion, realistic spatial distributions of neutron star birth sites traced by spiral arms and the Galactic disk, and selection effects that determine which sources would actually be detected by X-ray and radio surveys. Crucially, the code also follows each star’s physical evolution: its spin-down under magnetospheric torques computed from force-free magnetosphere models, and its magneto-thermal evolution, in which the coupled decay of the crustal and core magnetic field through Ohmic dissipation, Hall cascade and ambipolar diffusion simultaneously heats the star and reshapes its X-ray luminosity over millions of years.

A key ingredient is the assumed distribution of initial magnetic fields. The team adopted a bimodal model, a double log-normal distribution in which one component, centered near 10^13 gauss, corresponds to the radio pulsar population inferred from earlier pulsar population synthesis, while a second, broader component peaks at initial dipole fields of approximately 1 to 2.5 multiplied by 10^14 gauss and corresponds to magnetars. This bimodality reflects a growing realization that magnetars are not merely the extreme tail of a single continuous field distribution but may form through distinct physical channels, whether from exceptionally strong internal fossil fields inherited from their massive progenitor stars or from convective dynamo action in the first seconds after collapse. The relative normalization of the two components directly encodes the magnetar birth fraction, which the team varied to test which values reproduce the observed census.

The result is striking: magnetars peaking at those super-strong initial fields represent on average about fifty percent of the isolated neutron star population, roughly double the fraction of classical rotation-powered pulsars among the youngest detected sources and substantially larger than magnetar fractions of around ten percent or less previously inferred from earlier studies. When the authors compared simulated population counts of magnetars, rotation-powered neutron stars and X-ray-dim isolated neutron stars against the observed numbers for a grid of magnetar birth fractions, initial magnetar field strengths and supernova rates, only models with high magnetar fractions could simultaneously match all three classes. The best-fitting models also placed constraints on the environment in which these stars are born: reproducing the observed population requires a Galactic core-collapse supernova rate larger than two events per century.

That rate requirement is itself provocative. Traditional estimates of the Milky Way’s core-collapse supernova rate, often derived from observations of extragalactic supernovae scaled to the Galactic star formation rate, hover near one to two per century, and radio and gamma-ray surveys of supernova remnants have historically suggested lower values of a few tenths per century. The new analysis, calibrated against the census of supernova remnants with robust neutron star associations younger than ten thousand years within two kiloparsecs of the Sun, indicates that the true Galactic explosion rate must be higher than many of those inferences. A higher supernova rate simultaneously alleviates the tension of hosting so many very young, strongly magnetized neutron stars and implies that the Galaxy’s stellar graveyards are being restocked faster than conventional counts suggest.

The implications ripple outward to extragalactic transient astronomy. Magnetar-powered models are the leading interpretation for superluminous supernovae, whose peak luminosities can exceed ordinary supernovae by factors of a hundred, and for the extended emission of many gamma-ray bursts. If half of all neutron stars are born as magnetars, the raw supply of potential engines for these events is far larger than previously assumed, sharpening the question of why only a small subset of magnetars apparently powers such extreme displays. The elevated birth fraction also eases rate-budget tensions for fast radio burst populations, where the volumetric rate of bursts must be reconciled with the birth rate of their proposed progenitors. More broadly, recent work has shown that magnetars are evolutionarily intertwined with other neutron star classes: some central compact objects show magnetar-like outbursts, some young pulsars display magnetar-like flares, and transitional objects such as the fast-spinning magnetar Swift J1818.0-1607 blur the line between the classes. A population synthesis that treats all isolated neutron stars within a single magneto-thermal evolutionary framework, as this study does, is precisely the tool needed to capture those connections.

The team drew its observational backbone from publicly archived catalogues, including the ATNF Pulsar Catalog, the McGill Online Magnetar Catalog and comprehensive supernova remnant catalogues, with thermal luminosities for the seven X-ray-dim isolated neutron stars taken from the literature. The authors note that all observational data are publicly available through NASA and ESA archives, and the ML-Poppyns population-synthesis code has been released on GitHub, allowing the community to test and extend the results. As time-domain surveys such as the Vera C. Rubin Observatory begin to catch superluminous supernovae and fast radio bursts in unprecedented numbers, an accurate magnetar birth fraction becomes an essential ingredient for every population-level comparison. This study’s message is unambiguous: magnetars are not exotic outliers but a mainstream product of stellar death, and the Milky Way has been quietly forging them at a pace that our models of stellar collapse are only now catching up to.

The study’s reliance on a bimodal field distribution carries an important caveat: the inferred fraction depends on where the magnetar component peaks in field strength. Models with lower peak fields of roughly 10^14 gauss can accommodate somewhat different normalizations, which is why the authors report the result as an average across the plausible range rather than a single definitive number. This sensitivity mirrors a broader debate in the literature, where formation-rate estimates have historically ranged from a few percent to tens of percent depending on which observational class was used as a tracer.

The treatment of central compact objects deserves particular attention. These sources, found exclusively within young supernova remnants, have long resisted classification, and evidence that at least some of them harbor strong internal fields, including magnetar-like outbursts from otherwise quiescent objects, supports the paper’s decision to model them within the same magneto-thermal framework as magnetars. If a substantial share of central compact objects are hidden or incipient magnetars whose dipole fields have not yet surfaced, the effective magnetar fraction rises accordingly.

Independent constraints may soon emerge from magnetar giant flares, the rare gamma-ray flashes whose extragalactic detection rates provide a complementary census of the magnetar population. Combining flare statistics with the elevated birth fraction reported here will test whether the flaring subset is a small, biased sample or representative of the class. Similarly, the growing catalog of fast radio burst repeaters, now anchored by the Galactic magnetar SGR 1935+2154, offers an empirical cross-check that future population models can exploit.

Methodologically, the work illustrates the value of simulation-based inference in stellar population studies, replacing ad hoc forward modeling with a statistical comparison across parameter space. Extending the framework to binary neutron stars and to extragalactic populations with different metallicities and star formation histories would be a natural next step.

Subject of Research: The birth fraction of magnetars among neutron stars formed in core-collapse supernovae, constrained by population synthesis of the young Galactic neutron star population

Article Title: The magnetar fraction in core-collapse supernovae

Article References: Pardo-Araujo, C., Rea, N., Ronchi, M., & Graber, V. (2026). The magnetar fraction in core-collapse supernovae. Nature Astronomy. https://doi.org/10.1038/s41550-026-02915-5

Image Credits: AI Generated

DOI: 10.1038/s41550-026-02915-5

Keywords: magnetars, neutron stars, core-collapse supernovae, population synthesis, magneto-thermal evolution, fast radio bursts, superluminous supernovae, gamma-ray bursts, Galactic supernova rate, central compact objects, high-energy astrophysics, Nature Astronomy

Stars begin life fusing hydrogen into helium, releasing energy that counteracts gravitational collapse and maintains their stability. As hydrogen is exhausted, stars evolve through a sequence of fusion stages, each burning progressively heavier elements. This process carves the star into concentric shells: an outer hydrogen layer, followed inwardly by helium, carbon-oxygen, oxygen-neon-magnesium, and deeper layers rich in oxygen, silicon, and sulfur. This stratification culminates in the fusion of silicon and sulfur into iron-group elements, a pathway that ultimately heralds the star’s cataclysmic demise via core collapse, often manifesting as spectacular supernovae or formation of black holes.

Traditionally, direct evidence for these internal shells, especially the ones rich in silicon (Si) and sulfur (S), has been elusive. Most observed stellar explosions correspond to stars stripped only down to their helium or carbon-oxygen layers. Such stripped stars expose the evolutionary products of outer shells, but the deeper layers remain obscured by the stellar envelope or lost in explosive dynamics. This limitation has left a crucial gap in our empirical understanding of late-stage nucleosynthesis—the production of elements heavier than oxygen in the chaotic environment preceding core collapse.

The recent observations of SN 2021yfj mark a turning point. Astronomers have captured signals indicating the progenitor star had shed its outer hydrogen and helium layers, unveiling a massive shell where silicon and sulfur dominate. The implication is profound: this supernova originated from a star stripped down to its O/Si/S-rich inner shell, an unprecedented glimpse into advanced stellar evolution stages. The ejecta contain clear signatures of silicon, sulfur, and even argon, elements formed in highly energetic fusion reactions before the star’s final explosive death.

Exposing these inner layers before the explosion offers unique clues about the star’s evolutionary pathway and mass-loss mechanisms. The standard theory predicts that peeling back a star’s envelope to reveal such refractory, inner shells requires intense interactions or rare, violent mass-loss episodes shortly before collapse. The detection of a thick, circumstellar shell composed chiefly of Si and S material expelled immediately prior to the supernova suggests an atypical shedding process, potentially through pulsational instabilities or binary interactions, not commonly observed in massive star evolution.

Spectroscopic analyses of SN 2021yfj supplied decisive evidence for this deeply stripped progenitor. Early spectra revealed emission and absorption lines characteristic of silicon and sulfur ions, markedly different from typical Type Ib or Ic supernovae, where helium or carbon signatures prevail. The strength and velocity profiles of these lines indicate a dense, massive shell enveloping the star—a reservoir of freshly synthesized elements hurled outward before the star’s core collapsed.

This discovery extends our comprehension of nucleosynthesis and the diversity of supernova progenitors, challenging existing paradigms. While stellar evolutionary models have predicted layered interiors featuring silicon and sulfur shells, the direct detection of such material in the circumstellar environment confirms and refines these models. It provides a rare window into the final phases of massive star life, where fusion stages race towards the synthesis of the iron peak, shaping galactic chemical evolution.

Furthermore, the finding bears implications for understanding the mechanics behind different supernova types. Stripped-envelope supernovae—those lacking hydrogen and sometimes helium in their spectra—have long been linked to binary interactions or strong stellar winds removing outer layers. SN 2021yfj adds a novel category: a star exploded after extreme stripping that exposed and expelled its inner Si/S-rich strata. This challenges theorists to explain how such severe mass loss occurs naturally and what triggers it in the critical final years or months before core collapse.

Astrophysicists studying SN 2021yfj will likely investigate whether this mass loss was episodic, perhaps related to pulsational pair-instability or other advanced stellar instabilities causing violent outbursts. Alternatively, closely orbiting companions in binary systems might strip the progenitor’s outer layers during tight, late-stage interactions. Identifying and modeling these mechanisms could illuminate why such events are rare and how they influence the ultimate fate of massive stars.

Importantly, SN 2021yfj provides empirical evidence that enriches nucleosynthetic yields used in galactic chemical evolution studies. Knowing that silicon and sulfur can be ejected in circumstellar shells prior to explosion impacts predictions about elemental abundances traveling through the interstellar medium. This, in turn, affects interpretations of cosmic material cycling and the origins of elements essential to planet formation and life.

The achievement also underscores the power of multiwavelength observational campaigns in capturing transient phenomena. Coordinated spectroscopy and photometry, combined with theoretical modeling, enabled researchers to reconstruct the progenitor’s structure and mass-loss history despite the inherent challenges of studying distant, rapidly evolving supernovae. Such capabilities will be pivotal in identifying future rare events exposing even deeper layers, such as iron core material, pushing the boundaries of explosive stellar astrophysics.

Looking ahead, astronomers aim to monitor for similar stripped-envelope supernovae exhibiting unusual spectral features. Broader surveys may reveal whether SN 2021yfj represents an outlier or the first observed example of a subclass of stellar deaths previously hidden in observational biases. Improved modeling of mass-loss processes and nucleosynthesis will reshape how scientists interpret supernova progenitors and their explosive yields, informing our understanding of the cosmic origin story.

In essence, the discovery of SN 2021yfj’s Si/S-rich shell uncovers a hidden chapter in the lifecycle of massive stars, bridging theoretical predictions and observation. It elevates the field’s grasp of how massive stars craft intermediate-mass elements and spectacularly disperse them into space. By peeling back the layers of a star at the moment of death, astronomers reveal the intricate, layered forge that sustains the universe’s chemical diversity, redefining astrophysics and enriching humanity’s cosmic narrative.

Subject of Research: Advanced stages of nucleosynthesis and mass loss in massive stars revealed through a uniquely stripped progenitor supernova.

Article Title: Extremely stripped supernova reveals a silicon and sulfur formation site.

Article References:
Schulze, S., Gal-Yam, A., Dessart, L. et al. Extremely stripped supernova reveals a silicon and sulfur formation site. Nature 644, 634–639 (2025). https://doi.org/10.1038/s41586-025-09375-3

Image Credits: AI Generated

DOI: https://doi.org/10.1038/s41586-025-09375-3

For decades, astrophysicists have explored the connection between LGRBs and broad-lined type Ic supernovae – a subclass specifically tied to stripped-envelope Wolf–Rayet progenitors. The most energetic explosions launch highly relativistic jets that break through the stellar envelope, emitting intense gamma-ray radiation observable across cosmological distances. Yet not all jets succeed in emerging; some are “choked” within their host stars, yielding softer transients such as X-ray flashes or weaker, low-luminosity gamma-ray bursts. This gradation in jet success hints at a more nuanced interplay of progenitor properties, jet dynamics, and circumstellar environments than previously considered, motivating ongoing investigations into the continuum of relativistic outflows and transient behaviors.

Amid these efforts, a mysterious class of extragalactic fast X-ray transients has emerged, perplexing researchers because their rapid, luminous flares span timescales from mere seconds up to thousands of seconds. Their origins remain enigmatic, straddling theoretical models of jet physics and explosion mechanisms. Recent observations have failed to fit these phenomena neatly into existing schemes describing LGRBs or standard X-ray flashes, prompting proposals of alternative physical channels or progenitor conditions responsible for such transients.

A breakthrough in this field has now come with the detection of a notably bright X-ray transient, designated EP240414a, by the Einstein Probe—an advanced space observatory specialized in all-sky monitoring of X-ray emissions. Significantly, this transient coincides spatially and temporally with the type Ic broad-lined supernova SN 2024gsa, located at a cosmological redshift of 0.401. The association between the transient and the supernova provides a critical observational cornerstone to deepen our understanding of the diverse explosion scenarios linked to massive star deaths.

EP240414a’s X-ray light curve reveals an energy spectrum sharply distinct from classical LGRBs or their softer, low-luminosity cousins. The emission is extremely soft, peaking at energies less than 1.3 keV, which situates it in a spectral regime atypical for known relativistic jet-powered explosions. This softness coupled with the transient’s evolution rules out conventional high-energy jet breakout models, suggesting either substantially different jet properties or additional environmental interactions influencing the emission characteristics.

In a coordinated multiwavelength campaign following the initial X-ray detection, astronomers employed optical and radio telescopes to scrutinize the aftermath of the explosion. Observations uncovered the presence of a weak relativistic jet interacting with an extended circumstellar shell enveloping the progenitor star. This scenario contrasts with classical LGRBs where ultra-powerful jets penetrate the star’s envelope, but instead evokes a picture of a less powerful engine driving a successful, albeit relatively weak relativistic outflow that energizes the surrounding material.

The progenitor star implicated in SN 2024gsa and its transient, EP240414a, is believed to be a Wolf–Rayet star with notably reduced core angular momentum compared to traditional LGRB progenitors. This deficiency in rotation could account for the jet’s diminished power and the resulting observational signatures. Stellar rotation is a crucial parameter in magnetorotational core-collapse models that generate the conditions necessary for ultra-relativistic jets, meaning that even moderate variations can dramatically alter the explosion’s nature and the observable transient’s characteristics.

The supernova itself was located on the outskirts of a massive galaxy, a position suggestive of progenitor formation and evolution pathways differing from those in the star-forming regions typically producing classical LGRB progenitors. Environmental factors such as metallicity, binarity, and stellar feedback may have influenced the evolution of this Wolf–Rayet star and its final collapse. Thus, EP240414a and its supernova challenge astronomers to reconsider the diversity of explosion engines active in the universe and their dependence on progenitor and galactic environments.

From a theoretical perspective, the existence of such weak relativistic jets with successful but diminished breakout capabilities broadens the landscape of core-collapse end states. It indicates that there may be a continuum of jet powers governed by progenitor core spin and magnetic field properties rather than a binary classification of successful versus failed jets. Consequently, the gamma-ray and X-ray transient zoo may be more diverse and nuanced, including fast X-ray transients powered by weak jets rather than the high-luminosity events dominating the classical picture.

Moreover, the discovery of EP240414a highlights the essential role of all-sky monitoring instruments like the Einstein Probe in uncovering new transient populations. The ability to detect soft X-ray transients and coordinate multiwavelength follow-up observations is crucial for piecing together the complex interplay of jet physics, explosion dynamics, and circumstellar interactions. Such instruments open new discovery space by capturing events that would otherwise escape detection due to their intermediate luminosities and unusual spectral properties.

In terms of astrophysical implications, understanding weak relativistic jets bears significance beyond stellar death. These jets may contribute to cosmic ray acceleration, enrichment of the interstellar medium, and feedback processes that regulate star formation. Their observed interactions with circumstellar shells also shed light on mass-loss histories of massive stars, an area critical for reconstructing the final stages of stellar evolution.

Furthermore, the connection between weak jets and properties of progenitor angular momentum poses stringent tests for models of angular momentum transport and loss in massive stars. It underscores the importance of magnetohydrodynamic simulations and stellar evolution calculations that incorporate rotation, magnetic fields, and binary interactions to predict explosion outcomes and transient classifications accurately.

EP240414a thereby serves as a crucial piece in the puzzle of massive star explosions, opening pathways for future surveys to identify similar weak relativistic jet events. With improved observational capabilities, this new class of transients may become critical benchmarks for understanding jet launching mechanisms, progenitor diversity, and explosion energetics.

Scientifically, these findings demonstrate how nuanced the classification of cosmic transients has become, signaling a shift from broad categorizations toward a multidimensional parameter space capturing variations in jet power, progenitor structure, and environmental context. The binary distinction of LGRBs and failed jets is softened by discoveries like EP240414a, encouraging refinement of theoretical frameworks to incorporate intermediate cases.

Looking ahead, the synergy between transient detection facilities, wide-band follow-ups, and theoretical advances will illuminate whether weak relativistic jets are common endpoints for a significant fraction of Wolf–Rayet stars. Such understanding may bridge the gap between high-energy astrophysics, stellar evolution, and cosmology, enriching our knowledge of how massive stars influence and illuminate the universe.

The study of EP240414a and SN 2024gsa also exemplifies how serendipitous discoveries can reshape astrophysical paradigms. It is an invitation to remain vigilant for unconventional signatures that challenge current models and expand the landscape of known cosmic explosions. In this spirit, continued investment in sensitive all-sky X-ray monitors, rapid-response multiwavelength instrumentation, and theoretical modeling will drive the next leaps in revealing the lifecycle of the most massive stars.

In conclusion, the discovery of the fast X-ray transient EP240414a associated with the type Ic-BL supernova SN 2024gsa reveals a hidden population of Wolf–Rayet star explosions powered by weak yet successful relativistic jets. This new class of transients with softer X-ray spectra and intermediate jet powers challenges the classical understanding of LGRB progenitors and explosion mechanisms and highlights the complex interdependence of progenitor core rotation, jet dynamics, and circumstellar environments. Such advances promise to deepen our grasp of the most violent stellar deaths in the cosmos while unveiling new astrophysical processes shaping the universe.

Subject of Research:
Fast X-ray transients and weak relativistic jets associated with broad-lined type Ic supernovae originating from Wolf–Rayet stars.

Article Title:
A fast X-ray transient from a weak relativistic jet associated with a type Ic-BL supernova.

Article References:
Sun, H., Li, WX., Liu, LD. et al. A fast X-ray transient from a weak relativistic jet associated with a type Ic-BL supernova. Nat Astron (2025). https://doi.org/10.1038/s41550-025-02571-1

Image Credits:
AI Generated