Can a Neural Network Predict a Neutron Star's Radius Better Than a Handwritten Formula?

Neutron stars are extreme laboratories for matter above nuclear saturation density and for testing general relativity. Their bulk properties—mass, radius, moment of inertia, tidal deformability, and oscillation frequencies—depend on the nuclear equation of state, which is uncertain at high densities. That uncertainty makes it hard to use one measured property to predict another. A new paper by Christian J. Krüger (University of Tübingen) and Sebastian H. Völkel (Max Planck Institute for Gravitational Physics) combines simulation-based inference with universal relations to tackle that problem.

A Universal Shortcut, and Its Blind Spots

Astrophysicists have long tried to sidestep that uncertainty with universal relations. These are empirical, equation-of-state-insensitive links between a neutron star's bulk properties: its mass, radius, moment of inertia (how resistantly it spins up), tidal deformability (how easily it is stretched by a companion), and the frequencies at which it rings. If such a relation is accurate, measuring one property can reveal another without knowing the equation of state.

Yet universal relations have practical weaknesses. Constructing a new one is a matter of trial and error. Even after a fit, the fit's covariance matrix—the usual way of stating parameter errors—does not reliably describe the systematic uncertainty from the unknown equation of state. And the robustness of a relation depends on how thoroughly EOS uncertainties and biases were accounted for, including the size of the viable EOS set and the number of agnostic parametrizations considered.

Teaching a Network to Treat Ignorance as Noise

Christian J. Krüger of the University of Tübingen and Sebastian H. Völkel of the Max Planck Institute for Gravitational Physics propose a different route: simulation-based inference, or SBI. SBI is a machine-learning approach that learns a probability distribution directly from simulated data. It does not require scientists to write down a likelihood function, the mathematical rule that connects parameters to noisy observations. Instead, it is trained on many examples of how the data could look.

The authors considered 1,491 equation-of-state realizations based on four different parametrizations, keeping only those consistent with basic astrophysical constraints: the speed of sound stays causal, the maximum-mass model reaches at least 1.97 solar masses, the radius of a 1.6-solar-mass star exceeds 10.6 kilometers, the radius of a 1.4-solar-mass star lies between 11.5 and 13.5 kilometers, and its tidal deformability lies between 120 and 800.

For each equation of state, they randomly selected five non-rotating neutron-star models with masses of at least one solar mass. The quantities included mass, radius, moment of inertia, tidal deformability, the quadrupolar fundamental f-mode frequency, and the first pressure p1-mode frequency. Oscillation modes are the star's characteristic ring tones.

They then used neural posterior estimation to explore all 602 possible divisions of those six quantities into 'data' and 'observations' subsets. The key idea is to treat the spread caused by different equations of state—and the finite accuracy of universal relations—as intrinsic noise, which the authors call EOS noise. For a given division, the network learns a posterior distribution for the target quantity, not just a single best-fit value.

Known Relations Return, and a New One Appears

The framework recovered well-known universal relations, including links between f-mode frequency and moment of inertia, moment of inertia and tidal deformability, and f-mode frequency and tidal deformability. It also flagged a new relation for the radius, R, expressed as a function of mass M, f-mode frequency f, and first pressure-mode frequency p1. The authors then constructed an explicit five-term algebraic fit. In their tests, most deviations for SBI and the calibrated relation were below about 80 meters, which is well under one percent for a star with a radius near 12.7 kilometers.

When a Tiny Error Bar Is a Warning

The most striking comparison came from a single held-out test star with mass 2.33 solar masses, f-mode frequency 1.87 kHz, and p1-mode frequency 7.08 kHz. The true radius was 12.7 kilometers. The SBI posterior included that value within its 68 percent highest-density interval, a range only about 22 meters wide, and its mean was off by roughly 2 meters. The calibrated universal relation covered the true value within its 90 percent interval, with a 68 percent width of about 62 meters and a mean off by about 39 meters.

The uncalibrated universal relation looked far more precise: its 68 percent interval was only about 1.9 meters wide. But the true radius lay well outside that interval. That is the trap the paper highlights. A tiny error bar from a traditional fit can be an illusion, because the fit parameters do not reliably absorb the systematic spread from unknown nuclear physics.

Across 1,490 test points, SBI was slightly better than the calibrated relation and often narrower, with 68 percent widths as small as 18 meters. The calibrated relation usually had widths around 70 meters. The uncalibrated relation was narrower still but unreliable. Calibration plots showed that the SBI posteriors and the calibrated universal relation both delivered error bars close to the truth, while the uncalibrated relation badly underestimated its uncertainty.

Precision, Meet Honesty

SBI's advantage has a cost. The neural network optimizes more than 27,000 parameters, while the universal relation has only five free coefficients. The relation can be evaluated by hand; the network needs software and trained weights. But the relation can also be calibrated. The authors introduce an effective EOS-noise error, and a simpler residual-based method gives a standard deviation of about 37 meters, which produced a surprisingly good approximation to the systematic uncertainty.

The work is a methods paper, not a new observation. The neutron stars are simulated models, not detected neutron stars, and the analysis is limited to non-rotating stars. The authors note that extending it to rotating stars is technically possible with SBI, but generating rotating-star observables in large quantities can be harder, especially for oscillation modes. The paper says it would be interesting to quantify how measurement uncertainties from future detectors such as the Einstein Telescope or Cosmic Explorer would affect these applications; they are not included here. The paper is posted as an arXiv preprint.

Still, the result points to a broader lesson. In neutron-star science, precision without an honest account of systematic error can mislead. The new framework does not just shrink the error bars on a radius. It helps show which error bars deserve trust—and that may matter as much as the radius itself.