We report on Bayesian parameter estimation of the mass and equatorial radius of the millisecond pulsar PSRJ0030 +0451, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray spectraltiming event data. We perform relativistic ray-tracing of thermal emission from hot regions of the pulsar's surface. We assume two distinct hot regions based on two clear pulsed components in the phase-folded pulse-profile data; we explore a number of forms (morphologies and topologies) for each hot region, inferring their parameters in addition to the stellar mass and radius. For the family of models considered, the evidence (prior predictive probability of the data) strongly favors a model that permits both hot regions to be located in the same rotational hemisphere. Models wherein both hot regions are assumed to be simply connected circular single-temperature spots, in particular those where the spots are assumed to be reflection-symmetric with respect to the stellar origin, are strongly disfavored. For the inferred configuration, one hot region subtends an angular extent of only a few degrees (in spherical coordinates with origin at the stellar center) and we are insensitive to other structural details; the second hot region is far more azimuthally extended in the form of a narrow arc, thus requiring a larger number of parameters to describe. The inferred mass M and equatorial radius R eq are, respectively,-+ eq 2 0.010 0.008 is more tightly constrained; the credible interval bounds reported here are approximately the 16% and 84% quantiles in marginal posterior mass.

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 present first results from a LOFAR census of non-recycled pulsars. The census includes almost all such pulsars known (194 sources) at declinations Dec > 8• and Galactic latitudes |Gb| > 3• , regardless of their expected flux densities and scattering times. Each pulsar was observed for ≥20 min in the contiguous frequency range of 110-188 MHz. Full-Stokes data were recorded. We present the dispersion measures, flux densities, and calibrated total intensity profiles for the 158 pulsars detected in the sample. The median uncertainty in census dispersion measures (1.5 × 10 −3 pc cm −3 ) is ten times smaller, on average, than in the ATNF pulsar catalogue. We combined census flux densities with those in the literature and fitted the resulting broadband spectra with single or broken power-law functions. For 48 census pulsars such fits are being published for the first time. Typically, the choice between single and broken power-laws, as well as the location of the spectral break, were highly influenced by the spectral coverage of the available flux density measurements. In particular, the inclusion of measurements below 100 MHz appears essential for investigating the lowfrequency turnover in the spectra for most of the census pulsars. For several pulsars, we compared the spectral indices from different works and found the typical spread of values to be within 0.5-1.5, suggesting a prevailing underestimation of spectral index errors in the literature. The census observations yielded some unexpected individual source results, as we describe in the paper. Lastly, we will provide this unique sample of wide-band, low-frequency pulse profiles via the European Pulsar Network Database.