An official website of the United States government
The interiors of neutron stars reach densities and temperatures beyond the limits of terrestrial experiments, providing vital laboratories for probing nuclear physics. While the star's interior is not directly observable, its pressure and density determine the star's macroscopic structure which affects the spectra observed in telescopes. The relationship between the observations and the internal state is complex and partially intractable, presenting difficulties for inference. Previous work has focused on the regression from stellar spectra of parameters describing the internal state. We demonstrate a calculation of the full likelihood of the internal state parameters given observations, accomplished by replacing intractable elements with machine learning models trained on samples of simulated stars. Our machine-learning-derived likelihood allows us to performmaximum a posterioriestimation of the parameters of interest, as well as full scans. We demonstrate the technique by inferring stellar mass and radius from an individual stellar spectrum, as well as equation of state parameters from a set of spectra. Our results are more precise than pure regression models, reducing the width of the parameter residuals by 11.8% in the most realistic scenario. The neural networks will be released as a tool for fast simulation of neutron star properties and observed spectra.
more »
« less
Neural simulation-based inference of the neutron star equation of state directly from telescope spectra
Abstract Neutron stars provide a unique opportunity to study strongly interacting matter under extreme density conditions. The intricacies of matter inside neutron stars and their equation of state are not directly visible, but determine bulk properties, such as mass and radius, which affect the star's thermal X-ray emissions. However, the telescope spectra of these emissions are also affected by the stellar distance, hydrogen column, and effective surface temperature, which are not always well-constrained. Uncertainties on these nuisance parameters must be accounted for when making a robust estimation of the equation of state. In this study, we develop a novel methodology that, for the first time, can infer the full posterior distribution of both the equation of state and nuisance parameters directly from telescope observations. This method relies on the use of neural likelihood estimation, in which normalizing flows use samples of simulated telescope data to learn the likelihood of the neutron star spectra as a function of these parameters, coupled with Hamiltonian Monte Carlo methods to efficiently sample from the corresponding posterior distribution. Our approach surpasses the accuracy of previous methods, improves the interpretability of the results by providing access to the full posterior distribution, and naturally scales to a growing number of neutron star observations expected in the coming years.
more »
« less
- PAR ID:
- 10555447
- Publisher / Repository:
- Journal of Cosmology and Astroparticle Physics
- Date Published:
- Journal Name:
- Journal of Cosmology and Astroparticle Physics
- Volume:
- 2024
- Issue:
- 09
- ISSN:
- 1475-7516
- Page Range / eLocation ID:
- 009
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Gravitational-wave observations of binary neutron star mergers provide valuable information about neutron star structure and the equation of state of dense nuclear matter. Numerous methods have been proposed to analyze the population of observed neutron stars, and previous work has demonstrated the necessity of jointly fitting the astrophysical distribution and the equation of state in order to accurately constrain the equation of state. In this work, we introduce a new framework to simultaneously infer the distribution of binary neutron star masses and the nuclear equation of state using Gaussian mixture model density estimates, which mitigates some of the limitations previously used methods suffer from. Using our method, we reproduce previous projections for the expected precision of our joint mass distribution and equation-of-state inference with tens of observations. We also show that mismodeling the equation of state can bias our inference of the neutron star mass distribution. While we focus on neutron star masses and matter effects, our method is widely applicable to population inference problems.more » « less
-
Abstract: With recent advances in astronomical observations, major progress has been made in determining the pressure of neutron star matter at high density. This pressure is constrained by the neutron star deformability, determined from gravitational waves emitted in a neutron-star merger, and the mass-radii relation of two neutron stars, determined from a new X-ray observatory on the International Space Station. Previous studies have relied on nuclear theory calculations to constrain the equation of state at low density. Here we use a combination of constraints composed of three astronomical observations and twelve nuclear experimental constraints that extend over a wide range of densities. A Bayesian inference framework is then used to obtain a comprehensive nuclear equation of state. This data-centric result provides benchmarks for theoretical calculations and modeling of nuclear matter and neutron stars. Furthermore, it provides insights into the microscopic degrees of freedom of the nuclear matter equation of state and on the composition of neutron stars and their cooling via neutrino radiation.more » « less
-
Neutron stars provide a window into the properties of dense nuclear matter. Several recent observational and theoretical developments provide powerful constraints on their structure and internal composition. Among these are the first observed binary neutron star merger, GW170817, whose gravitational radiation was accompanied by electromagnetic radiation from a short γ-ray burst and an optical afterglow believed to be due to the radioactive decay of newly minted heavy r-process nuclei. These observations give important constraints on the radii of typical neutron stars and on the upper limit to the neutron star maximum mass and complement recent pulsar observations that established a lower limit. Pulse-profile observations by the Neutron Star Interior Composition Explorer (NICER) X-ray telescope provide an independent, consistent measure of the neutron star radius. Theoretical many-body studies of neutron matter reinforce these estimates of neutron star radii. Studies using parameterized dense matter equations of state (EOSs) reveal several EOS-independent relations connecting global neutron star properties.more » « less
-
We present two recent parametrizations of the equation of state (FSU2R and FSU2H models) that reproduce the properties of nuclear matter and finite nuclei, fulfill constraints on high-density matter stemming from heavy-ion collisions, produce 2 M_Sun neutron stars, and generate neutron star radii below 13 km. Making use of these equations of state, cooling simulations for isolated neutron stars are performed. We find that two of the models studied, FSU2R (with nucleons) and, in particular, FSU2H (with nucleons and hyperons), show very good agreement with cooling observations, even without including nucleon pairing. This indicates that cooling observations are compatible with an equation of state that produces a soft nuclear symmetry energy and, thus, generates small neutron star radii. Nevertheless, both schemes produce cold isolated neutron stars with masses above 1.8 M_Sun.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1088/1475-7516/2024/09/009
