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)

Magnetars are the strongest magnets in the present universe and the combination of extreme magnetic field, gravity and density makes them unique laboratories to probe current physical theories (from quantum electrodynamics to general relativity) in the strong field limit. Magnetars are observed as peculiar, burst-active x-ray pulsars, the anomalous x-ray pulsars (AXPs) and the soft gamma repeaters (SGRs); the latter emitted also three 'giant flares', extremely powerful events during which luminosities can reach up to 10(47) erg s(-1) for about one second. The last five years have witnessed an explosion in magnetar research which has led, among other things, to the discovery of transient, or 'outbursting', and 'low-field' magnetars. Substantial progress has been made also on the theoretical side. Quite detailed models for explaining the magnetars' persistent x-ray emission, the properties of the bursts, the flux evolution in transient sources have been developed and confronted with observations. New insight on neutron star asteroseismology has been gained through improved models of magnetar oscillations. The long-debated issue of magnetic field decay in neutron stars has been addressed, and its importance recognized in relation to the evolution of magnetars and to the links among magnetars and other families of isolated neutron stars. The aim of this paper is to present a comprehensive overview in which the observational results are discussed in the light of the most up-to-date theoretical models and their implications. This addresses not only the particular case of magnetar sources, but the more fundamental issue of how physics in strong magnetic fields can be constrained by the observations of these unique sources.