Previous XMM-Newton observations of the thermally emitting isolated neutron star RX J1605.3+3249 provided a candidate for a shallow periodic signal and evidence of a fast spin down, which suggested a high dipolar magnetic field and an evolution from a magnetar. We obtained a large programme with XMM-Newton to confirm its candidate timing solution, understand the energy-dependent amplitude of the modulation, and investigate the spectral features of the source. We performed extensive high-resolution and broadband periodicity searches in the new observations, using the combined photons of the three EPIC cameras and allowing for moderate changes of pulsed fraction and the optimal energy range for detection. A deep $4\sigma$ upper limit of $1.33(6)\%$ for modulations in the relevant frequency range conservatively rules out the candidate period previously reported. Blind searches revealed no other periodic signal above the $1.5\%$ level $(3\sigma$) in any of the four new observations. While theoretical models fall short at physically describing the complex energy distribution of the source, best-fit X-ray spectral parameters are obtained for a fully or partially ionized neutron star hydrogen atmosphere model with $B=10^{13}$ G, modified by a broad Gaussian absorption line at $385\pm10$ eV. The deep limits from the timing analysis disfavour equally well-fit double temperature blackbody models where both the star surface and small hotspots contribute to the X-ray flux of the source. We identified a low significance ($1\sigma$) temporal trend on the parameters of the source in the analysis of RGS data dating back to 2002, which may be explained by unaccounted calibration issues and spectral model uncertainties. The new dataset also shows no evidence of the previously reported narrow absorption feature at $\sim570$ eV, whose possible transient nature disfavours an atmospheric origin.

The Astronomical Roentgen Telescope X-ray Concentrator (ART-XC) is a hard X-ray telescope with energy response up to 30 keV, to be launched on board the Spectrum Roentgen Gamma (SRG) spacecraft in 2018. ART-XC consists of seven identical co-aligned mirror modules. Each mirror assembly is coupled with a CdTe double-sided strip (DSS) focal-plane detector. Eight X-ray mirror modules (seven flight and one spare units) for ART-XC were developed and fabricated at the Marshall Space Flight Center (MSFC), NASA, USA. We present results of testing procedures performed with an X-ray beam facility at MSFC to calibrate the point spread function (PSF) of the mirror modules. The shape of the PSF was measured with a high-resolution CCD camera installed in the focal plane with defocusing of 7 mm, as required by the ART-XC design. For each module, we performed a parametrization of the PSF at various angular distances Theta. We used a King function to approximate the radial profile of the near on-axis PSF (Theta<9 arcmin) and an ellipse fitting procedure to describe the morphology of the far off-axis angular response (9<Theta<24 arcmin). We found a good agreement between the seven ART-XC flight mirror modules at the level of 10%. The on-axis angular resolution of the ART-XC optics varies between 27 and 33 arcsec (half-power diameter), except for the spare module.

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.

The group of thermally emitting isolated neutron stars (INSs) known as the "Magnificent Seven" (M7) is unique among the various neutron star populations. Crustal heating by means of magnetic field decay and an evolutionary link with magnetars may explain why these objects rotate more slowly and have higher thermal luminosities and magnetic field intensities than standard pulsars of similar age. The third brightest INS, RX J1605.3+3249, is the only object amidst the seven still lacking a detected periodicity. We observed the source with the XMM-Newton Observatory for 60 ks aiming at unveiling the neutron star rotation rate and investigating its spectrum in detail. A periodic signal at P=3.387864(16) s, most likely the neutron star spin period, is detected at the 4-sigma confidence level. The coherent combination of the new data with a past XMM-Newton EPIC-pn observation of the source constrains the pulsar spin-down rate at the 2-sigma confidence level, implying a dipolar magnetic field of B~7.4e13 G. If confirmed, RX J1605.3+3249 would be the neutron star with the highest dipolar field amongst the M7. The spectrum of the source shows evidence of a cool blackbody component, as well as for the presence of two broad absorption features. Furthermore, high-resolution spectroscopy with the RGS cameras confirms the presence of a narrow absorption feature at energy 0.57 keV in the co-added spectrum of the source, also seen in other thermally emitting isolated neutron stars. Phase-resolved spectroscopy, as well as a dedicated observing campaign aimed at determining a timing solution, will give invaluable constraints on the neutron star geometry and will allow one to confirm the high value of spin down, which would place the source closer to a magnetar than any other M7 INS.

We report on the discovery of gamma-ray pulsations from five millisecond pulsars (MSPs) using the Fermi Large Area Telescope (LAT) and timing ephemerides provided by various radio observatories. We also present confirmation of the gamma-ray pulsations from a sixth source, PSR J2051-0827. Five of these six MSPs are in binary systems: PSRs J1713+0747, J1741+1351, J1600-3053 and the two black widow binary pulsars PSRs J0610-2100 and 2051-0827. The only isolated MSP is the nearby PSR J1024-0719, which is also known to emit X-rays. We present X-ray observations in the direction of PSRs J1600-3053 and J2051-0827. While the latter is firmly detected, we an only give upper limits for the X-ray flux of the former. There are no dedicated X-ray observations available for the other 3 objects. The MSPs mentioned above, together with most of the MSPs detected by Fermi, are used to put together a sample of 30 gamma-ray MSPs which is used to study the morphology and phase connection of radio and gamma-ray pulse profiles. We show that MSPs with pulsed gamma-ray emission which is phase aligned with the radio emission present the steepest radio spectra and the largest magnetic fields at the light cylinder among all MSPs. As well, we also observe a trend towards very low, or undetectable, radio linear polarisation levels. These properties could be attributed to caustic radio emission produced at a range of different altitudes in the magnetosphere. We note that most of these characteristics are also observed in the Crab pulsar, the only other radio pulsar known to exhibit phase-aligned radio and gamma-ray emission.

Since the launch of the Fermi satellite, the number of known gamma-ray pulsars has increased tenfold. Most gamma-ray detected pulsars are young and energetic, and many are associated with TeV sources. PSR J1357-6429 is a high spin-down power pulsar (Edot = 3.1 * 10^36 erg/s), discovered during the Parkes multibeam survey of the Galactic plane, with significant timing noise typical of very young pulsars. In the very-high-energy domain, H.E.S.S. has reported the detection of the extended source HESS J1356-645 (intrinsic Gaussian width of 12') whose centroid lies 7' from PSR J1357-6429. Using a rotational ephemeris obtained with 74 observations made with the Parkes telescope at 1.4 GHz, we phase-fold more than two years of gamma-ray data acquired by Fermi-LAT as well as those collected with XMM-Newton, and perform gamma-ray spectral modeling. Significant gamma and X-ray pulsations are detected from PSR J1357-6429. The light curve in both bands shows one broad peak. Gamma-ray spectral analysis of the pulsed emission suggests that it is well described by a simple power-law of index 1.5 +/- 0.3stat +/- 0.3syst with an exponential cut-off at 0.8 +/- 0.3stat +/- 0.3syst GeV and an integral photon flux above 100 MeV of (6.5 +/- 1.6stat +/- 2.3syst) * 10^-8 cm^-2 s^-1. The X-ray spectra obtained from the new data provide results consistent with those reported by Zavlin (2007). Upper limits on the gamma-ray emission from its potential pulsar wind nebula (PWN) are also reported. Assuming a distance of 2.4 kpc, the Fermi LAT energy flux yields a gamma-ray luminosity for PSR J1357-6429 of L_gamma = (2.13 +/- 0.25stat +/- 0.83syst) * 10^34 erg/s, consistent with an L_gamma \propto sqrt(Edot) relationship. The Fermi non-detection of the pulsar wind nebula associated with HESS J1356-645 provides new constraints on the electron population responsible for the extended TeV emission.

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.

We present a new study of the X-ray spectral properties of the Crab Pulsar. The superb angular resolution of the Chandra X-ray Observatory enables distinguishing the pulsar from the surrounding nebulosity. Analysis of the spectrum as a function of pulse phase allows the least-biased measure of interstellar X-ray extinction due primarily to photoelectric absorption and secondarily to scattering by dust grains in the direction of the Crab Nebula. We modify previous findings that the line-of-sight to the Crab is under-abundant in oxygen and provide measurements with improved accuracy and less bias. Using the abundances and cross sections from Wilms, Allen & McCray (2000) we find [O/H] = $(5.28 \pm 0.28)\times10^{-4}$ ($4.9 \times10^{-4}$ is solar abundance). We also measure for the first time the impact of scattering of flux out of the image by interstellar grains. We find $\tau_{\rm scat} = 0.147 \pm 0.043$. Analysis of the spectrum as a function of pulse phase also measures the X-ray spectral index even at pulse minimum --- albeit with increasing statistical uncertainty. The spectral variations are, by and large, consistent with a sinusoidal variation. The only significant variation from the sinusoid occurs over the same phase range as some rather abrupt behavior in the optical polarization magnitude and position angle. We compare these spectral variations to those observed in Gamma-rays and conclude that our measurements are both a challenge and a guide to future modeling and will thus eventually help us understand pair cascade processes in pulsar magnetospheres. The data were also used to set new, and less biased, upper limits to the surface temperature of the neutron star for different models of the neutron star atmosphere.

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.

Motivated by a paper (Kirsch et al. 2005) on possible use of the Crab Nebula as a standard candle for calibrating X-ray response functions, we examine consequences of intrinsic departures from a single (absorbed) power law upon such calibrations. We limit our analyses to three more modern X-ray instruments-the ROSAT/PSPC, the RXTE/PCA, and the XMM-Newton/EPIC-pn (burst mode). The results indicate a need to refine two of the three response functions studied. We are also able to distinguish between two current theoretical models for the system spectrum.

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

The remarkable astrometric capabilities of Chandra offer the possibility to measure proper motions of X-ray sources with an unprecedented accuracy in this wavelength range. We recently completed a proper motion survey of three of the seven thermally emitting radio-quiet isolated neutron stars (INSs) discovered in the ROSAT all-sky survey. These INSs (RX J0420.0-5022, RX J0806.4-4123, and RX J1308.6+2127) either lack an optical counterpart or have one so faint that ground based or space born optical observations push the current possibilities of the instrumentation to the limit. Pairs of ACIS observations were acquired 3 to 5 years apart to measure the displacement of the sources on the X-ray sky using as reference the background of extragalactic or remote Galactic X-ray sources. We derive 2 sigma upper limits of 123 mas/yr and 86 mas/yr on the proper motion of RX J0420.0-5022 and RX J0806.4-4123, respectively. RX J1308.6+2127 exhibits a very significant displacement (~ 9 sigma) yielding mu = 220 +/- 25 mas/yr, the second fastest measured among all ROSAT discovered INSs. The source is probably moving away rapidly from the Galactic plane at a speed which precludes any significant accretion of matter from the interstellar medium. Its transverse velocity of ~ 740 (d/700pc) km/s might be the largest of all ROSAT INSs and its corresponding spatial velocity stands among the fastest recorded for neutron stars. RX J1308.6+2127 is thus a middle-aged (age ~ 1 My) high velocity cooling neutron star. We investigate its possible origin in nearby OB associations or from a field OB star. In most cases, the flight time from birth place appears significantly shorter than the characteristic age derived from spin down rate. The distribution in transverse velocity of ROSAT INSs is not statistically different from that of normal radio pulsars.

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.

The unprecedented spatial resolution of the Chandra observatory opens the possibility to detect with relatively high accuracy proper motions at X-ray wavelengths. We have conducted an astrometric study of three of the "Magnificent Seven", the thermally emitting radio quiet isolated neutron stars (INSs) discovered by ROSAT. These three INSs (RX J0420.0-5022, RX J0806.4-4123 and RX J1308.6+2127/RBS 1223) either lack an optical counterpart or have one too faint to be used for astrometric purposes. We obtained ACIS observations 3 to 5 years apart to constrain or measure the displacement of the sources on the X-ray sky using as reference the background of extragalactic or remote galactic X-ray sources. Upper limits of 138 mas/yr and 76 mas/yr on the proper motion of RX J0420.0-5022 and RX J0806.4-4123, respectively, have already been presented in Motch et al. (2007). Here we report the very significant measurement (~ 10 sigma) of the proper motion of the third INS of our program, RX J1308.6+2127/RBS1223. Comparing observations obtained in 2002 and 2007 reveals a displacement of 1.1 arcsec implying a yearly proper motion of 223 mas, the second fastest measured for the ROSAT discovered INSs. The source is rapidly moving away from the galactic plane at a speed which precludes any significant accretion of matter from the interstellar medium. Its transverse velocity of ~ 740 (d/700pc) km/s might be the largest of the "Magnificent Seven" and among the fastest recorded for neutron stars. RX J1308.6+2127/RBS1223 is thus a young high velocity cooling neutron star. The source may have its origin in the closest part of the Scutum OB2 association about 0.8 Myr ago, an age consistent with that expected from cooling curves, but significantly younger than inferred from pulse timing measurements (1.5 Myr).

The first short Chandra and XMM-Newton observations of the young and energetic pulsar J1357-6429 provided strong indications of a tail-like pulsar-wind nebula associated with this object, as well as pulsations of its X-ray flux with a pulsed fraction above 50% and a thermal component dominating at lower photon energies (below 2 keV). The elongated nebula is very compact in size and may be interpreted as evidence for a pulsar jet. The thermal radiation is most plausibly emitted from the entire neutron star surface of a 10 km radius and a 1.0+/-0.1 K temperature, covered with a magnetized hydrogen atmosphere. At higher energies the pulsar's emission is of a nonthermal (magnetospheric) origin, with a power-law spectrum of a photon index of 1.1-1.3. This makes the X-ray properties of PSR J1357-6429 very similar to those of the youngest pulsars J1119-6127 and Vela with a detected thermal radiation.

The possibility for direct investigation of thermal emission from isolated neutron stars was opened about 25 years ago with the launch of the first X-ray observatory, Einstein. A significant contribution to this study was provided by ROSAT in 1990's. The outstanding capabilities of the currently operating observatories, Chandra and XMM-Newton, have greatly increased the potential to observe and analyze thermal radiation from the neutron star surfaces. Confronting observational data with theoretical models of thermal emission, presumably formed in neutron star atmospheres, allows one to infer the surface temperatures, magnetic fields, chemical composition, and neutron star masses and radii. This information, supplemented with model equations of state and neutron star cooling models, provides an opportunity to understand the fundamental properties of the superdense matter in the neutron star interiors. I review the current status and most important results obtained from modeling neutron star thermal emission and present selected Chandra and XMM-Newton results on thermal radiation from various types of these objects: ordinary radio pulsars with ages ranging from about 2 kyr to 20 Myr (J1119-6127, Vela, B1706-44, J0538+2817, B2334+61, B0656+14, B1055-52, Geminga, B0950+08, J2043+2740), millisecond pulsars (J0030+0451, J2124-3358, J1024-0719, J0437-4715), putative pulsars (CXOU J061705.3+222127, RX J0007.0+7302), central compact objects in supernova remnats (in particular, 1E 1207.4-5209), and isolated radio-quiet neutron stars.

We report observations with the Chandra X-ray Observatory of the single, cool, magnetic white dwarf GD 356. For consistent comparison with other X-ray observations of single white dwarfs, we also re-analyzed archival ROSAT data for GD 356 (GJ 1205), G 99-47 (GR 290 = V1201 Ori), GD 90, G 195-19 (EG250 = GJ 339.1), and WD 2316+123 and archival Chandra data for LHS 1038 (GJ 1004) and GD 358 (V777 Her). Our Chandra observation detected no X rays from GD 356, setting the most restrictive upper limit to the X-ray luminosity from any cool white dwarf -- L_{X} < 6.0 x 10^{25} ergs/s, at 99.7% confidence, for a 1-keV thermal-bremsstrahlung spectrum. The corresponding limit to the electron density is n_{0} < 4.4 x 10^{11} cm^{-3}. Our re-analysis of the archival data confirmed the non-detections reported by the original investigators. We discuss the implications of our and prior observations on models for coronal emission from white dwarfs. For magnetic white dwarfs, we emphasize the more stringent constraints imposed by cyclotron radiation. In addition, we describe (in an appendix) a statistical methodology for detecting a source and for constraining the strength of a source, which applies even when the number of source or background events is small.

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.

Millisecond pulsars represent an evolutionarily distinct group among rotation-powered pulsars. Outside the radio band, the soft X-ray range ($\sim 0.1$--10 keV) is most suitable for studying radiative mechanisms operating in these fascinating objects. X-ray observations revealed diverse properties of emission from millisecond pulsars. For the most of them, the bulk of radiation is of a thermal origin, emitted from small spots (polar caps) on the neutron star surface heated by relativistic particles produced in pulsar acceleration zones. On the other hand, a few other very fast rotating pulsars exhibit almost pure nonthermal emission generated, most probably, in pulsar magnetospheres. There are also examples of nonthermal emission detected from X-ray nebulae powered by millisecond pulsars, as well as from pulsar winds shocked in binary systems with millisecond pulsars as companions. These and other most important results obtained from X-ray observations of millisecond pulsars are reviewed in this paper, as well as results from the search for millisecond pulsations in X-ray flux of the radio-quite neutron star RX J1856.5-3754.

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.

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.

I present an analysis of the XMM-Newton observations of four millisecond pulsars, J0437-4715, J2124-3358, J1024-0719, and J0034-0534. The new data provide strong evidence of thermal emission in the X-ray flux detected from the first three objects. This thermal component is best interpreted as radiation from pulsar polar caps covered with a nonmagnetic hydrogen atmosphere. A nonthermal power-law component, dominating at energies E>3 keV, can also be present in the detected X-ray emission. For PSR J0437-4715, the timing analysis reveals that the shape and pulsed fraction of the pulsar light curves are energy dependent. This, together with the results obtained from the phase-resolved spectroscopy, supports the two-component (thermal plus nonthermal) interpretation of the pulsar's X-ray radiation. Highly significant pulsations have been found in the X-ray flux of PSRs J2124-3358 and J1024-0719. For PSR J0034-0534, a possible X-ray counterpart of the radio pulsar has been suggested. The inferred properties of the detected thermal emission are compared with predictions of radio pulsar models.

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 obtained deep optical imaging of the thermally emitting X-ray bright and radio-quiet isolated neutron star RX J1605.3+3249 with the Subaru telescope in 1999 and 2003. Together with archival HST images acquired in 2001 these data reveal a proper motion of mu = 144.5 +/- 13.2 mas/yr. This implies a relatively high spatial velocity and indicates that the star is unlikely to be re-heated by accretion of matter from the interstellar medium. Assuming that RX J1605.3+3249 is a young (10^5-10^6 yr) cooling neutron star, its apparent trajectory is consistent with a birth in the nearby Sco OB2 OB association at a location close to that derived for RX J1856.5-3754 and perhaps also to that of RX J0720.4-3125. This suggests that the X-ray bright part of ROSAT-discovered isolated neutron stars is dominated by the production of the Sco OB2 complex which is the closest OB association and a part of the Gould belt. The B and R magnitudes of the faint optical counterpart did not vary from 1999 to 2003 at B = 27.22 +/- 0.10. Its B-R colour index of +0.32 +/- 0.17 is significantly redder than that of other isolated neutron stars and the optical flux lies a factor 11.5 above the extrapolation of the X-ray blackbody-like spectrum. The red optical colour reveals the presence of an additional emitting component in the optical regime over the main neutron star thermal emission. We also discovered a small elongated Halpha nebula approximately centered on the neutron star and aligned with the direction of motion. The width of the nebula is unresolved and smaller than ~ 0.4" for a length of about 1". The shape of the Balmer emitting nebula is very different from those seen close to other neutron stars and should be confirmed by follow-up observations. We shortly discuss the possible mechanisms which could give rise to such a geometry.

High quantum efficiency over a broad spectral range is one of the main properties of the EPIC pn camera on-board XMM-Newton. The quantum efficiency rises from ~75% at 0.2 keV to ~100% at 1 keV, stays close to 100% until 8 keV, and is still ~90% at 10 keV. The EPIC pn camera is attached to an X-ray telescope which has the highest collecting area currently available, in particular at low energies (more than 1400 cm2 between 0.1 and 2.0 keV). Thus, this instrument is very sensitive to the low-energy X-ray emission. However, X-ray data at energies below ~0.2 keV are considerably affected by detector effects, which become more and more important towards the lowest transmitted energies. In addition to that, pixels which have received incorrect offsets during the calculation of the offset map at the beginning of each observation, show up as bright patches in low-energy images. Here we describe a method which is not only capable of suppressing the contaminations found at low energies, but which also improves the data quality throughout the whole EPIC pn spectral range. This method is then applied to data from the Vela supernova remnant.

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 report on the analysis of new X-ray data obtained with XMM-Newton and Chandra from two ROSAT-discovered X-ray dim isolated neutron stars. RX J0806.4-4123 was observed with XMM-Newton in April 2003, 2.5 years after the first observation. The EPIC-pn data confirm that this object is an X-ray pulsar with 11.371 s neutron star spin period. The X-ray spectrum is consistent with absorbed black-body emission with a temperature kT = 96 eV and NH = 4x10^19 cm^-2 without significant changes between the two observations. Four XMM-Newton observations of RX J0420.0-5022 between December 2002 and July 2003 did not confirm the 22.7 s pulsations originally indicated in ROSAT data, but clearly reveal a 3.453 s period. A fit to the X-ray spectrum using an absorbed black-body model yields kT = 45 eV, the lowest value found from the small group of XDINs and NH = 1.0x10^20 cm^-2. Including a broad absorption line improves the quality of the spectral fits considerably for both objects and may indicate the presence of absorption features similar to those reported from RBS1223, RX J1605.3+3249 and RX J0720.4-3125. For both targets we derive accurate X-ray positions from the Chandra data and present an optical counterpart candidate for RX J0420.0-5022 with B = 26.6+/-0.3 mag from VLT imaging.

The XMM-Newton spectra of the isolated neutron star RX J0720.4-3125 show deviations from a Planckian energy distribution below 400 eV, similar to the spectra of RBS1223, another long-period X-ray pulsar, as reported recently by Haberl et al. (2003). For a Gaussian-shaped absorption line we derive an energy of 271 eV and an equivalent width of -40 eV from the phase-averaged spectra of RX J0720.4-3125. We investigate the spectral variations seen in hardness ratios as function of pulse phase and find that they are best described by changes in the depth of the absorption line. The line equivalent width changes between -31 eV around intensity maximum of the pulse and -58 eV at the declining part of the pulse. Small variations (<20 eV) of the line energy with pulse phase may still be caused by statistical fluctuations. On the other hand, the black-body temperature varies significantly by 2.5 eV (statistical 90% errors typically 0.7 eV) reaching the highest value at pulse maximum. One possible interpretation for the absorption line is cyclotron resonance scattering of protons in a magnetic field with B about 5x10^13 G. This field strength is compatible with estimates inferred from recent spin down measurements of the pulsar.

We present a reanalysis of the X-ray data for RX J0720.4-3125 presented in Zane et al. (2002), using more data recently available from XMM-Newton and Chandra. This analysis also corrects the Rosat data used in that paper to the TDB time system, incorporates the revised XMM-Newton barycentric correction available since then, and corrects the definition of the instantaneous period in the maximum likelihood periodogramme search. However, we are now unable to find a single coherent period that is consistent with all Rosat, Chandra and XMM-Newton datasets. From an analysis of the separate datasets, we have derived limits on the period change of Pdot = 1.4+/-0.6x10^-13 s/s at 99% confidence level. This is stronger than the value presented in Zane et al. (2002), but sufficiently similar that their scientific conclusions remain unchanged. We examine the implications in more detail, and find that RX J0720.4-3125 can have been born as a magnetar provided that it has a young age of ~10^4 yr. A more conservative interpretation is that the field strength has remained relatively unchanged at just over 10^13 G, over the ~10^6 yr lifetime of the star.

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.

We discuss recent Chandra and XMM-Newton observations of the bright isolated neutron star RX J1856.5-3754 and suggest that the absence of any line features is due to effects of a high magnetic field strength (~10^13 G). Using different models for the temperature distribution across the neutron star surface assuming blackbody emission to fit the optical and X-ray spectrum and we derive a conservative lower limit of the "apparent" neutron star radius of 16.5 km x (d/117 pc). This corresponds to the radius for the "true" (de-redshifted) radius of 14 km for a 1.4 Msun neutron star, indicating a stiff equation of state at high densities. A comparison of the result with mass-radius diagrams shows that quark stars and neutron stars with quark matter cores can be ruled out with high confidence.

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.

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

ESO 4m class telescope and VLT deep imaging of the isolated neutron star RX J0720.4-3125 reveals a proper motion of mu = 97 +/-12 mas/yr and a blue U-B color index. We show that a neutron star atmosphere model modified to account for a limited amount of hydrogen on the star's surface can well represent both the optical and X-ray data without invoking any additional components. The large proper motion almost completely excludes the possibility that accretion from interstellar medium is the powering mechanism of the X-ray emission. It also implies that the proposed spin down is entirely due to magnetic dipole losses. RX J0720.4-3125 is thus a very likely middle aged cooling neutron star. Its overall properties are quite similar to some of the long period radio pulsars recently discovered, giving further support to the idea that RX J0720.4-3125 may be a pulsar whose narrow radio beam does not cross the Earth.

We present results of the analysis of data collected in 57-ks XMM-Newton and 505-ks Chandra observations of the nearby ~= 120 pc isolated neutron star RX J1856.5-3754. We confirm most of the statements made by Burwitz et al. (2001) who discussed the original 55-ks Chandra data. Detailed spectral analysis of the combined X-ray and optical data rules out the currently available nonmagnetic light and heavy element neutron star atmosphere (LTE) models with hydrogen, helium, iron and solar compositions. We find that strongly magnetized atmosphere models also are unable to represent the data. The X-ray and optical data show no spectral features and are best fitted with a two-component blackbody model with kT_{bb,X}^{\infty} ~= 63.5 eV and R_{bb,X}^{\infty} ~= 4.4(d/120 pc) km for the hot X-ray emitting region, and kT_{bb,opt}^{\infty}<33 eV and R_{bb,opt}^{\infty}>17(d/120 pc) km for the rest of the neutron star surface responsible for the optical flux. The large number of counts collected with XMM-Newton allows us to reduce the upper limit on periodic variation in the X-ray range down to 1.3% (at a 2\sigma confidence level) in the 10^{-3}-50 Hz frequency range. In an attempt to explain this small variability, we discuss an one-component model with kT_{bb}^{\infty} ~= 63 eV and R_{bb}^{\infty} ~= 12.3(d/120 pc) km. This model requires a low radiative efficiency in the X-ray domain, which may be expected if the neutron star has a condensed matter surface.

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.

The isolated neutron star RX J0806.4-4123 was observed with XMM-Newton in November 2000. The data from the three EPIC instruments allowed us (i) to derive an improved X-ray position to an accuracy of 2-3'', (ii) to accumulate the first medium-resolution soft X-ray spectra of high statistical quality and (iii) to find a candidate for the neutron star rotation period. Although this period of 11.3714 s is formally detected at a 3.5-\sigma level in the EPIC-pn data, the similar pulse profiles deduced from all three EPIC instruments increase the confidence that the period is real. The pulsed fraction of ~6% would then be the weakest X-ray flux modulation detected from dim isolated neutron stars. We fitted the X-ray spectra with blackbody and neutron star atmosphere models and discuss the results with respect to the brightness limit placed by optical images. The reduced size of the error circle on the X-ray position should allow deeper searches for an optical counterpart.

The pn camera of EPIC is inherently robust against radiation damage effects. Nevertheless, the knowledge of charge transfer inefficiency (CTI) of the pn camera is crucial for obtaining the correct energy scale. We describe detailed in-orbit monitoring of this effect, utilizing the internal calibration source. We find that during the first two years in orbit the CTI increased by about 4% for Al-K and by about 7% for Mn-K. The increases in the CTI are well within expectations, with no measurable effect on the energy resolution.

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.

One of the three cameras behind the X-ray telescopes on board XMM-Newton is equipped with a pn CCD detector. Here the current status of the pn spectral response calibration is summarized. Several parameters describing the spectral response of the EPIC-pn detector show spatial dependencies. E.g. spectral resolution and single- to double-pixel event ratios are functions of the charge transfer losses and therefore depend on the distance to the read-out node of the CCD (the RAWY detector coordinate).

We present a combined analysis of XMM-Newton, Chandra and Rosat observations of the isolated neutron star RX J0720.4-3125, spanning a total period of \sim 7 years. We develop a maximum likelihood periodogramme based on \Delta C statistic and maximum likelihood method, which are appropriate for sparse event lists. As an "a posteriori" check, we have folded a further BeppoSAX dataset with the period predicted at the time of that observation, finding that the phase is consistent. The value of the spin down rate, here measured for the first time, is \approx 10^{-14} s/s and can not be explained in terms of propeller or torque from a fossil disk. When interpreted in terms of dipolar losses, it gives a magnetic field of B \approx 10^{13} G, making also implausible that the source is accreting from the underdense surroundings. We discuss the implications of this measure for the different mechanisms that have been suggested to explain the X-ray emission. We conclude that the observed properties are more compatible with a scenario in which the source is \approx 10^6 yrs old, and its magnetic field has not changed substantially over the lifetime.

We present a combined analysis of XMM-Newton, Chandra and Rosat observations of the isolated neutron star RXJ0720.4-3125, spanning a total period of \sim 7 years. We develop a maximum likelihood periodogramme for our analysis based on the \Delta C-statistic and the maximum likelihood method, which are appropriate for the treatment of sparse event lists. Our results have been checked "a posteriori" by folding a further BeppoSAX dataset with the period predicted at the time of that observation: the phase is found to be consistent. The study of the spin history and the measure of the spin-down rate is of extreme importance in discriminating between the possible mechanisms suggested for the nature of the X-ray emission. The value of \dot P, here measured for the first time, is \approx 10^{-14} s/s. This value can not be explained in terms of torque from a fossil disk. When interpreted in terms of dipolar losses, it gives a magnetic field of B \approx 10^{13} G, making also implausible that the source is accreting from the underdense surroundings. On the other hand, we also find unlikely that the field decayed from a much larger value (B\approx 10^{15} G, as expected for a magnetar powered by dissipation of a superstrong field) since this scenario predicts a source age of \approx 10^4 yrs, too young to match the observed X-ray luminosity. The observed properties are more compatible with a scenario in which the source is \approx 10^6 yrs old, and its magnetic field has not changed substantially over the lifetime.

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 the Chandra LETGS X-ray spectrum of the nearby (~60 pc) neutron star RX J1856.5-3754. Detailed spectral analysis of the combined X-ray and optical data rules out the nonmagnetic neutron star atmosphere models with hydrogen, helium, iron and solar compositions. We also conclude that strongly magnetized atmosphere models are unable to represent the data. The data can be explained with a two-component blackbody model. The harder component with temperature of kT_bb~63 eV and a radius R_bb~2.2 km of the emitting region well fits the X-ray data and can be interpreted as radiation from a hot region on the star's surface.

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

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.

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.

We present the XMM-Newton first light image, taken in January 2000 with the EPIC pn camera during the instrument's commissioning phase, when XMM-Newton was pointing towards the Large Magellanic Cloud (LMC). The field is rich in different kinds of X-ray sources: point sources, supernova remnants (SNRs) and diffuse X-ray emission from LMC interstellar gas. The observations are of unprecedented sensitivity, reaching a few 10^32 erg/s for point sources in the LMC. We describe how these data sets were analysed and discuss some of the spectroscopic results. For the SNR N157B the power law spectrum is clearly steeper than previously determined from ROSAT and ASCA data. The existence of a significant thermal component is evident and suggests that N157B is not a Crab-like but a composite SNR. Most puzzling is the spectrum of the LMC hot interstellar medium, which indicates a significant overabundance of Ne and Mg of a few times solar.

We present the high resolution spectrum of the isolated neutron star RXJ0720.4-3125, obtained with the Reflection Grating Spectrometer on XMM-Newton, complemented with the broad band spectrum observed with the EPIC PN camera. The spectrum appears smooth, with no evidence for strong photospheric absorption or emission features. We briefly discuss the implications of our failure to detect structure in the spectrum.

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

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.

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 use models of magnetized hydrogen atmospheres to fit the soft X-ray thermal spectrum of the Vela pulsar obtained by ROSAT. The distance and hydrogen column density required by our fits are in good agreement with other independents estimates and can be considered as an argument for the presence of hydrogen at the surface of this pulsar. The low temperature obtained, T = 7.85 +/- 0.25 10^5 K, is below the predictions of the `standard' model of neutron star cooling and strengthens previous claims for the necessity of the occurrence of fast neutrino cooling in this neutron star.

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.

We present a model to explain the decrease in the amplitude of the pulse profile with increasing energy observed in Geminga's soft X-ray surface thermal emission. We assume the presence of plates surrounded by a surface with very distinct physical properties: these two regions emit spectra of very distinct shapes which present a crossover, the warm plates emitting a softer spectrum than the colder surrounding surface. The strongly pulsed emission from the plates dominates at low energy while the surroundings emission dominates at high energy, producing naturally a strong decrease in the pulsed fraction. In our illustrative example the plates are assumed to be magnetized while the rest of the surface is field free. This plate structure may be seen as a schematic representation of a continuous but very nonuniform distribution of the surface magnetic field or as a quasi realistic structure induced by past tectonic activity on Geminga.