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 -+ 10 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. Unified Astronomy Thesaurus concepts: Millisecond pulsars (1062); Rotation powered pulsars (1408); Pulsars (1306); Radio pulsars (1353); X-ray astronomy (1810); Neutron stars (1108)

We report results from continued timing observations of PSR J0740+6620, a high-mass, 2.8 ms radio pulsar in orbit with a likely ultracool white dwarf companion. Our data set consists of combined pulse arrival-time measurements made with the 100 m Green Bank Telescope and the Canadian Hydrogen Intensity Mapping Experiment telescope. We explore the significance of timing-based phenomena arising from general relativistic dynamics and variations in pulse dispersion. When using various statistical methods, we find that combining ∼1.5 yr of additional, high-cadence timing data with previous measurements confirms and improves on previous estimates of relativistic effects within the PSR J0740+6620 system, with the pulsar mass m p = 2.08 − 0.07 + 0.07 M ⊙ (68.3% credibility) determined by the relativistic Shapiro time delay. For the first time, we measure secular variation in the orbital period and argue that this effect arises from apparent acceleration due to significant transverse motion. After incorporating contributions from Galactic differential rotation and off-plane acceleration in the Galactic potential, we obtain a model-dependent distance of d = 1.14 − 0.15 + 0.17 kpc (68.3% credibility). This improved distance confirms the ultracool nature of the white dwarf companion determined from recent optical observations. We discuss the prospects for future observations with next-generation facilities, which will likely improve the precision on m p for J0740+6620 by an order of magnitude within the next few years.

A nanohertz-frequency stochastic gravitational-wave background can potentially be detected through the precise timing of an array of millisecond pulsars. This background produces low-frequency noise in the pulse arrival times that would have a characteristic spectrum common to all pulsars and a well-defined spatial correlation. Recently the North American Nanohertz Observatory for Gravitational Waves collaboration (NANOGrav) found evidence for the common-spectrum component in their 12.5 yr data set. Here we report on a search for the background using the second data release of the Parkes Pulsar Timing Array. If we are forced to choose between the two NANOGrav models—one with a common-spectrum process and one without—we find strong support for the common-spectrum process. However, in this paper, we consider the possibility that the analysis suffers from model misspecification. In particular, we present simulated data sets that contain noise with distinctive spectra but show strong evidence for a common-spectrum process under the standard assumptions. The Parkes data show no significant evidence for, or against, the spatially correlated Hellings–Downs signature of the gravitational-wave background. Assuming we did observe the process underlying the spatially uncorrelated component of the background, we infer its amplitude to be A = 2.2 − 0.3 + 0.4 × 10 − 15 in units of gravitational-wave strain at a frequency of 1 yr−1. Extensions and combinations of existing and new data sets will improve the prospects of identifying spatial correlations that are necessary to claim a detection of the gravitational-wave background.

We present results from the search for a stochastic gravitational-wave background (GWB) as predicted by the theory of General Relativity using six radio millisecond pulsars from the Data Release 2 (DR2) of the European Pulsar Timing Array (EPTA) covering a timespan up to 24 yr. A GWB manifests itself as a long-term low-frequency stochastic signal common to all pulsars, a common red signal (CRS), with the characteristic Hellings-Downs (HD) spatial correlation. Our analysis is performed with two independent pipelines, ENTERPRISE, and TEMPONEST+FORTYTWO, which produce consistent results. A search for a CRS with simultaneous estimation of its spatial correlations yields spectral properties compatible with theoretical GWB predictions, but does not result in the required measurement of the HD correlation, as required for GWB detection. Further Bayesian model comparison between different types of CRSs, including a GWB, finds the most favoured model to be the common uncorrelated red noise described by a power law with $A = 5.13_{-2.73}^{+4.20} \times 10^{-15}$ and $\gamma = 3.78_{-0.59}^{+0.69}$ (95 per cent credible regions). Fixing the spectral index to γ = 13/3 as expected from the GWB by circular, inspiralling supermassive black hole binaries results in an amplitude of $A =2.95_{-0.72}^{+0.89} \times 10^{-15}$. We implement three different models, BAYESEPHEM, LINIMOSS, and EPHEMGP, to address possible Solar system ephemeris (SSE) systematics and conclude that our results may only marginally depend on these effects. This work builds on the methods and models from the studies on the EPTA DR1. We show that under the same analysis framework the results remain consistent after the data set extension.