A new Bayesian method is proposed for inferring the equation of state of neutron stars from mass and radius observations and neutron matter theory. Unlike traditional parametrizations in the pressure–energy density space, here prior distributions are placed directly in the mass–radius space, which directly characterizes the observables. Exact inverse approximations are used to map back to the equation of state. A systematic comparison is carried out with inferences from standard parametrizations, quantifying systematic uncertainties stemming from prior assumptions in both approaches. It is shown that accounting for prior uncertainties is essential in all Bayesian methods. The proposed approach more fully covers the physically permissible mass–radius space, especially for small-radius configurations, while also boosting computational efficiency and greatly reducing dependence on the choice of priors. Direct parametrization in the space of observables offers a robust and efficient alternative to traditional methods.
Inside a neutron star, matter is compressed to unimaginable densities. Previously, scientists speculated about its exact makeup, cycling through possible equations linking pressure and density. The new method tosses out guesswork: it directly ties the measurable mass and radius to the internal structure.
It’s like a watermelon: by weight and girth, you can judge ripeness without cracking it open. But with neutron stars, it’s even wilder: the heavier the star, the smaller it can be in diameter—thanks to monstrous gravity. Now, by measuring the mass and radius of actual pulsars (rapidly spinning neutron stars) or during supernova explosions, the method quickly reconstructs the “stuffing.”
In the future, this approach will not only help astrophysicists understand the birth of neutron stars, but also peek into the laws of matter under extreme conditions that can’t be recreated on Earth.
🎯 Collisions of neutron stars are the only known places where precious metals like gold and platinum are born.