This catalog summarizes 117 high-confidence 0.1 GeV gamma-ray pulsar detections using three years of data acquired by the Large Area Telescope (LAT) on the Fermi satellite. Half are neutron stars discovered using LAT data through periodicity searches in gamma-ray and radio data around LAT unassociated source positions. The 117 pulsars are evenly divided into three groups: millisecond pulsars, young radio-loud pulsars, and young radio-quiet pulsars. We characterize the pulse profiles and energy spectra and derive luminosities when distance information exists. Spectral analysis of the off-peak phase intervals indicates probable pulsar wind nebula emission for four pulsars, and off-peak magnetospheric emission for several young and millisecond pulsars. We compare the gammaray properties with those in the radio, optical, and X-ray bands. We provide flux limits for pulsars with no observed gamma-ray emission, highlighting a small number of gamma-faint, radio-loud pulsars. The large, varied gamma-ray pulsar sample constrains emission models. Fermi's selection biases complement those of radio surveys, enhancing comparisons with predicted population distributions.

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)

Context. Accretion-powered X-ray pulsars exhibit significant variability of the cyclotron resonance scattering feature (CRSF) centroid energy on pulse-to-pulse timescales, and also on much longer timescales. Two types of spectral variability are observed. For sources in group 1, the CRSF energy is negatively correlated with the variable source luminosity, and for sources in group 2, the opposite behavior is observed. The physical basis for this bimodal behavior is currently not well understood. Aims. We explore the hypothesis that the accretion dynamics in the group 1 sources is dominated by radiation pressure near the stellar surface, and that Coulomb interactions decelerate the gas to rest in the group 2 sources. Methods. We derive a new expression for the critical luminosity, L crit , such that radiation pressure decelerates the matter to rest in sources with X-ray luminosity L X > L crit . The formula for L crit is based on a simple physical model for the structure of the accretion column in luminous X-ray pulsars that takes into account radiative deceleration, the energy dependence of the cyclotron cross section, the thermodynamics of the accreting gas, the dipole structure of the pulsar magnetosphere, and the diffusive escape of radiation through the column walls. We show that for typical neutron star parameters, L crit = 1.5 × 10 37 B 16/15 12 erg s −1 , where B 12 is the surface magnetic field strength in units of 10 12 G. Results. The formula for the critical luminosity is evaluated for five sources, using the maximum value of the CRSF centroid energy to estimate the surface magnetic field strength B 12 . The results confirm that the group 1 sources are supercritical (L X > L crit ) and the group 2 sources are subcritical (L X < L crit ), although the situation is less clear for those highly variable sources that cross over the line L X = L crit . We also explain the variation of the CRSF energy with luminosity as a consequence of the variation of the characteristic emission height. The sign of this dependence is opposite in the supercritical and subcritical cases, hence creating the observed bimodal behavior. Conclusions. We have developed a new model for the critical luminosity in accretion-powered X-ray pulsars that explains the bimodal dependence of the CRSF centroid energy on the X-ray luminosity L X . Our model provides a physical basis for the observed variation of the CRSF energy as a function of L X for both the group 1 (supercritical) and the group 2 (subcritical) sources as a result of the variation of the emission height in the column.

We develop a new theoretical model for the spectral formation process in accretion-powered X-ray pulsars based on a detailed treatment of the bulk and thermal Comptonization occurring in the accreting, shocked gas. A rigorous eigenfunction expansion method is employed to obtain the analytical solution for the Green's function describing the scattering of radiation injected into the column from a monochromatic source located at an arbitrary height above the stellar surface. The emergent spectrum is calculated by convolving the Green's function with source terms corresponding to bremsstrahlung, cyclotron, and blackbody emission. The energization of the photons in the shock, combined with cyclotron absorption, naturally produces an X-ray spectrum with a relatively flat continuum shape and a high-energy quasi-exponential cutoff. We demonstrate that the new theory successfully reproduces the phaseaveraged spectra of the bright pulsars Her X-1, LMC X-4, and Cen X-3. In these luminous sources, it is shown that the emergent spectra are dominated by Comptonized bremsstrahlung emission.