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   
                    
                            
                            Determination of the equation of state from nuclear experiments and neutron star observations
                        
                    
    
            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   
        
    
                            - Award ID(s):
- 2209145
- PAR ID:
- 10519252
- Publisher / Repository:
- Nature Astronomy
- Date Published:
- Journal Name:
- Nature Astronomy
- Edition / Version:
- 1
- Volume:
- 8
- Issue:
- 3
- ISSN:
- 2397-3366
- Page Range / eLocation ID:
- 328 to 336
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
- 
            
- 
            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
- 
            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
- 
            Abstract Over the past decade, an abundance of information from neutron-star observations, nuclear experiments and theory has transformed our efforts to elucidate the properties of dense matter. However, at high densities relevant to the cores of neutron stars, substantial uncertainty about the dense matter equation of state (EoS) remains. In this work, we present a semiparametric equation of state framework aimed at better integrating knowledge across these domains in astrophysical inference. We use a Meta-model and realistic crust at low densities, and Gaussian Process extensions at high densities. Comparisons between our semiparametric framework to fully nonparametric EoS representations show that imposing nuclear theoretical and experimental constraints through the Meta-model up to nuclear saturation density results in constraints on the pressure up to twice nuclear saturation density. We also show that our Gaussian Process trained on EoS models with nucleonic, hyperonic, and quark compositions extends the range of EoS explored at high density compared to a piecewise polytropic extension schema, under the requirements of causality of matter and of supporting the existence of heavy pulsars. We find that maximum TOV masses above $$3.2 M_{\odot}$$ can be supported by causal EoS compatible with nuclear constraints at low densities. We then combine information from existing observations of heavy pulsar masses, gravitational waves emitted from binary neutron star mergers, and X-ray pulse profile modeling of millisecond pulsars within a Bayesian inference scheme using our semiparametric EoS prior. With information from all public NICER pulsars (including PSR J0030$$+$$0451, PSR J0740$$+$$6620, PSR J0437-4715, and PSR J0614-3329), we find an astrophysically favored pressure at two times nuclear saturation density of $$P(2\rho_{\rm nuc}) = 1.98^{+2.13}_{-1.08}\times10^{34}$$ dyn/cm$$^{2}$$, a radius of a $$1.4 M_{\odot}$$ neutron star value of $$R_{1.4} = 11.4^{+0.98}_{-0.60}$$\;km, and $$M_{\rm max} = 2.31_{-0.23}^{+0.35} M_{\odot}$$ at the 90\% credible level.more » « less
- 
            Born in the aftermath of core-collapse supernovae, neutron stars contain matter under extraordinary conditions of density and temperature that are difficult to reproduce in the laboratory. In recent years, neutron star observations have begun to yield novel insights into the nature of strongly interacting matter in the high-density regime where current theoretical models are challenged. At the same time, chiral effective field theory has developed into a powerful framework to study nuclear matter properties with quantified uncertainties in the moderate-density regime for modeling neutron stars. In this article, we review recent developments in chiral effective field theory and focus on many-body perturbation theory as a computationally efficient tool for calculating the properties of hot and dense nuclear matter. We also demonstrate how effective field theory enables statistically meaningful comparisons among nuclear theory predictions, nuclear experiments, and observational constraints on the nuclear equation of state.more » « less
 An official website of the United States government
An official website of the United States government 
				
			 
					 
					
 
                                    