Supernova fallback disks around neutron stars have been discussed to influence the evolution of the diverse neutron star populations. Slowly rotating neutron stars are most promising to find such disks. Searching for the cold and warm debris of old f
allback disks, we carried out Herschel PACS (70 $mu$m, 160 $mu$m) and Spitzer IRAC (3.6 $mu$m, 4.5 $mu$m) observations of eight slowly rotating ($Papprox 3 - 11$ s) nearby ($<1$ kpc) isolated neutron stars. Herschel detected 160 $mu$m emission ($>5sigma$) at locations consistent with the positions of the neutron stars RX J0806.4-4123 and RX J2143.0+0654. No other significant infrared emission was detected from the eight neutron stars. We estimate probabilities of 63%, 33% and 3% that, respectively, none, one, or both Herschel PACS 160 $mu$m detections are unrelated excess sources due to background source confusion or an interstellar cirrus. If the 160 $mu$m emission is indeed related to cold (10 K to 22 K) dust around the neutron stars, this dust is absorbing and re-emitting $sim 10$% to $sim 20$% of the neutron stars X-rays. Such high efficiencies would be at least three orders of magnitude larger than the efficiencies of debris disks around nondegenerate stars. While thin dusty disks around the neutron stars can be excluded as counterparts of the 160 $mu$m emission, dusty asteroid belts constitute a viable option.
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 $Papprox 48.6$ ms, and its derivative $dot{P}approx 8.6times 10^{-14}$, twic
e 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 pulsars spin-down power $dot{E}=3.0 times 10^{37}$ erg s$^{-1}$ and magnetic field $B=2.1times 10^{12}$ G. The pulsar spectrum is well modeled with a hard power-law (PL) model (photon index $Gamma = 0.9pm0.1 $, hydrogen column density $n_H = (2.3pm0.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.7pm0.7$), poorly constrained because of contamination in the EPIC-pn timing mode data. The pulsars 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 2times 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 compone
nts, 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.
We present 35 ks Chandra ACIS observations of the 42 Myr old radio pulsar PSR B1451-68. A point source is detected 0.32 +/- 0.73 from the expected radio pulsar position. It has ~200 counts in the 0.3-8 keV energy range. We identify this point source
as the X-ray counterpart of the radio pulsar. PSR B1451-68 is located close to a 2MASS point source, for which we derive 7% as the upper limit on the flux contribution to the measured pulsar X-ray flux. The pulsar spectrum can be described by either a power-law model with photon index Gamma=2.4 (+0.4/-0.3) and a unrealistically high absorbing column density N(H)= (2.5 (+1.2/-1.3)) * 10^(21) cm^-2, or by a combination of a kT=0.35 (+0.12/-0.07) keV blackbody and a Gamma = 1.4 +/- 0.5 power-law component for N(H)[DM]= 2.6 * 10^(20) cm^-2, estimated from the pulsar dispersion measure. At the parallactic, Lutz-Kelker bias corrected distance of 480 pc, the non-thermal X-ray luminosities in the 0.3-8 keV energy band are either Lx(nonth)= (11.3 +/- 1.7) * 10^(29) erg/s or Lx(nonth)= (5.9 (+4.9/-5.0)) * 10^(29) erg/s, respectively. This corresponds to non-thermal X-ray efficiencies of either eta(nonth)= Lx(nonth) / (dE/dt) ~ 0.005 or 0.003, respectively.
All the neutron star (NS) atmosphere models published so far have been calculated in the cold plasma approximation, which neglects the relativistic effects in the radiative processes, such as cyclotron emission/absorption at harmonics of cyclotron fr
equency. Here we present new NS atmosphere models which include such effects. We calculate a set of models for effective temperatures T_eff =1-3 MK and magnetic fields B sim 10^{10}-10^{11} G, typical for the so-called central compact objects (CCOs) in supernova remnants, for which the electron cyclotron energy E_{c,e} and its first harmonics are in the observable soft X-ray range. Although the relativistic parameters, such as kT_eff /(m_e c^2) and E_{c,e} /(m_e c^2), are very small for CCOs, the relativistic effects substantially change the emergent spectra at the cyclotron resonances, E approx sE_{c,e} (s=1, 2,...). Although the cyclotron absorption features can form in a cold plasma due to the quantum oscillations of the free-free opacity, the shape and depth of these features change substantially if the relativistic effects are included. In particular, the features acquire deep Doppler cores, in which the angular distribution of the emergent intensity is quite different from that in the cold plasma approximation. The relative contributions of the Doppler cores to the equivalent widths of the features grow with increasing the quantization parameter b_eff = E_{c,e}/kT_eff and harmonic number s. The total equivalent widths of the features can reach sim 150-250 eV; they increase with growing b_eff and are smaller for higher harmonics.
Synchrotron radiation of ultra-relativistic particles accelerated in a pulsar wind nebula may dominate its spectrum up to gamma-ray energies. Because of the short cooling time of the gamma-ray emitting electrons, the gamma-ray emission zone is in the
immediate vicinity of the acceleration site. The particle acceleration likely occurs at the termination shock of the relativistic striped wind, where multiple forced magnetic field reconnections provide strong magnetic fluctuations facilitating Fermi acceleration processes. The acceleration mechanisms imply the presence of stochastic magnetic fields in the particle acceleration region, which cause stochastic variability of the synchrotron emission. This variability is particularly strong in the steep gamma-ray tail of the spectrum, where modest fluctuations of the magnetic field lead to strong flares of spectral flux. In particular, stochastic variations of magnetic field, which may lead to quasi-cyclic gamma-ray flares, can be produced by the relativistic cyclotron ion instability at the termination shock. Our model calculations of the spectral and temporal evolution of synchrotron emission in the spectral cut-off regime demonstrate that the intermittent magnetic field concentrations dominate the gamma-ray emission from highest energy electrons and provide fast, strong variability even for a quasi-steady distribution of radiating particles. The simulated light curves and spectra can explain the very strong gamma-ray flares observed in the Crab nebula and the lack of strong variations at other wavelengths. The model predicts high polarization in the flare phase, which can be tested with future polarimetry observations.
Candidate supernova remnants G23.5+0.1 and G25.5+0.0 were observed by XMM-Newton in the course of a snap-shot survey of plerionic and composite SNRs in the Galactic plane. In the field of G23.5+0.1, we detected an extended source, ~3 in diameter, whi
ch we tentatively interpret as a pulsar-wind nebula (PWN) of the middle-aged radio pulsar B1830-08. Our analysis suggests an association between PSR B1830-08 and the surrounding diffuse radio emission. If the radio emission is due to the SNR, then the pulsar must be significantly younger than its characteristic age. Alternatively, the radio emission may come from a relic PWN. In the field of G25.5+0.0, which contains the extended TeV source HESS J1837-069, we detected the recently discovered young high-energy pulsar J1838-0655 embedded in a PWN with extent of 1.3. We also detected another PWN candidate (AX J1837.3-0652) with an extent of 2 and unabsorbed luminosity L_(2-10 keV) ~ 4 x 10^33 erg/s at d=7 kpc. The third X-ray source, located within the extent of the HESS J1837-069, has a peculiar extended radio counterpart, possibly a radio galaxy with a double nucleus or a microquasar. We did not find any evidence of the SNR emission in the G25.5+0.0 field. We provide detailed multiwavelength analysis and identifications of other field sources and discuss robustness of the G25.5+0.0 and G23.5+0.1 classifications as SNRs. (abstract abridged)
Previous observations of the middle-aged pulsar Geminga with XMM-Newton and Chandra have shown an unusual pulsar wind nebula (PWN), with a 20 long central (axial) tail directed opposite to the pulsars proper motion and two 2 long, bent lateral (outer
) tails. Here we report on a deeper (78 ks) Chandra observation and a few additional XMM-Newton observations of the Geminga PWN. The new Chandra observation has shown that the axial tail, which includes up to three brighter blobs, extends at least 50 (i.e., 0.06 d_{250} pc) from the pulsar. It also allowed us to image the patchy outer tails and the emission in the immediate vicinity of the pulsar with high resolution. The PWN luminosity, L_{0.3-8 keV} ~ 3times 10^{29} d_{250}^2 erg/s, is lower than the pulsars magnetospheric luminosity by a factor of 10. The spectra of the PWN elements are rather hard (photon index ~ 1). Comparing the two Chandra images, we found evidence of PWN variability, including possible motion of the blobs along the axial tail. The X-ray PWN is the synchrotron radiation from relativistic particles of the pulsar wind; its morphology is connected with the supersonic motion of Geminga. We speculate that the outer tails are either (1) a sky projection of the limb-brightened boundary of a shell formed in the region of contact discontinuity, where the wind bulk flow is decelerated by shear instability, or (2) polar outflows from the pulsar bent by the ram pressure from the ISM. In the former case, the axial tail may be a jet emanating along the pulsars spin axis, perhaps aligned with the direction of motion. In the latter case, the axial tail may be the shocked pulsar wind collimated by the ram pressure.
Some isolated neutron stars show harmonically spaced absorption features in their thermal soft X-ray spectra. The interpretation of the features as a cyclotron line and its harmonics has been suggested, but the usual explanation of the harmonics as c
aused by relativistic effects fails because the relativistic corrections are extremely small in this case. We suggest that the features correspond to the peaks in the energy dependence of the free-free opacity in a quantizing magnetic field, known as quantum oscillations. The peaks arise when the transitions to new Landau levels become allowed with increasing the photon energy; they are strongly enhanced by the square-root singularities in the phase-space density of quantum states in the case when the free (non-quantized) motion is effectively one-dimensional. To explore observable properties of these quantum oscillations, we calculate models of hydrogen neutron star atmospheres with B sim 10^{10} - 10^{11} G (i.e., electron cyclotron energy E_{c,e} = 0.1 - 1 keV) and T_{eff} = 1 - 3 MK. Such conditions are thought to be typical for the so-called central compact objects in supernova remnants, such as 1E 1207.4-5209 in PKS 1209-51/52. We show that observable features at the electron cyclotron harmonics form at moderately large values of the quantization parameter, b_{eff} = E_{c,e}/kT_{eff} = 0.5 - 20. The equivalent widths of the features can reach 100 - 200 eV; they grow with increasing b_{eff} and are lower for higher harmonics.
PSR J0108-1431 is a nearby, 170 Myr old, very faint radio pulsar near the pulsar death line in the P-Pdot diagram. We observed the pulsar field with the Chandra X-ray Observatory and detected a point source (53 counts in a 30 ks exposure, energy flux
(9+/-2)times 10^{-15} ergs cm^{-2} s^{-1} in the 0.3-8 keV band) close to the radio pulsar position. Based on the large X-ray/optical flux ratio at the X-ray source position, we conclude that the source is the X-ray counterpart of PSR J0108-1431.The pulsar spectrum can be described by a power-law model with photon index Gamma approx 2.2 and luminosity L_{0.3-8 keV} sim 2times 10^{28} d_{130}^2 ergs s^{-1}, or by a blackbody model with the temperature kTapprox 0.28 keV and bolometric luminosity L_{bol} sim 1.3times 10^{28} d_{130}^2 ergs s^{-1}, for a plausible hydrogen column density NH = 7.3times 10^{19} cm^{-2} (d_{130}=d/130 pc). The pulsar converts sim 0.4% of its spin-down power into the X-ray luminosity, i.e., its X-ray efficiency is higher than for most younger pulsars. From the comparison of the X-ray position with the previously measured radio positions, we estimated the pulsar proper motion of 0.2 arcsec yr^{-1} (V_perp sim 130 d_{130} km s^{-1}), in the south-southeast direction.
