We report on the results of Chandra X-ray observations of the $\gamma$-ray radio-quiet pulsar J1932+1916. We confirm the previous detection of the pulsar counterpart and its pulsar wind nebula in X-rays from low spatial resolution data obtained by Suzaku and Swift. The Chandra data with much better spatial resolution resolved the fine structure of the nebula in the pulsar vicinity, which can be interpreted as jets bent by the ram pressure. The size of this compact part is about 30'', and it is surrounded by a weak asymmetric diffusive emission extended up to $\approx$ 2.8'. The X-ray spectra of the pulsar, the compact and extended parts of the nebula can be described by the power-law models with photon indexes of 0.44, 2.34 and 2.14, respectively. The actual pulsar's X-ray flux is several times weaker than it was obtained with Suzaku and Swift where it was dominated by the unresolved nebula. We discuss possible associations of J1932+1916 with the nearby supernova remnant G54.4-0.3.

We have performed a series of Chandra X-ray Observatory observations of the gamma-ray binary LS 2883, which is comprised of a young pulsar (PSR B1259-63) orbiting a massive Be star with a period of 1236.7 days. The system was observed in 5 epochs, spanning a range from 352 to 1175 days after the periastron passage on 2014 May 4. The observations confirmed the recurrent nature of the high-speed ejecta that appear as an extended X-ray structure (clump) moving away from the binary. Compared to the results of the previous monitoring campaign (between the 2010 and 2014 periastron passages), this time we find evidence suggesting that the clump is accelerated to a projected velocity $ v_{\perp}\approx0.15c$ with an acceleration $a_{\perp}=47\pm2$ cm s$^{-2}$ (for uniformly accelerated motion), assuming that it was launched near periastron passage. The observed X-ray properties of the clump are consistent with synchrotron emission from pulsar wind particles accelerated at the interface between the pulsar wind and the clump. We have also performed contemporaneous observations with the {\sl Hubble Space Telescope}, which are used to set an upper limit on the optical flux of the extended emission.

We consider acceleration of leptons up to GeV-TeV energies in the bow shock wind nebula of PSR J0437-4715 and their subsequent diffusion through the interstellar magnetic fields. The leptons accelerated at the pulsar wind termination surface are injected into re-acceleration in colliding shock flows. Modelled spectra of synchrotron emission from the accelerated electrons and positrons are consistent with the far-ultraviolet and X-ray observations of the nebula carried out with the Hubble Space Telescope and Chandra X-ray Observatory. These observations are employed to constrain the absolute fluxes of relativistic leptons, which are escaping from the nebula and eventually reaching the Solar System after energy-dependent diffusion through the local interstellar medium accompanied by synchrotron and Compton losses. It is shown that accelerated leptons from the nebula of PSR J0437-4715 can be responsible both for the enhancement of the positron fraction above a few GeV detected by PAMELA and AMS-02 spectrometers and for the TeV range lepton fluxes observed with H.E.S.S., VERITAS, Fermi, CALET, and DAMPE.

We report non-detections of the $\sim 3\times 10^8$ yr old, slow, isolated, rotation-powered pulsar PSR J2144$-$3933 in observations with the Hubble Space Telescope in one optical band (F475X) and two far-ultraviolet bands (F125LP and F140LP), yielding upper bounds $F_{\rm F475X}< 22.7$ nJy, $F_{\rm F125LP}< 5.9$ nJy, $F_{\rm F140LP}< 19.5$ nJy, at the pivot wavelengths 4940 \AA, 1438 \AA\ and 1528 \AA, respectively. Assuming a blackbody spectrum, we deduce a conservative upper bound on the surface (unredshifted) temperature of the pulsar of $T<42,000$ K. This makes PSR~J2144--3933 the coldest known neutron star, allowing us to study thermal evolution models of old neutron stars. This temperature is consistent with models with either direct or modified Urca reactions including rotochemical heating, and, considering frictional heating from the motion of neutron vortex lines, it puts an upper bound on the excess angular momentum in the neutron superfluid, $J<10^{44}\,\mathrm{erg\,s}$.

Observing the famous high-mass, eccentric X-ray and gamma-ray binary PSR B1259-63/LS 2883 with Chandra, we detected X-ray emitting clumps moving from the binary with speeds of about 0.1 of the speed of light, possibly with acceleration. The clumps are being ejected at least once per binary period, 3.4 years, presumably around binary periastrons. The power-law spectra of the clumps can be interpreted as synchrotron emission of relativistic electrons. Here we report the results of 8 observations of the clumps in 2011-2017 (two binary cycles) and discuss possible interpretations of this unique phenomenon.

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$.

Following up on a faint detection of a near-infrared (NIR) source at the position of the X-ray thermal isolated neutron star RX J0806.4-4123, we present new Hubble Space Telescope observations in the H-band. The NIR source is unambiguously detected with a Vega magnitude of 23.7 +/- 0.2 (flux density of 0.40 +/- 0.06 microJy at lambda =1.54 microm. The source position is coincident with the neutron star position, and the implied NIR flux is strongly in excess of what one would expect from an extrapolation of the optical-UV spectrum of RX J0806.4-4123. The NIR source is extended with a size of at least 0.8arcsec and shows some asymmetry. The conservative upper limit on the flux contribution of a point source is 50%. Emission from gas and dust in the ambient diffuse interstellar medium can be excluded as cause for the extended emission. The source parameters are consistent with an interpretation as either the first NIR-only detected pulsar wind nebula or the first resolved disk around an isolated neutron star.

PSR J1809-1917 is a young ($\tau=51$ kyr) energetic ($\dot{E}=1.8\times10^{36}$ erg s$^{-1}$) radio pulsar powering a pulsar wind nebula (PWN). We report on the results of three Chandra X-ray Observatory observations which show that the PWN consists of a small ($\sim 20''$) bright compact nebula (CN) and faint extended emission seen up to $2'$ from the pulsar. The CN is elongated in the northeast-southwest direction and exhibits morphological and flux variability on a timescale of a few months. We also find evidence of small arcsecond-scale jets extending from the pulsar along the same direction, and exhibiting a hard power-law (PL) spectrum with photon index $\Gamma_{\rm jet}=1.2\pm0.1$. The more extended emission and CN share the same symmetry axis, which is also aligned with the direction toward the TeV $\gamma$-ray source HESS J1809--193, supporting their association. The spectrum of the extended nebula (EN) fits an absorbed PL with about the same slope as that of the CN, $\Gamma_{\rm CN}\approx\Gamma_{\rm EN}=1.55\pm0.09$; no spectral changes across the EN's 2 pc extent are seen. The total PWN 0.5-8 keV luminosity is $L_{\rm PWN}\approx 9\times10^{32}$ erg s$^{-1}$, about half of which is due to the EN.

The Central Compact Object (CCO) in the Cassiopeia A supernova remnant is most likely a very young ($\approx 300$ yr) neutron star. If a previously reported decrease of its surface temperature by 4% in 10 years could be confirmed, it would have profound theoretical implications for neutron star physics. However, the temperature decrease was inferred from Chandra ACIS data affected by instrumental effects which could cause time-dependent spectral distortions. Employing a different instrument setup which minimizes spectral distortions, our 2006 and 2012 Chandra spectra of the CCO did not show a statistically significant temperature decrease. Here, we present additional observations from 2015 taken in the same instrument mode. During the time span of 8.5 years, we detect no significant temperature decrease, using either carbon or hydrogen atmosphere models in the X-ray spectral fits. Our conservative $3\sigma$ upper limits correspond to $<3.3$% and $<2.4$% temperature decrease in 10 years for carbon atmosphere model fits with varying or constant values of the absorbing hydrogen column density, respectively. The recently revised model for the ACIS filter contaminant has a strong effect on the fit results, reducing the significance of the previously reported temperature and flux changes. We expect that a further improved contaminant model and longer time coverage can significantly lower the upper limits in the future.

The young energetic pulsar J1747-2958 ($\tau=26$ kyr, $\dot{E}=2.5\times 10^{36}$ erg s$^{-1}$) powers the Mouse pulsar wind nebula (PWN), famous for its spectacular tail spanning 45$"$ in X-rays and 12$'$ in radio ($d\sim5$ kpc). We present the results of Chandra observations of the PWN and the analysis of archival lower-frequency data. The Chandra HRC image reveals a point-like source at the pulsar position, $\approx$1$'$ behind the bow shock apex of the PWN. The flattened appearance of the compact nebula is consistent with an equatorial outflow deformed by the ram pressure, implying that the angle between the pulsar's spin axis and line of sight is $\sim70^\circ$ (in agreement with the radio and $\gamma$-ray pulse profiles). The spatially-resolved spectroscopy with Chandra ACIS shows that the power-law (PL) spectrum steepens from $\Gamma=1.65\pm0.02$ to $3.0\pm0.1$ over the 45$'$ extent of the X-ray tail. However, the tail's X-ray spectrum integrated over its 45$"$ length fits a single absorbed PL with $\Gamma=2.09\pm0.03$. We also found the Mouse PWN in 150 MHz GMRT data, and a possible counterpart in 24 ${\rm \mu}$m Spitzer data. The multiwavelength data suggest that, at low frequencies, the spectrum of the X-ray-emitting tail region can be described by a broken PL with at least one turnover between radio and X-rays. This is consistent with synchrotron cooling of electrons injected at the termination shock (with an SED slope of 2.2) in an equipartition magnetic field $B\sim200$ ${\rm \mu G}$ and a bulk flow speed $v\sim4000$ km s$^{-1}$.

We discuss the observational properties of pulsar wind nebulae (PWNe) linking them to the injected (at the termination shock) electron spectral energy distribution and parameters of pulsar magnetospheres. In particular, we (1) present spatially-resolved Chandra ACIS spectral maps of twelve PWNe and measure the slopes of the uncooled PWN spectra just downstream of the termination shock obtained from these maps, (2) consider the connections between PWN morphologies and predictions of the magnetospheric emission models and (3) discuss the limits on the maximum energies of particles in PWNe from X-ray and gamma-ray observations.

We report on a Hubble Space Telescope detection of the nearby, old pulsar B0950+08 ($d\simeq 262$ pc, spin-down age 17.5 Myr) in two far-ultraviolet (FUV) bands. We measured the mean flux densities $\bar{f}_\nu = 109\pm 6$ nJy and $83\pm 14$ nJy in the F125LP and F140LP filters (pivot wavelengths 1438 and 1528 \AA). Using the FUV data together with previously obtained optical-UV data, we conclude that the optical-FUV spectrum consists of two components -- a nonthermal (presumably magnetospheric) power-law spectrum ($f_\nu\propto \nu^\alpha$) with slope $\alpha\sim -1.2$ and a thermal spectrum emitted from the bulk of the neutron star surface with a temperature in the range of $(1-3)\times 10^5$ K, depending on interstellar extinction and neutron star radius. These temperatures are much higher than predicted by neutron star cooling models for such an old pulsar, which means that some heating mechanisms operate in neutron stars. A plausible mechanism responsible for the high temperature of PSR B0950+08 is the interaction of vortex lines of the faster rotating neutron superfluid with the slower rotating normal matter in the inner neutron star crust (vortex creep heating).

We review multiwavelength properties of pulsar wind nebulae (PWNe) created by supersonically moving pulsars and the effects of pulsar motion on the PWN morphologies and the ambient medium. Supersonic pulsar wind nebulae (SPWNe) are characterized by bow-shaped shocks around the pulsar and/or cometary tails filled with the shocked pulsar wind. In the past several years significant advances in SPWN studies have been made in deep observations with the Chandra and XMM-Newton X-ray Observatories as well as the Hubble Space Telescope. In particular, these observations have revealed very diverse SPWN morphologies in the pulsar vicinity, different spectral behaviors of long pulsar tails, the presence of puzzling outflows misaligned with the pulsar velocity, and far-UV bow shocks. Here we review the current observational status focusing on recent developments and their implications.

A charged particle in a magnetic field possesses discrete energy levels associated with particle's rotation around the field lines. The radiative transitions between these levels are the well-known cyclotron transitions. We show that a bound complex of particles with a non-zero net charge displays analogous transitions between the states of confined motion of the entire complex in the field. The latter {\em bound-ion cyclotron transitions} are affected by a coupling between the collective and internal motions of the complex and, as a result, differ from the transitions of a "reference" bare ion with the same mass and charge. We analyze the cyclotron transitions for complex ions by including the coupling within a rigorous quantum approach. Particular attention is paid to comparison of the transition energies and oscillator strengths to those of the bare ion. Selection rules based on integrals of collective motion are derived for the bound-ion cyclotron transitions analytically, and the perturbation and coupled-channel approaches are developed to study the transitions quantitatively. Representative examples are considered and discussed for positive and negative atomic and cluster ions.

We review observations of several classes of neutron-star-powered outflows: pulsar-wind nebulae (PWNe) inside shell supernova remnants (SNRs), PWNe interacting directly with interstellar medium (ISM), and magnetar-powered outflows. We describe radio, X-ray, and gamma-ray observations of PWNe, focusing first on integrated spectral-energy distributions (SEDs) and global spectral properties. High-resolution X-ray imaging of PWNe shows a bewildering array of morphologies, with jets, trails, and other structures. Several of the 23 so far identified magnetars show evidence for continuous or sporadic emission of material, sometimes associated with giant flares, and a few possible "magnetar-wind nebulae" have been recently identified.

We report the results of three Chandra observations covering most of the extent of the TeV gamma-ray source HESS J1616-508 and a search for a lower energy counterpart to this source. We detect 56 X-ray sources, of which 37 have counterparts at lower frequencies, including a young massive star cluster, but none of them appears to be a particularly promising counterpart to the TeV source. The brightest X-ray source, CXOU J161423.4-505738 with a flux of 5x10^{-13} c.g.s. in the 0.5-7.0 keV band, has a hard spectrum that is well fit by a power-law model with a photon index of 0.2+/-0.3 and is a likely intermediate polar CV candidate. No counterparts of this source were detected at other wavelengths. CVs are not known to produce extended TeV emission, and the source is also largely offset (19') from HESS J1616-508, making them unlikely to be associated. We have also set an upper limit on the X-ray flux of PSR J1614-5048 in the 0.5-8 keV band (F<5x10^{-15} c.g.s. at a 90% confidence level). This makes PSR J1614-5048 one of the least X-ray efficient pulsars known, with an X-ray efficiency {\eta}=L_x/Edot < 2x10^{-5} in the 0.5-8 keV band. We find no evidence supporting the association between the pulsar and the TeV source. We rule out a number of X-ray sources as possible counterparts to the TeV emission and do not find a plausible counterpart among the other sources. Lastly, we discuss the possible relation of PSR J1617-5055 to HESS J1616-508 in light of the new observations.

We present a detailed analysis of the morphological properties of the Crab Nebula across the electromagnetic spectrum based on new and previous high-quality data from radio to X-rays. In the radio range we obtained an image of the entire nebula at 3 GHz with subarcsecond angular resolution using the VLA (NRAO) and an image at 100 GHz of the central region using the ALMA array. Simultaneously with the VLA observations we performed HST WFPC3 near infrared and Chandra X-ray observations of the central region of the nebula. In addition we produced a new UV image of the Crab nebula at 291 nm by co-adding 75 individual exposures of the Optical-UV Monitor on board XMM-Newton. The high-angular resolution and high-dynamic range radio image at 3 GHz allowed us to improve the detection and characterization of peculiar morphological features including arches with foot brightening and intercrossed loop-like structures, likely originating in plasma confined to magnetic field lines. Based on the new radio image, we carried out a detailed multiwavelength correlation. In the central area, the comparison of the almost simultaneous images confirms that the wisps in the three spectral ranges do not generally coincide in location, the radio emission being the most discordant, which is suggestive of the existence of two different synchrotron components. The X-ray pulsar jet does not have a radio counterpart. Instead, another jet-like feature is seen in radio, though with different curvature and starting point.

We observed a nearby millisecond pulsar J2124-3358 with the Hubble Space Telescope in broad far-UV (FUV) and optical filters. The pulsar is detected in both bands with fluxes F(1250-2000 A)= (2.5+/-0.3)x10^-16 erg/s/cm^2 and F(3800-6000 A)=(6.4+/-0.4)x10^-17 erg/s/cm^2, which correspond to luminosities of ~5.8x10^27 and 1.4x10^27 erg/s, for d=410 pc and E(B-V)=0.03. The optical-FUV spectrum can be described by a power-law model, f_nu~nu^alpha, with slope alpha=0.18-0.48 for a conservative range of color excess, E(B-V)=0.01-0.08. Since a spectral flux rising with frequency is unusual for pulsar magnetospheric emission in this frequency range, it is possible that the spectrum is predominantly magnetospheric (power law with alpha<0) in the optical while it is dominated by thermal emission from the neutron star surface in the FUV. For a neutron star radius of 12 km, the surface temperature would be between 0.5x10^5 and 2.1x10^5 K, for alpha ranging from -1 to 0, E(B-V)=0.01-0.08, and d=340-500 pc. In addition to the pulsar, the FUV images reveal extended emission spatially coincident with the known Halpha bow shock, making PSR J2124-3358 the second pulsar (after PSR J0437-4715) with a bow shock detected in FUV.

We report on six new Chandra observations of the Geminga pulsar wind nebula (PWN). The PWN consists of three distinct elongated structures - two $\approx 0.2 d_{250}$ pc long lateral tails and a segmented axial tail of $\approx 0.05 d_{250}$ pc length, where $d_{250}=d/(250 {\rm pc})$. The photon indices of the power law spectra of the lateral tails, $\Gamma \approx 1$, are significantly harder than those of the pulsar ($\Gamma \approx 1.5$) and the axial tail ($\Gamma \approx 1.6$). There is no significant diffuse X-ray emission between the lateral tails -- the ratio of the X-ray surface brightness between the south tail and this sky area is at least 12. The lateral tails apparently connect directly to the pulsar and show indication of moving footpoints. The axial tail comprises time-variable emission blobs. However, there is no evidence for constant or decelerated outward motion of these blobs. Different physical models are consistent with the observed morphology and spectra of the Geminga PWN. In one scenario, the lateral tails could represent an azimuthally asymmetric shell whose hard emission is caused by the Fermi acceleration mechanism of colliding winds. In another scenario, the lateral tails could be luminous, bent polar outflows, while the blobs in the axial tail could represent a crushed torus. In a resemblance to planetary magnetotails, the blobs of the axial tail might also represent short-lived plasmoids which are formed by magnetic field reconnection in the relativistic plasma of the pulsar wind tail.

We report on Chandra X-ray Observatory (CXO) observations of the pulsar wind nebula (PWN) associated with PSR B0355+54 (eight observations with a 395 ks total exposure, performed over an 8 month period). We investigated the spatial and spectral properties of the emission coincident with the pulsar, compact nebula (CN), and extended tail. We find that the CN morphology can be interpreted in a way that suggests a small angle between the pulsar spin axis and our line-of-sight, as inferred from the radio data. On larger scales, emission from the 7' (2 pc) tail is clearly seen. We also found hints of two faint extensions nearly orthogonal to the direction of the pulsar's proper motion. The spectrum extracted at the pulsar position can be described with an absorbed power-law + blackbody model. The nonthermal component can be attributed to magnetospheric emission, while the thermal component can be attributed to emission from either a hot spot (e.g., a polar cap) or the entire neutron star surface. Surprisingly, the spectrum of the tail shows only a slight hint of cooling with increasing distance from the pulsar. This implies either a low magnetic field with fast flow speed, or particle re-acceleration within the tail. We estimate physical properties of the PWN and compare the morphologies of the CN and the extended tail with those of other bow shock PWNe observed with long CXO exposures.

Pulsars traveling at supersonic speeds are often accompanied by cometary bow shocks seen in Halpha. We report on the first detection of a pulsar bow shock in the far-ultraviolet (FUV). We detected it in FUV images of the nearest millisecond pulsar J0437-4715 obtained with the Hubble Space Telescope. The images reveal a bow-like structure positionally coincident with part of the previously detected Halpha bow shock, with an apex at 10" ahead of the moving pulsar. Its FUV luminosity, L(1250-2000 A) ~ 5x10^28 erg/s, exceeds the Halpha luminosity from the same area by a factor of 10. The FUV emission could be produced by the shocked ISM matter or, less likely, by relativistic pulsar wind electrons confined by strong magnetic field fluctuations in the bow shock. In addition, in the FUV images we found a puzzling extended (~3" in size) structure overlapping with the limb of the bow shock. If related to the bow shock, it could be produced by an inhomogeneity in the ambient medium or an instability in the bow shock. We also report on a previously undetected X-ray emission extending for about 5" ahead of the pulsar, possibly a pulsar wind nebula created by shocked pulsar wind, with a luminosity L(0.5-8 keV) ~ 3x10^28 erg/s.

PSR J1509-5850 is a middle-aged pulsar with the period P ~ 89 ms, spin-down power Edot = 5.1 x 10^35 erg/s, at a distance of about 3.8 kpc. We report on deep Chandra X-ray Observatory observations of this pulsar and its pulsar wind nebula (PWN). In addition to the previously detected tail extending up to 7' southwest from the pulsar (the southern outflow), the deep images reveal a similarly long, faint diffuse emission stretched toward the north (the northern outflow) and the fine structure of the compact nebula (CN) in the pulsar vicinity. The CN is resolved into two lateral tails and one axial tail pointing southwest (a morphology remarkably similar to that of the Geminga PWN), which supports the assumption that the pulsar moves towards the northeast. The luminosities of the southern and northern outflows are about 1 x 10^33 and 4 x 10^32 erg/s, respectively. The spectra extracted from four regions of the southern outflow do not show any softening with increasing distance from the pulsar. The lack of synchrotron cooling suggests a high flow speed or in-situ acceleration of particles. The spectra extracted from two regions of the northern outflow show a hint of softening with distance from the pulsar, which may indicate slower particle propagation. We speculate that the northern outflow is associated with particle leakage from the bow shock apex into the ISM, while the southern outflow represents the tail of the shocked pulsar wind behind the moving pulsar. We estimate the physical parameters of the observed outflows and compare the J1509-5850 PWN with PWNe of other supersonically moving pulsars.

We report on the analysis of two deep XMM-Newton observations of the magnetar Swift J1834.9-0846 and its surrounding extended emission taken in March 2014 and October 2014, 2.5 and 3.1 years after the source went into outburst. The magnetar is only weakly detected in the first observation with an absorption corrected flux $F_{\rm 0.5-10 keV}\approx4\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$, and a $3\sigma$ upper limit during the second observation of about $3\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$. This flux level is more than 3 orders of magnitude lower than the flux measured at the outburst onset on September 2011. The extended emission, centered at the magnetar position and elongated towards the south-west, is clearly seen in both observations; it is best fit by a highly absorbed power-law (PL), with a hydrogen column density of $N_{\rm H}=8.0\times10^{22}$ cm$^{-2}$ and PL photon index $\Gamma=2.2\pm0.2$. Its flux is constant between the two observations at $F_{\rm 0.5-10 keV}=1.3\times10^{-12}$ erg s$^{-1}$ cm$^{-2}$. We find no statistically significant changes in the spectral shape or the flux of this extended emission over a period of 9 years from 2005 to 2014. These new results strongly support the extended emission nature as a wind nebula and firmly establish Swift J1834.9-0846 as the first magnetar to show a surrounding wind nebula. Further, our results imply that such nebulae are no longer exclusive to rotation-powered pulsars and narrow the gap between these two sub-populations of isolated neutron stars. The size and spectrum of the nebula are compatible with those of pulsar-wind nebulae but its radiative efficiency $\eta_{\rm X}=L_{\rm X}/\dot{E}\approx0.1$ is markedly high, possibly pointing to an additional wind component in Swift J1834.9-0846.

PSR B0656+14 is a middle-aged pulsar with a characteristic age $\tau_c=110$ kyr and spin-down power $\dot{E}= 3.8\times 10^{34}$ erg s$^{-1}$. Using Chandra data, we searched for a pulsar wind nebula (PWN) and found evidence of extended emission in a 3.5-15 arcsec annulus around the pulsar, with a luminosity $L_{\rm 0.5-8\,keV}^{\rm ext} \sim 8\times 10^{28}$ erg s$^{-1}$ (at the distance of 288 pc), which is a fraction of $\sim 0.05$ of the non-thermal pulsar luminosity. If the extended emission is mostly due to a PWN, its X-ray efficiency, $\eta_{\rm pwn} = L_{\rm 0.5-8\,keV}^{\rm ext}/\dot{E} \sim 2\times 10^{-6}$, is lower than those of most other known PWNe but similar to that of the middle-aged Geminga pulsar. The small radial extent and nearly round shape of the putative PWN can be explained if the pulsar is receding (or approaching) in the direction close to the line of sight. The very soft spectrum of the extended emission ($\Gamma\sim 8$), much softer than those of typical PWNe, could be explained by a contribution from a faint dust scattering halo, which may dominate in the outer part of the extended emission.

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.

The nearby, middle-aged PSR B1055-52 has many properties in common with the Geminga pulsar. Motivated by the Geminga's enigmatic and prominent pulsar wind nebula (PWN), we searched for extended emission around PSR B1055-52 with Chandra ACIS. For an energy range 0.3-1 keV, we found a 4 sigma flux enhancement in a 4.9-20 arcsec annulus around the pulsar. There is a slight asymmetry in the emission close, 1.5-4 arcsec, to the pulsar. The excess emission has a luminosity of about 10^{29} erg s^{-1} in an energy range 0.3-8 keV for a distance of 350 pc. Overall, the faint extended emission around PSR B1055-52 is consistent with a PWN of an aligned rotator moving away from us along the line of sight with supersonic velocity, but a contribution from a dust scattering halo cannot be excluded. Comparing the properties of other nearby, middle-aged pulsars, we suggest that the geometry -- the orientations of rotation axis, magnetic field axis, and the sight-line -- is the deciding factor for a pulsar to show a prominent PWN. For PSR B1055-52, we also report on a flux decrease of at least 30% between the 2000 XMM-Newton and our 2012 Chandra observation. We tentatively attribute this flux decrease to a cross-calibration problem, but further investigations of the pulsar are required to exclude actual intrinsic flux changes.

Deep Chandra ACIS observations of the region around the putative pulsar, CXOU J061705.3+222127, in the supernova remnant IC443 reveal an ~5$^{\prime\prime}$-radius ring-like structure surrounding the pulsar and a jet-like feature oriented roughly north-south across the ring and through the pulsar's location at 06$^{\rm h}$17$^{\rm m}$5.200$^{\rm s}$ +22$^{\circ}$21$^{\prime}$27.52$^{\prime\prime}$ (J2000.0 coordinates). The observations further confirm that (1) the spectrum and flux of the central object are consistent with a rotation-powered pulsar, (2) the non-thermal spectrum and morphology of the surrounding nebula are consistent with a pulsar wind and, (3) the spectrum at greater distances is consistent with thermal emission from the supernova remnant. The cometary shape of the nebula, suggesting motion towards the southwest, appears to be subsonic: There is no evidence either spectrally or morphologically for a bow shock or contact discontinuity; the nearly circular ring is not distorted by motion through the ambient medium; and the shape near the apex of the nebula is narrow. Comparing this observation with previous observations of the same target, we set a 99% confidence upper limit to the proper motion of CXOU J061705.3+222127 to be less than 44 mas/yr (310 km/s for a distance of 1.5 kpc), with the best-fit (but not statistically significant) projected direction toward the west.

We present results from two {\it Chandra} observations of the 16.6 d X-ray/$\gamma$-ray high-mass binary 1FGL J1018.6$-$5856 located at the center of the supernova remnant G284.3$-$1.8. The binary spectra, separated by 0.25 in binary phase, are fit with an absorbed power-law model with $\Gamma\approx1.7-1.8$ for both observations (the flux during the second observation is a factor of 1.7 smaller). In the high-resolution ACIS-I image we found a hint of extended emission $\approx2''$--$3''$ southeast of the binary, significant at the $3\sigma$ level. Binary evolution codes reproduce the system's observed properties with two massive stars with an initial 18-day period, undergoing mass transfer and leaving behind a heavy $\approx2$ $\msun$ neutron star. The initial mass of the progenitor star in this scenario is 27 $\pm 4$ $\msun$. \emph{Chandra} and {\it XMM-Newton} images of the remnant show it has a relatively low X-ray surface brightness. The two brightest regions of extended X-ray emission, with luminosities $\sim10^{33}$ erg s$^{-1}$ for $d=5$ kpc, lie in the northern and western portions and show significantly different spectra. The northern patch is consistent with shocked ISM, with a low temperature and long ionization timescale. However, the western patch is dominated by ejecta, and shows significantly enhanced Mg content relative to other ejecta products. The abundance ratios inferred resemble those from the Large Magellanic Cloud remnant N49B. To our knowledge, this is only the second case of such Mg-rich ejecta found in an SNR. Nucleosynthesis models for core-collapse SNe predict Mg-rich ejecta from very massive progenitors of $>25$ $\msun$.

We present the analysis of the Chandra X-ray Observatory observations of the eccentric gamma-ray binary PSR B1259-63/LS 2883. The analysis shows that the extended X-ray feature seen in previous observations is still moving away from the binary with an average projected velocity of about 0.07c and shows a hint of acceleration. The spectrum of the feature appears to be hard (photon index of 0.8) with no sign of softening compared to previously measured values. We interpret it as a clump of plasma ejected from the binary through the interaction of the pulsar with the decretion disk of the O-star around periastron passage. We suggest that the clump is moving in the unshocked relativistic pulsar wind (PW), which can accelerate the clump. Its X-ray emission can be interpreted as synchrotron radiation of the PW shocked by the collision with the clump.

We obtained six observations of PSR J1741-2054 using the $Chandra$ ACIS-S detector totaling $\sim$300 ks. By registering this new epoch of observations to an archival observation taken 3.2 years earlier using X-ray point sources in the field of view, we have measured the pulsar proper motion at $\mu =109 \pm 10 {\rm mas\ yr}^{-1}$ in a direction consistent with the symmetry axis of the observed H$\alpha$ nebula. We investigated the inferred past trajectory of the pulsar but find no compelling association with OB associations in which the progenitor may have originated. We confirm previous measurements of the pulsar spectrum as an absorbed power law with photon index $\Gamma$=2.68$\pm$0.04, plus a blackbody with an emission radius of (4.5$^{+3.2}_{-2.5})d_{0.38}$ km, for a DM-estimated distance of $0.38d_{0.38}$ kpc and a temperature of $61.7\pm3.0$ eV. Emission from the compact nebula is well described by an absorbed power law model with a photon index of $\Gamma$ = 1.67$\pm$0.06, while the diffuse emission seen as a trail extending northeast of the pulsar shows no evidence of synchrotron cooling. We also applied image deconvolution techniques to search for small-scale structures in the immediate vicinity of the pulsar, but found no conclusive evidence for such structures.

PSR J2021+3651 is a 17 kyr old rotation powered pulsar detected in the radio, X-rays, and $\gamma$-rays. It powers a torus-like pulsar wind nebula with jets, dubbed the Dragonfly, which is very similar to that of the Vela pulsar. The Dragonfly is likely associated with the extended TeV source VER J2019+368 and extended radio emission. We conducted first deep optical observations with the GTC in the Sloan $r'$ band to search for optical counterparts of the pulsar and its nebula. No counterparts were detected down to $r'\gtrsim27.2$ and $\gtrsim24.8$ for the point-like pulsar and the compact X-ray nebula, respectively. We also reanalyzed Chandra archival X-ray data taking into account an interstellar extinction--distance relation, constructed by us for the Dragonfly line of sight using the red-clump stars as standard candles. This allowed us to constrain the distance to the pulsar, $D=1.8^{+1.7}_{-1.4}$ kpc at 90% confidence. It is much smaller than the dispersion measure distance of $\sim$12 kpc but compatible with a $\gamma$-ray "pseudo-distance" of 1 kpc. Based on that and the optical upper limits, we conclude that PSR J2021+3651, similar to the Vela pulsar, is a very inefficient nonthermal emitter in the optical and X-rays, while its $\gamma$-ray efficiency is consistent with an average efficiency for $\gamma$-pulsars of similar age. Our optical flux upper limit for the pulsar is consistent with the long-wavelength extrapolation of its X-ray spectrum while the nebula flux upper limit does not constrain the respective extrapolation.

HESS J1809-193 is an extended TeV gamma-ray source in the Galactic Plane. Multiwavelength observations of the HESS J1809-193 field reveal a complex picture. We present results from three CXO and two Suzaku observations of a region in the northeastern outskirts of HESS J1809-193, where enhanced TeV emission has been reported. Our analysis also includes GeV gamma-ray and radio data. One of the X-ray sources in the field is the X-ray binary XTE J1810-189, for which we present the outburst history from multiple observatories and confirm that XTE J1810-189 is a strongly variable type I X-ray burster, which can hardly produce TeV emission. We investigate whether there is any connection between the possible TeV extension of HESS J1809-193 and the sources seen at lower energies. We find that another X-ray binary candidate, Suzaku J1811-1900, and a radio supernova remnant, SNR G11.4-0.1, can hardly be responsible for the putative TeV emission. Our multiwavelength classification of fainter X-ray point sources also does not produce a plausible candidate. We conclude that the northeast extension of HESS J1809-193, if confirmed by deeper observations, can be considered as a dark accelerator - a TeV source without visible counterpart at lower energies.

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 fallback 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 ($P\approx 3 - 11$ s) nearby ($<1$ kpc) isolated neutron stars. Herschel detected 160 $\mu$m emission ($>5\sigma$) 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 $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 detection of the double pulsar system J0737-3039 in the far-UV with the ACS/SBC detector aboard HST. We measured the energy flux F = 4.5+/-1.0e-17 erg cm-2s-1 in the 1250-1550 AA band, which corresponds to the extinction-corrected luminosity L~1.5e28 erg s-1 for the distance d=1.1 kpc and a plausible reddening E(B-V)=0.1. If the detected emission comes from the entire surface of one of the neutron stars with a 13 km radius, the surface blackbody temperature is in the range T~2-5e5 K for a reasonable range of interstellar extinction. Such a temperature requires an internal heating mechanism to operate in old neutron stars, or it might be explained by heating of the surface of the less energetic Pulsar B by the relativistic wind of Pulsar A. If the far-UV emission is non-thermal (e.g., produced in the magnetosphere of Pulsar A), its spectrum exhibits a break between the UV and X-rays.

We present observations of the eccentric gamma-ray binary B1259-63/LS2883 with the Chandra X-ray Observatory. The images reveal a variable, extended about 4'', or about 1000 times the binary orbit size) structure, which appears to be moving away from the binary with the velocity of 0.05 of the speed of light. The observed emission is interpreted as synchrotron radiation from relativistic particles supplied by the pulsar. However, the fast motion through the circumbinary medium would require the emitting cloud to be loaded with a large amount of baryonic matter. Alternatively, the extended emission can be interpreted as a variable extrabinary shock in the pulsar wind outflow launched near binary apastron. The resolved variable X-ray nebula provides an opportunity to probe pulsar winds and their interaction with stellar winds in a previously inaccessible way.

The accuracy of measuring the basic parameters of neutron stars is limited in particular by uncertainties in chemical composition of their atmospheres. For example, atmospheres of thermally - emitting neutron stars in supernova remnants might have exotic chemical compositions, and for one of them, the neutron star in CasA, a pure carbon atmosphere has recently been suggested by Ho & Heinke (2009). To test such a composition for other similar sources, a publicly available detailed grid of carbon model atmosphere spectra is needed. We have computed such a grid using the standard LTE approximation and assuming that the magnetic field does not exceed 10^8 G. The opacities and pressure ionization effects are calculated using the Opacity Project approach. We describe the properties of our models and investigate the impact of the adopted assumptions and approximations on the emergent spectra.

To examine the previously claimed fast cooling of the Central Compact Object (CCO) in the Cas A supernova remnant (SNR), we analyzed two Chandra observations of this CCO, taken in a setup minimizing instrumental spectral distortions. We fit the two CCO X-ray spectra from 2006 and 2012 with hydrogen and carbon neutron star atmosphere models. The temperature and flux changes in the 5.5 years between the two epochs depend on the adopted constraints on the fitting parameters and the uncertainties of the effective area calibrations. If we allow a change of the equivalent emitting region size, R_Em, the effective temperature remains essentially the same. If R_Em is held constant, the best-fit temperature change is negative, but its statistical significance ranges from 0.8sigma to 2.5sigma, depending on the model. If we assume that the optical depth of the ACIS filter contaminant in 2012 was +/-10% different from its default calibration value, the significance of the temperature drop becomes 0.8sigma to 3.1sigma, for the carbon atmospheres with constant R_Em. Thus, we do not see a statistically significant temperature drop in our data, but the involved uncertainties are too large to firmly exclude the previously reported fast cooling. Our analysis indicate a decrease of 4%-6% (1.9-2.9sigma significance) for the absorbed flux in the energy range 0.6-6keV between 2006 and 2012, most prominent in the 1.4-1.8 keV energy range. It could be caused by unaccounted changes of the detector response or contributions from unresolved SNR material along the line of sight to the CCO.

We review the multiwavelength properties of the few known gamma-ray binaries, focusing on extended emission recently resolved with Chandra. We discuss the implications of these findings for the nature of compact objects and for physical processes operating in these systems.

We review properties of Galactic VHE sources detected at TeV energies. The number of associations between the VHE sources and pulsars has grown in recent years, making pulsar-wind nebulae the dominant population, although there is still a substantial number of VHE sources which remain to be identified. Among the latter there are several "dark" sources which do not have plausible counterparts at any other wavelengths. In this review we compile and compare the TeV and X-ray properties of pulsar wind nebulae (PWNe), PWN candidates, and unidentified TeV sources.

We have studied the fascinating dynamics of the nearby Vela pulsar's nebula in a campaign comprising eleven 40ks observations with Chandra X-ray Observatory (CXO). The deepest yet images revealed the shape, structure, and motion of the 2-arcminute-long pulsar jet. We find that the jet's shape and dynamics are remarkably consistent with that of a steadily turning helix projected on the sky. We discuss possible implications of our results, including free precession of the neutron star and MHD instability scenarios.

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.

We report on radio and X-ray observations of PSR 1832+0029, a 533-ms radio pulsar discovered in the Parkes Multibeam Pulsar Survey. From radio observations taken with the Parkes, Lovell and Arecibo telescopes, we show that this pulsar exhibits two spindown states akin to PSRs B1931+24 reported by Kramer et al. and J1841-0500 reported by Camilo et al. Unlike PSR B1931+24, which switches between "on" and "off" states on a 30--40 day time-scale, PSR 1832+0029 is similar to PSR J1841-0500 in that it spends a much longer period of time in the off-state. So far, we have fully sampled two off-states. The first one lasted between 560 and 640 days and the second one lasted between 810 and 835 days. From our radio timing observations, the ratio of on/off spindown rates is $1.77 \pm 0.03$. Chandra observations carried out during both the on- and off-states of this pulsar failed to detect any emission. Our results challenge but do not rule out models involving accretion onto the neutron star from a low-mass stellar companion. In spite of the small number of intermittent pulsars currently known, difficulties in discovering them and in quantifying their behavior imply that their total population could be substantial.

The vast majority of known non-accreting neutron stars (NSs) are rotation-powered radio and/or gamma-ray pulsars. So far, their multiwavelength spectra have all been described satisfactorily by thermal and non-thermal continuum models, with no spectral lines. Spectral features have, however, been found in a handful of exotic NSs and thought to be a manifestation of their unique traits. Here we report the detection of absorption features in the X-ray spectrum of an ordinary rotation-powered radio pulsar, J1740+1000. Our findings bridge the gap between the spectra of pulsars and other, more exotic, NSs, suggesting that the features are more common in the NS spectra than they have been thought so far.

We report on the analysis of two XMM-Newton observations of the recently discovered soft gamma repeater Swift J1834.9-0846, taken in September 2005 and one month after the source went into outburst on 2011 August 7. We performed timing and spectral analyses on the point source as well as on the extended emission. We find that the source period is consistent with an extrapolation of the Chandra ephemeris reported earlier and the spectral properties remained constant. The source luminosity decreased to a level of 1.6x10^34 erg s^-1 following a decay trend of $\propto t^{-0.5}$. Our spatial analysis of the source environment revealed the presence of two extended emission regions around the source. The first (Region A) is a symmetric ring around the point source, starting at 25arcsec and extending to ~50arcsec. We argue that Region A is a dust scattering halo. The second (Region B) has an asymmetrical shape extending between 50arcsec and 150arcsec, and is detected both in the pre- and post-outburst data. We argue that this region is a possible magnetar wind nebula (MWN). The X-ray efficiency of the MWN with respect to the rotation energy loss is substantially higher than those of rotation powered pulsars: $\eta_{\rm X}\equiv L_{\rm MWN,0.5-8 keV}/\dot{E}_{\rm rot}\approx0.7$. The higher efficiency points to a different energy source for the MWN of Swift J1834.9-0846, most likely bursting activity of the magnetar, powered by its high magnetic field, B=1.4x10^14 G.

The number of plausible associations of extended VHE (TeV) sources with pulsars has been steadily growing, suggesting that many of these sources are pulsar wind nebulae (PWNe). Here we overview the recent progress in X-ray and TeV observations of PWNe and summarize their properties.

PSR J1357-6429 is a Vela-like radio pulsar that has been recently detected with Chandra and Fermi, which, like Vela, powers a compact X-ray pulsar wind nebula and X-ray-radio plerion associated with an extended TeV source. We present our deep optical observations with the Very Large Telescope to search for an optical counterpart of the pulsar and its nebula. We detected a point-like source in V, R, and I bands whose position is in agreement with the X-ray position of the pulsar, and whose colours are distinct from those of ordinary stars. The tentative optical luminosity and efficiency of the source are similar to those of the Vela pulsar, which also supports the optical identification. However, the source spectrum is unusually steep, with a spectral index of about 5, which is not typical of optical pulsars. The source offset from the radio position of PSR J1357-6429, which is in line with the corresponding offset of the X-ray position, implies the pulsar transverse velocity of 1600-2000 km/s at the distance of 2-2.5 kpc, making it the fastest moving pulsar known.

(Abridged) PSR J1357-6429 is a Vela-like radio pulsar that has been recently detected in X-rays and gamma-rays. It powers a compact tail-like X-ray pulsar wind nebula and X-ray-radio plerion associated with an extended TeV source HESS J1356-645. We present our deep optical observations with the Very Large Telescope to search for an optical counterpart of the pulsar and its nebula. We detect a point-like source in V, R, and I bands whose position is within the 1-sigma error circle of the X-ray position of the pulsar, and whose colours are distinct from those of ordinary stars. We consider it as a candidate optical counterpart of the pulsar. If it is indeed the counterpart, its 5-sigma offset from the radio pulsar position, measured about 9 yr earlier, implies that the transverse velocity of the pulsar is in the range of 1600--2000 km s^{-1} at the distance of 2--2.5 kpc, making it the fastest moving pulsar known. The direction of the estimated proper motion coincides with the extension of the pulsar's X-ray tail, suggesting that this is a jet. The tentative optical luminosity and efficiency of the pulsar are similar to those of the Vela pulsar, which also supports the optical identification. However, the candidate undergoes an unusually steep dereddened flux increase towards the infrared with a spectral index of about 5, that is not typical of optical pulsars. It implies a strong double-knee spectral break in the pulsar emission between the optical and X-rays. The reasons for the spectral steepness are unclear. It may be caused by a nebula knot projected onto the jet and strongly overlapping with the pulsar, as observed for the Crab, where the knot has a significantly steeper spectrum than the pulsar. We find no other signs of the pulsar nebula in the optical. Alternatively, the detected source may be a faint AGN, that has not yet been seen at other wavelengths.

Recent X-ray observations of the supernova remnant IC443 interacting with molecular clouds have shown the presence of a new population of hard X-ray sources related to the remnant itself, which has been interpreted in terms of fast ejecta fragment propagating inside the dense environment. Prompted by these studies, we have obtained a deep {\sl XMM-Newton} observation of the supernova remnant (SNR) Kes 69, which also shows signs of shock-cloud interaction. We report on the detection of 18 hard X-ray sources in the field of Kes 69, a significant excess of the expected galactic source population in the field, spatially correlated with CO emission from the cloud in the remnant environment. The spectra of 3 of the 18 sources can be described as hard power laws with photon index <2 plus line emission associated to K-shell transitions. We discuss the two most promising scenarios for the interpretation of the sources, namely fast ejecta fragments (as in IC443) and cataclysmic variables. While most of the observational evidences are consistent with the former interpretation, we cannot rule out the latter.

Using high-quality Hubble Space Telescope observations, we construct the near infra-red (NIR) to far ultra-violet (FUV) spectral energy distribution (SED) of PSR B0656+14. The SED is non-monotonic. Fitting it with a simple combination of a Rayleigh-Jeans spectrum (UV) and non-thermal power-law (optical/NIR) leaves significant residuals, strongly hinting at one or more spectral features. We consider various models (combination of continuum components, and absorption/emission lines) with possible interpretations, and place them in the context of the broader spectral energy distribution. Surprisingly, the extrapolation of the best-fit X-ray spectral model roughly match the NIR-FUV data, and the power-law component is also consistent with the gamma-ray fluxes. We compare the multiwavelength SED of B0656+14 with those of other optical, X-ray and gamma-ray detected pulsars, and notice that a simple power-law spectrum crudely accounts for most of the non-thermal emission.

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.

Taking advantage of the high sensitivity of the Chandra Advanced CCD Imaging Spectrometer, we have conducted a snap-shot survey of pulsars previously undetected in X-rays. We detected 12 pulsars and established deep flux limits for 11 pulsars. Using these new results, we revisit the relationship between the X-ray luminosity, L_psr_x, and spin-down power, Edot. We find that the obtained limits further increase the extremely large spread in the non-thermal X-ray efficiencies, eta_psr_x=L_psr_x/Edot, with some of them being now below 1e-5. Such a spread cannot be explained by poorly known distances or by beaming of pulsar radiation. We also find evidence of a break in the dependence of L_psr_x on Edot, such that pulsars become more X-ray efficient at Edot<~ 1e34-1e35 erg/s. We examine the relationship between the gamma-ray luminosity, L_psr_g, and Edot, which exhibits a smaller scatter compared to that in X-rays. This confirms that the very large spread in the X-ray efficiencies cannot be explained just by beaming because the gamma-ray emission is generally expected to be beamed stronger than the X-ray emission. Intriguingly, there is also an indication of a break in the L_psr_g(Edot) dependence at Edot~1e35 erg/s, with lower-Edot pulsars becoming less gamma-ray efficient. We also examine the distance-independent L_psr_f/L_psr_x ratio as a function of Edot for a sample of gamma-ray pulsars observed by Chandra and find that it peaks at Edot~1e35 erg/s, showing that the breaks cannot originate from poorly measured distances. We discuss the implications of our findings for existing models of magnetospheric emission and venues for further exploration.

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 frequency. 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.

We present X-ray observations of the new transient magnetar Swift J1834.9-0846, discovered with Swift BAT on 2011 August 7. The data were obtained with Swift, RXTE, CXO, and XMM-Newton both before and after the outburst. Timing analysis reveals singe peak pulsations with a period of 2.4823 s and an unusually high pulsed fraction, 85+/-10%. Using the RXTE and CXO data, we estimated the period derivative, dot{P}=8\times 10^{-12} s/s, and confirmed the high magnetic field of the source, B=1.4\times 10^{14} G. The decay of the persistent X-ray flux, spanning 48 days, is consistent with a power law, t^{-0.5}. In the CXO/ACIS image, we find that the highly absorbed point source is surrounded by extended emission, which most likely is a dust scattering halo. Swift J1834.9-0846 is located near the center of the radio supernova remnant W41 and TeV source HESS J1834-087. An association with W41 would imply a source distance of about 4 kpc; however, any relation to the HESS source remains unclear, given the presence of several other candidate counterparts for the latter source in the field. Our search for an IR counterpart of Swift J1834.9-0846 revealed no source down to K_s=19.5 within the 0.6' CXO error circle.

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.

We present extensive spectral and photometric observations of the recycled pulsar/white-dwarf binary containing PSR J0437-4715, which we analyzed together with archival X-ray and gamma-ray data, to obtain the complete mid-infrared to gamma-ray spectrum. We first fit each part of the spectrum separately, and then the whole multi-wavelength spectrum. We find that the optical-infrared part of the spectrum is well fit by a cool white dwarf atmosphere model with pure hydrogen composition. The model atmosphere (Teff = 3950pm150K, log g=6.98pm0.15, R_WD=(1.9pm0.2)e9 cm) fits our spectral data remarkably well for the known mass and distance (M=0.25pm0.02Msun, d=156.3pm1.3pc), yielding the white dwarf age (tau=6.0pm0.5Gyr). In the UV, we find a spectral shape consistent with thermal emission from the bulk of the neutron star surface, with surface temperature between 1.25e5 and 3.5e5K. The temperature of the thermal spectrum suggests that some heating mechanism operates throughout the life of the neutron star. The temperature distribution on the neutron star surface is non-uniform. In the X-rays, we confirm the presence of a high-energy tail which is consistent with a continuation of the cut-off power-law component (Gamma=1.56pm0.01, Ecut=1.1pm0.2GeV) that is seen in gamma-rays and perhaps even extends to the near-UV.

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, which 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)

We observed the young pulsar J1357--6429 with the {\it Chandra} and {\it XMM-Newton} observatories. The pulsar spectrum fits well a combination of absorbed power-law model ($\Gamma=1.7\pm0.6$) and blackbody model ($kT=140^{+60}_{-40}$ eV, $R\sim2$ km at the distance of 2.5 kpc). Strong pulsations with pulsed fraction of $42%\pm5%$, apparently associated with the thermal component, were detected in 0.3--1.1 keV. Surprisingly, pulsed fraction at higher energies, 1.1--10 keV, appears to be smaller, $23%\pm4%$. The small emitting area of the thermal component either corresponds to a hotter fraction of the neutron star (NS) surface or indicates inapplicability of the simplistic blackbody description. The X-ray images also reveal a pulsar-wind nebula (PWN) with complex, asymmetric morphology comprised of a brighter, compact PWN surrounded by the fainter, much more extended PWN whose spectral slopes are $\Gamma=1.3\pm0.3$ and $\Gamma=1.7\pm0.2$, respectively. The extended PWN with the observed flux of $\sim7.5\times10^{-13}$ erg s$^{-1}$ cm$^{-2}$ is a factor of 10 more luminous then the compact PWN. The pulsar and its PWN are located close to the center of the extended TeV source HESS J1356--645, which strongly suggests that the VHE emission is powered by electrons injected by the pulsar long ago. The X-ray to TeV flux ratio, $\sim0.1$, is similar to those of other relic PWNe. We found no other viable candidates to power the TeV source. A region of diffuse radio emission, offset from the pulsar toward the center of the TeV source, could be synchrotron emission from the same relic PWN rather than from the supernova remnant.

We report broad-band Hubble Space Telescope imaging of the field of soft gamma-ray repeater SGR 0418+5729 with ACS/WFC and WFC3/IR. Observing in two wide filters F606W and F110W, we find no counterpart within the positional error circle derived from Chandra observations, to limiting magnitudes mF606W>28.6, mF110W>27.4 (Vega system), equivalent to reddening-corrected luminosity limits LF606W<5e28, LF110W<6e28 erg s-1 for a distance d=2 kpc, at 3sig confidence. This, in turn, imposes lower limits on the contemporaneous X-ray/optical flux ratio of 1100 and X-ray/near-infra-red flux ratio of 1000. We derive an upper limit on the temperature and/or size of any fall-back disk around the magnetar. We also compare the detection limits with observations of other magnetars.

We observed the double neutron star binary (DNSB) containing PSR J1537+1155 (also known as B1534+12) with the Chandra X-ray Observatory. This is one of the two DNSBs detected in X-rays and the only one where a hint of variability with orbital phase was found (in the previous Chandra observation). Our follow-up observation supports the earlier result: the distribution of photon arrival times with orbital phase again shows a deficit around apastron. The significance of the deficit in the combined dataset exceeds 99%. Such an orbital light-curve suggests that the X-ray emission is seen only when neutron star B passes through the equatorial pulsar wind of neutron star A. We describe statistical tests we used to determine the significance of the deficit, and conclusions that can be drawn from its existence, such as interaction of the pulsar wind with the neutron star companion. We also provide better constrained spectral model parameters obtained from the joint spectral fits to the data from both observations. A power-law successfully fits the data, with best-fit photon index Gamma=3.1pm0.4 and unabsorbed flux f=(3.2pm0.8)x10^-15 erg s-1cm-2 (0.3-8 keV range).

This paper reports the search for quiescent low-mass X-ray binaries (qLMXBs) in the globular cluster (GC) NGC 6553 using an XMM-Newton observation designed specifically for that purpose. We spectrally identify one candidate qLMXB in the core of the cluster, based on the consistency of the spectrum with a neutron star H-atmosphere model at the distance of NGC 6553. Specifically, the best-fit radius found using the three XMM European Photon Imaging Camera spectra is R_NS=6.3(+2.3)(-0.8) km (for M_NS=1.4 Msun) and the best-fit temperature is kTeff=136 (+21)(-34) eV. Both physical parameters are in accordance with typical values of previously identified qLMXBs in GC and in the field, i.e., R_NS~5-20 km and kTeff~50-150 eV. A power-law (PL) component with a photon index Gamma=2.1(+0.5)(-0.8) is also required for the spectral fit and contributes to ~33% of the total flux of the X-ray source. A detailed analysis supports the hypothesis that the PL component originates from nearby sources in the core, unresolved with XMM. The analysis of an archived Chandra observation provides marginal additional support to the stated hypothesis. Finally, a catalog of all the sources detected within the XMM field of view is presented here.

High-magnetic-field pulsars represent an important class of objects for studying the relationship between magnetars and radio pulsars. Here we report on four Chandra observations of the high-magnetic-field pulsar J1718-3718 ($B=7.4\times10^{13}$\,G) taken in 2009 as well as on a re-analysis of 2002 Chandra observations of the region. We also report an improved radio position for this pulsar based on ATCA observations. We detect X-ray pulsations at the pulsar's period in the 2009 data, with pulsed fraction 52%$\pm$13% in the 0.8--2.0 keV band. We find that the X-ray pulse is aligned with the radio pulse. The data from 2002 and 2009 show consistent spectra and fluxes; a merged overall spectrum is well fit by a blackbody of temperature $186^{+19}_{-18}$\,eV, slightly higher than predicted by standard cooling models, however, the best-fit neutron star atmosphere model is consistent with standard cooling. We find the bolometric luminosity $L^{\infty}_{\rm bb}=4^{+5}_{-2} \times10^{32}$\,erg\,s$^{-1}\sim0.3\dot{E}$, for a distance of 4.5 kpc. We compile measurements of the temperatures of all X-ray detected high-$B$ pulsars as well as those of low-$B$ radio pulsars and find evidence for the former being on average hotter than the latter.

High-resolution Chandra observations of Tycho's SNR have revealed several sets of quasi-steady, high-emissivity, nearly-parallel X-ray stripes in some localized regions of the SNR. These stripes are most likely the result of cosmic-ray (CR) generated magnetic turbulence at the SNR blast wave. However, for the amazingly regular pattern of these stripes to appear requires the simultaneous action of a number of shock-plasma phenomena and is not predicted by most models of magnetic field amplification. A consistent explanation of these stripes yields information on the complex nonlinear plasma processes connecting efficient CR acceleration and magnetic field fluctuations in strong collisionless shocks. The nonlinear diffusive shock acceleration (NL-DSA) model described here, which includes magnetic field amplification from a cosmic-ray current driven instability, does predict stripes consistent with the synchrotron emission observations of Tycho's SNR. We argue that the local ambient mean magnetic field geometry determines the orientation of the stripes and therefore it can be reconstructed with the high resolution X-ray imaging. The estimated maximum energy of the CR protons responsible for the stripes is $\sim 10^{15}$\,eV. Furthermore, the model predicts that a specific X-ray polarization pattern, with a polarized fraction ~ 50%, accompanies the stripes, which can be tested with future X-ray polarimeter missions.

LS 5039 is a high-mass binary with a period of 4 days, containing a compact object and an O star, one of the few high-mass binaries detected in gamma-rays. Our Chandra ACIS observation of LS 5039 provided a high-significance (~10sigma) detection of extended emission clearly visible for up to 1' from the point source. The spectrum of this emission can be described by an absorbed power-law model with photon index Gamma=1.9pm0.3, somewhat softer than the point source spectrum Gamma=1.44pm0.07, with the same absorption, N_H=(6.4pm0.6)e21 /cm2. The observed 0.5-8 keV flux of the extended emission is 8.8e-14 erg/s/cm2, or 5% of the point source flux; the latter is a factor of ~2 lower than the lowest flux detected so far. Fainter extended emission with comparable flux and a softer (Gamma~3) spectrum is detected at even greater radii (up to 2'). Two possible interpretations of the extended emission are a dust scattering halo and a synchrotron nebula powered by energetic particles escaping the binary. We discuss both of these scenarios and favor the nebula interpretation, although some dust contribution is possible. We have also found transient sources located within a narrow stripe south of LS 5039. We discuss the likelihood of these sources to be related to LS 5039.

PSR B1259-63 is a middle-aged radio pulsar (P=48 ms, tau=330 kyr, Edot=8.3*10^{35} erg/s) in an eccentric binary (P_orb =3.4 yr, e=0.87) with a high-mass Be companion, SS 2883. We observed the binary near apastron with the Chandra ACIS detector on 2009 May 14 for 28 ks. In addition to the previously studied pointlike source at the pulsar's position, we detected extended emission on the south-southwest side of this source. The pointlike source spectrum can be described by the absorbed power-law model with the hydrogen column density N_H = (2.5+/-0.6)*10^{21} cm^{-2}, photon index Gamma = 1.6+/-0.1, and luminosity L_{0.5-8 keV} = 1.3*10^{33} d_3^2 erg/s, where d_3 is the distance scaled to 3 kpc. This emission likely includes an unresolved part of the pulsar wind nebula (PWN) created by the colliding winds from the pulsar and the Be companion, and a contribution from the pulsar magnetosphere. The extended emission apparently consists of two components. The highly significant compact component looks like a southward extension of the pointlike source image, seen up to about 4 arcsec from the pulsar position. Its spectrum has about the same slope as the pointlike source spectrum, while its luminosity is a factor of 10 lower. We also detected an elongated feature extended ~15 arcsec southwest of the pulsar, but significance of this detection is marginal. We tentatively interpret the resolved compact PWN component as a shocked pulsar wind blown out of the binary by the wind of the Be component, while the elongated component could be a pulsar jet.

Chandra ACIS observed the field of the extended TeV source HESS J1834-087 for 47 ks. A previous XMM-Newton EPIC observation of the same field revealed a point-like source (XMMU J183435.3-084443) and an offset region of faint extended emission. In the low-resolution, binned EPIC images the two appear to be connected. However, the high-resolution Chandra ACIS images do not support the alleged connection. Instead, in these images XMMU J183435.3-084443 is resolved into a point source, CXOU~J183434.9--084443 (L=2.5 x 10e33 ergs/s, for a distance of 4 kpc; photon index 1.1), and a compact (~20'') nebula with an isotropic morphology and a softer spectrum (L=4.2 x10e33 ergs/s, photon index 2.7). The nature of the nebula is uncertain. We discuss a dust scattering halo and a pulsar-wind nebula as possible interpretations. Based on our analysis of the X-ray data, we re-evaluate the previously suggested interpretations of HESS J1834-087 and discuss a possible connection to the Fermi LAT source 1FGL J1834.3-0842c. We also obtained an upper limit of 3 x 10e-14 ergs/s cm^2 on the unabsorbed flux of the SGR J1833--0832 (in quiescence), which happened to be in the ACIS field of view.

The X-ray spectra of some magnetized isolated neutron stars (NSs) show absorption features with equivalent widths (EWs) of 50 - 200 eV, whose nature is not yet well known. To explain the prominent absorption features in the soft X-ray spectra of the highly magnetized (B ~ 10^{14} G) X-ray dim isolated NSs (XDINSs), we theoretically investigate different NS local surface models, including naked condensed iron surfaces and partially ionized hydrogen model atmospheres, with semi-infinite and thin atmospheres above the condensed surface. We also developed a code for computing light curves and integral emergent spectra of magnetized neutron stars with various temperature and magnetic field distributions over the NS surface. We compare the general properties of the computed and observed light curves and integral spectra for XDINS RBS\,1223 and conclude that the observations can be explained by a thin hydrogen atmosphere above the condensed iron surface, while the presence of a strong toroidal magnetic field component on the XDINS surface is unlikely. We suggest that the harmonically spaced absorption features in the soft X-ray spectrum of the central compact object (CCO) 1E 1207.4-5209 (hereafter 1E 1207) correspond to peaks in the energy dependence of the free-free opacity in a quantizing magnetic field, known as quantum oscillations. To explore observable properties of these quantum oscillations, we calculate models of hydrogen NS atmospheres with B ~ 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 1E 1207. We show that observable features at the electron cyclotron harmonics with EWs \approx 100 - 200 eV can arise due to these quantum oscillations.

We analyzed XMM-Newton observations of two center-filled supernova remnants, G27.8+0.6 and G28.8+1.5, to search for pulsars/pulsar wind nebulae (PWNe) and other types of neutron stars associated with these remnants. We discovered a PWN candidate within the extent of the centrally-peaked radio emission of G27.8+0.6. The X-ray morphology of the PWN candidate and its position with respect to the host supernova remnant suggest that the alleged pulsar might be moving away from the supernova remnant's center with a transverse velocity of around 100-200 km/s. The majority of the detected X-ray point sources in both fields are classified as main-sequence stars, based on their bright optical counterparts and relatively soft X-ray spectra. The remaining medium-hard and hard sources are most probably either AGNs or cataclysmic variables, although we cannot completely rule out the possibility that some of these sources are neutron stars.

PSRB1055-52 is a middle-aged (~535 kyr) radio, X-ray, and gamma-ray pulsar showing X-ray thermal emission from the neutron star (NS) surface. A candidate optical counterpart to PSRB1055-52 was proposed by Mignani and coworkers based on Hubble Space Telescope (HST) observations performed in 1996, in one spectral band only. We report on HST observations of this field carried out in 2008, in four spectral bands. The astrometric and photometric analyses of these data confirm the identification of the proposed candidate as the pulsar's optical counterpart. Similarly to other middle-aged pulsars, its optical-UV spectrum can be described by the sum of a power-law (PLO) component, presumably emitted from the pulsar magnetosphere, and a Rayleigh-Jeans (RJ) component emitted from the NS surface. The spectral index of the PLO component, alpha_O=1.05+/-0.34, is larger than for other pulsars with optical counterparts. The RJ component, with the brightness temperature TO=(0.66+/-0.10) d_350**2 R_O,13**-2 MK (where d_350 and R_O,13 are the distance to the pulsar in units of 350 pc and the radius of the emitting area in units of 13 km), shows a factor of 4 excess with respect to the extrapolation of the X-ray thermal component into the UV-optical. This hints that the RJ component is emitted from a larger, colder area, and suggests that the distance to the pulsar is smaller than previously thought. From the absolute astrometry of the HST images we measured the pulsar coordinates with a position accuracy of 0.15". From the comparison with previous observations we measured the pulsar proper motion, mu = 42+/-5 mas/yr, which corresponds to a transverse velocity V_t = (70+/-8) d_350 km/s.

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 pulsar's 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} ~ 3\times 10^{29} d_{250}^2 erg/s, is lower than the pulsar's 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 pulsar's 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 caused 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.

Pulsars are known to be efficient accelerators that produce copious amounts of relativistic particles and inject them into the Galactic medium. The radiation emitted by such a pulsar wind can be seen from radio through gamma-rays as a pulsar-wind nebula (PWN). Here we overview and summarize recent progress in X-ray and TeV observations of PWNe.

We present results of a recent Chandra X-ray Observatory observation of the central compact object (CCO) in the supernova remnant Cassiopeia A. This observation was obtained in an instrumental configuration that combines a high spatial resolution with a minimum spectral distortion, and it allowed us to search for pulsations with periods longer than 0.68 s. We found no evidence of extended emission associated with the CCO, nor statistically significant pulsations (the 3-sigma upper limit on pulsed fraction is about 16%). The fits of the CCO spectrum with the power-law model yield a large photon index, Gamma\approx 5, and a hydrogen column density larger than that obtained from the SNR spectra. The fits with the blackbody model are statistically unacceptable. Better fits are provided by hydrogen or helium neutron star atmosphere models, with the best-fit effective temperature kT_{eff}^\infty \approx 0.2 keV, but they require a small star's radius, R = 4 - 5.5 km, and a low mass, M < 0.8 M_sol. A neutron star cannot have so small radius and mass, but the observed emission might emerge from an atmosphere of a strange quark star. More likely, the CCO could be a neutron star with a nonuniform surface temperature and a low surface magnetic field (the so-called anti-magnetar), similar to three other CCOs for which upper limits on period derivative have been established. The bolometric luminosity, L_{bol}^\infty \sim 6\times 10^{33} erg s^{-1}, estimated from the fits with the hydrogen atmosphere models, is consistent with the standard neutron star cooling for the CCO age of 330 yr. The origin of the surface temperature nonuniformity remains to be understood; it might be caused by anisotropic heat conduction in the neutron star crust with very strong toroidal magnetic fields.

This paper presents the analysis of candidate quiescent low mass xray binarie (qLMXBs) observed during a short Chandra/ACIS observation of the globular cluster (GC) NGC 6304. Two out of the three candidate qLMXBs of this cluster, XMMU 171433-292747 and XMMU 171421-292917, lie within the field of view. This permits comparison with the discovery observation of these sources. The one in the GC core -- XMMU 171433-292747 -- is spatially resolved into two separate X-ray sources, one of which is consistent with a pure H-atmosphere qLMXB, and the other is an X-ray power-law spectrum source. These two spectral components separately account for those observed from XMMU 171433-292747 in its discovery observation. We find that the observed flux and spectral parameters of the H-atmosphere spectral components are consistent with the previous observation, as expected from a qLMXB powered by deep crustal heating. XMMU 171421-292917 also has neutron star atmosphere spectral parameters consistent with those in the XMM-Newton observation and the observed flux has decreased by a factor 0.54^{+0.30}_{-0.24}.

Pulsar winds shocked in the ambient medium produce spectacular nebulae observable from the radio through gamma-rays. The shape and the spectrum of a pulsar wind nebula (PWN) depend on the angular distribution, magnetization and energy spectrum of the wind streaming from the pulsar magnetosphere, as well as on the pulsar velocity and the properties of the ambient medium. The advent of Chandra, with its unprecedented angular resolution and high sensitivity, has allowed us not only to detect many new PWNe, but also study their spatial and spectral structure and dynamics, which has significantly advanced our understanding of these objects. Here we overview recent observational results on PWNe, with emphasis on Chandra observations.

PSR J1906+0746 is a 112-kyr-old radio pulsar in a tight relativistic binary with a compact high-mass companion, at the distance of about 5 kpc. We observed this unique system with the Chandra ACIS detector for 31.6 ks. Surprisingly, not a single photon was detected within the 3" radius from the J1906+0746 radio position. For a plausible range of hydrogen column densities, n_H=(0.5-1)\times10^{22} cm^{-2}, the nondetection corresponds to the 90% upper limit of (3-5)\times10^{30} erg s^{-1} on the unabsorbed 0.5-8 keV luminosity for the power-law model with Gamma=1.0-2.0, and ~10^{32} erg s^{-1} on the bolometric luminosity of the thermal emission from the NS surface. The inferred limits are the lowest known for pulsars with spin-down properties similar to those of PSR J1906+0746. We have also tentatively detected a puzzling extended structure which looks like a tilted ring with a radius of 1.6' centered on the pulsar. The measured 0.5-8 keV flux of the feature, 3.1\times10^{-14} erg cm^{-2} s^{-1}, implies an unabsorbed luminosity of 1.2\times10^{32} erg s^{-1} (4.5\times10^{-4} of the pulsar's spin-down power). Although all conventional interpretations of the ring appear to be problematic, the pulsar-wind nebula with an unusually underluminous pulsar remains the most viable interpretation.

The cores of neutron stars harbor the highest matter densities known to occur in nature, up to several times the densities in atomic nuclei. Similarly, magnetic field strengths can exceed the strongest fields generated in terrestrial laboratories by ten orders of magnitude. Hyperon-dominated matter, deconfined quark matter, superfluidity, even superconductivity are predicted in neutron stars. Similarly, quantum electrodynamics predicts that in strong magnetic fields the vacuum becomes birefringent. The properties of matter under such conditions is governed by Quantum Chromodynamics (QCD) and Quantum Electrodynamics (QED), and the close study of the properties of neutron stars offers the unique opportunity to test and explore the richness of QCD and QED in a regime that is utterly beyond the reach of terrestrial experiments. Experimentally, this is almost virgin territory.

We report the search for low-mass X-ray binaries in quiescence (qLMXBs) in the globular cluster NGC 6304 using XMM observations. We present the spectral analysis leading to the identification of three candidate qLMXBs within the field of this globular cluster (GC), each consistent with the X-ray spectral properties of previously identified qLMXBs in the field and in other globular clusters -- specifically, with a hydrogen atmosphere neutron star with radius between 5--20\km. One (source 4, with R=11.7^{+8.3}_{-0.4} (D/5.97 kpc) km and kT_eff=117^{+59}_{-44} eV) is located within one core radius (r_c) of the centre of NGC 6304. This candidate also presents a spectral power-law component contributing 49 per cent of the 0.5-10 keV flux. A second one (source 9 with R=15.3^{+11.2}_{-6.5} (D/5.97 kpc) km and kT_eff=100^{+24}_{-19} eV) is found well outside the optical core (at 32 r_c) but still within the tidal radius. From spatial coincidence, we identify a bright 2MASS infrared counterpart which, at the distance of NGC 6304, seems to be a post-asymptotic giant branch star. The third qLMXB (source 5 with R=23^{+38}_{-14} (D/5.97 kpc) km and kT_eff=70^{+28}_{-20} eV) is a low signal-to-noise candidate for which we also identify from spatial coincidence a bright 2MASS infrared counterpart, with 99.916 per cent confidence. Three qLMXBs from this GC is marginally consistent with that expected from the encounter rate of NGC 6304. We also report a low signal-to-noise source with an unusually hard photon index (\alpha=-2.0^{+1.2}_{-2.2}). Finally, we present an updated catalogue of the X-ray sources lying in the field of NGC 6304, and compare this with the previous catalogue compiled from ROSAT observations.

We present results from our X-ray data analysis of the SNR G330.2+1.0 and its CCO, CXOU J160103.1--513353 (J1601). Using our XMM-Newton and Chandra observations, we find that the X-ray spectrum of J1601 can be described by neutron star atmosphere models (T ~ 2.5--3.7 MK). Assuming the distance of d ~ 5 kpc for J1601 as estimated for SNR G330.2+1.0, a small emission region of R ~ 1--2 km is implied. X-ray pulsations previously suggested by Chandra are not confirmed by the XMM-Newton data, and are likely not real. However, our timing analysis of the XMM-Newton data is limited by poor photon statistics, and thus pulsations with a relatively low amplitude (i.e., an intrinsic pulsed-fraction < 40%) cannot be ruled out. Our results indicate that J1601 is a CCO similar to that in the Cassiopeia A SNR.X-ray emission from SNR G330.2+1.0 is dominated by power law continuum (Gamma ~ 2.1--2.5) which primarily originates from thin filaments along the boundary shell. This X-ray spectrum implies synchrotron radiation from shock-accelerated electrons with an exponential roll-off frequency ~ 2--3 x 10^17 Hz. For the measured widths of the X-ray filaments (D ~ 0.3 pc) and the estimated shock velocity (v_s ~ a few x 10^3 km s^-1), a downstream magnetic field B ~ 10--50 $\mu$G is derived. The estimated maximum electron energy E_max ~ 27--38 TeV suggests that G330.2+1.0 is a candidate TeV gamma-ray source. We detect faint thermal X-ray emission in G330.2+1.0. We estimate a low preshock density n_0 ~ 0.1 cm^-3, which suggests a dominant contribution from an inverse Compton mechanism (than the proton-proton collision) to the prospective gamma-ray emission. Follow-up deep radio, X-ray, and gamma-ray observations will be essential to reveal the details of the shock parameters and the nature of particle accelerations in this SNR.

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 2\times 10^{28} d_{130}^2 ergs s^{-1}, or by a blackbody model with the temperature kT\approx 0.28 keV and bolometric luminosity L_{bol} \sim 1.3\times 10^{28} d_{130}^2 ergs s^{-1}, for a plausible hydrogen column density NH = 7.3\times 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.

We report on two Chandra observations of the 3-Myr pulsar B1929+10, which reveal a faint compact (~9"x5") nebula elongated in the direction perpendicular to the pulsar's proper motion, two patchy wings, and a possible short (~3") jet emerging from the pulsar. In addition, we detect a tail extending up to at least 4' in the direction opposite to the pulsar's proper motion, aligned with the 15'-long tail detected in ROSAT and XMM-Newton observations. The overall morphology of the nebula suggests that the shocked pulsar wind is confined by the ram pressure due to the pulsar's supersonic speed. The shape of the compact nebula in the immediate vicinity of the pulsar seems to be consistent with the current MHD models. However, since these models do not account yet for the change of the flow velocity at larger distances from the pulsar, they are not able to constrain the extent of the long pulsar tail. The luminosity of the whole nebula as seen by Chandra is ~10^30 ergs/s in the 0.3-8 keV band, for the distance of 361 pc. Using the Chandra and XMM-Newton data, we found that the pulsar spectrum is comprised of non-thermal (magnetospheric) and thermal components. The non-thermal component can be described by a power-law model with photon index ~1.7 and luminosity 1.7x10^30 ergs/s in the 0.3-10 keV band. The blackbody fit for the thermal component, which presumably emerges from hot polar caps, gives the temperature kT~0.3 keV and projected emitting area 3x10^3 m^2, corresponding to the bolometric luminosity ~(1-2)x10^30 ergs/s.

The multi-wavelength study of old (>100 Myr) radio pulsars holds the key to understanding the long-term evolution of neutron stars, including the advanced stages of neutron star cooling and the evolution of the magnetosphere. Optical/UV observations are particularly useful for such studies because they allow one to explore both thermal and non-thermal emission processes. In particular, studying the optical/UV emission constrains temperature of the bulk of the neutron star surface, too cold to be measured in X-ray observations.Aim of this work is to identify the optical counterpart of the very old (166 Myr) radio pulsar J0108-1431. We have re-analyzed our original VLT observations (Mignani et al. 2003), where a very faint object was tentatively detected close to the radio position, near the edge of a field galaxy. We found that the backward extrapolation of the PSR J0108-1431 proper motion recently measured by CHANDRA(Pavlov et al. 2008) nicely fits the position of this object. Based on that, we propose it as a viable candidate for the optical counterpart to PSR J0108-1431. The object fluxes (U =26.4+/-0.3; B =27.9; V >27.8) are consistent with a thermal spectrum with a brightness temperature of 9X10^4 K (for R = 13 km at a distance of 130 pc), emitted from the bulk of the neutron star surface. New optical observations are required to confirm the optical identification of PSR J0108-1431 and measure its spectrum.

We observed the young energetic pulsar J1617-5055 with the Chandra ACIS detector for 60 ks. In addition to the pulsar, the X-ray images show a faint pulsar-wind nebula (PWN) seen up to ~1' from the pulsar. Deconvolution and reconstruction of the image reveal a brighter compact PWN component of ~1'' size, possibly with a jet-torus morphology. The PWN spectrum fits an absorbed power-law (PL) model with the photon index Gamma=1.5. The total PWN luminosity in the 0.5-8 keV band, L_{pwn}=3.2\times10^{33} ergs s^{-1} for d=6.5 kpc, is a fraction of 2\times 10^{-4} of the pulsar's spin-down power and a fraction of 0.2 of the pulsar's X-ray luminosity, which is a factor of 20 lower than one would expect from an average empirical relation found for a sample of PWNe observed with Chandra. The pulsar's spectrum can be described by an absorbed PL with n_{H}=3.5\times10^{22} cm^{-2} and Gamma=1.1, harder than any other pulsar spectrum reliably measured in the soft X-ray band. This non-thermal emission is 50% pulsed showing one peak per period. We have also investigated a possible connection between J1617 and the extended TeV source HESS J1616-508 whose center is located about 10' west of the pulsar. We find no preferential extension of the X-ray PWN toward the TeV source. Therefore, the Chandra data do not provide conclusive evidence for PSR J1617-5055 and HESS J1616-508 association. We have also analyzed archival X-ray, radio, and IR data on the HESS J1616-508 region and found traces of diffuse emission (resembling a shell in the radio) coinciding with the central part of HESS J1616-508. We speculate that the TeV source may be multiple, with most of the emission coming from an unknown SNR or a star-forming region, while some fraction of the TeV emission still may be attributed to the J1617 PWN.

Chandra and XMM-Newton resolved extremely long tails behind two middle-aged pulsars, J1509-5850 and J1740+1000. The tail of PSR J1509-5850 is discernible up to 5.6' from the pulsar (6.5 pc at a distance of 4 kpc), with a flux of 2*10^{-13} erg s^{-1} cm^{-2} in 0.5-8 keV. The tail spectrum fits an absorbed power-law (PL) model with the photon index of 2.3\pm0.2, corresponding to the 0.5-8 keV luminosity of 1*10^{33} ergs s^{-1}, for n_H= 2.1*10^{22} cm^{-2}. The tail of PSR J1740+1000 is firmly detected up to 5' (2 pc at a 1.4 kpc distance), with a flux of 6*10^{-14} ergs cm^{-2} s^{-1} in 0.4-10 keV. The PL fit yields photon index of 1.4-1.5 and n_H=1*10^{21} cm^{-2}. The large extent of the tails suggests that the bulk flow in the tails starts as mildly relativistic downstream of the termination shock, and then gradually decelerates. Within the observed extent of the J1509-5850 tail, the average flow speed exceeds 5,000 km s^{-1}, and the equipartition magnetic field is a few times 10^{-5} G. For the J1740+1000 tail, the equipartition field is a factor of a few lower. The harder spectrum of the J1740+1000 tail implies either less efficient cooling or a harder spectrum of injected electrons. For the high-latitude PSR J1740+1000, the orientation of the tail on the sky shows that the pulsar is moving toward the Galactic plane, which means that it was born from a halo-star progenitor. The comparison between the J1509 and J1740 tails and the X-ray tails of other pulsars shows that the X-ray radiation efficiency correlates poorly with the pulsar spin-down luminosity or age. The X-ray efficiencies of the ram-pressure confined pulsar wind nebulae (PWNe) are systematically higher than those of PWNe around slowly moving pulsars with similar spin-down parameters.

The nature of the extended hard X-ray source XMMU J061804.3+222732 and its surroundings is investigated using XMM-Newton, Chandra, and Spitzer observations. This source is located in an interaction region of the IC 443 supernova remnant with a neighboring molecular cloud. The X-ray emission consists of a number of bright clumps embedded in an extended structured non-thermal X-ray nebula larger than 30" in size. Some clumps show evidence for line emission at ~1.9 keV and ~3.7 keV at the 99% confidence level. Large-scale diffuse radio emission of IC 443 passes over the source region, with an enhancement near the source. An IR source of about 14" x 7" size is prominent in the 24 um, 70 um, and 2.2 um bands, adjacent to a putative Si K-shell X-ray line emission region. The observed IR/X-ray morphology and spectra are consistent with those expected for J/C-type shocks of different velocities driven by fragmented supernova ejecta colliding with the dense medium of a molecular cloud. The IR emission of the source detected by Spitzer can be attributed to both continuum emission from an HII region created by the ejecta fragment and line emission excited by shocks. This source region in IC 443 may be an example of a rather numerous population of hard X-ray/IR sources created by supernova explosions in the dense environment of star-forming regions. Alternative Galactic and extragalactic interpretations of the observed source are also discussed.

We report on an observation of the Vela-like pulsar B1823-13 and its synchrotron nebula with Chandra.The pulsar's spectrum fits a power-law model with a photon index Gamma_PSR=2.4 for the plausible hydrogen column density n_H=10^{22} cm^{-2}, corresponding to the luminosity L_PSR=8*10^{31} ergs s^{-1} in the 0.5-8 keV band, at a distance of 4 kpc. The pulsar radiation likely includes magnetospheric and thermal components, but they cannot be reliably separated because of the small number of counts detected and strong interstellar absorption. The pulsar is surrounded by a compact, 25''x 10'', pulsar wind nebula (PWN) elongated in the east-west direction, which includes a brighter inner component, 7''x 3'', elongated in the northeast-southwest direction. The slope of the compact PWN spectrum is Gamma_comp=1.3, and the 0.5-8 keV luminosity is L_comp~3*10^{32} ergs s^{-1}. The compact PWN is surrounded by asymmetric diffuse emission (extended PWN) seen up to at least 2.4' south of the pulsar, with a softer spectrum (Gamma_ext=1.9), and the 0.5-8 keV luminosity L_ext~10^{33}-10^{34} ergs s^{-1}. We also measured the pulsar's proper motion using archival VLA data: \mu_\alpha=23.0+/-2.5 mas yr^{-1}, \mu_\delta=-3.9+/-3.3 mas yr^{-1}, which corresponds to the transverse velocity v_perp=440 km s^{-1}. The direction of the proper motion is approximately parallel to the elongation of the compact PWN, but it is nearly perpendicular to that of the extended PWN and to the direction towards the center of the bright VHE gamma-ray source HESS J1825-137, which is likely powered by PSR B1823-13.

We detected X-ray emission from the 50-kyr-old pulsar J1809-1917 and resolved its pulsar wind nebula (PWN) with Chandra. The pulsar spectrum fits PL+BB model with the photon index of 1.2 and the BB temperature of 2 MK for n_{H}=0.7\times 10^{22} cm^{-2}. The luminosities are(4\pm 1)\times 10^{31} ergs s^{-1} for the PL component (in the 0.5-8 keV band) and ~1\times 10^{32} ergs s^{-1} for the BB component (bolometric) at a plausible distance of 3.5 kpc. The bright inner PWN component of a 3''\times12'' size is elongated in the north-south direction, with the pulsar close to its south end. This component is immersed in a larger (20''\times40''), similarly elongated outer PWN component of lower surface brightness. The elongated shape of the compact PWN can be explained by the ram pressure confinement of the pulsar wind due to the supersonic motion of the pulsar. The PWN spectrum fits a PL model with photon index of 1.4\pm0.1 and 0.5-8 keV luminosity of 4\times10^{32} ergs s^{-1}. The compact PWN appears to be inside a large-scale (~4'\times4') emission more extended to the south of the pulsar, i.e. in the direction of the alleged pulsar motion. To explain the extended X-ray emission ahead of the moving pulsar, one has to invoke strong intrinsic anisotropy of the pulsar wind or assume that this emission comes from a relic PWN swept by the asymmetrical reverse SNR shock. The pulsar and its PWN are located within the extent of the unidentified TeV source HESS J1809-193 whose brightest part is offset by ~8' to the south of the pulsar, i.e. in the same direction as the large-scale X-ray emission. Although the association between J1809-1917 and HESS J1809-193 is plausible, an alternative source of relativistic electrons powering HESS J1809-193 might be the serendipitously discovered X-ray source CXOU J180940.7-192544.

We report the detection of the millisecond pulsar B1257+12 with the Chandra X-ray Observatory. In a 20 ks exposure we detected 25 photons from the pulsar, with energies between 0.4 and 2.0 keV, corresponding to the flux F_X=(4.4+/- 0.9)*10^{-15} ergs s^{-1} cm^{-2} in this energy range. The X-ray spectrum can be described by a power-law model with photon index Gamma = 2.8 and luminosity L_X \approx 2.5*10^{29} ergs s^{-1} in the 0.3--8 keV band, for a plausible distance of 500 pc and hydrogen column density N_H=3*10^{20} cm^{-2}. Alternatively, the spectrum can be fitted by a blackbody model with kT ~ 0.22 keV and projected emitting area ~2000 m^2. If the thermal X-rays are emitted from two symmetric polar caps, the bolometric luminosity of the two caps is 2 L_bol ~ 3*10^{29} ergs s^{-1}. We compared our results with the data on other 30 millisecond pulsars observed in X-rays and found that the apparent X-ray efficiency of PSR B1257+12, L_X/Edot ~ 3*10^{-5} for d=500 pc, is lower than those of most of millisecond pulsars. This might be explained by an unfavorable orientation of the X-ray pulsar beam if the radiation is magnetospheric, or by strong asymmetry of polar caps if the radiation is thermal (e.g., one of the polar caps is much brighter than the other and remains invisible for most part of the pulsar period). Alternatively, it could be attributed to absorption of X-rays in circumpulsar matter, such as a flaring debris disk left over after formation of the planetary system around the pulsar.

RX J0720.4-3125 is an isolated neutron star that, uniquely in its class, has shown changes in its thermal X-ray spectrum. We use new spectra taken with Chandra's Low Energy Transmission Grating Spectrometer, as well as archival observations, to try to understand the timescale and nature of these changes. We construct lightcurves, which show both small, slow variations on a timescale of years, and a larger event that occurred more quickly, within half a year. From timing, we find evidence for a `glitch' coincident with this larger event, with a fractional increase in spin frequency of 5x10^{-8}. We compare the `before' and `after' spectra with those from RX J1308.6+2127, an isolated neutron star with similar temperature and magnetic field strength, but with a much stronger absorption feature in its spectrum. We find that the `after' spectrum can be represented remarkably well by the superposition of the `before' spectrum, scaled by two thirds, and the spectrum of RX J1308.6+2127, thus suggesting that the event affected approximately one third of the surface. We speculate the event reflects a change in surface composition caused by, e.g., an accretion episode.

The field of the extended TeV source HESS J1804-216 was serendipitously observed with the Chandra ACIS detector. The data reveal several X-ray sources within the bright part of HESS J1804-216. The brightest of these objects, CXOU J180432.4-214009, which has been also detected with Swift and Suzaku, is consistent with being a point-like source, with the 0.3-7 keV flux of (1.7\pm0.2)\times10^{-13} ergs s^{-1} cm^{-2}. Its hard and strongly absorbed spectrum can be fitted by the absorbed power-law model with the best-fit photon index of 0.45 and hydrogen column density of 4\times10^{22} cm^{-2}, both with large uncertainties due to the strong correlation between these parameters. A search for pulsations resulted in a 106 s period candidate, which however has a low significance of 97.9%. We found no infrared-optical counterparts for this source. The second brightest source, CXOU J180441.9-214224, which has been detected with Suzaku, is either extended or multiple, with the flux of about 10^{-13} ergs cm^{-2} s^{-1}. We found a nearby M dwarf within the X-ray source extension, which could contribute a fraction of the observed X-ray flux. The remaining sources are very faint (<3\times 10^{-14} ergs cm^{-2} s^{-1}), and at least some of them are likely associated with nearby stars. Although one or both of the two brighter X-ray sources could be faint accreting binaries or remote pulsars with pulsar wind nebulae (hence possible TeV sources), their relation to HESS J1804-216 remains elusive. The possibility that HESS J1804-216 is powered by the relativistic wind from the young pulsar B1800-21, located at a distance of about 10 pc from the TeV source, still remains a more plausible option.

We detected X-ray emission from PSR B1800-21 and its synchrotron nebula with the Chandra X-ray Observatory. The pulsar's observed flux is (1.4+/-0.2) 10^{-14} ergs cm^{-2} s^{-1} in the 1-6 keV band. The spectrum can be described by a two-component PL+BB model, suggesting a mixture of thermal and magnetospheric emission. For a plausible hydrogen column density n_{H}=1.4 10^{22} cm^{-2}, the PL component has a slope Gamma=1.4+/-0.6 and a luminosity L_{psr}^{nonth}=4 10^{31}(d/4 kpc)^2 ergs s^{-1}. The properties of the thermal component (kT=0.1-0.3 keV, L^{bol}=10^{31}-10^{33} ergs s^{-1}) are very poorly constrained because of the strong interstellar absorption. The compact, 7''\times4'', inner pulsar-wind nebula (PWN), elongated perpendicular to the pulsar's proper motion, is immersed in a fainter asymmetric emission. The observed flux of the PWN is (5.5+/-0.6) 10^{-14} ergs cm^{-2} s^{-1} in the 1-8 keV band. The PWN spectrum fits by a PL model with Gamma=1.6+/-0.3, L=1.6 10^{32} (d/4 kpc})^2 ergs s^{-1}. The shape of the inner PWN suggests that the pulsar moves subsonically and X-ray emission emerges from a torus associated with the termination shock in the equatorial pulsar wind. The inferred PWN-pulsar properties (e.g., the PWN X-ray efficiency, L_{pwn}/\dot{E}~10^{-4}; the luminosity ratio, L_{pwn}/L_{psr}^{nonth}=4; the pulsar wind pressure at the termination shock, p_s=10^{-9} ergs cm^{-3}) are very similar to those of other subsonically moving Vela-like objects detected with Chandra (L_{pwn}/\dot{E}=10^{-4.5}-10^{-3.5}, L_{pwn}/L_{psr}^{nonth}~5, p_s=10^{-10}-10^{-8} ergs cm^{-1}).

We review the state of the art for measuring the X-ray polarization of neutron stars. We discuss how valuable precision measurements of the degree and position angle of polarization as a function of energy and, where relevant, of pulse phase, would provide deeper insight into the details of the emission mechanisms. We then review the current state of instrumentation and its potential for obtaining relevant data. Finally, we conclude our discussion with some opinions as to future directions.

We report on the discovery of a pointlike source (CXOU J160103.1$-$513353) at the center of a Galactic supernova remnant (SNR) G330.2+1.0 with {\it Chandra X-Ray Observatory}. The X-ray spectrum fits a black-body (BB) model with $kT$ $\sim$ 0.49 keV, implying a small emission region of $R$ $\sim$ 0.4 km at the distance of 5 kpc. The estimated X-ray luminosity is $L_X$ $\sim$ 1 $\times$ 10$^{33}$ ergs s$^{-1}$ in the 1 $-$ 10 keV band. A power law model may also fit the observed spectrum, but the fit results in a very large photon index, $\Gamma$ $\sim$ 5. We find no counterparts at other wavelengths. The X-ray emission was steady over the $\sim$13 hr observation period, showing no variability. While we find marginal evidence for X-ray pulsations ($P$ $\approx$ 7.5 s), the presence of a pulsar at the position of this object is not conclusive with the current data, requiring an independent confirmation. These results are generally consistent with an interpretation of this object as a Central Compact Object associated with SNR G330.2+1.0.

We found a very likely counterpart to the recently discovered hard X-ray source IGR J2018+4043 in the multi-wavelength observations of the source field. The source, originally discovered in the 20-40 keV band, is now confidently detected also in the 40-80 keV band, with a flux of (1.4 +/- 0.4) x 10(-11) erg cm(-2) s(-1). A 5 ks Swift observation of the IGR J2018+4043 field revealed a hard point-like source with the observed 0.5-10 keV flux of 3.4(+0.7)(-0.8) x 10(-12) erg cm(-2) s(-1) (90% confidence level) at alpha = 20h18m38.55s, delta = +40d41m00.4s (with a 4.2" uncertainty). The combined Swift-INTEGRAL spectrum can be described by an absorbed power-law model with photon index gamma = 1.3 +/- 0.2 and N_H = 6.1(+3.2)(-2.2) x 10(22) cm(-2). In archival optical and infrared data we found a slightly extended and highly absorbed object at the Swift source position. There is also an extended VLA 1.4 GHz source peaked at a beam-width distance from the optical and X-ray positions. The observed morphology and multiwavelength spectra of IGR J2018+4043 are consistent with those expected for an obscured accreting object, i.e. an AGN or a Galactic X-ray binary. The identification suggests possible connection of IGR J2018+4043 to the bright gamma-ray source GEV J2020+4023 (3EG J2020+4017) detected by COS B and CGRO EGRET in the gamma-Cygni SNR field.

We present the UV spectroscopy and timing of three nearby pulsars (Vela, B0656+14 and Geminga) recently observed with the Space Telescope Imaging Spectrograph. We also review the optical and X-ray properties of these pulsars and establish their connection with the UV properties. We show that the multiwavelengths properties of neutron stars (NSs) vary significantly within the sample of middle-aged pulsars. Even larger differences are found between the thermal components of Ge-minga and B0656+14 as compared to those of radio-quiet isolated NSs. These differences could be attributed to different properties of the NS surface layers.

Unique among neutron stars, 1E 1207.4-5209 is an X-ray pulsar with a spin period of 424 ms that contains at least two strong absorption features in its energy spectrum. This neutron star has been identified as a member of the radio-quiet compact central objects in supernova remnants. It has been found that 1E 1207.4-5209 is not spinning down monotonically suggesting that this neutron star undergoes strong, frequent glitches, contains a fall-back disk, or possess a binary companion. Here, we report on a sequence of seven XMM-Newton observations of 1E 1207.4-5209 performed during a 40 day window in June/July 2005. Due to unanticipated variance in the phase measurements beyond the statistical uncertainties, we could not identify a unique phase-coherent timing solution. The three most probable timing solutions give frequency time derivatives of +0.9, -2.6, and +1.6 X 10^(-12) Hz/s (listed in descending order of significance). We conclude that the local frequency derivative during our XMM-Newton observing campaign differs from the long-term spin-down rate by more than an order of magnitude, effectively ruling out glitch models for 1E 1207.4-5209. If the long-term spin frequency variations are caused by timing noise, the strength of the timing noise in 1E 1207.4-5209 is much stronger than in other pulsars with similar period derivatives. Therefore, it is highly unlikely that the spin variations are caused by the same physical process that causes timing noise in other isolated pulsars. The most plausible scenario for the observed spin irregularities is the presence of a binary companion to 1E 1207.4-5209. We identified a family of orbital solutions that are consistent with our phase-connected timing solution, archival frequency measurements, and constraints on the companions mass imposed by deep IR and optical observations.

We present a brief review of Chandra X-ray Observatory observations of neutron stars. The outstanding spatial and spectral resolution of this great observatory have allowed for observations of unprecedented clarity and accuracy. Many of these observations have provided new insights into neutron star physics. We present an admittedly biased and overly brief overview of these observations, highlighting some new discoveries made possible by the Observatory's unique capabilities. We also include our analysis of recent multiwavelength observations of the putative pulsar and its pulsar-wind nebula in the IC 443 SNR.

We report the detection of the double neutron star binary (DNSB) B1534+12 (= J1537+1155) with the Chandra X-ray Observatory. This DNSB (orbital period 10.1 hr) consists of the millisecond (recycled) pulsar J1537+1155A (P_A=37.9 ms) and a neutron star not detected in the radio. After the remarkable double pulsar binary J0737-3039, it is the only other DNSB detected in X-rays. We measured the flux of (2.2\pm 0.6)\times10^{-15} ergs s^{-1} cm^{-2} in the 0.3-6 keV band. The small number of collected counts allows only crude estimates of spectral parameters. The power-law fit yields the photon index of 3.2\pm 0.5 and the unabsorbed 0.2-10 keV luminosity L_X=6\times10^{29} ergs s^{-1} = 3\times 10^{-4}Edot_A, where Edot_A is the spin-down power of J1537+1155A. Alternatively, the spectrum can be fitted by a blackbody model with T = 2.2 MK and the projected emitting area of ~ 5\times 10^3 m^2. The distribution of photon arrival times over binary orbital phase shows a deficit of X-ray emission around apastron, which suggests that the emission is caused by interaction of the relativistic wind from J1537+1155A with its neutron star companion. We also reanalyzed the Chandra and XMM-Newton observations of J0737-3039 and found that its X-ray spectrum is similar to the spectrum of B1534+12, and its X-ray luminosity is about the same fraction of Edot_A, which suggests similar X-ray emission mechanisms. However, the X-ray emission from J0737-3039 does not show orbital phase dependence. This difference can be explained by the smaller eccentricity of J0737-3039 or a smaller misalignment between the equatorial plane of the millisecond pulsar and the orbital plane of the binary.

The superb spatial resolution of Chandra has allowed us to detect a 20''-long tail behind the Geminga pulsar, with a hard spectrum (photon index 1.0+/-0.2) and a luminosity (1.3+/-0.2) 10^{29} ergs/s in the 0.5 - 8 keV band, for an assumed distance of 200 pc. The tail could be either a pulsar jet, confined by a toroidal magnetic field of about 100 microGauss, or it can be associated with the shocked relativistic wind behind the supersonically moving pulsar confined by the ram pressure of the oncoming interstellar medium. We also detected an arc-like structure 5'' - 7'' ahead of the pulsar, extended perpendicular to the tail, with a factor of 3 lower luminosity. We see a 3-sigma enhancement in the Chandra image apparently connecting the arc with the southern outer tail that has been possibly detected with XMM-Newton. The observed structures imply that the Geminga's pulsar wind is intrinsically anisotropic.

We present an X-ray tour of diffuse emission in the 30 Doradus star-forming complex in the Large Magellanic Cloud using high-spatial-resolution X-ray images and spatially-resolved spectra obtained with the Advanced CCD Imaging Spectrometer aboard the Chandra X-ray Observatory. The dominant X-ray feature of the 30 Doradus nebula is the intricate network of diffuse emission generated by interacting stellar winds and supernovae working together to create vast superbubbles filled with hot plasma. We construct maps of the region showing variations in plasma temperature (T = 3--9 million degrees), absorption (N_H = 1--6 x 10^{21} cm^{-2}), and absorption-corrected X-ray surface brightness (S_X = 3--126 x 10^{31} ergs s^{-1} pc^{-2}). Enhanced images reveal the pulsar wind nebula in the composite supernova remnant N157B and the Chandra data show spectral evolution from non-thermal synchrotron emission in the N157B core to a thermal plasma in its outer regions. In a companion paper we show that R136, the central massive star cluster, is resolved at the arcsecond level into almost 100 X-ray sources. Through X-ray studies of 30 Doradus the complete life cycle of such a massive stellar cluster can be revealed.

We detected a nearby (d=360 pc), old (5 Myr) pulsar B1133+16 with Chandra. The observed pulsar's flux is $(0.8\pm 0.2)\times 10^{-14}$ ergs cm$^{-2}$ s$^{-1}$ in the 0.5-8 keV band. Because of the small number of counts detected, the spectrum can be described by various models. A power-law fit of the spectrum gives a photon index $\Gamma = 2.5$ and an isotropic luminosity of $1.4\times 10^{29}$ ergs s$^{-1}$ in the 0.5-8 keV band, which is about $1.6\times 10^{-3}$ of the spin-down power $\dot{E}$. The spectrum can also be fitted by a blackbody model with a temperature of $\approx 2.8$ MK and a projected emitting area of $\sim 500$ m$^2$, possibly a hot polar cap. The X-ray properties of PSR B1133+16 are similar to those of other old pulsars observed in X-rays, particularly the drifting pulsar B0943+10.

Near-infrared images of the center of the Cas A supernova remnant taken with the NICMOS 2 camera on HST using the F110W and F160W filters (~J & ~H bands) have magnitude limits of 26.2 and 24.6 respectively, but reveal no sources within a 1.2" radius of the Chandra detected central X-ray point source (XPS). The NICMOS data, taken together with broadband optical magnitude limits (R ~ 28 mag) obtained from a deep STIS CCD exposure taken with a clear filter (50CCD), indicate that the XPS luminosities are very low in the optical/NIR bands (e.g., L_H < 3 x E29 erg/s) with no optical, J, or H band counterpart to the XPS easily detectable by HST. We discuss the nature of the Cas A central compact object based upon these new near-infrared and optical flux limits.

It is generally believed that magnetic fields of some neutron stars, the so-called magnetars, are enormously strong, up to 10^{14} - 10^{15} G. Recent investigations have shown that the atmospheres of magnetars are possibly composed of helium. We calculate the structure and bound-bound radiative transitions of the He^+ ion in superstrong fields, including the effects caused by the coupling of the ion's internal degrees of freedom to its center-of-mass motion. We show that He^+ in superstrong magnetic fields can produce spectral lines with energies of up to about 3 keV, and it may be responsible for absorption features detected recently in the soft X-ray spectra of several radio-quiet isolated neutron stars. Quantization of the ion's motion across a magnetic field results in a fine structure of spectral lines, with a typical spacing of tens electron-volts in magnetar-scale fields. It also gives rise to ion cyclotron transitions, whose energies and oscillator strengths depend on the state of the bound ion.

A Chandra observation of the isolated hard X-ray source XMMU J061804.3+222732, located in the region of apparent interaction of the supernova remnant IC 443 with a molecular cloud, resolved the complex structure of the source in a few bright clumps embedded in an extended emission of a ~ 30 arcsec size. The X-ray spectra of the clumps and the extended emission are dominated by a hard power-law component with a photon index of 1.2--1.4. In addition, we see some indications of an optically thin thermal plasma of a ~ 0.3 keV temperature. The observed X-ray morphology and spectra are consistent with those expected for an isolated supernova ejecta fragment interacting with a dense ambient medium. A possible alternative interpretation is a pulsar wind nebula associated with either IC 443 or another SNR, G189.6+3.3.

Radio pulsar subpulse drifting has been interpreted as rotation of sub-beams (sparks) of pair plasma produced by intermittent breakdowns of an inner vacuum gap above the pulsar polar cap. This model also predicts strong thermal X-ray emission from the polar cap caused by inflowing particles created in spark discharges. We have observed the best-studied drifting pulsar B0943+10 with XMM-Newton and detected a point source coincident with the radio pulsar position. Its spectrum could be fitted with a thermal blackbody model, although a power-law model is also acceptable. The thermal fit gives a bolometric luminosity L_bol ~ 5 x 10^{28} erg/s and a surface area A ~ 10^3 (T/3MK)^{-4} m^2, much smaller than the conventional polar cap area, 6 x 10^4 m^2. Such thermal radiation can be interpreted as emitted from footprints of sparks drifting in an inner gap of a height h ~ 0.1 - 0.2 r_pc, where r_pc is the polar cap radius. However, the original vacuum gap model by Ruderman and Sutherland requires some modification to reconcile the X-ray and radio data.

We describe time-tagged observations of the Vela pulsar (PSR B0833-45) with the HST STIS MAMA-UV detectors. Using an optimal extraction technique, we obtain a NUV light curve and crude FUV phase-resolved spectra. The pulse is dominated by complex non-thermal emission at all phases. We obtain approximate spectral indices for the brightest pulse components, relate these to pulsations observed in other wavebands, and constrain the thermal surface emission from the neutron star. An upper limit on the Rayleigh-Jeans component at UV pulse minimum suggests a rather low mean surface temperature T_eff^\infty <4.6 x 10^5 (d_{300}/R_{14.2})^2 K. If confirmed by improved phase-resolved spectroscopy, this observation has significant impact on our understanding of neutron star interiors.

We observed the gamma-ray pulsar Geminga with the HST STIS/MAMA detectors to measure the Geminga's UV spectrum and pulsations. The slope of the far-UV (FUV) spectrum is close to that of a Rayleigh-Jeans spectrum, suggesting that the FUV radiation is dominated by thermal emission from the neutron star (NS) surface. The measured FUV flux, F_FUV=(3.7+/-0.2)x10^-15 ergs cm^-2 s^-1 in 1155-1702 A band, corresponds to a brightness temperature T_RJ=(0.3-0.4)(d_200/R_13)^2 MK. The soft thermal component of the Geminga's X-ray spectrum measured with the XMM-Newton observatory corresponds to a temperature T_s=0.49+/-0.01 MK and radius R_s = (12.9+/-1.0)d_200 km. Unlike other NSs detected in the UV-optical, the FUV spectrum of Geminga lies below the extrapolation of the soft thermal component, which might be associated with the Geminga's very low temperature. Surprisingly, the thermal FUV radiation is strongly pulsed, showing a narrow dip at a phase close to that of a broader minimum of the soft X-ray light curve. The strong pulsations might be attributed to partial occultations of the thermal UV radiation by regions of the magnetosphere filled with electron/positron plasma. In contrast with the FUV spectrum, the near-IR (NIR) through near-UV (NUV) spectrum of Geminga is clearly nonthermal. It can be described by a power-law model, F_nu \propto nu^{-alpha} with a slope alpha = 0.43+/-0.15, close to alpha = 0.56+/-0.24 of the hard X-ray (E>2.5 keV) magnetospheric component. The observed NIR-optical-NUV spectrum is marginally consistent with the extrapolation of the X-ray magnetospheric spectrum. The NUV pulsations, however, do not show a clear correlation with the hard X-ray pulsations.

We present the results of deep, high-resolution, multi-band optical observations of the field of the young (~ 5,000 yrs) 16 ms X-ray pulsar PSR J0537$-$6910 performed with the Advanced Camera for Surveys (ACS) aboard the Hubble Space Telescope (HST). Although a few new potential counterparts have been detected within or close to the revised Chandra X-ray error circle (~1 arcsec) of the pulsar, only two of them (with magnitudes m_{814W} ~ 23.9 and m_{814W} ~ 24.2) show indications of a peculiar spectrum which could be related to optical emission from the pulsar. This might be true also for a third, fainter, candidate detected only in one filter (with magnitude m_{814W} ~ 26.7). If either of the two brighter candidates is indeed the actual counterpart, the optical output of PSR J0537-6910 would make it similar to young Crab-like pulsars. If not, it would mean that PSR J0537-6910 is significantly underluminous with respect to all pulsars detected in the optical

We present the timing and spectral analyses of theXMM-newton data on the 17-Myr-old, nearby radio pulsar B0950+08. This observation revealed pulsations of the X-ray flux of the pulsar at its radio period. The pulse shape and pulsed fraction are apparently different at lower and higher energies of the observed 0.2-10 keV energy range, which suggests that the radiation cannot be explained by a single emission mechanism. The X-ray spectrum of the pulsar can be fitted with a power-law model with a photon index about 1.75 and an (isotropic) luminosity about 9.8e29 erg/s in the 0.2-10 keV. Better fits are obtained with two-component, power-law plus thermal, models with index of 1.30 and 9.7e29 erg/s for the power-law component that presumably originates from the pulsar's magnetosphere. The thermal component, dominating at E>0.7 keV, can be interpreted as radiation from heated polar caps on the neutron star surface covered with a hydrogen atmosphere. The inferred effective temperature, radius, and bolometric luminosity of the polar caps are about 1 MK, 250 m, and 3e29 erg/s. Optical through X-ray nonthermal spectrum of the pulsar can be described as a single power-law with index 1.3-1.4 for the two-component X-ray fit. The ratio of the nonthermal X-ray (1-10 keV) luminosity to the nonthermal optical (4000-9000 \AA) luminosity is within the range of 1e2-1e3 observed for younger pulsars, which suggests that the magnetospheric X-ray and optical emissions are powered by the same mechanism in all pulsars. An upper limit on the temperature of the bulk of the neutron star surface, inferred from the optical and X-ray data, is about 0.15 MK. We also analyze X-ray observations of several other old pulsars, B2224+65, J2043+2740, B0628-28, B1813-36, B1929+10, and B0823+26.

We present spectral analysis of the Crab Nebula obtained with the {\it Chandra} X-ray observatory. The X-ray spectrum is characterized by a power-law whose index varies across the nebula. The variation can be discussed in terms of the particle injection from the pulsar in two different directions: the equatorial plane containing the torus and the symmetry axis along the jet. In the equatorial plane, spectra within the torus are the hardest, with photon index $\alpha \approx 1.9$, and are almost independent of the surface brightness. At the periphery of the torus, the spectrum gradually softens in the outer, lower surface brightness regions, up to $\alpha \approx 3.0$. This indicates that synchrotron losses become significant to X-ray emitting particles at the outer boundary of the torus. We discuss the nature of the torus, incorporating information from observations at other wavelengths. Spectral variations are also seen within the southern jet. The core of the jet is the hardest with $\alpha \approx$ 2.0, and the outer sheath surrounding the core becomes softer with $\alpha$ up to 2.5 at the outermost part. Based on the similarity between the spectra of the jet core and the torus, we suggest that the electron spectra of the particles injected from the pulsar are also similar in these two different directions. The brightness ratio between the near and far sides of the torus can be explained by Doppler boosting and relativistic aberration; however, the observed ratio cannot be derived from the standard weakly magnetized pulsar wind model. We also found a site where an optical filament comprised of supernova ejecta is absorbing the soft X-ray emission ($<$ 2 keV).

The X-ray source 1E 1207.4-5209 is a compact central object in the G296.5+10.0 supernova remnant. Its spin period of 424 ms, discovered with the Chandra X-ray Observatory, suggests that it is a neutron star. The X-ray spectrum of this radio-quiet pulsar shows at least two absorption lines, first spectral features discovered in radiation from an isolated neutron star. Here we report the results of timing analysis of Chandra and XMM-Newton observations of this source showing a non-monotonous behavior of its period. We discuss three hypotheses which may explain the observational result. The first one assumes that 1E 1207.$-5209 is a glitching pulsar, with frequency jumps of \Delta f > 5 \muHz occurring every 1-2 years. The second hypothesis explains the deviations from a steady spin-down as due to accretion, with accretion rate varying from \sim 10^{13} to >10^{16} g s^{-1}, from a disk possibly formed from ejecta produced in the supernova explosion. Finally, the period variations could be explained assuming that the pulsar is in a wide binary system with a long period, P_orb \sim 0.2-6 yr, and a low-mass companion, M_2 < 0.3 M_\odot.

X-ray observations of middle-aged pulsars allow one to study nonthermal radiation from pulsar magnetospheres and thermal radiation from neutron star (NS) surfaces. In particular, from the analysis of thermal radiation one can infer the surface temperatures and radii of NSs, which is important for investigating evolution of these objects and constraining the equation of state of the superdense matter in the NS interiors. Here we present results of XMM observations of three middle-aged pulsars, J0538+2817, B0656+14 and J0633+1746 (Geminga), and briefly discuss mechanisms of their X-ray emission.

There are point-like sources in central regions of several supernova remnants which have not been detected outside the X-ray range. The X-ray spectra of these Central Compact Objects (CCOs) have thermal components with blackbody temperatures of 0.2-0.5 keV and characteristic sizes of 0.3-3 km. Most likely, the CCOs are neutron stars born in supernova explosions. We overview their observational properties, emphasizing the Chandra data, and compare them with magnetars.

We present year-scale morphological variations of the Crab Nebula revealed by Chandra X-ray observatory. Observations have been performed about every 1.7 years over the 3 years from its launch. The variations are clearly recognized at two sites: the torus and the southern jet. The torus, which had been steadily expanding until 1.7 years ago, now appears to have shrunk in the latest observation. Additionally, the circular structures seen to the northeast of the torus have decayed into several arcs. On the other hand, the southern jet shows the growth of its overall kinked-structure. We discuss the nature of these variations in terms of the pulsar wind nebula (PWN) mechanism.

We observed PSR J0437-4715 with the FUV-MAMA detector of the Hubble Space Telescope Imaging Spectrometer (STIS) to measure the pulsar's spectrum and pulsations. For the first time, UV emission from a millisecond pulsar is detected. The measured flux, $(2.0\pm 0.2)\times 10^{-15}$ erg s$^{-1}$ cm$^{-2}$ in the 1150-1700 \AA range, corresponds to the luminosity $L_{\rm FUV}=(4.7\pm 0.5)\times 10^{27}$ erg s$^{-1}$, for the distance of 140 pc and negligible interstellar extinction. The shape of the observed spectrum suggests thermal emission from the neutron star surface with a surprisingly high temperature of about $1\times 10^5$ K, above the upper limit on the surface temperature of the younger ``ordinary'' pulsar J0108-1431. For the few-Gyr-old J0437-4715, such a temperature requires a heating mechanism to operate. The spectrum of J0437-4715 shows marginal evidence of an emission line at 1372 \AA, which might be a gravitationally redshifted Zeeman component of the Hydrogen Ly$\alpha$ line in a magnetic field $\sim 7\times 10^8$ G. No pulsations are detected, with a $3\sigma$ upper limit of 50% on pulsed fraction.

A series of 13 Chandra observations provided the deepest images of the Vela PWN yet available. In addition to the fine structure of the inner PWN features, a much larger and fainter asymmetric X-ray nebula emerges in the summed images. The shape of this outer PWN is similar to that of the radio PWN. We also present the spectral map of the Vela PWN that reveals a shell of soft emission surrounding the inner PWN and an extended region of harder emission south-west of the pulsar. This may indicate that the outer jet is supplying particles to this region.

Chandra observations of the Vela pulsar-wind nebula (PWN) have revealed a jet in the direction of the pulsar's proper motion, and a counter-jet in the opposite direction, embedded in diffuse nebular emission. The jet consists of a bright, 8''-long inner jet, between the pulsar and the outer arc, and a dim, curved outer jet that extends up to 100'' in approximately the same direction. From the analysis of thirteen Chandra observations spread over about 2.5 years we found that this outer jet shows particularly strong variability, changing its shape and brightness. We observed bright blobs in the outer jet moving away from the pulsar with apparent speeds (0.3-0.6)c and fading on time-scales of days to weeks. The spectrum of the outer jet fits a power-law model with a photon index of 1.3\pm0.1. The X-ray emission of the outer jet can be interpreted as synchrotron radiation of ultrarelativistic electrons/positrons. This interpretation allows one to estimate the magnetic field, ~100 microGauss, maximum energy of X-ray emitting electrons, ~2\times 10^{14} eV, and energy injection rate, ~8\times 10^{33} erg/s, for the outer jet. In the summed PWN image we see a dim, 2'-long outer counter-jet, which also shows a power-law spectrum with photon ined of 1.2-1.5. Southwest of the jet/counter-jet an extended region of diffuse emission is seen. Relativistic particles responsible for this radiation are apparently supplied by the outer jet.

Observations of the Vela pulsar region with the Chandra X-ray observatory have revealed the fine structure of its synchrotron pulsar-wind nebula (PWN), which showed an overall similarity with the Crab PWN. However, contrary to the Crab, no firm detection of the Vela PWN in optical has been reported yet. To search for the optical counterpart of the X-ray PWN, we analyzed deep optical observations performed with different telescopes. We compared the optical images with those obtained with the Chandra ACIS to search for extended emission patterns which could be identified as counterparts of the X-ray nebula elements. Although some features are seen in the optical images, we find no correlation with the X-ray structure. Thus, we conclude that the diffuse optical emission is more likely associated with filaments in the host Vela SNR. The derived upper limits on the optical flux from the PWN are compatibile, within the uncertainties, with the values expected on the basis of the extrapolations of the X-ray data.

Observations of the Vela pulsar-wind nebula (PWN) with the Chandra X-ray Observatory have revealed a complex, variable PWN structure, including inner and outer arcs, a jet in the direction of the pulsar's proper motion, and a counter-jet in the opposite direction, embedded in diffuse nebular emission. From the analysis of thirteen Chandra observations spread over about 2.5 years we found that this outer jet shows particularly strong variability, changing its shape and brightness. We observed bright blobs in the outer jet moving away from the pulsar with apparent speeds (0.3-0.6)c and fading on time-scales of days to weeks. The spectrum of the outer jet fits a power-law model with a photon index of 1.3. The X-ray emission of the outer jet can be interpreted as synchrotron radiation of ultrarelativistic electrons/positrons. This interpretation allows one to estimate the magnetic field, ~100 microGauss, maximum energy of X-ray emitting electrons, $2 \times 10^{14}$ eV, and energy injection rate, $8\times 10^{33}$ erg/s, for the outer jet. In the summed PWN image, we see a faint, strongly bent, extension of the outer jet. The more extreme bends closer to the pulsar, as well as the apparent side motions of the outer jet, can be associated with kink instabilities of a magnetically confined, pinched jet flow. Another feature found in the summed image is a dim, 2'-long outer counter-jet, which also shows a power-law spectrum with photon index of 1.2-1.5. Southwest of the jet/counter-jet (i.e., approximately perpendicular to the direction of pulsar's proper motion), an extended region of diffuse emission is seen. Relativistic particles responsible for this radiation are apparently supplied by the outer jet.

Observations of thermal radiation from neutron stars allow one to measure the surface temperatures and confront them with cooling scenarios. Detection of gravitationally redshifted spectral lines can yield the mass-to-radius ratio. In the few cases when the distance is known, one can measure the neutron star radius, which is particularly important to constrain the equation of state of the superdense matter in the neutron star interiors. Finally, one can infer the chemical composition of the neutron star surface layers, which provides information about formation of neutron stars and their interaction with the environments. We present the observational results on thermal radiation from active pulsars and radio-quiet neutron stars, with emphasis on the results obtained with the Chandra X-ray Observatory and discuss some implications of these results. In particular, we present a new interpretation of the thermal emission detected in the broad energy range (from optical to X-rays) from the puzzling object RX J1856.5-3754

We report here the results of the first Chandra X-Ray Observatory observations of the globular cluster M28 (NGC 6626). 46 X-ray sources are detected, of which 12 lie within one core radius of the center. We show that the apparently extended X-ray core emission seen with the ROSAT HRI is due to the superposition of multiple discrete sources for which we determine the X-ray luminosity function down to a limit of about 6xE30 erg/s. For the first time the unconfused phase-averaged X-ray spectrum of the 3.05-ms pulsar B1821--24 is measured and found to be best described by a power law with photon index ~ 1.2. Marginal evidence of an emission line centered at 3.3 keV in the pulsar spectrum is found, which could be interpreted as cyclotron emission from a corona above the pulsar's polar cap if the the magnetic field is strongly different from a centered dipole. The unabsorbed pulsar flux in the 0.5--8.0 keV band is ~3.5xE-13 ergs/s/cm^2. Spectral analysis of the 5 brightest unidentified sources is presented. Based on the spectral parameters of the brightest of these sources, we suggest that it is a transiently accreting neutron star in a low-mass X-ray binary, in quiescence. Fitting its spectrum with a hydrogen neutron star atmosphere model yields the effective temperature T_eff^\infty = 90^{+30}_{-10} eV and the radius R_NS^\infty = 14.5^{+6.9}_{-3.8} km. In addition to the resolved sources, we detect fainter, unresolved X-ray emission from the central core of M28. Using the Chandra-derived positions, we also report on the result of searching archival Hubble Space Telescope data for possible optical counterparts.

New results of HST and VLT observations of isolated pulsars and their environments are presented. We present the first deep optical observations of the nearby pulsar PSR J0108-1431, the identification evidence of PSR 1929+10 based on the proper motion measurement of its optical counterpart and a deep investigation at optical wavelengths of the X-ray PWN around the Vela pulsar.

We present the results of first deep optical observations of the field of the old ($\sim 10^8$ yr), nearby, isolated pulsar J0108-1431, in an attempt to detect its optical counterpart. The observations were performed using the FORS1 instrument at the focus of the European Southern Observatory Antu Telescope of the VLT. Observations with the Australia Telescope Compact Array (ATCA) were made to determine an accurate position for the radio pulsar at the current epoch. The imaging data, obtained in the V, B, and U passbands reveal no counterpart at the revised radio position down to $V \simeq 28$, $B\simeq 28.6$ and $U \simeq 26.4$. For a distance of 130 pc, estimated from the pulsar's dispersion measure, our constraints on the optical flux put an upper limit of $T=4.5\times 10^4$ K for the surface temperature of the neutron star, assuming a stellar radius $R_\infty=13$ km. Our new radio position allows us to place an upper limit on the pulsar proper motion of 82 mas yr$^{-1}$ which, for d=130 pc, implies a transverse velocity $\la 50$ km sec$^{-1}$.

CONTENTS The Milky Way -- Pulsars and Isolated Neutron Stars / Introduction: Historical Overview / Physics and Astrophysics of Isolated Neutron Stars / Rotation-powered Pulsars: The Magnetic Braking Model / High-energy Emission Models / Magnetospheric Emission Models / Thermal Evolution of Neutron Stars / Photospheric Emission from Cooling Neutron Stars / The Current Picture of High-Energy Emission Properties of Isolated Neutron Stars / Young Neutron Stars in Supernova Remnants / Crab-like Pulsars / Vela-like pulsars / Radio-silent Neutron Stars in Supernova Remnants / Anomalous X-ray Pulsars and Soft Gamma-Ray Repeaters / Thermal Emission from Middle-Aged Pulsars / Old Nearby Radio Pulsars / Isolated Radio-quiet Neutron Stars / Recycled Millisecond Pulsars / Impressive Achievements and Great Expectations / References

We observed the compact central object CXOU J085201.4--461753 in the supernova remnant G266.2--1.2 (RX J0852.0--4622) with the Chandra ACIS detector in timing mode. The spectrum of this object can be described by a blackbody model with the temperature kT=404 eV and radius of the emitting region R=0.28 km, at a distance of 1 kpc. Power-law and thermal plasma models do not fit the source spectrum. The spectrum shows a marginally significant feature at 1.68 keV. Search for periodicity yields two candidate periods, about 301 ms and 33 ms, both significant at a 2.1 sigma level; the corresponding pulsed fractions are 13% and 9%, respectively. We find no evidence for long-term variability of the source flux, nor do we find extended emission around the central object. We suggest that CXOU J085201.4--461753 is similar to CXOU J232327.9+584842, the central source of the supernova remnant Cas A. It could be either a neutron star with a low or regular magnetic field, slowly accreting from a fossil disk, or, more likely, an isolated neutron star with a superstrong magnetic field. In either case, a conservative upper limit on surface temperature of a 10 km radius neutron star is about 90 eV, which suggests accelerated cooling for a reasonable age of a few thousand years.

We observed 1E 1207.4--5209, a neutron star in the center of the supernova remnant PKS 1209--51/52, with the ACIS detector aboard the Chandra X-ray observatory and detected two absorption features in the source spectrum. The features are centered near 0.7 keV and 1.4 keV, their equivalent widths are about 0.1 keV. We discuss various possible interpretations of the absorption features and exclude some of them. A likely interpretation is that the features are associated with atomic transitions of once-ionized helium in the neutron star atmosphere with a strong magnetic field. The first clear detection of absorption features in the spectrum of an isolated neutron star provides an opportunity to measure the mass-to-radius ratio and constrain the equation of state of the superdense matter.

The outstanding capabilities of the Chandra X-ray observatory have greatly increased our potential to observe and analyze thermal radiation from the surfaces of neutron stars (NSs). Such observations allow one to measure the surface temperatures and confront them with the predictions of the NS cooling models. Detection of gravitationally redshifted spectral lines can yield the NS mass-to-radius ratio. In rare cases when the distance is known, one can measure the NS radius, which is particularly important to constrain the equation of state of the superdense matter in the NS interiors. Finally, one can infer the chemical composition of the NS surface layers, which provides information about formation of NSs and their interaction with the environment. We overview the recent Chandra results on the thermal radiation from various types of NSs -- active pulsars, young radio-quiet neutron stars in supernova remnants, old radio-silent ``dim'' neutron stars -- and discuss their implications.

Models of thermal emission of neutron stars, presumably formed in their atmospheres, are needed to infer the surface temperatures, magnetic fields, chemical composition, and neutron star masses and radii from the observational data. This information, supplemented with model equations of state and neutron star cooling models, is expected to move us further in understanding the fundamental properties of the superdense matter in the neutron star interiors. The neutron star atmospheres are very different from those of usual stars due to the immense gravity and huge magnetic fields. In this presentation we review the current status of the neutron star atmosphere modeling and present most important results.

We obtained four Chandra/ACIS-S observations beginning two weeks after the end of the November 2000 outburst of the neutron star (NS) transient Aql X-1. Over the five month span in quiescence, the X-ray spectra are consistent with thermal emission from a NS with a pure hydrogen photosphere and R_{\infty}=15.9+{0.8}-{2.9} (d/5 kpc) km at the optically implied X-ray column density. We also detect a hard power-law tail during two of the four observations. The intensity of Aql X-1 first decreased by 50+/-4% over three months, then increased by 35+/-5% in one month, and then remained constant (<6% change) over the last month. These variations in the first two observations cannot be explained by a change in the power-law spectral component, nor in the X-ray column density. Presuming that R_{\infty} is not variable and a pure hydrogen atmosphere, the long-term changes can only be explained by variations in the NS effective temperature, from kT_{eff, \infty}=130+3-5 eV, down to 113+3-4 eV, finally increasing to 118+9-4 eV for the final two observations. During one of these observations, we observe two phenomena which were previously suggested as indicators of quiescent accretion onto the NS: short-timescale (<1e4 sec) variability (at 32+8-6% rms), and a possible absorption feature near 0.5 keV. The possible absorption feature can potentially be explained as due to a time-variable response in the ACIS detector. Even so, such a feature has not been detected previously from a NS, and if confirmed and identified, can be exploited for simultaneous measurements of the photospheric redshift and NS radius.

Second Chandra observation of 1E 1207.4-5209, the central source of the supernova remnant PKS 1209-51/52, allowed us to confirm the previously detected period of 424 ms and, assuming a uniform spin-down, estimate the period derivative, Pdot (0.7-3)e-14 s/s. The corresponding characteristic age of the pulsar, P/2Pdot~200-900 kyr, is much larger than the estimated age of the SNR, \~7 kyr. The values of the spin-down luminosity, Edot (0.4-1.6)e34 erg/s, and conventional magnetic field, B (2-4)e12 G, are typical for a middle-aged radio pulsar, although no manifestations of pulsar activity have been observed. If 1E 1207.4-5209 is indeed the neutron star formed in the same supernova explosion that created PKS 1209-51/52, such a discrepancy in ages could be explained either by a long initial period, close to its current value, or, less likely, by a very large braking index of the pulsar. Alternatively, the pulsar could be a foreground object unrelated to the supernova remnant, but the probability of such a coincidence is very low.

We report on spectral and timing observations of the nearest millisecond pulsar J0437-4715 with Chandra. The pulsar spectrum, detected up to 7 keV, cannot be described by a simple one-component model. We suggest that it consists of two components, a nonthermal power-law spectrum generated in the pulsar magnetosphere, with a photon index about 2, and a thermal spectrum emitted by heated polar caps, with a temperature decreasing outwards from 2 MK to 0.5 MK. The lack of spectral features in the thermal component suggests that the neutron star surface is covered by a hydrogen (or helium) atmosphere. The timing analysis shows one X-ray pulse per period, with a pulsed fraction of about 40% and the peak at the same pulse phase as the radio peak. No synchrotron pulsar-wind nebula is seen in the X-rays.

X-ray images of some young supernova remnants show bright point sources which have not been detected in radio, optical and gamma-ray bands. Despite the similarity of the X-ray spectra of these objects, they show a variety of temporal properties. Most likely, they are neutron stars whose properties (spin periods? magnetic fields? environments?) are different from those of radio and/or gamma-ray pulsars. We present an overview of observational results on several objects of this class -- the central sources of Cassiopeia A, RX J0852-4622, RCW 103, Puppis A, and PKS 1209-51/52 -- with emphasis on the recent Chandra observations.

We report the results of the spectral and timing analysis of observations of the Vela pulsar with the Chandra X-ray Observatory. The spectrum shows no statistically significant spectral lines in the observed 0.25--8.0 keV band. It consists of two distinct continuum components. The softer component can be modeled as either a magnetic hydrogen atmosphere spectrum with kT = 59 +- 3 eV, R = 15.5 +- 1.5 km, or a standard blackbody with kT = 129 +- 4 eV, R = 2.1 +- 0.2 km (the radii are for a distance of 250 pc). The harder component, modeled as a power-law spectrum, gives photon indices depending on the model adopted for the soft component: gamma = 1.5 +- 0.3 for the magnetic atmosphere soft component, and gamma = 2.7 +- 0.4 for the blackbody soft component. Timing analysis shows three peaks in the pulse profile, separated by about 0.3 in phase. Energy-resolved timing provides evidence for pulse profile variation with energy. The higher energy (E > 1.8 keV) profile shows significantly higher pulsed fraction.

We present a series of monitoring observations of the Crab Nebula with the Chandra X-ray Observatory, focusing on the temporal evolution of the structure. This series of 8 observations, spanning a period of approximately six months, shows the dynamic nature of the inner X-ray structures. We detected outward moving ``wisps'' from the recently discovered inner ring seen in optical observations. We also find that the inner ring itself shows temporal variations in structure. The torus also appears to be expanding. Such temporal variations generally match the canonical scenario that an expanding synchrotron nebula injected from the pulsar is confined by the supernova ejecta.

G21.5-0.9 is a plerionic supernova remnant (SNR) whose properties hint at the presence of a pulsar -- yet no pulsations have been found at any wavelength. Early observations with Chandra led to the discovery of an extended component, making the SNR at least twice as big as originally thought. Our analysis indicates that this low-surface brightness extended component is non-thermal. We perform a spatially resolved spectroscopy using a 72 ksec exposure with ACIS-S, and find no evidence of line emission. The 5-arcminute diameter remnant is well fitted with a power law with a photon index steepening from 1.5 (central core) to 2.7 (extended component) and a column density N_H=2.2 10^22 cm^-2. Using a 76 ksec exposure with the HRC, we derive an upper limit of 16% on the pulsed fraction from a putative pulsar. We also infer the parameters of the `hidden pulsar' in G21.5-0.9 in comparison with the Crab and 3C 58 pulsars. This remnant remains unique and intriguing since it is, to date, the only candidate whose size is bigger in X-rays than in the radio.

We have implemented a Monte Carlo algorithm to model and predict the response of various kinds of CCDs to X-ray photons and minimally-ionizing particles and have applied this model to the CCDs in the Chandra X-ray Observatory's Advanced CCD Imaging Spectrometer. This algorithm draws on empirical results and predicts the response of all basic types of X-ray CCD devices. It relies on new solutions of the diffusion equation, including recombination, to predict the radial charge cloud distribution in field-free regions of CCDs. By adjusting the size of the charge clouds, we can reproduce the event grade distribution seen in calibration data. Using a model of the channel stops developed here and an insightful treatment of the insulating layer under the gate structure developed at MIT, we are able to reproduce all notable features in ACIS calibration spectra. The simulator is used to reproduce ground and flight calibration data from ACIS, thus confirming its fidelity. It can then be used for a variety of calibration tasks, such as generating spectral response matrices for spectral fitting of astrophysical sources, quantum efficiency estimation, and modeling of photon pile-up.

The central region of the recently discovered supernova remnant RX J0852.0-4622 was observed with the ACIS detector aboard the Chandra X-ray Observatory. We found only one relatively bright source, about 4' north of the SNR center, with a flux of $\sim 2\times 10^{-12}$ erg s$^{-1}$ cm$^{-2}$ in the 0.5--10 keV band. The position of this point-like source, CXOU J085201.4-461753, rules out its association with the two bright stars in the field, HD 76060 and Wray 16-30. Observations of the field with the CTIO 0.9-m telescope show a star ($R\approx 17$, $B\approx 19$) at about 2\farcs4 from the nominal X-ray position. We consider association of this star with the X-ray source unlikely and estimate a limiting magnitude of the optical counterpart as $B \ge 22.5$ and $R \ge 21.0$. Based on the X-ray-to-optical flux ratio, we argue that the X-ray source is likely the compact remnant of the supernova explosion that created the RX J0852.0-4622 SNR. The observed X-ray spectrum of the source is softer than spectra of magnetospheric radiation of rotation-powered pulsars, but it is harder than spectra of cooling neutron stars emitting thermal radiation from the entire surface, similar to the central compact source of the Cas A SNR. We suggest that CXOU J085201.4-461753 belongs to the growing family of radio-quiet compact central sources, presumably neutron stars, recently discovered in a number of SNRs.

(Abridged). The type-I X-ray bursting low mass X-ray binary KS 1731-260 was recently detected for the first time in quiescence by Wijnands et al., following an approximately 13 yr outburst which ended in Feb 2001. Unlike all other known transient neutron stars, the duration of this recent outburst is as long as the thermal diffusion time of the crust. The large amount of heat deposited by reactions in the crust will have heated the crust to temperatures much higher than the equilibrium core temperature. As a result, the thermal luminosity currently observed from the neutron star is dominated not by the core, but by the crust. Moreover, the level and the time evolution of quiescent luminosity is determined mostly by the amount of heat deposited in the crust during the most recent outburst. Using estimates of the outburst mass accretion rate, our calculations of the quiescent flux immediately following the end of the outburst agree with the observed quiescent flux to within a factor of a few. We present simulations of the evolution of the quiescent lightcurve for different scenarios of the crust microphysics, and demonstrate that monitoring observations (with currently flying instruments) spanning from 1--30 yr can measure the crust cooling timescale and the total amount of heat stored in the crust. These quantities have not been directly measured for any neutron star.

We present the analysis of archival X-ray observations of the supernova remnant (SNR) G21.5-0.9. Based on its morphology and spectral properties, G21.5-0.9 has been classified as a Crab-like SNR. In their early analysis of the CHANDRA calibration data, Slane et al. (2000) discovered a low-surface-brightness, extended emission. They interpreted this component as the blast wave formed in the supernova (SN) explosion. In this paper, we present the CHANDRA analysis using a total exposure of ~150 ksec. We also include ROSAT and ASCA observations. Our analysis indicates that the extended emission is non-thermal -- a result in agreement with XMM observations. The entire remnant of radius ~ 2'.5 is best fitted with a power law model with a photon index steepening away from the center. The total unabsorbed flux in the 0.5-10 keV is 1.1E-10 erg/cm2/s with an 85% contribution from the 40" radius inner core. Timing analysis of the High-Resolution Camera (HRC) data failed to detect any pulsations. We put a 16% upper limit on the pulsed fraction. We derive the physical parameters of the putative pulsar and compare them with those of other plerions (such as the Crab and 3C 58). G21.5-0.9 remains the only plerion whose size in X-rays is bigger than in the radio. Deep radio observations will address this puzzle.

Using the Chandra/ACIS-I detector, we have identified an X-ray source (CXOU 132619.7-472910.8) in the globular cluster NGC 5139 with a thermal spectrum identical to that observed from transiently accreting neutron stars in quiescence. The absence of intensity variability on timescales as short as 4 seconds (< 25% rms variability) and as long as 5 years (<50% variability) supports the identification of this source as a neutron star, most likely maintained at a high effective temperature (approximately 1e6 K) by transient accretion from a binary companion. The ability to spectrally identify quiescent neutron stars in globular clusters (where the distance and interstellar column densities are known) opens up new opportunities for precision neutron star radius measurements.

We report on the quiescent spectrum measured with Chandra/ACIS-S of the transient, type-I X-ray bursting neutron star Aql X-1, immediately following an accretion outburst. The neutron star radius, assuming a pure hydrogen atmosphere and hard power-law spectrum, is $R_\infty$=13.4{+5}{-4} (d/5 \kpc) km. Based on the historical outburst record of RXTE/ASM, the quiescent luminosity is consistent with that predicted by Brown, Bildsten and Rutledge from deep crustal heating, lending support to this theory for providing a minimum quiescent luminosity of transient neutron stars. While not required by the data, the hard power-law component can account for 18+/-8% of the 0.5-10 keV thermal flux. Short-timescale intensity variability during this observation is less than 15% rms (3 sigma; 0.0001-1 Hz, 0.2-8 keV). Comparison between the Chandra spectrum and three X-ray spectral observations made between Oct 1992 and Oct 1996 find all spectra consistent with a pure H atmosphere, but with temperatures ranging from 145--168 eV, spanning a factor of 1.87+/-0.21 in observed flux. The source of variability in the quiescent luminosity on long timescales (greater than years) remains a puzzle. If from accretion, then it remains to be explained why the quiescent accretion rate provides a luminosity so nearly equal to that from deep crustal heating.

The observations of the pulsar-wind nebula (PWN) around the Vela pulsar with the Advanced CCD Imaging Spectrometer aboard the Chandra X-ray Observatory, taken on 2000 April 30 and November 30, reveal its complex morphology reminiscent of that of the Crab PWN. Comparison of the two observations shows changes up to 30% in the surface brightness of the PWN features. Some of the PWN elements show appreciable shifts, up to a few arcseconds (about 10^{16} cm), and/or spectral changes. To elucidate the nature of the observed variations, further monitoring of the Vela PWN is needed.

We report the results of the spectral analysis of two observations of the Vela pulsar with the Chandra X-ray observatory. The spectrum of the pulsar does not show statistically significant spectral lines in the observed 0.25-8.0 keV band. Similar to middle-aged pulsars with detected thermal emission, the spectrum consists of two distinct components. The softer component can be modeled as a magnetic hydrogen atmosphere spectrum - for the pulsar magnetic field $B=3\times 10^{12}$ G and neutron star mass $M=1.4 M_\odot$ and radius $R^\infty =13$ km, we obtain $\tef^\infty =0.68\pm 0.03$ MK, $L_{\rm bol}^\infty = (2.6\pm 0.2)\times 10^{32}$ erg s$^{-1}$, $d=210\pm 20$ pc (the effective temperature, bolometric luminosity, and radius are as measured by a distant observer). The effective temperature is lower than that predicted by standard neutron star cooling models. A standard blackbody fit gives $T^\infty =1.49\pm 0.04$ MK, $L_{\rm bol}^\infty=(1.5\pm 0.4)\times 10^{32} d_{250}^2$ erg s$^{-1}$ ($d_{250}$ is the distance in units of 250 pc); the blackbody temperature corresponds to a radius, $R^\infty =(2.1\pm 0.2) d_{250}$ km, much smaller than realistic neutron star radii. The harder component can be modeled as a power-law spectrum, with parameters depending on the model adopted for the soft component - $\gamma=1.5\pm 0.3$, $L_x=(1.5\pm 0.4)\times 10^{31} d_{250}^2$ erg s$^{-1}$ and $\gamma=2.7\pm 0.4$, $L_x=(4.2\pm 0.6)\times 10^{31} d_{250}^2$ erg s$^{-1}$ for the hydrogen atmosphere and blackbody soft component, respectively ($\gamma$ is the photon index, $L_x$ is the luminosity in the 0.2--8 keV band). The extrapolation of the power-law component of the former fit towards lower energies matches the optical flux at $\gamma\simeq 1.35$--1.45.

We report on spectral and intensity variability analysis from a Chandra/ACIS-S observation of the transient, type-I X-ray bursting low-mass X-ray binary Cen X-4. The quiescent X-ray spectrum during this observation is statistically identical to one observed previously with Beppo/SAX, and close, but not identical, to one observed previously with ASCA. The X-ray spectrum is best described as a pure Hydrogen atmosphere thermal spectrum plus a power-law component that dominates the spectrum above 2 keV. The best-fit radius of the neutron star is r=12.9+/-2.6 (d/1.2 kpc) km if the interstellar absorption is fixed at the value implied by the optical reddening. Allowing the interstellar absorption to be a free parameter yields r=19+45-10 (d/1.2 kpc) km (90% confidence). The thermal spectrum from the neutron star surface is inconsistent with a solar metallicity. We find a 3sigma upper-limit of root-mean-square variability <18% (0.2-2.0 keV; 0.0001-1 Hz) during the observation. On the other hand, the 0.5-10.0 keV luminosity decreased by 40+/-8% in the 4.9 years between the Asca and Chandra observations. This variability can be attributed to the power-law component. Moreover, we limit the variation in thermal temperature to <10% over these 4.9 years. The stability of the thermal temperature and emission area radius supports the interpretation that the quiescent thermal emission is due to the hot neutron star core.

HST observations have contributed significantly to our knowledge on the behaviour of Isolated Neutron Stars (INSs) as optical emitters. First, HST has been instrumental both to discover new optical counterparts (PSR B1055-52, PSR B1929+10, PSR B0950+08) and to confirm proposed identifications (PSR B0656+14). Second, HST multicolor photometry provided useful information to characterize the optical emission mechanism(s) at work in middle-aged INSs like PSR B0656+14 and Geminga. Last, but not least, the superior angular resolution of the HST allowed both to perform a very accurate morphological study of the plerionic environments of young pulsars (e.g. the Crab and PSR B0540-69) and to perform very accurate astrometric measurements yielding proper motions (Crab, Vela, Geminga, PSR B0656+14) and parallaxes (Geminga).

We present the results of the broad-band photometry of the nearby middle-aged radio pulsar PSR B0656+14 and its neighborhood obtained with the 6-meter telescope of the SAO RAS and with the Hubble Space Telescope. The broad-band spectral flux $F_\nu$ of the pulsar decreases with increasing frequency in the near-IR range and increases with frequency in the near-UV range. The increase towards UV can be naturally interpreted as the Rayleigh-Jeans tail of the soft thermal component of the X-ray spectrum emitted from the surface of the cooling neutron star. Continuation of the power-law component, which dominates in the high-energy tail of the X-ray spectrum, to the IR-optical-UV frequencies is consistent with the observed fluxes. This suggests that the non-thermal pulsar radiation may be of the same origin in a broad frequency range from IR to hard X-rays. We also studied 4 objects detected in the pulsar's 5" neighborhood.

The central source of the supernova remnant PKS 1209-52 was observed with the Advanced CCD Imaging Spectrometer aboard Chandra X-ray observatory on 2000 January 6-7. The use of the Continuos Clocking mode allowed us to perform the timing analysis of the data with time resolution of 2.85 ms and to find a period P=0.42412927+/-2.3e-7 s. The detection of this short period proves that the source is a neutron star. It may be either an active pulsar with unfavorably directed radio beam or a truly radio-silent neutron star whose X-ray pulsations are caused by a nonuniform distribution of surface temperature. To infer the actual properties of this neutron star, the period derivative should be measured.

The central pointlike X-ray source of the Cas A supernova remnant was discovered in the Chandra First Light Observation and found later in the archival ROSAT and Einstein images. The analysis of these data does not show statistically significant variability of the source. The power-law fit yields the photon index 2.6-4.1, and luminosity (2-60)e34 erg/s, for d=3.4 kpc. The power-law index is higher, and the luminosity lower, than those observed fromvery young pulsars. One can fit the spectrum equally well with a blackbody model with T=6-8 MK, R=0.2-0.5 km, L=(1.4-1.9)e33 erg/s. The inferred radii are too small, and the temperatures too high, for the radiationcould be interpreted as emitted from the whole surface of a uniformly heated neutron star. Fits with the neutron star atmosphere models increase the radius and reduce the temperature, but these parameters are still substantially different from those expected for a young neutron star. One cannot exclude, however, that the observed emission originates from hot spots on a cooler neutron star surface. Because of strong interstellar absorption, the possible low-temperature component gives a small contribution to the observed spectrum; an upper limit on the (gravitationally redshifted) surface temperature is < 1.9-2.3 MK. Amongst several possible interpretations, we favor a model of a strongly magnetized neutron star with magnetically confined hydrogen or helium polar caps on a cooler iron surface. Alternatively, the observed radiation may be interpreted as emitted by a compact object (more likely, a black hole) accreting from a fossil disk or from a late-type dwarf in a close binary.

Since the opacity of a magnetized plasma depends on polarization of radiation, the radiation emergent from atmospheres of neutron stars with strong magnetic fields is expected to be strongly polarized. The degree of linear polarization, typically ~10-30%, depends on photon energy, effective temperature and magnetic field. The spectrum of polarization is more sensitive to the magnetic field than the spectrum of intensity. Both the degree of polarization and the position angle vary with the neutron star rotation period so that the shape of polarization pulse profiles depends on the orientation of the rotational and magnetic axes. Moreover, as the polarization is substantially modified by the general relativistic effects, observations of polarization of X-ray radiation from isolated neutron stars provide a new method for evaluating the mass-to-radius ratio of these objects, which is particularly important for elucidating the properties of the superdense matter in the neutron star interiors.

We fit hydrogen atmosphere models to the X-ray data for four neutron stars (three from a previous paper, plus 4U 2129+47) and six black hole candidates (A0620-00, GS 2000+25, GS 1124-68, GS 2023+33, GRO J1655-40, and GRO J0422+32). While the neutron stars are similar in their intrinsic X-ray spectra (similar effective temperatures and emission area radii ~10 km), the spectra of two black hole candidates are significantly different, and the spectra of the remaining four are consistent with a very large parameter space that includes the neutron stars. The spectral differences between the neutron stars and black hole candidates favors the interpretation that the quiescent neutron star emission is predominantly thermal emission from the neutron star surface. Our work suggests that an X-ray spectral comparison in quiescence provides an additional means for distinguishing between neutron stars and black holes. The faint X-ray sources in globular clusters are also a class of objects which can be investigated in this manner.

There are several supernova remnants which contain unresolved X-ray sources close to their centers, presumably radio-quiet neutron stars. To prove that these objects are indeed neutron stars, to understand the origin of their X-ray radiation, and to explain why they are radio-quiet, one should know their periods and period derivatives. We searched for pulsations of the X-ray flux from the radio-quiet neutron star candidate RX J0822-4300 near the center of the Puppis A supernova remnant observed with the ROSAT PSPC and HRI. A standard timing analysis of the separate PSPC and HRI data sets does not allow one to detect the periodicity unequivocally. However, a thorough analysis of the two observations separated by 4.56 yr enabled us to find a statistically significant period $P\simeq 75.3$ ms and its derivative $\dot{P}\simeq 1.49\times 10^{-13}$ s s$^{-1}$. The corresponding characteristic parameters of the neutron star, age $\tau=P/(2\dot{P})=8.0$ kyr, magnetic field $B=3.4\times 10^{12}$ G, and rotational energy loss $\dot{E}=1.4\times 10^{37}$ erg s$^{-1}$, are typical for young radio pulsars. Since the X-ray radiation has a thermal-like spectrum, its pulsations may be due to a nonuniform temperature distribution over the neutron star surface caused by anisotropy of the heat conduction in the strongly magnetized crust.

We re-analyze the available X-ray spectral data of the type I bursting neutron star transients Aql X-1, Cen X-4, and 4U 1608-522 using realistic hydrogen atmosphere models. Previous spectral fits assumed a blackbody spectrum; because the free-free dominated photospheric opacity decreases with increasing frequency, blackbody spectral fits overestimate the effective temperature and underestimate, by as much as two orders of magnitude, the emitting area. Hydrogen atmosphere spectral models, when fit to the available observational data, imply systematically larger emission area radii, consistent with the canonical 10 km radius of a neutron star. This suggests that a substantial fraction of the quiescent luminosity is thermal emission from the surface of the neutron star. The magnitude of the equivalent hydrogen column density toward these systems, however, presents a considerable systematic uncertainty, which can only be eliminated by high signal-to-noise X-ray spectral measurements (e.g., with AXAF or XMM) which would permit simultaneous determination of the equivalent hydrogen column density, emission area, and thermal temperature.

We re-analyzed soft X-ray data collected with the ROSAT and ASCA observatories on a candidate neutron star (NS) near the center of the supernova remnant PKS 1209-52. We fitted the observed spectra with NS atmosphere models. The hydrogen atmosphere fits yield more realistic parameters of the NS and the intervening hydrogen column than the traditional blackbody fit. In particular, for a NS of mass 1.4 M_\odot and radius 10 km, we obtained the NS surface temperature T_{eff}=(1.4-1.9)*10^6 K and distance d=1.6-3.3 kpc versus T=(4.2-4.6)*10^6 K and (implausible) d=11-13 kpc for the blackbody fit, at a 90% confidence level. Our fits suggest that the surface magnetic field is either very weak, B < 10^{10} G, or it exceeds 2*10^{12} G. The hydrogen column density inferred from the atmosphere fits, n_H=(0.7-2.2)*10^{21} cm^{-2}, agrees fairly well with independent estimates obtained from UV observations of nearby stars, radio data, and X-ray spectrum of the shell of the supernova remnant, whereas the blackbody and power-law fits give considerably lower and greater values, n_H=(0.2-0.4)*10^{21} and (5.2-7.0)*10^{21} cm^{-2}, respectively. The inferred NS surface temperature is consistent with standard NS cooling models.

Properties of X-ray radiation emitted from the polar caps of a radio pulsar depend not only on the cap temperature, size, and position, but also on the surface chemical composition, magnetic field, and neutron star's mass and radius. Fitting the spectra and the light curves with neutron star atmosphere models enables one to infer these parameters. As an example, we present here results obtained from the analysis of the pulsed X-ray radiation of a nearby millisecond pulsar J0437-4715. In particular, we show that stringent constraints on the mass-to-radius ratio can be obtained if orientations of the magnetic and rotation axes are known, e.g., from the radio polarization data.

We observed the middle-aged radio pulsar 0656+14 with a CCD detector at the 6-m telescope. Broadband BVRI images show the following magnitudes of the pulsar counterpart: B = 24.85 (+0.19, -0.16), V = 24.90 (+0.16, -0.14), R = 24.52 (+0.12, -0.11), I = 23.81 (+0.27, -0.21). We fitted the UV-optical (space + ground-based) data with a two-component model which combines a power law (non-thermal component) with a thermal spectrum emitted by the neutron star surface. The power law component, with the energy power-law index \alpha=1.5 (+1.1,-1.2), dominates in the observed range. Constraints on the thermal component correspond to the Rayleigh-Jeans parameter G\equiv T_6(R_{10}/d_{500})^2=4.1 (+2.1, -4.1), where $T=10^6T_6$ K is the brightness temperature, R_\infty = 10 R_{10}$ km is the neutron star radius as seen by a distant observer, and $d=500 d_{500}$ pc is the distance. The shape of the optical-UV spectrum of PSR 0656+14 differs considerably from those observed from other pulsars: the middle-aged Geminga and young Crab, Vela.

We show that the soft X-ray spectra and light curves observed with the ROSAT and EUVE from the closest known millisecond pulsar J0437--4715 can be interpreted as thermal radiation from two hot polar caps whose emitting layers (atmospheres) are comprised of hydrogen. The simplest model yields a uniform temperature of $(0.8-0.9)\times 10^6 K$ within a cap radius of $0.7-0.9 km$. The spectral fits indicate that the temperature may be nonuniformly distributed along the cap surface. The distribution can be approximated by a central core heated up to $(1-2)\times 10^6 K$ within a radius of $0.2-0.4 km$, surrounded by a colder rim with temperatures $(3-5)\times 10^5 K$ extending out to $2-6 km$. The polar cap interpretation implies low column densities, $(1-3)\times 10^{19} cm^{-2}$, and a high degree of ionization, $> 20$%, of the interstellar hydrogen towards the pulsar. The inferred bolometric luminosity of the polar caps, $(1.0-1.6)\times 10^{30} erg s^{-1}$, is in excellent agreement with the predictions of the slot-gap model of radio pulsars developed by Arons and his coworkers. Similar polar cap radiation should be emitted by other millisecond pulsars, although in some of them (e.~g., PSR B1821--24) the soft X-ray flux is dominated by the nonthermal radiation from pulsar magnetospheres.

The properties of X-ray radiation from the polar caps predicted by the radio pulsar models depend on the surface chemical composition, magnetic field and star's mass and radius as well as on the cap temperature, size and position. Fitting the radiation spectra and light curves with the neutron star atmosphere models enables one to infer these parameters. We present here results obtained from the analysis of the soft X-ray radiation of PSR J0437-4715. In particular, with the aid of radio polarization data, we put constraints on the pulsar mass-to-radius ratio.

The strong magnetic fields (B ~ 10^{12} - 10^{13} G) characteristic of neutron stars make all the properties of an atom strongly dependent on the transverse component K_\perp of its generalized momentum. In particular, the photoionization process is modified substantially: (i) threshold energies are decreased as compared with those for an atom at rest, (ii) cross section values are changed significantly, and (iii) selection rules valid for atoms at rest are violated by the motion so that new photoionization channels become allowed. To calculate the photoionization cross sections, we, for the first time, employ exact numerical treatment of both initial and final atomic states. This enables us to take into account the quasi-bound (autoionizing) atomic states as well as coupling of different ionization channels. We extend the previous consideration, restricted to the so-called centered states corresponding to relatively small values of K_\perp, to arbitrary states of atomic motion. We fold the cross sections with the thermal distribution of atoms over K. For typical temperatures of neutron star atmospheres, the averaged cross sections differ substantially from those of atoms at rest. In particular, the photoionization edges are strongly broadened by the thermal motion of atoms; this "magnetic broadening" exceeds the usual Doppler broadening by orders of magnitude. The decentered states of the atoms give rise to the low-energy component of the photoionization cross section. This new component grows significantly with increasing temperature above 10^{5.5} K and decreasing density below 1 g/cm^3, i.e., for the conditions expected in atmospheres of middle-aged neutron stars.

We show that an analysis of multiwavelength observations of isolated neutron stars based on neutron star atmosphere models can be used not only to evaluate the neutron star effective temperature, but also to determine chemical composition of its surface. To demonstrate how this new method can be applied to a specific object, we chose the old isolated neutron star candidate RXJ1856.5-3754, whose soft X-ray radiation has been studied recently by Walter et al. (1996). We fitted the soft X-ray spectrum of this object with neutron star atmosphere models of different chemical compositions and used these fits to calculate the source spectrum over a broad wavelength range. We showed, in particular, that the optical/UV flux expected from this object depends drastically on the composition of its surface. In particular, the neutron star covered with hydrogen would be 5-6 magnitudes brighter than the neutron star with an iron surface. The object should also be observable with $EUVE$; the EUV flux is expected to be almost twice higher for the iron surface than for the hydrogen one. Thus, multiwavelength observations of this object would enable one to examine, for the first time, chemical composition of the neutron star surface. The method proposed can be applied to other nearby isolated neutron stars.

High magnetic fields in neutron stars, B ~ 10^{11} - 10^{13} G, substantially modify the properties of atoms and their interaction with radiation. In particular, the photoionization cross section becomes anisotropic and strongly polarization dependent. In a number of previous works based on the adiabatic approximation the conclusion was drawn that the transverse cross section vanishes for frequencies smaller than the electron cyclotron frequency. In other works (which employed a different form of the interaction operator) appreciable finite values were obtained. An adequate interpretation of the neutron star thermal-like radiation requires a resolution of this controversy. In this work the atomic wave functions for both discrete and continuum states are calculated by solving the coupled channel equations allowing the admixture between different Landau levels, which provides much higher accuracy than the adiabatic approximation. This enables to resolve the above contradiction in favour of the finite transverse cross sections. The non-adiabatic treatment of the continuum includes coupling between closed and open channels, which leads to the autoionization of quasi-bound energy levels associated with the electron cyclotron excitations and gives rise to Beutler - Fano resonances of the photoionization cross section. Autoionization widths of these quasi-bound levels are calculated and compared with the radiative widths. The results are important for investigations of the radiation emergent from the surface layers of neutron stars.

We present a detailed investigation of neutron star atmospheres with low magnetic fields, $B < 10^8 -10^{10}$ G, which do not affect opacities and equation of state of the atmospheric matter. We compute the atmospheric structure, emergent spectral fluxes and specific intensities for hydrogen, helium and iron atmospheres in a wide domain of effective temperatures and gravitational accelerations expected for neutron stars. The iron atmospheres are computed with the opacities and equations of state from the OPAL opacity library. We show that the model atmosphere spectra are substantially different from the blackbody spectra. For light element atmospheres, the flux is greater than the blackbody flux, and the spectrum is harder, at high photon energies, whereas at low energies the spectral flux follows the Rayleigh-Jeans law with a (surface) temperature lower than the effective temperature. The spectra of iron atmospheres display prominent spectral features in the soft X-ray range. The emergent specific intensity is anisotropic, with the anisotropy depending on energy. These properties of the atmospheric radiation should be taken into account for the proper interpretation of the thermal component of the neutron star radiation detectable at X-ray through UV energies. In particular, fitting of the thermal component with the blackbody model may result in substantially incorrect parameters of the radiating regions. We discuss an application of our atmosphere models to the nearby millisecond pulsar PSR J0437-4715.

We investigate the properties of the normal modes in a strongly magnetized hydrogen gas for conditions expected in atmospheres of isolated neutron stars. We derive and compute the polarizations and absorption coefficients of the normal modes from the polarizabily tensor using both analytical approximations and numerical calculations. We find that the spectral features and anisotropy associated with the bound-bound and bound-free transitions in the magnetized hydrogen are manifested in the polarization characteristics, which affects substantially the spectral and angular dependences of the absorption coefficients of the normal modes and the transfer of radiation in neutron star atmospheres.

We present first results of modeling convective photospheres of neutron stars. We show that in photospheres composed of the light elements convection arises only at relatively low effective temperatures (< 3-5*10^4 K), whereas in the case of iron compositon it arises at < 3*10^ K. Convection changes the depth dependence of the photosphere temperature and the shapes of the emergent spectra. Thus, it should be taken into account for the proper interpretation of EUV/soft-X-ray observations of the thermal radiation from neutron stars.