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 Star Radii, Deformabilities, and Moments of Inertia from Experimental and Ab Initio Theory Constraints of the 208Pb Neutron Skin Thickness
                        
                    
    
            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   
        
    
                            - Award ID(s):
- 1652199
- PAR ID:
- 10465255
- Date Published:
- Journal Name:
- Galaxies
- Volume:
- 10
- Issue:
- 5
- ISSN:
- 2075-4434
- Page Range / eLocation ID:
- 99
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
- 
            
- 
            We combine equation of state of dense matter up to twice nuclear saturation density (nsat = 0.16 fm−3 ) obtained using chiral effective field theory (χEFT), and recent observations of neutron stars to gain insights about the high-density matter encountered in their cores. A key element in our study is the recent Bayesian analysis of correlated EFT truncation errors based on order-byorder calculations up to next-to-next-to-next-to-leading order in the χEFT expansion. We refine the bounds on the maximum mass imposed by causality at high densities, and provide stringent limits on the maximum and minimum radii of ∼ 1.4 M and ∼ 2.0 M stars. Including χEFT predictions from nsat to 2 nsat reduces the permitted ranges of the radius of a 1.4 M star, R1.4, by ∼ 3.5 km. If observations indicate R1.4 < 11.2 km, our study implies that either the squared speed of sound c 2 s > 1/2 for densities above 2 nsat, or that χEFT breaks down below 2 nsat. We also comment on the nature of the secondary compact object in GW190814 with mass ' 2.6 M , and discuss the implications of massive neutron stars > 2.1 M (2.6 M ) in future radio and gravitational-wave searches. Some form of strongly interacting matter with c 2 s > 0.35 (0.55) must be realized in the cores of such massive neutron stars. In the absence of phase transitions below 2 nsat, the small tidal deformability inferred from GW170817 lends support for the relatively small pressure predicted by χEFT for the baryon density nB in the range 1−2 nsat. Together they imply that the rapid stiffening required to support a high maximum mass should occur only when nB & 1.5 − 1.8 nsat.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
- 
            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.05+0.11−0.06𝑀⊙ (median and 90% symmetric credible interval), and merging neutron star binaries, 𝑀pop,GW =1.85+0.39−0.16𝑀⊙, are consistent with the maximum mass of nonrotating neutron stars set by nuclear physics, 𝑀TOV =2.28+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
- 
            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
 An official website of the United States government
An official website of the United States government 
				
			 
					 
					
 
                                    