The intellectual foundations of the subject were laid in the 1960s and 1970s. Salpeter and, independently, Zeldovich and Novikov recognized that radiation pressure from an accreting nucleus could push surrounding gas outward, limiting growth and luminosity. The balance between outward radiation force and inward gravity defines the Eddington luminosity, which scales linearly with mass. In 1973, Shakura and Sunyaev supplied the master framework that still governs the field: at high accretion rates, radiation pressure inflates the geometrically thin disc until, beyond a characteristic radius, the disc half-thickness becomes comparable to its radius. Inside this so-called spherisation radius, the disc can no longer hold itself together, and matter is blown off its surface in a powerful, radiation-driven wind. The self-regulating mechanism proposed by Shakura and Sunyaev, in which the outflow carries away excess mass while photons are effectively trapped in the optically thick inflow, remains the backbone of every modern super-critical accretion model.

Two further theoretical innovations shaped the modern understanding. First, Begelman showed in 1979 that in optically thick flows, photons become trapped and are advected inward with the gas, so that much of the accretion energy is swallowed by the black hole rather than radiated. Second, Abramowicz, Czerny, Lasota and colleagues developed the slim disc model in 1988, a solution in which radial advection of heat is an unavoidable cooling channel and the flow remains thermally stable even far above the Eddington rate. The signature prediction is subtle but profound: because of photon trapping and wind-driven mass loss, the bolometric luminosity grows only logarithmically with accretion rate, roughly as the Eddington luminosity multiplied by one plus the natural logarithm of the dimensionless accretion rate. A black hole fed a thousand times its Eddington supply therefore shines only a few times brighter than the classical limit, while the rest of the inflowing mass is flung back into space through the wind.

The review also traces the stranger corners of the theory, including the famous Polish doughnut, an elegant three-dimensional analytic solution for very high mass accretion rates in which a thick, low-viscosity, optically opaque torus forms around the hole. Such tori concentrate their emission into a narrow polar funnel and can geometrically collimate radiation and jets, an idea originally proposed by Lynden-Bell in 1978 to explain both active galactic nuclei and the enigmatic Galactic source SS 433. Later work showed that strong advective cooling thins and dims these doughnuts considerably, and the modern consensus favors a hybrid picture: a radiation-pressure-dominated, advective inner disc sheathed in an optically thick wind, with mass loss concentrated near the spherisation radius. For magnetized neutron stars, the auto-regulation picture acquires an extra layer, since the magnetosphere truncates the inner disc, and the maximum accretion rate onto the stellar surface depends on the magnetic dipole moment as well as the Eddington luminosity, scaling with the dipole moment to the four-ninths power.

The observational revolution arrived with the recognition that ULXs in nearby galaxies are genuine super-Eddington accretors rather than hidden intermediate-mass black holes. The decisive twist came when pulsations were discovered in several ULXs, proving that at least some of these extreme sources are powered by neutron stars with solid surfaces and strong magnetic fields. Accommodating such sources demanded new physics: magnetospheric truncation of the disc, columnar accretion along field lines onto the magnetic poles, and a reduced effective scattering cross-section in strong magnetic fields that allows the star to radiate far above the canonical Eddington limit. Accretion columns form above the surface, shocks settle within them, and radiation escapes largely through the sides of the column in a fan-beamed pattern rather than as pencil beams, consistent with the moderate pulse fractions observed in ultraluminous X-ray pulsars.

Numerical simulation has transformed the field from a collection of one-dimensional analytic models into a genuinely multi-dimensional science. Pioneering two-dimensional radiation-hydrodynamic simulations by Ohsuga and colleagues in 2005 demonstrated self-consistently how a geometrically and optically thick, radiation-pressure-dominated disc forms and launches outflows without assuming any disc configuration in advance. Photons are visibly trapped and dragged into the hole with the gas, while radiatively driven winds emerge from the disc surface with mildly collimated radiation along the rotation axis. Subsequent radiation-magnetohydrodynamic simulations removed the artificial alpha-viscosity prescription and incorporated magnetic turbulence directly, revealing radiatively accelerated, magnetically collimated jets. General relativistic versions of these calculations showed that black hole spin dramatically raises the energy conversion efficiency, from roughly five percent for a non-spinning hole to about thirty-three percent for a spin parameter of 0.9, and up to 140 percent in the magnetically arrested disc regime, where the Blandford-Znajek mechanism extracts rotational energy to power powerful jets. Simulations in this state even suggest that super-Eddington discs in a magnetically arrested configuration can spin their black holes down over time.

The simulations also predict structural features that observations can test. Disk winds fragment into clumpy gas clouds through Rayleigh-Taylor instabilities, and such clumpy, structured winds may explain both the X-ray variability of ULXs and the multiple absorption lines recently detected by the XRISM satellite in ultrafast outflows from a distant quasar. Large-domain simulations show that outflows are launched across the entire region within the photon trapping radius, with the highest mass-loss rates occurring not closest to the black hole but somewhat farther out, and with failed outflows that stall and fall back near the trapping radius. The outflow mechanical power inferred for ULXs, in the range of 10^39 to 10^41 erg per second, is comfortably consistent with the energetics of the vast bubble nebulae inflated around some of these sources. Even the puzzling X-ray weakness of the so-called Little Red Dots in the early Universe has been explained using simulation-based spectra of mildly super-Eddington, slowly spinning black holes viewed at moderate inclinations.

Observationally, the broadband X-ray spectra of ULXs now split naturally into components that map onto the theory: a soft, outflow-modified disc component peaking near the spherisation radius that violates the standard luminosity-temperature relation expected of thin discs, and a harder component that is comparatively insensitive to accretion rate. Because the thick wind obscures the innermost regions from most viewing angles, the apparent luminosity of a super-critical source depends strongly on inclination, giving rise to a geometric unification model in which face-on systems appear ultraluminous while edge-on ones look comparatively modest or even supersoft. Resonant absorption lines resolved in high-resolution X-ray spectra have confirmed powerful winds across the ULX population, and related ultrafast outflows have now been detected in tidal disruption events and quasi-periodic eruption sources, suggesting a common physics spanning eight orders of magnitude in accretor mass. Polarization measurements with IXPE have even revealed the predicted funnel geometry in the Galactic source Cygnus X-3, while NuSTAR phase-resolved spectroscopy showed that SS 433, viewed more face-on, would radiate at super-Eddington levels.

Substantial puzzles remain. Whether the inner flow is magnetically arrested or not, how much mass actually reaches the compact object, and how advection competes with wind-driven mass loss all remain contested, and different feeding prescriptions in simulations yield divergent answers. Timing features such as quasi-periodic oscillations in ULXs, possibly produced by Lense-Thirring precession of the tilted super-Eddington disc, still lack definitive confirmation, and recent work proposes that quasi-periodic eruptions in galactic nuclei may be the high-mass cousins of the same precessing flows. The path forward, the review argues, lies in the tight coupling of longer-duration general relativistic radiation-magnetohydrodynamic simulations running on modern GPU-based codes with the next generation of instruments: XRISM’s high-resolution spectroscopy of structured winds, the Vera Rubin Observatory’s harvest of tidal disruption events, NewAthena’s sensitivity to faint pulsating neutron star ULXs, and the SKA’s ability to chart the jets. Together these efforts promise to convert super-Eddington accretion from a beautiful theoretical curiosity into a precision tool for understanding how black holes grow, how neutron stars survive impossible feeding rates, and how radiation-hungry monsters across the cosmos regulate themselves.

Subject of Research: Super-Eddington accretion onto stellar-mass black holes and neutron stars

Article Title: Super-Critical Accretion onto Stellar Mass Black Holes and Neutron Stars: A Short Review

Article References: Middleton, M., Lipunova, G., Ohsuga, K., & Abramowicz, M. (2026). Super-Critical Accretion onto Stellar Mass Black Holes and Neutron Stars: A Short Review. Space Science Reviews, 222(6), Article 71. https://doi.org/10.1007/s11214-026-01318-2

Image Credits: AI Generated

DOI: 10.1007/s11214-026-01318-2

Keywords: super-Eddington accretion, Eddington limit, ultraluminous X-ray sources, black holes, neutron stars, accretion discs, disk winds, photon trapping, slim disc model, GR-RMHD simulations, ULX pulsars, Space Science Reviews

The research, carried out by Khushbu Zala and Sreemoyee Sarkar of SVKM’s NMIMS University in Mumbai, focuses on a phenomenon known as the r-mode instability. R-modes are non-radial oscillation modes in a rotating star that couple to gravity: as the stellar fluid sloshes back and forth, it emits gravitational radiation that carries away angular momentum. In the absence of any counteracting effect, these modes would grow without bound, spinning the star down dramatically in a runaway process. Yet astronomers observe that many millisecond pulsars—stars rotating hundreds of times per second—remain remarkably stable. Something inside them must be damping the oscillations, and the leading candidates are the two forms of viscosity that govern how dense matter dissipates energy: shear viscosity, which resists the sliding of adjacent fluid layers, and bulk viscosity, which dissipates energy when the fluid is periodically compressed and expanded.

What makes the new work distinctive is its treatment of a star whose core contains a mixed phase, a region where hadronic and quark matter coexist in thermodynamic equilibrium. Below a transition density of about 2.23 times nuclear saturation density, the outer layer is modeled as ordinary baryonic matter described by the relativistic mean-field equation of state known as DDME2. Above that threshold, the researchers assume a hybrid interior in which quarks and hadrons mingle. The problem, as the authors emphasize, is inherently ill-posed: the equation of state at supranuclear densities is poorly constrained, quantum chromodynamics cannot yet be solved reliably in the non-perturbative regime relevant to stellar cores, and the transport coefficients—the viscosities that control damping—are even harder to pin down from first principles.

To cut through these uncertainties, the team turned to Bayesian inference, a statistical framework that updates prior physical knowledge with observational data to produce posterior distributions for unknown parameters. The prior distributions for the key dimensionless coefficients—labeled S-tilde, V-tilde, W-tilde and J-tilde, which encode the shear and bulk viscous response of the mixed phase along with the star’s equilibrium structure—were drawn from existing calculations for neutron stars, strange stars and hybrid stars. The likelihood function combined two independent observational constraints: mass-radius measurements from NASA’s Neutron Star Interior Composition Explorer, or NICER, mission, and the spin-frequency and temperature observations of neutron stars in low-mass X-ray binaries. Using the UltraNest nested sampling algorithm, the researchers explored the parameter space efficiently, focusing on regions of high likelihood and estimating the Bayesian evidence for their model.

The formalism at the heart of the analysis describes how the amplitude of an r-mode evolves in time as an exponential whose decay constant is set by the competition between three timescales: gravitational radiation, which drives the instability, and shear and bulk viscosity, which suppress it. At low temperatures, shear viscosity dominates and stabilizes the star; at high temperatures, bulk viscosity takes over. Between these regimes lies a window where damping is least effective and the star is most vulnerable. The minimum of the instability curve—the lowest spin frequency at which the mode can grow—is therefore an exquisitely sensitive probe of the microphysics inside the star, including the equation of state and the weak-interaction processes, such as the direct Urca reaction, that generate bulk viscosity.

Applying this framework to two hybrid star configurations of 1.5 and 1.75 solar masses, the team obtained concrete estimates for the dissipation timescales. The shear viscous damping time came out as approximately 4.99 times ten to the eight, multiplied by the temperature to the five-thirds power, in seconds, while the bulk viscous timescale follows a more complex dependence on both temperature and spin, scaling inversely with the square of the angular velocity. From the inferred coefficients, the researchers calculated the minima of the instability curves: the critical angular velocity reaches its lowest value of about 451.87 hertz at a temperature of 0.259 megaelectronvolts for the 1.5 solar mass star, and 517.47 hertz at 0.234 megaelectronvolts for the 1.75 solar mass star. Normalized to the Kepler frequency, the maximum spin rate a star can sustain before mass shedding, these minima correspond to ratios of roughly 0.069 and 0.071.

Crucially, the resulting instability window does more than produce elegant numbers—it matches what astronomers actually see. When the team compared their inferred instability curves with the observed spin frequencies and temperatures of real millisecond pulsars, they found that the enhanced viscous dissipation from the mixed hadron-quark phase provides sufficient damping to explain the stability of several well-known objects. Among them are the accreting low-mass X-ray binary sources XTE J0929-314 and XTE J1807-294, and the radio millisecond pulsars J0437-4715 and J2124-3358. These stars all spin faster than 100 hertz, placing them squarely in the frequency range where r-mode physics matters, and all of them sit safely outside the region where the instability would grow—precisely as the two-layer hybrid model predicts.

The posterior distributions themselves carried informative structure. The corner plots of the inferred parameters showed moderate correlations among the shear and bulk viscous coefficients, reflecting the coupled role of the two viscosities in setting the instability boundary, while correlations involving the equilibrium parameter remained comparatively weak. Notably, the posterior contours for the more massive 1.74 solar mass configuration were narrower and more tightly localized than those for the lighter star, indicating that frequency-temperature observations constrain the dense-core dissipation properties more stringently in heavier hybrid stars. The team also found that the position of the instability minimum is remarkably robust: it barely shifts when the equilibrium and viscous parameters vary across their full credible intervals, suggesting that the result is not an artifact of statistical noise.

The broader significance of the work lies in its demonstration that r-mode phenomenology, combined with modern statistical inference, can serve as a practical observational tool for probing phase transitions in ultra-dense matter. Because pulsar rotational frequencies and their time derivatives are among the most precisely measured quantities in all of astrophysics, and because NICER continues to deliver mass-radius constraints, the approach links macroscopic observables directly to microscopic transport physics. If the inferred viscous properties of the mixed phase continue to align with observations, it would strengthen the case that some neutron stars genuinely harbor deconfined quark matter in their cores—a question that has remained open since the earliest theoretical speculations about quark stars.

The authors outline several directions for extending the framework. A more realistic three-layer stellar model could better capture the stratification of a hybrid star’s interior, and generalizing the formalism to derive quantitative constraints on the shear and bulk viscosities of each individual layer would sharpen the physical picture. Perhaps most ambitiously, they aim to characterize the nature of the hadron-quark phase transition itself—determining whether it is first order, second order, or a smooth crossover—by performing statistical inference on transport coefficients constrained by the gravitational-wave signatures that r-mode oscillations generate. As gravitational-wave detectors grow more sensitive, the faint hum of a wobbling hybrid star may one day confirm what this Bayesian analysis already hints at: that the universe’s most extreme matter hides its secrets in the way it dissipates motion.

Subject of Research: Bayesian inference of viscous dissipation timescales governing r-mode instability in hybrid stars with hadron-quark mixed phases

Article Title: Modelling dissipative dynamics of r-mode instability in hybrid stars

Article References: Zala, K., & Sarkar, S. (2026). Modelling dissipative dynamics of r-mode instability in hybrid stars. The European Physical Journal C, 86(9), Article 1064. https://doi.org/10.1140/epjc/s10052-026-16199-6

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16199-6

Keywords: hybrid stars, r-mode instability, neutron stars, bulk viscosity, shear viscosity, Bayesian inference, quark matter, millisecond pulsars, NICER, low-mass X-ray binaries, gravitational waves, dense matter

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

The central challenge the authors tackle is deceptively simple to state: neutron stars are born hot and magnetized, and every observable we measure — their X-ray glow, their spin periods, their sudden outbursts — is shaped by how heat and magnetic field decay and interact over time. Solving this problem requires numerically integrating two coupled families of equations: the heat transfer equation governing the cooling of the stellar interior, and the induction equation describing the evolution of the magnetic field. Both depend sensitively on microphysics such as thermal and electrical conductivities, neutrino emission rates and the behavior of superfluid and superconducting particles at densities found nowhere else in nature.

The cooling story begins dramatically. A newborn neutron star starts at more than ten billion kelvin, transparent only to neutrinos after roughly a minute, and shrinks from about a hundred kilometers to its final radius of ten to fourteen kilometers. Within the first year, its core — containing about ninety-nine percent of the mass — becomes nearly isothermal thanks to enormous thermal conductivity. Long-term cooling then proceeds over hundreds of thousands of years, with temperature gradients confined to the thin crust and envelope. Crucially, the absence of internal convection means that neutron stars cannot operate the kind of self-sustained dynamos that generate magnetic fields in planets and ordinary stars, so their fields must originate elsewhere and be preserved or dissipated over time.

One of the most consequential recent developments highlighted in the review concerns enhanced neutrino cooling. For decades, theorists have debated whether neutron stars cool via ‘minimal’ channels involving modified Urca processes or whether some stars undergo dramatically faster ‘direct Urca’ cooling, possibly triggered by exotic matter such as hyperons or quarks. The Vela pulsar’s anomalously low temperature hinted at such processes, but definitive evidence was elusive. A recent analysis of three young, nearby and extremely cold neutron stars now shows that their properties require enhanced cooling, allowing researchers for the first time to place meaningful constraints on the dense-matter equation of state directly from thermal observations.

The same cooling physics reaches into particle physics beyond the Standard Model. Because hypothetical light particles such as axions would drain energy from a neutron star’s interior, cooling measurements can bound their properties. Studies cited in the review use the thermal behavior of young neutron stars to constrain the QCD axion mass, and more recent work incorporating envelope structure has tightened those limits further. Meanwhile, the famous case of the neutron star in the Cassiopeia A supernova remnant — once thought to show a rapid temperature decline signaling the onset of neutron superfluidity — illustrates the difficulty of the enterprise. The inferred cooling rate has shrunk from about four percent per decade to between roughly 1.6 and 2.2 percent as data and detector calibrations improved, and recent critical reanalyses question whether the decline is robust at all.

On the magnetic side, the review lays out the three dominant interior processes. Ohmic dissipation, the slow diffusion of currents through the highly conductive crustal lattice, acts on timescales of hundreds of thousands to millions of years. The Hall drift, a nonlinear advection of magnetic field by the electron fluid, becomes dominant in magnetars and can transfer magnetic energy from large scales to small scales, accelerating dissipation. Ambipolar diffusion, the coupled drift of charged particles relative to the neutron superfluid, may dominate in the cores of strongly magnetized stars during their first hundred thousand years, though its exact role in superfluid and superconducting conditions remains actively contested. Newer ingredients — including the chiral magnetic effect, in which a tiny imbalance between left- and right-handed electrons can amplify fields, and crustal failure mechanisms ranging from brittle starquakes to plastic flow — are given detailed treatment, with the authors cautioning that some widely used failure criteria may overestimate how often the crust actually cracks.

Numerical methodology occupies a substantial portion of the review, and for good reason: the equations are stiff, nonlinear and span many orders of magnitude in the relevant coefficients. Older simulations relied on spectral methods elegant in accuracy but fragile near the discontinuities that the Hall term inevitably produces. Recent codes have shifted toward finite-difference and finite-volume schemes borrowed from high-resolution shock-capturing techniques in computational fluid dynamics. The new three-dimensional code MATINS, developed by members of the review team, goes further by adopting a cubed-sphere grid that sidesteps the numerical pathologies of the polar axis, while incorporating realistic microphysics throughout the crust. Benchmark tests against analytical solutions confirm accuracy at the sub-percent level.

What happens when such codes are initialized with realistic, turbulent magnetic fields inherited from the proto-neutron star dynamo phase is striking. Fully three-dimensional simulations starting from complex fields dominated by small and intermediate scales show that the tangled structure persists for hundreds of thousands of years, with the Hall term continuously feeding energy from large scales while Ohmic dissipation erodes the small scales. These turbulent initial conditions naturally reproduce the observed properties of central compact objects and the so-called low-field magnetars — stars that erupt in magnetar-like bursts despite possessing relatively weak large-scale dipoles. However, the simulations still struggle to generate the ultra-strong, dominant dipole that classical magnetars require, leaving the origin of magnetar-strength dipolar fields an open question. One provocative possibility explored in recent work is that the chiral magnetic effect could grow the dipolar component to magnetar strengths within fifty to a hundred years of birth, driven by an extraordinarily small but persistent chiral imbalance sustained by the star’s magnetic helicity.

The final piece of the puzzle is rotation. Spin period and its derivative are the most precisely measured neutron star observables, and they encode the electromagnetic torque exerted by the magnetosphere. The review emphasizes that the familiar textbook formula for inferring magnetic fields from timing data, based on a vacuum dipole, systematically overestimates the true field strength. Plasma-filled force-free magnetospheres exert torques even on aligned rotators, and general relativistic effects near the stellar surface further amplify the spin-down luminosity. A properly relativistic formula yields dipolar field estimates roughly a factor of four lower than the classical expression, potentially reshaping how the entire neutron star population is classified — including how many objects truly qualify as magnetars.

As instruments such as the Square Kilometre Array Observatory prepare to detect many thousands of new pulsars and X-ray observatories continue to monitor magnetar outbursts, the theoretical framework assembled in this review provides the essential bridge between microphysics and observation. The authors’ outlook is clear: the future lies in fully three-dimensional, self-consistent magneto-thermal simulations coupled to dynamical magnetospheres, with machine-learning solvers such as physics-informed neural networks emerging as a computationally efficient complement to traditional schemes. What emerges from the review as a whole is a picture of neutron stars not as static lighthouses but as dynamic, evolving systems in which magnetism, heat and rotation are inseparably entangled — a picture that is only now becoming possible to compute.

Beyond the theoretical machinery, the review situates itself within a rapidly growing observational census. Nearly four thousand rotation-powered radio pulsars are now catalogued, a number set to multiply as the Square Kilometre Array Observatory comes online, while high-energy instruments have identified several hundred gamma-ray pulsars and roughly a hundred X-ray sources. Only a few dozen neutron stars show the soft, thermal surface emission in X-rays that directly probes interior cooling, making each detection disproportionately valuable for constraining the underlying microphysics.

Magnetars, the extreme end of this population, remain scarce — a few dozen confirmed candidates — but their energetic output is extraordinary. Their persistent X-ray luminosities of 10^33 to 10^35 erg per second routinely exceed their rotational energy loss rates by orders of magnitude, implying that magnetic energy, not spin-down, powers their emission. The most violent events, giant flares, release up to 10^46 erg in under a second; only three have ever been recorded. This magnetic reservoir interpretation, first advanced by Thompson and Duncan in the mid-1990s, now underpins most theoretical modeling of magnetar activity.

Perhaps the most significant conceptual shift the review documents is the blurring of historical subclass boundaries. Low-field magnetars — nominally ordinary pulsars that occasionally erupt in bursts — and some central compact objects demonstrate that a strong surface dipole is not the decisive ingredient for magnetar-like behavior. Instead, attention has turned to the hidden internal field architecture: how magnetic energy is partitioned between toroidal and poloidal components and across spatial scales, where electrical currents reside, and how helicity is transported outward into the magnetosphere. These questions set the agenda for the three-dimensional simulations and future code development efforts the review outlines.

Subject of Research: Magneto-thermal and rotational evolution models of isolated neutron stars

Article Title: Magnetic, thermal and rotational evolution of isolated neutron stars

Article References: Pons, J. A., Dehman, C., & Viganò, D. (2026). Magnetic, thermal and rotational evolution of isolated neutron stars. Living Reviews in Computational Astrophysics, 12(1), Article 5. https://doi.org/10.1007/s41115-026-00028-4

Image Credits: AI Generated

DOI: 10.1007/s41115-026-00028-4

Keywords: neutron stars, magnetars, pulsars, magnetic field decay, neutron star cooling, neutrino emission, Hall drift, Ohmic dissipation, ambipolar diffusion, magneto-thermal simulations, dense matter equation of state, numerical astrophysics

The sheer density of neutron stars renders them some of the most extreme objects known to science. Imagine an object with a mass greater than our Sun packed into a sphere no larger than a city. This incredible compression leads to physics far removed from our everyday experiences, where gravitational forces dominate to an extent that protons and electrons are crushed together to form neutrons. However, the conventional models describing these cosmic behemoths, while incredibly successful, may not encompass the full picture. Modern cosmological theories, particularly those attempting to unify gravity with quantum mechanics, often invoke the existence of extra spatial dimensions beyond the three we perceive. Braneworld theories, a prominent example of such frameworks, suggest that our universe might be embedded within a higher-dimensional reality, with our everyday forces confined to our three-dimensional brane. The research presented here daringly proposes that the immense gravitational pull and the exotic matter configurations within neutron stars could be sensitive to subtle influences from these hypothetical extra dimensions, specifically through a property known as brane tension.

Brane tension, in this context, refers to the inherent energy density of the brane itself. Think of it as a stretching force that holds the brane together. If our universe is a brane within a larger bulk, then this tension would be a fundamental property of our cosmic existence. The idea is that phenomena occurring on our brane, especially those involving extreme densities and energies like those found in neutron stars, might interact with or be affected by this fundamental tension. This interaction could manifest as deviations from the predictions of standard general relativity, offering a potential avenue for observational verification of these speculative, yet deeply compelling, theories about the architecture of spacetime. The intricate interplay between the immense gravity of neutron stars and the fundamental properties of our cosmic membrane could therefore provide a unique laboratory for probing the very nature of reality.

The brilliance of the research lies in its innovative approach to constraining these theoretical ideas. Instead of relying solely on abstract mathematical models, Murshid and his team have ingeniously sought to utilize observational data from actual neutron stars. By analyzing the properties of these pulsars, such as their mass, radius, and the emitted radiation, physicists can infer the internal structure and the equation of state that governs the matter within them. The equation of state describes how pressure changes with density, a critical factor in understanding the stability and behavior of neutron stars. The theoretical models that incorporate braneworld effects predict subtly different equations of state compared to those rooted in traditional four-dimensional spacetime. It is precisely these predicted differences that the researchers aimed to detect through careful analysis of observational data.

The process of constraining brane tension involves a meticulous comparison between theoretical predictions and actual astronomical observations. The researchers developed sophisticated models that incorporate the influence of brane tension on the internal structure and observable properties of neutron stars. These models predict specific correlations between the mass and radius of a neutron star, or how its surface behaves under extreme conditions. Any deviations from the predictions made by standard general relativity, when fed into these braneworld models, could then be attributed to the presence and magnitude of brane tension. It’s akin to searching for a faint whisper of a different physics regime amidst the colossal roar of a neutron star’s gravitational field, a testament to the precision of modern astrophysics.

The data used in this study likely comprises a curated collection of precise measurements from radio telescopes and X-ray observatories, focusing on neutron stars with well-determined masses and radii. These crucial parameters allow theorists to test various equations of state. For instance, if a neutron star’s observed mass and radius suggest a stiffer equation of state than predicted by standard models, this could be an indirect signal of braneworld effects. The strength of the braneworld influence, and thus the effective brane tension, would then be inferred from how well these braneworld models can reproduce the observed properties. The challenge lies in disentangling these subtle braneworld effects from other astrophysical uncertainties and systematic errors in the observations, a task demanding immense computational power and rigorous statistical analysis.

The research highlights the power of astrophysical objects like neutron stars as natural laboratories for testing the limits of our physical theories. While particle accelerators on Earth can probe energies up to a certain point, the extreme conditions within neutron stars—densities reaching nuclear saturation and gravitational fields far exceeding anything we can replicate—provide a unique opportunity to explore physics at energy scales far beyond our current experimental reach. By observing neutron stars, we are, in essence, performing experiments on the fundamental laws of nature under conditions that have not existed on Earth since the earliest moments of the universe. This paper represents a significant step in leveraging these cosmic laboratories to probe the exotic realms of extra dimensions and brane theories.

The implications of finding a non-zero brane tension could be profound. It would provide strong empirical support for braneworld scenarios, suggesting that our universe is indeed embedded in a richer, higher-dimensional landscape. This discovery would have far-reaching consequences for our understanding of gravity, cosmology, and potentially even the origin of mass itself. It could offer new insights into dark matter and dark energy, two of the most significant mysteries in modern cosmology, by providing a new framework within which to formulate theoretical explanations. The notion that the properties of everyday objects are influenced by the very structure of spacetime is a concept that sparks the imagination and pushes the boundaries of scientific inquiry ever further.

The methodology employed by Murshid, Moneer, Zotos, and their collaborators involves the meticulous construction and refinement of theoretical models that describe neutron stars within the context of braneworld scenarios. These models incorporate the effects of the extra dimensions and the inherent tension of our brane on the equilibrium structure and the dynamical behavior of neutron star matter. By considering various possible values of brane tension, the researchers can predict how the mass-radius relationship of neutron stars, or their vibrational modes, might deviate from predictions made by standard general relativity. These precisely calculated deviations are then compared with the actual observational data, allowing the team to place stringent constraints on the allowed values of brane tension.

The paper’s findings offer tangible results in the form of numerical constraints on the magnitude of this theoretical brane tension. While the exact values remain under intense scrutiny and may evolve with further data, the study asserts that observational data from neutron stars can indeed limit the possible range for this fundamental cosmic parameter. This is a critical achievement because it moves the concept of braneworlds from purely theoretical speculation towards experimentally verifiable physics. Such constraints are vital for guiding future theoretical developments and for identifying which braneworld models are most consistent with our observed universe, marking a significant step in empirical physics.

The potential for these findings to be a “viral” scientific discovery stems from their ability to capture the public’s imagination. The idea that our universe is a “brane” in a larger reality, and that the exotic objects like neutron stars can reveal secrets about this hidden architecture, is a narrative that resonates deeply. It taps into humanity’s innate curiosity about the unknown and our place in the cosmos. If these results hold up to further scrutiny and are corroborated by other studies, they could usher in a new era of cosmological and astrophysical research, inspiring widespread public interest and potentially leading to a re-evaluation of our fundamental understanding of reality.

Furthermore, the research is not a static conclusion but rather an invitation for more extensive investigation. The team emphasizes the need for more precise observational data and the continued development of sophisticated theoretical models to further refine the constraints on brane tension. As new generations of telescopes and detectors come online, promising unprecedented accuracy in astronomical measurements, the opportunities to test these braneworld scenarios will only increase. This ongoing interplay between theory and observation is the engine that drives scientific progress, and this work has provided a powerful new direction for that engine to pursue.

The visual representation provided with the research, depicting a neutron star with an ethereal glow, hints at the profound nature of the forces at play. While the image itself might be an artistic rendering, it serves as a powerful reminder of the vast cosmic phenomena that scientists are striving to understand. The immense gravitational fields and the extreme densities within neutron stars are not just abstract concepts; they are tangible, observable realities that hold clues to the deepest mysteries of the universe, including its potential higher dimensions and the fundamental tension of its very fabric. The research signifies a triumph of human curiosity and ingenuity, pushing the boundaries of our knowledge into the most extreme and fascinating corners of existence. The convergence of abstract theoretical physics with the raw, observable data from the cosmos has never been more compelling.

Subject of Research: Investigating the properties of neutron stars to constrain theoretical models of extra spatial dimensions, specifically focusing on the concept of brane tension in braneworld scenarios.

Article Title: Braneworld neutron stars: constraining brane tension with observational data.

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-025-14561-8

Keywords**: Neutron stars, braneworlds, brane tension, general relativity, astrophysics, cosmology, extra dimensions, equation of state, observational constraints.