Title: Interplay of astrophysics and nuclear physics in determining the properties of neutron stars
Neutron star properties depend on both nuclear physics and astrophysical processes, and thus observations of neutron stars offer constraints on both large-scale astrophysics and the behavior of cold, dense matter. In this study, we use astronomical data to jointly infer the universal equation of state of dense matter along with two distinct astrophysical populations: Galactic neutron stars observed electromagnetically and merging neutron stars in binaries observed with gravitational waves. We place constraints on neutron star properties and quantify the extent to which they are attributable to macrophysics or microphysics. We confirm previous results indicating that the Galactic and merging neutron stars have distinct mass distributions. The inferred maximum mass of both Galactic neutron stars, 𝑀pop,EM=2.0⁢5+0.11−0.06⁢𝑀⊙ (median and 90% symmetric credible interval), and merging neutron star binaries, 𝑀pop,GW =1.8⁢5+0.39−0.16⁢𝑀⊙, are consistent with the maximum mass of nonrotating neutron stars set by nuclear physics, 𝑀TOV =2.2⁢8+0.41−0.21⁢𝑀⊙. The radius of a 1.4⁢𝑀⊙ neutron star is 12.2+0.8−0.9  km, consistent with, though ∼20% tighter than, previous results using an identical equation of state model. Even though observed Galactic and merging neutron stars originate from populations with distinct properties, there is currently no evidence that astrophysical processes cannot produce neutron stars up to the maximum value imposed by nuclear physics.  more » « less
Award ID(s):
2207758
PAR ID:
10613700
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
American Physical Society
Date Published:
Journal Name:
Physical Review D
Volume:
111
Issue:
2
ISSN:
2470-0010
Subject(s) / Keyword(s):
neutron stars neutron star properties merging neutron stars gravitational waves
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. 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
  3. Abstract Neutron stars may experience differential rotation on short, dynamical timescales following extreme astrophysical events like binary neutron star mergers. In this work, the masses and radii of differentially rotating neutron star models are computed. We employ a set of equations of states for dense hypernuclear and ‐admixed‐hypernuclear matter obtained within the framework of CDF theory in the relativistic Hartree‐Fock (RHF) approximation. Results are shown for varying meson‐ couplings, or equivalently the ‐potential in nuclear matter. A comparison of our results with those obtained for nonrotating stars shows that the maximum mass difference between differentially rotating and static stars is independent of the underlying particle composition of the star. We further find that the decrease in the radii and increase in the maximum masses of stellar models when ‐isobars are added to hyperonuclear matter (as initially observed for static and uniformly rotating stars) persist also in the case of differentially rotating neutron stars. 
    more » « less
  4. In numerical simulations of binary neutron star systems, the equation of state of the dense neutron star matter is an important factor in determining both the physical realism and the numerical accuracy of the simulations. Some equations of state used in simulations are 𝐶2 or smoother in the pressure/density relationship function, such as a polytropic equation of state, but may not have the flexibility to model stars or remnants of different masses while keeping their radii within known astrophysical constraints. Other equations of state, such as tabular or piecewise polytropic, may be flexible enough to model additional physics and multiple stars’ masses and radii within known constraints, but are not as smooth, resulting in additional numerical error. We will study in this paper a recently developed family of equation of state, using a spectral expansion with sufficient free parameters to allow for a larger flexibility than current polytropic equations of state, and with sufficient smoothness to reduce numerical errors compared to tabulated or piecewise polytropic equations of state. We perform simulations at three mass ratios with a common chirp mass, using two distinct spectral equations of state, and at multiple numerical resolutions. We evaluate the gravitational waves produced from these simulations, comparing the phase error between resolutions and equations of state, as well as with respect to analytical models. From our simulations, we estimate that the phase difference at the merger for binaries with a dimensionless weighted tidal deformability difference greater than Δ⁢Λ≈55 can be captured by the spectral Einstein code for these equations of state. 
    more » « less
  5. Recent experimental and ab initio theory investigations of the 208Pb neutron skin thickness have the potential to inform the neutron star equation of state. In particular, the strong correlation between the 208Pb neutron skin thickness and the pressure of neutron matter at normal nuclear densities leads to modified predictions for the radii, tidal deformabilities, and moments of inertia of typical 1.4M⊙ neutron stars. In the present work, we study the relative impact of these recent analyses of the 208Pb neutron skin thickness on bulk properties of neutron stars within a Bayesian statistical analysis. Two models for the equation of state prior are employed in order to highlight the role of the highly uncertain high-density equation of state. From our combined Bayesian analysis of nuclear theory, nuclear experiment, and observational constraints on the dense matter equation of state, we find at the 90% credibility level R1.4=12.36−0.73+0.38 km for the radius of a 1.4M⊙ neutron star, R2.0=11.96−0.71+0.94 km for the radius of a 2.0M⊙ neutron star, Λ1.4=440−144+103 for the tidal deformability of a 1.4M⊙ neutron star, and I1.338=1.425−0.146+0.074×1045gcm2 for the moment of inertia of PSR J0737-3039A whose mass is 1.338M⊙. 
    more » « less