⚛️ nuclear theory

Extending multi-messenger constraints on neutron star matter through the inclusion of direct Urca cooling

This paper presents a Bayesian study that incorporates direct Urca cooling constraints into neutron star equation of state inference, revealing that while these composition-sensitive observations have a minor impact on macroscopic properties like mass and radius when combined with other data, they significantly refine our understanding of the proton fraction within neutron stars.

Original authors: M. Dio Danarianto, Guilherme Grams, Thibeau Wouters, Tim Dietrich

Published 2026-09-25

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: M. Dio Danarianto, Guilherme Grams, Thibeau Wouters, Tim Dietrich

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Deep in the cosmos, hidden within the crushing hearts of dead stars, lies a form of matter so dense that a single teaspoon would weigh as much as a mountain. This is the realm of neutron stars, the collapsed cores of massive stars that have exploded in supernovae. For decades, scientists have tried to understand the rules governing this extreme material, but the densities involved are far beyond what we can recreate in laboratories on Earth. While we can study atomic nuclei in particle accelerators, the conditions inside a neutron star are a different beast entirely, pushing matter to limits where our current theories of physics begin to fray. To solve this puzzle, researchers have turned to the stars themselves, using them as natural laboratories. By observing how these stars behave—how heavy they are, how big they are, and how they cool down over time—scientists hope to piece together the recipe for the densest matter in the universe.

A key piece of this puzzle involves how heat escapes from a neutron star. When a star is young, it is incredibly hot, but it must eventually cool down. Most stars cool slowly, like a cup of coffee left on a table. However, some neutron stars cool much faster, shedding their heat at a breakneck pace. This rapid cooling happens because of a specific process called the direct Urca reaction. Imagine the star's core as a crowded room where neutrons, protons, and electrons are packed together. In a normal star, these particles interact gently, but if the crowd becomes dense enough and the number of protons reaches a certain threshold, a rapid exchange occurs: a neutron turns into a proton and an electron, shooting out a ghostly particle called a neutrino that carries the heat away instantly. This process only kicks in if the star is heavy enough and the mix of particles is just right. If we can figure out which stars are cooling this fast, we can learn exactly how much proton-rich matter exists inside them.

In a new study, researchers set out to combine these cooling observations with other ways of measuring neutron stars to get a clearer picture of their internal composition. They used a sophisticated computer framework to build thousands of possible models for the matter inside these stars. These models were tested against three different types of real-world data. First, they looked at the theoretical limits of nuclear physics, which tell us how matter behaves at lower densities. Second, they included measurements from telescopes and gravitational wave detectors that tell us the mass and size of various neutron stars. Third, and most importantly for this study, they added the information from the rapidly cooling stars. The team specifically looked at objects like the Cassiopeia A supernova remnant and several pulsars, which are known to be cooling quickly, and compared them to a slowly cooling star to see where the transition happens.

The researchers found that adding the cooling data to their models did not drastically change the overall size or maximum weight of the stars they predicted. The models still pointed to stars with radii around 12 kilometers and maximum masses just over two times the mass of our Sun, consistent with what other observations suggest. However, the cooling data did reveal something crucial about the ingredients inside the star. The models that included the rapid cooling constraint showed that the stars must contain a higher fraction of protons in their cores than models that ignored cooling. This higher proton content is exactly what is needed to trigger the fast cooling process. In essence, the cooling observations acted as a filter, ruling out models where the core was too "neutron-heavy" and keeping only those where the proton mix was sufficient to allow the rapid heat loss.

This work highlights the power of combining different types of cosmic clues. While the cooling data alone could not pin down the exact size of the stars, it provided a unique window into their chemical makeup that mass and size measurements alone could not offer. The study suggests that the stars cooling rapidly are likely around 1.4 times the mass of the Sun, and for these stars to cool that fast, their cores must be rich in protons. The findings also show that the specific way scientists model the transition from known nuclear physics to the unknown high-density regions can influence the results, but the core conclusion about the proton fraction remains robust across different modeling choices. By weaving together the story of how these stars cool with how they weigh and how they deform, scientists are slowly closing in on the true nature of matter at its most extreme, turning the distant, silent hearts of dead stars into a readable text for the laws of physics.

Technical Summary: Extending Multi-Messenger Constraints on Neutron Star Matter through the Inclusion of Direct Urca Cooling

Problem Statement
The properties of strongly interacting matter at extreme densities remain a fundamental open question in modern physics. While terrestrial experiments and theoretical methods like chiral effective field theory ( $\chi$ EFT) provide constraints up to approximately twice nuclear saturation density ( $2n_{\text{sat}}$ ), and perturbative quantum chromodynamics applies at extremely high densities ( $\gtrsim 40n_{\text{sat}}$ ), the intermediate regime remains theoretically uncertain and difficult to probe directly. Neutron stars (NSs) offer a unique laboratory for this density region. Current multi-messenger approaches typically rely on macroscopic observables such as mass, radius, and tidal deformability. However, these observables do not directly constrain the microscopic composition of the star, specifically the proton fraction ( $Y_p$ ). The direct Urca (dUrca) process, a rapid neutrino emission mechanism, is highly sensitive to composition, occurring only when the proton fraction exceeds a specific threshold ( $Y_{\text{DU}}$ ). Integrating cooling observations, which indicate the presence or absence of dUrca, into equation of state (EoS) inference requires a self-consistent treatment of $\beta$ -equilibrium and charge neutrality to determine the proton fraction, a capability often lacking in standard macroscopic EoS analyses.

Methodology
The authors perform a Bayesian inference study to constrain the NS EoS by explicitly incorporating dUrca cooling constraints alongside traditional macroscopic and nuclear data.

Framework Extension: The study utilizes and extends the jester framework. A key modification is the implementation of a self-consistent treatment of $\beta$ -equilibrium and charge neutrality within nucleonic models. This allows for the calculation of the proton fraction ( $Y_p$ ) and the determination of the dUrca onset density ( $n_{\text{DU}}$ ) in the presence of electrons and muons, moving beyond previous approximations.
EoS Construction: Two classes of EoS are employed: a metamodel with a sound-speed extension (CSE) and extended-Skyrme energy-density functionals. The nucleonic description is valid up to a breakdown density ( $n_{\text{break}}$ ). Two scenarios for $n_{\text{break}}$ are tested: (1) a random sampling between $1n_{\text{sat}}$ and $4n_{\text{sat}}$ , and (2) a causal limit where the model breaks down when the speed of sound reaches zero or one. Above $n_{\text{break}}$ , an agnostic six-segment sound-speed extension is used up to $25n_{\text{sat}}$ .
Likelihood Construction: The total likelihood combines three components:
1. Nuclear Constraints: Data from $\chi$ EFT.
2. Astrophysical Constraints: Masses and radii from NICER observations (PSR J0030+0451, J0740+6620, etc.), heavy pulsar masses (PSR J1614-2230, J0348+0432), and tidal deformability from GW170817.
3. dUrca Constraints: A likelihood term derived from cooling observations. Rapidly cooling sources (SAX J1808.4−3658, Cassiopeia A, PSR B2334+61, Vela pulsar) are assumed to have active dUrca, implying their masses exceed the threshold mass $M_{\text{DU}}$ where $n_c \ge n_{\text{DU}}$ . Slowly cooling sources (HESS J1731−347) are assumed to lack active dUrca.
Handling Model Breakdown: If the dUrca threshold occurs within the nucleonic branch ( $n_{\text{DU}} < n_{\text{break}}$ ), the likelihood is evaluated directly. If the threshold falls within the agnostic CSE extension ( $n_{\text{DU}} > n_{\text{break}}$ ), the likelihood is marginalized over the unknown density where dUrca would occur, as the microscopic composition cannot be deduced in the CSE region.
Sampling: Posterior distributions are sampled using a sequential Monte Carlo algorithm with 4000 particles, accelerated via JAX on GPUs.

Key Results

Impact on Macroscopic Properties: The inclusion of dUrca constraints has a minor impact on macroscopic NS properties when combined with strong astrophysical and nuclear constraints. The posterior distributions for the radius of a $1.4M_\odot$ star ( $R_{1.4}$ ) and the maximum mass ( $M_{\text{TOV}}$ ) remain largely consistent with analyses using only $\chi$ EFT and astrophysical data. The dUrca constraint tends to favor stiffer EoSs (larger radii), but this effect is subdominant to the constraints provided by tidal deformability and mass measurements.
Impact on Composition: In contrast to macroscopic properties, the dUrca information provides a significant constraint on the composition of canonical-mass NSs. The combined posterior ( $\chi$ EFT + Astro + dUrca) concentrates the proton fraction at the center of a $1.4M_\odot$ star around $Y_p \approx 0.15$ . This is notably higher than values derived from laboratory symmetry-energy data alone (which suggest $Y_p \approx 0.07\text{--}0.09$ and disfavor dUrca onset below $1.4M_\odot$ ).
Model Dependence: The results are qualitatively consistent across the metamodel and extended-Skyrme constructions. However, the extended-Skyrme functional retains the nucleonic hypothesis up to higher densities than the metamodel, resulting in a larger fraction of the posterior providing a well-defined composition for canonical-mass stars. The metamodel results are more sensitive to the choice of $n_{\text{break}}$ (random vs. causal limit) compared to the Skyrme results.
Tidal Deformability and Symmetry Energy: The dUrca constraint yields a broad distribution for the tidal deformability ( $\Lambda_{1.4}$ ) and the slope of the symmetry energy ( $L_{\text{sym}}$ ), indicating that current cooling data adds little precision to these macroscopic parameters compared to gravitational wave and NICER data.

Significance and Claims
The paper claims that incorporating composition-sensitive observables, specifically dUrca cooling constraints, into multi-messenger inference is a viable step toward a more comprehensive treatment of the NS EoS. The primary contribution is demonstrating that while cooling data may not drastically alter macroscopic inferences (mass/radius) in the presence of other strong constraints, it provides complementary information regarding the microscopic composition (proton fraction) of neutron stars.

The authors emphasize that their results highlight the potential of using cooling observations to constrain the proton fraction, suggesting that rapid cooling near $1.4M_\odot$ implies a higher proton fraction than previously inferred from nuclear data alone. They note that this conclusion relies on the assumption that the observed rapid cooling is driven by nucleonic dUrca; if the cooling is due to exotic degrees of freedom or suppressed by superfluidity, the interpretation would differ. The study serves as a proof of principle for systematically including composition information in Bayesian EoS inference, moving beyond purely macroscopic observables.

More like this