PSR J0108$-$1431 is an old pulsar where the X-ray emission is expected to have a thermal component from the polar cap and a non-thermal component from the magnetosphere. Although the phase-integrated spectra are fit best with a single non-thermal component modeled with a power-law (PL) of photon index $\Gamma=2.9$, the X-ray pulse profiles do show the presence of phase-separated thermal and non-thermal components. The spectrum extracted from half the rotational phase away from the X-ray peak fits well with either a single blackbody (BB) or a neutron star atmosphere (NA) model, whereas, the spectrum from the rest of the phase range is dominated by a PL. From Bayesian analysis, the estimated BB area is smaller than the expected polar cap area for a dipolar magnetic field with a probability of 86% whereas the area estimate from the NA model is larger with a probability of 80%. Due to the ambiguity in the thermal emission model, the polar cap area cannot be reliably estimated and hence cannot be used to understand the nature of the surface magnetic field. Instead, we can infer the presence of multipolar magnetic field from the misalignment between the pulsar's thermal X-ray peak and the radio emission peak. For J0108$-$1431, we estimated a phase-offset $\Delta\phi > 0.1$ between the thermal polar cap emission peak and the radio emission peak and argue that this is best explained by the presence of a multipolar surface magnetic field.

We report on the energy-resolved timing and phase-resolved spectral analysis of X-ray emission from PSR J0659+1414 observed with XMM-Newton and NuSTAR. We find that the new data rule out the previously suggested model of the phase-dependent spectrum as a three-component (2 blackbodies + power-law) continuum, which shows large residuals between $0.3-0.7$ keV. Fitting neutron star atmosphere models or several blackbodies to the spectrum does not provide a better description of the spectrum, and requires spectral model components with unrealistically large emission region sizes. The fits improve significantly if we add a phase-dependent absorption feature with central energy $0.5-0.6$ keV and equivalent width up to $\approx 50$ eV. We detected the feature for about half of the pulse cycle. Energy-resolved pulse profiles support the description of the spectrum with a three-component continuum and an absorption component. The absorption feature could be interpreted as an electron cyclotron line originating in the pulsar magnetosphere and broadened by the non-uniformity of the magnetic field along the line of sight. The significant phase-variability in the thermal emission from the entire stellar surface may indicate multi-polar magnetic fields and a non-uniform temperature distribution. The strongly pulsed non-thermal spectral component detected with NuSTAR in the $3-20$ keV range is well fit by a power-law model with a photon index $\Gamma=1.5\pm0.2$.

We present the results of detailed X-ray analysis of two black-widow pulsars (BWPs), J1446-4701 and J1311-3430. PSR J1446-4701 is a BWP with orbital parameters near the median values of the sample of known BWPs. Its X-ray emission detected by $XMM-Newton$ is well characterized by a soft power-law (PL) spectrum (photon index $\Gamma \approx 3$), and it shows no significant orbital modulations. In view of a lack of radio eclipses and an optical non-detection, the system most likely has a low orbital inclination. PSR J1311-3430 is an extreme BWP with a very compact orbit and the lowest minimum mass companion. Our $Chandra$ data confirm the hard, $\Gamma \approx 1.3$, emission seen in previous observations. Through phase-restricted spectral analysis, we found a hint ($\sim 2.6 \sigma$) of spectral hardening around pulsar inferior conjunction. We also provide a uniform analysis of the 12 BWPs observed with $Chandra$ and compare their X-ray properties. Pulsars with soft, $\Gamma > 2.5$, emission seem to have lower than average X-ray and $\gamma$-ray luminosities. We do not, however, see any other prominent correlation between the pulsar's X-ray emission characteristics and any of its other properties. The contribution of the intra-binary shock to the total X-ray emission, if any, is not discernible in this sample of pulsars with shallow observations.

We report on $XMM-Newton$ EPIC observations of the young pulsar J2022+3842, with a characteristic age of 8.9 kyr. We detected X-ray pulsations and found the pulsation period $P\approx 48.6$ ms, and its derivative $\dot{P}\approx 8.6\times 10^{-14}$, twice larger than the previously reported values. The pulsar exhibits two very narrow (FWHM $\sim 1.2$ ms) X-ray pulses each rotation, separated by $\approx 0.48$ of the period, with a pulsed fraction of $\approx 0.8$. Using the correct values of $P$ and $\dot{P}$, we calculate the pulsar's spin-down power $\dot{E}=3.0 \times 10^{37}$ erg s$^{-1}$ and magnetic field $B=2.1\times 10^{12}$ G. The pulsar spectrum is well modeled with a hard power-law (PL) model (photon index $\Gamma = 0.9\pm0.1 $, hydrogen column density $n_H = (2.3\pm0.3) \times 10^{22}\,{\rm cm}^{-2}$). We detect a weak off-pulse emission which can be modeled with a softer PL ($\Gamma \approx 1.7\pm0.7$), poorly constrained because of contamination in the EPIC-pn timing mode data. The pulsar's X-ray efficiency in the $0.5-8$ keV energy band, $\eta_{\rm PSR}= L_{\rm PSR}/\dot{E} = 2 \times 10^{-4} (D/10\,{\rm kpc})^2$, is similar to those of other pulsars. The $XMM-Newton$ observation did not detect extended emission around the pulsar. Our re-analysis of $Chandra$ X-ray observatory archival data shows a hard, $\Gamma \approx 0.9 \pm 0.5$, spectrum and a low efficiency, $\eta_{\rm PWN}\sim 2\times 10^{-5} (D/10\,{\rm kpc})^2$, for the compact pulsar wind nebula, unresolved in the $XMM-Newton$ images.

We report on an X-ray observation of the 166 Myr old radio pulsar J0108-1431 with XMM-Newton. The X-ray spectrum can be described by a power-law model with a relatively steep photon index Gamma~3 or by a combination of thermal and non-thermal components, e.g., a power-law component with fixed photon index Gamma~2 plus a blackbody component with a temperature of kT=0.11 keV. The two-component model appears more reasonable considering different estimates for the hydrogen column density. The non-thermal X-ray efficiency in the single power-law model is eta^PL (1-10 keV) = L^PL (1-10 keV) / Edot ~ 0.003, higher than in most other X-ray detected pulsars. In the case of the combined model, the non-thermal and thermal X-ray efficiencies are even higher, eta^PL (1-10 keV) ~ eta^bb ~ 0.006. We detected X-ray pulsations at the radio period of P=0.808s with significance of 7sigma. The pulse shape in the folded X-ray lightcurve (0.15-2 keV) is asymmetric, with statistically significant contributions from up to 5 leading harmonics. Pulse profiles at two different energy ranges differ slightly: the profile is asymmetric at low energies, 0.15-1 keV, while at higher energies, 1-2 keV, it has a nearly sinusodial shape. The radio pulse peak leads the 0.15-2 keV X-ray pulse peak by delta phi = 0.06 +/- 0.03.