A new way to study neutron stars using simulations and universal patterns
Simulation-based inference helps to deal with the uncertainty about neutron star matter.
If one were to describe neutron star properties in terms of simple rules, universal relations act as shortcuts linking quantities like mass, radius, oscillation frequencies or others almost independently of the uncertain nuclear equation of state (EOS). However, their predictive power is limited by the difficulty of identifying optimal parameter combinations and, in particular, by reliably quantifying systematic uncertainties, an issue that is especially relevant for gravitational-wave observations. In this study, researchers from the University of Tübingen and the AEI introduce a new approach based on simulation-based inference (SBI), a machine-learning method that learns directly from simulated data and treats variations across EOS models as intrinsic “EOS noise.” This allows for a systematic exploration of correlations, recovers known universal relations, and reveals a new relation linking the neutron star radius to its mass and two oscillation frequencies. The study shows that SBI can capture universal relations, and often, given a sufficiently large data sample, even outperform them while providing reliable uncertainty estimates out of the box. In contrast, standard universal relations do not come with a measure of their systematic uncertainty and a naive approach tends to considerably underestimate it. The work demonstrates that SBI offers a powerful complementary perspective, turning universal relations into part of a broader probabilistic framework for precision neutron star physics.
Paper abstract
In this work, we propose a novel approach for identifying, constructing, and validating precise and accurate universal relations for neutron star bulk quantities. A central element is simulation-based inference (SBI), which we adopt to treat uncertainties due to the unknown nuclear equation of state (EOS) as intrinsic non-trivial noise. By assembling a large set of bulk properties of non-rotating neutron stars across multiple state-of-the-art EOS models, we are able to systematically explore universal relations in high-dimensional parameter spaces. Our framework further identifies the most promising parameter combinations, enabling a more focused and traditional construction of explicit universal relations. At the same time, SBI does not rely on explicit relations; instead, it directly provides predictive distributions together with a quantitative measure of systematic uncertainties, which are not captured by conventional approaches. As an example, we report a new universal relation that allows us to obtain the radius as a function of mass, fundamental mode, and one pressure mode. Our analysis shows that SBI can surpass the predictive power of this universal relation while also mitigating systematic errors. Finally, we demonstrate how universal relations can be further calibrated to mitigate systematic errors accurately.
