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

At the heart of this investigation lies the concept of the equation of state (EoS) of dense nuclear matter. This is not merely a theoretical construct; it is the rulebook that dictates how matter behaves under immense pressure and density. For neutron stars, this equation of state is paramount in determining their maximum possible mass. Imagine trying to compress a substance indefinitely; at some point, the internal forces resisting compression will overpower the external force. For neutron stars, this internal resistance is governed by the complex interplay of nuclear forces, including the strong nuclear force that binds protons and neutrons together, and potentially more exotic phenomena like the presence of hyperons or quark matter at even higher densities. The EoS essentially maps out the pressure experienced by the matter within a neutron star as a function of its density. A stiffer EoS, meaning matter strongly resists compression, will allow for more massive neutron stars, while a softer EoS will lead to a lower maximum mass. The challenge is that the EoS is not directly observable, and its form at the densities found inside neutron stars is still a subject of intense theoretical debate.

The recent observations of PSR J0740+6620 have provided an unprecedentedly accurate measurement of its mass, placing it at an astonishing 2.14 solar masses. This is not just another data point; it is a crucial anchor for theoretical models. Finding a neutron star with such a substantial mass strongly suggests that the nuclear matter within it is remarkably incompressible at these extreme densities, hinting at a “stiff” equation of state. If neutron stars could only exist up to a certain mass, and we then observe one that surpasses the previously accepted theoretical limits, it forces a reevaluation of our understanding of nuclear interactions at these densities. This observation serves as a powerful constraint, ruling out many theoretical EoS models that predict a lower maximum mass. The sheer existence of such massive neutron stars, packed into such compact volumes, is a testament to the extraordinary strength and complexity of the forces at play beyond the realm of everyday experience.

Adding another layer of complexity and observational power to this puzzle is the detection of gravitational waves, particularly the event GW190814. This event involved the merger of two compact objects, one of which was definitively identified as a neutron star with a mass around 1.4 solar masses. The other object’s mass, however, was a tantalizing enigma, falling into a mass gap between typical neutron stars and known black holes, estimated to be around 23 solar masses. While the precise nature of this companion is still debated – it could be an extremely massive neutron star or a low-mass black hole – the gravitational wave signal provides invaluable information about the inspiral and merger process. The way these objects orbit each other and distort spacetime as they merge leaves an imprint on the gravitational waves that can be used to infer their masses and radii. The properties of the neutron star involved, particularly its tidal deformability during the inspiral, as imprinted in the gravitational waveform, offer a complementary probe of the nuclear EoS.

The European Physical Journal C study by Zhou and colleagues meticulously stitches together these observational threads with theoretical calculations. They explore a range of modern nuclear EoS models, each representing different assumptions about the behavior of nuclear matter under extreme conditions. These models are then tested against the stringent constraints provided by the mass of PSR J0740+6620 and the information gleaned from the GW190814 merger. The interplay between these two distinct observational messengers is critical. While the mass of PSR J0740+6620 directly probes the maximum possible mass, and thus the stiffness of the EoS at its upper limit, the gravitational wave data from GW190814, particularly concerning tidal effects during the inspiral, provides information about the EoS at slightly lower, but still extremely high, densities.

The correlations explored in the paper highlight a profound link: the maximum mass of a neutron star is not an isolated parameter. It is intrinsically tied to other fundamental nuclear matter properties, such as the nuclear saturation density, the symmetry energy, and the pressure at supranuclear densities. The symmetry energy, in particular, describes how the energy of nuclear matter changes with the ratio of neutrons to protons. This is a key ingredient in nuclear physics, and its value at high densities has significant consequences for neutron star structure and maximum mass. A higher symmetry energy generally leads to a stiffer EoS and thus potentially more massive neutron stars. The study investigates how different theoretical assumptions about these properties translate into predictions for the maximum mass and tidal deformability, and then quantitatively assesses how well these predictions match the observed data.

The findings of Zhou et al. are poised to send ripples through the astrophysics community. By analyzing the detailed correlations between maximum mass and various nuclear matter properties, and critically evaluating them against the precise constraints from PSR J0740+6620 and GW190814, the researchers have managed to narrow down the viable parameter space for theoretical nuclear EoS models. This is a significant step forward in our understanding of the fundamental forces that govern matter at densities far exceeding those found in terrestrial laboratories or even within atomic nuclei. The study provides compelling evidence that favors certain nuclear physics models over others, bringing us closer to a unified and accurate description of ultradense matter. This rigorous comparison of theory with observation is the engine of scientific progress, turning abstract theories into physically grounded realities.

This research also has profound implications for our understanding of potential exotic matter phases within neutron stars. At densities exceeding approximately twice the nuclear saturation density, it is theoretically possible that neutrons and protons could “dissolve” into a soup of quarks and gluons, forming quark matter or strange quark matter. The presence and properties of such phases would dramatically alter the equation of state, potentially leading to a softening that could limit the maximum neutron star mass. The observational constraints from PSR J0740+6620 and GW190814 are crucial in determining whether these exotic phases are permitted under realistic astrophysical conditions. If the maximum mass is indeed as high as indicated, it suggests that if quark matter exists, it either does not significantly soften the EoS or it forms at even higher densities than previously thought, or perhaps the neutron star is masquerading as something even stranger.

The implications of this work extend beyond neutron stars themselves, potentially shedding light on the enigmatic nature of dark matter. While not directly addressed in this specific study, the fundamental properties of matter at extreme densities are deeply intertwined with our understanding of the universe’s composition. Theories that seek to explain dark matter often involve new particles and interactions that could manifest themselves in the structure and evolution of dense objects like neutron stars. By refining our understanding of the EoS and the limits of nuclear matter, this research helps to constrain broader cosmological models and the fundamental physics that underpins them. It’s a testament to how advancements in one field of physics can illuminate seemingly unrelated areas of inquiry.

The precision of modern astrophysical observations is a key driver of these advances. The ability to measure neutron star masses with such accuracy, and to detect gravitational waves from their mergers, provides a level of detail previously unimaginable. These observations act as empirical lighthouses, guiding theorists through the vast landscape of possible models and theories. The synergy between cutting-edge observational facilities like advanced radio telescopes and gravitational wave detectors, and sophisticated theoretical frameworks, is what allows us to probe the universe’s most extreme phenomena and unlock its deepest secrets. The ongoing improvements in these observational capabilities promise even more exciting discoveries in the years to come.

The specific correlations examined in the paper are subtle but crucial. For instance, the study quantifies how the neutron star radius evolves with its mass, and how tidal deformability, a measure of how much an object is stretched by an external gravitational field, changes with compactness. These are all directly related to the underlying equation of state. A stiffer EoS leads to more compact, less deformable neutron stars with potentially higher maximum masses. By mapping out these relationships and comparing them to the observational data, Zhou and colleagues can effectively “sound out” the interior of neutron stars, probing densities and pressures that are otherwise inaccessible. This process of inferring fundamental properties from macroscopic behavior is a hallmark of scientific investigation.

The publication in the European Physical Journal C signifies the broad impact and acceptance of this research within the physics community. The rigorous peer-review process ensures that the methodology is sound, the calculations are accurate, and the conclusions are well-supported by the evidence. This kind of detailed theoretical work, grounded in solid observational constraints, is essential for building a reliable picture of the fundamental physics governing the universe. It’s through such diligent scientific contributions that our understanding of the cosmos is steadily advanced, moving us from speculation to well-founded knowledge.

Looking ahead, this research opens up new avenues for exploration. Future gravitational wave observations, potentially involving mergers of even more massive neutron stars or neutron star-black hole binaries, will provide even tighter constraints on the EoS. Similarly, continued observations of isolated pulsars like PSR J0740+6620, especially those with precisely measured masses and radii, will further refine our understanding of these extreme objects. The quest to fully unravel the mysteries of dense nuclear matter is far from over, but this study marks a significant milestone in our journey, bringing us closer to understanding the ultimate fate of matter in the universe and the fundamental forces that shape it. The continued interplay between theory and observation will undoubtedly lead to further paradigm shifts in our comprehension of the cosmos.

Subject of Research: The relationship between the maximum mass of neutron stars and the fundamental properties of nuclear matter, constrained by astronomical observations.

Article Title: Correlations between maximum mass of neutron stars and the nuclear matter properties and the constraints from PSR J0740+6620 and GW190814.

Article References: Zhou, M., Liu, H.M., Zheng, H. et al. Correlations between maximum mass of neutron stars and the nuclear matter properties and the constraints from PSR J0740+6620 and GW190814. Eur. Phys. J. C 85, 825 (2025). https://doi.org/10.1140/epjc/s10052-025-14557-4

Image Credits: AI Generated

DOI: https://doi.org/10.1140/epjc/s10052-025-14557-4

Keywords: Neutron stars, maximum mass, equation of state, nuclear matter, PSR J0740+6620, GW190814, gravitational waves, dense matter, nuclear physics, astrophysics