Neutron stars are not only of astrophysical interest, but are also of great interest to nuclear physicists because their attributes can be used to determine the properties of the dense matter in their cores. One of the most informative approaches for determining the equation of state (EoS) of this dense matter is to measure both a star’s equatorial circumferential radius R _e and its gravitational mass M. Here we report estimates of the mass and radius of the isolated 205.53 Hz millisecond pulsar PSR J0030+0451 obtained using a Bayesian inference approach to analyze its energy-dependent thermal X-ray waveform, which was observed using the Neutron Star Interior Composition Explorer (NICER). This approach is thought to be less subject to systematic errors than other approaches for estimating neutron star radii. We explored a variety of emission patterns on the stellar surface. Our best-fit model has three oval, uniform-temperature emitting spots and provides an excellent description of the pulse waveform observed using NICER. The radius and mass estimates given by this model are km and (68%). The independent analysis reported in the companion paper by Riley et al. explores different emitting spot models, but finds spot shapes and locations and estimates of R _e and M that are consistent with those found in this work. We show that our measurements of R _e and M for PSR J0030+0451 improve the astrophysical constraints on the EoS of cold, catalyzed matter above nuclear saturation density.

PSR J0740+6620 has a gravitational mass of 2.08 ± 0.07 M _⊙, which is the highest reliably determined mass of any neutron star. As a result, a measurement of its radius will provide unique insight into the properties of neutron star core matter at high densities. Here we report a radius measurement based on fits of rotating hot spot patterns to Neutron Star Interior Composition Explorer (NICER) and X-ray Multi-Mirror (XMM-Newton) X-ray observations. We find that the equatorial circumferential radius of PSR J0740+6620 is 13.7 − 1.5 + 2.6 km (68%). We apply our measurement, combined with the previous NICER mass and radius measurement of PSR J0030+0451, the masses of two other ∼2 M _⊙ pulsars, and the tidal deformability constraints from two gravitational wave events, to three different frameworks for equation-of-state modeling, and find consistent results at ∼1.5–5 times nuclear saturation density. For a given framework, when all measurements are included, the radius of a 1.4 M _⊙ neutron star is known to ±4% (68% credibility) and the radius of a 2.08 M _⊙ neutron star is known to ±5%. The full radius range that spans the ±1σ credible intervals of all the radius estimates in the three frameworks is 12.45 ± 0.65 km for a 1.4 M _⊙ neutron star and 12.35 ± 0.75 km for a 2.08 M _⊙ neutron star.

We report on Bayesian estimation of the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray Timing Instrument event data. We condition on informative pulsar mass, distance, and orbital inclination priors derived from the joint North American Nanohertz Observatory for Gravitational Waves and Canadian Hydrogen Intensity Mapping Experiment/Pulsar wideband radio timing measurements of Fonseca et al. We use XMM-Newton European Photon Imaging Camera spectroscopic event data to inform our X-ray likelihood function. The prior support of the pulsar radius is truncated at 16 km to ensure coverage of current dense matter models. We assume conservative priors on instrument calibration uncertainty. We constrain the equatorial radius and mass of PSR J0740+6620 to be 12.39 − 0.98 + 1.30 km and 2.072 − 0.066 + 0.067 M _⊙ respectively, each reported as the posterior credible interval bounded by the 16% and 84% quantiles, conditional on surface hot regions that are non-overlapping spherical caps of fully ionized hydrogen atmosphere with uniform effective temperature; a posteriori, the temperature is log 10 ( T [ K ] ) = 5.99 − 0.06 + 0.05 for each hot region. All software for the X-ray modeling framework is open-source and all data, model, and sample information is publicly available, including analysis notebooks and model modules in the Python language. Our marginal likelihood function of mass and equatorial radius is proportional to the marginal joint posterior density of those parameters (within the prior support) and can thus be computed from the posterior samples.

The structure of neutron stars is considered from theoretical and observational perspectives. We demonstrate an important aspect of neutron star structure: the neutron star radius is primarily determined by the behavior of the pressure of matter in the vicinity of nuclear matter equilibrium density. In the event that extreme softening does not occur at these densities, the radius is virtually independent of the mass and is determined by the magnitude of the pressure. For equations of state with extreme softening, or those that are self-bound, the radius is more sensitive to the mass. Our results show that in the absence of extreme softening, a measurement of the radius of a neutron star more accurate than about 1 km will usefully constrain the equation of state. We also show that the pressure near nuclear matter density is primarily a function of the density dependence of the nuclear symmetry energy, while the nuclear incompressibility and skewness parameters play secondary roles. In addition, we show that the moment of inertia and the binding energy of neutron stars, for a large class of equations of state, are nearly universal functions of the star's compactness. These features can be understood by considering two analytic, yet realistic, solutions of Einstein's equations, due, respectively, to Buchdahl and Tolman. We deduce useful approximations for the fraction of the moment of inertia residing in the crust, which is a function of the stellar compactness and, in addition, the presssure at the core-crust interface.Comment: 51 pages, 10 figures, submitted to Ap

Neutron stars are some of the densest manifestations of massive objects in the universe. They are ideal astrophysical laboratories for testing theories of dense matter physics and provide connections among nuclear physics, particle physics, and astrophysics. Neutron stars may exhibit conditions and phenomena not observed elsewhere, such as hyperon-dominated matter, deconfined quark matter, superfluidity and superconductivity with critical temperatures near 10 ¹⁰ kelvin, opaqueness to neutrinos, and magnetic fields in excess of 10 ¹³ Gauss. Here, we describe the formation, structure, internal composition, and evolution of neutron stars. Observations that include studies of pulsars in binary systems, thermal emission from isolated neutron stars, glitches from pulsars, and quasi-periodic oscillations from accreting neutron stars provide information about neutron star masses, radii, temperatures, ages, and internal compositions.