Ages and thermal luminosities of neutron stars, inferred from observations, can be interpreted with the aid of the neutron star cooling theory to gain information on the properties of superdense matter in neutron-star interiors. We present a survey of estimated ages, surface temperatures and thermal luminosities of middle-aged neutron stars with relatively weak or moderately strong magnetic fields, which can be useful for these purposes. The catalogue includes results selected from the literature, supplemented with new results of spectral analysis of a few cooling neutron stars. The data are compared with the theory. We show that overall agreement of theoretical cooling curves with observations improves substantially for models where neutron superfluidity in stellar core is weak.

We present observations of the pulsar-wind nebula (PWN) region ofSNR 0540-69.3. The observations were made with the Atacama Compact Array (ACA) in Bands 4 and 6. We also add radio observations from the Australia Compact Array (ATCA) at 3 cm. For 1.449 - 233.50 GHz we obtain a synchrotron spectrum $F_{\nu} \propto \nu^{-\alpha_{\nu}}$, with the spectral index $\alpha_{\nu} = 0.17\pm{0.02}$. To conclude how this joins the synchrotron spectrum at higher frequencies we include hitherto unpublished AKARI mid-infrared data, and evaluate published data in the ultraviolet (UV), optical and infrared (IR). In particular, some broad-band filter data in the optical must be discarded from our analysis due to contamination by spectral line emission. For the UV/IR part of the synchrotron spectrum, we arrive at $\alpha_{\nu} = 0.87^{+0.08}_{-0.10}$. There is room for $2.5\times10^{-3}$ solar masses of dust with temperature $\sim 55$ K if there are dual breaks in the synchrotron spectrum, one around $\sim 9\times10^{10}$ Hz, and another at $\sim 2\times10^{13}$ Hz. The spectral index then changes at $\sim 9\times10^{10}$ Hz from $\alpha_{\nu} = 0.14\pm0.07$ in the radio, to $\alpha_{\nu} = 0.35^{-0.07}_{+0.05}$ in the millimetre to far-IR range. The ACA Band 6 data marginally resolves the PWN. In particular, the strong emission 1.5" south-west of the pulsar, seen at other wavelengths, and resolved in the 3-cm data with its 0.8" spatial resolution, is also strong in the millimeter range. The ACA data clearly reveal the supernova remnant shell 20-35 arcsec west of the pulsar, and for the shell we derive $\alpha_{\nu} = 0.64\pm{0.05}$ for the range $8.6-145$~GHz.

We analyse archival optical data on the binary companion to the millisecond pulsar \j1045\ obtained with the Very Large Telescope. A possible optical counterpart candidate %source position coincides with the pulsar coordinates is found at the pulsar position. Its magnitudes are $V>26.4$, $R=25.7(2)$ and $I=25.4(2)$. The data are compared with white dwarf evolutionary models. Depending on still poorly known distance to the pulsar, two alternative interpretations are possible. For the radio timing parallax distance, the candidate can be a very cool and old white dwarf with the hydrogen atmosphere and the temperature of $< 3000$ K and the age of $> 5$ Gyr. The dispersion measure distance suggests a younger ($\sim$ 2 Gyr) and hotter ($\sim$ 6000 K) white dwarf, whose atmosphere composition remains unknown.

We report results of XMM-Newton observations of a $\gamma$-ray pulsar J0633+0632 and its wind nebula. We reveal, for the first time, pulsations of the pulsar X-ray emission with a single sinusoidal pulse-profile and a pulsed fraction of $23 \, \pm \, 6$ per cent in the 0.3--2 keV band. We confirm previous Chandra findings that the pulsar X-ray spectrum consists of thermal and non-thermal components. However, we do not find the absorption feature that was previously detected at about 0.8 keV. Thanks to the greater sensitivity of XMM-Newton, we get stronger constraints on spectral model parameters compared to previous studies. The thermal component can be equally well described by either blackbody or neutron star atmosphere models, implying that this emission is coming from either hot pulsar polar caps with a temperature of about 120 eV or from the colder bulk of the neutron star surface with a temperature of about 50 eV. In the latter case, the pulsar appears to be one of the coolest among other neutron stars of similar ages with estimated surface temperatures. We discuss cooling scenarios relevant to this neutron star. Using an interstellar absorption--distance relation, we also constrain the distance to the pulsar to the range of 0.7--2 kpc. Besides the pulsar and its compact nebula, we detect regions of weak large-scale diffuse non-thermal emission in the pulsar field and discuss their possible nature.

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 measured the proper motion of a $\gamma$-ray radio-quiet pulsar J0633+0632 using Chandra observations performed in 2009 and 2017. The measured proper motion is $53\pm15$~mas~yr$^{-1}$. We found that the proper motion direction does not follow the extension of the J0633+0632 pulsar wind nebula. The J0633+0632 pulsar wind nebula therefore can be a jet-like feature or a misaligned outflow. We also discuss a possible birth cite of the pulsar.

We report on multi-band photometric observations of four binary millisecond pulsars with the Gran Telescopio Canarias. The observations led to detection of binary companions to PSRs J1630+3734, J1741+1351 and J2042+0246 in the Sloan g', r' and i' bands. Their magnitudes in the r' band are $\approx$24.4, 24.4 and 24.0, respectively. We also set a 3$\sigma$ upper limit on the brightness of the PSR J0557+1550 companion in the r' band of $\approx$25.6 mag. Combining the optical data with the radio timing measurements and white dwarf cooling models, we show that the detected companions are cool low-mass white dwarfs with temperatures and ages in the respective ranges of (4-7)$\times 10^3$ K and 2-5 Gyr. All the detected white dwarfs are found to likely have either pure hydrogen or mixed helium-hydrogen atmospheres.

We present optical time-resolved multi-band photometry of the black widow binary millisecond pulsar J2052+1219 using direct-imaging observations with the 2.1m telescope of Observatorio Astronomico Nacional San Pedro Martir, Mexico (OAN-SPM). The observations revealed a variable optical source whose position and periodicity P = 2.752h coincide with the pulsar coordinates and the orbital period obtained from radio timing. This allowed us to identify it with the binary companion of the pulsar. We reproduce light curves of the source modelling the companion heating by the pulsar and accounting for the system parameters obtained from the radio data. As a result, we independently estimate the distance to the system of 3.94(16) kpc, which agrees with the dispersion measure distance. The companion star size is 0.12-0.15 Rsun, close to filling its Roche lobe. It has a surface temperature difference of about 3000 K between the side facing the pulsar and the back side. We summarise characteristics of all black widow systems studied in the optical and compare them with the PSR J2052+1219 parameters derived from our observations.

We have analyzed archival XMM-Newton and Chandra observations of the gamma-ray radio-quiet pulsar J1826-1256 and its pulsar wind nebula. The pulsar spectrum can be described by a power-law model with a photon index $\Gamma\approx1$. We find that the nebular spectrum softens with increasing distance from the pulsar, implying synchrotron cooling. The empirical interstellar absorption-distance relation gives a distance of $\approx3.5$ kpc to J1826-1256. We also discuss the nebula geometry and association between the pulsar, the very high energy source HESS J1826-130, the supernova remnant candidate G18.45-0.42 and the open star cluster Bica 3.

The bright radio pulsar B1727$-$47 with a characteristic age of 80 kyr is among the first pulsars discovered 50 yr ago. Using its regular timing observations and interferometric positions at three epochs, we measured, for the first time, the pulsar proper motion of $151 \pm 19$ mas yr$^{-1}$. At the dispersion measure distance of $\gtrsim 2.7$ kpc, this would suggest a record transverse velocity of the pulsar $\gtrsim 1900$ km s$^{-1}$. However, a backward extrapolation of the pulsar track to its birth epoch points remarkably close to the center of the evolved nearby supernova remnant RCW 114, which suggests genuine association of the two objects. In this case, the pulsar is substantially closer ($\sim 0.6$ kpc) and younger ($\sim 50$ kyr), and its velocity ($\sim400$ km s$^{-1}$) is compatible with the observed pulsar velocity distribution. We also identified two new glitches of the pulsar. We discuss implications of our results on the pulsar and remnant properties.

We report detection of the binary companion to the millisecond pulsar J1741+1351 with the Gran Telescopio Canarias. The optical source position coincides with the pulsar coordinates and its magnitudes are g' = 24.84(5), r' = 24.38(4) and i' = 24.17(4). Comparison of the data with the white dwarf evolutionary models shows that the source can be a He-core white dwarf with a temperature of $\approx 6000$ K and a mass of $\approx 0.2$ M$_{\odot}$. The latter is in excellent agreement with the companion mass obtained from the radio timing solution for PSR J1741+1351.

We report the detection of the likely companion of the binary millisecond pulsar J0740+6620 with the Gran Telescopio Canarias in the $r'$ and $i'$ bands. The position of the detected starlike source coincides with the pulsar coordinates within the 1$\sigma$ uncertainty of $\approx$ 0.2 arcsec. Its magnitudes are $r'=26.51\pm0.17$ and $i'=25.49\pm0.15$. Comparing the data with the white dwarf cooling tracks suggests that it can be an ultracool helium-atmosphere white dwarf with the temperature $\lesssim$ 3500 K and cooling age $\gtrsim$ 5 Gyr. The age is consistent with the pulsar characteristic age corrected for kinematic effects. This is the reddest source among known white dwarf companions of millisecond pulsars. Detection of the source in other bands would be useful to clarify its properties and nature.

We report detection of the binary companion to the millisecond pulsar J2302+4442 based on the deep observations performed with the Gran Telescopio Canarias. The observations revealed an optical source with $r'$ = 23.33$\pm$0.02 and $i'$=23.08$\pm$0.02, whose position coincides with the pulsar radio position. By comparing the source colour and magnitudes with the white dwarf cooling predictions, we found that it likely represents a He or CO-core white dwarf and estimated its mass of 0.52$^{+0.25}_{-0.19}$ M$_{\odot}$ and effective temperature of 6300$^{+1000}_{-800}$ K. Combining our results with the radio timing measurements, we set constraints on the binary system inclination angle $i=73^{+6}_{-5}$ degrees and the pulsar mass $m_p=3.1^{+2.7}_{-2.0}$ M$_{\odot}$.

Using the SDSS and Pan-STARRS1 survey data, we found a likely companion of the recently discovered binary $\gamma$-ray radio-loud millisecond pulsar J0621+2514. Its visual brightness is about 22 mag. The broad band magnitudes and colours suggest that this is a white dwarf. Comparing the data with various white dwarfs evolutionary tracks, we found that it likely belongs to a class of He-core white dwarfs with a temperature of about 10 000 K and a mass of $\lesssim 0.5$M$_\odot$. For a thin hydrogen envelope of the white dwarf its cooling age is $\lesssim0.5$ Gyr which is smaller than the pulsar characteristic age of 1.8 Gyr. This may indicate that the pulsar age is overestimated. Otherwise, this may be explained by the presence of a thick hydrogen envelope or a low metallicity of the white dwarf progenitor.

Using archival X-ray data we have found point-like X-ray counterpart candidates positionally coincident with six $\gamma$-ray pulsars discovered recently in the Fermi Gamma-ray Space Telescope data by the Einstein@Home project. The candidates for PSRs J0002$+$6216, J0554$+$3107, J1844$-$0346 and J1105$-$6037 are detected with Swift, and those for PSRs J0359$+$5414 and J2017$+$3625 are detected with Chandra. Despite a low count statistics for some candidates, assuming plausible constraints on the absorbing column density towards the pulsars, we show that X-ray spectral properties for all of them are consistent with those observed for other pulsars. J0359$+$5414 is the most reliably identified object. We detect a nebula around it, whose spectrum and extent suggest that this is a pulsar wind nebula powered by the pulsar. Associations of J0002$+$6216 and J1844$-$0346 with supernova remnants CTB 1 and G28.6$-$0.1 are proposed.

We report preliminary results of the analysis of the proper motion of the bright radio pulsar B1727$-$47. Using archival Parkes timing data, as well as original and archival ATCA interferometry observations, we, for the first time, constrain the pulsar proper motion at the level of 148$\pm$11 mas yr$^{-1}$. The backward extrapolation of the proper motion vector to the pulsar birth epoch points at the center of the Galactic supernova remnant RCW 114 suggesting the genuine association between the two objects. We discuss the implications of the association and argue that the distance to the system is less than $1$ kpc. This value is at least two times lower than the dispersion measure distance estimates. This suggests that the existing Galaxy electron density models are incomplete in the direction to the pulsar.

We present the analysis of the archival Suzaku and Swift X-ray observations of the young $\gamma$-ray pulsar J1932+1916 field. The data revealed a point-like object at the $\gamma$-ray position of the pulsar and diffuse X-ray emission around it. Spectra of the point-like source and diffuse emission are well-described by absorbed power-law models with spectral parameters typical for pulsar plus pulsar wind nebula systems. Therefore we suggest that Suzaku and Swift detected the X-ray counterpart of PSR J1932+1916. Assuming this interpretation, we constrain the distance to the pulsar in the range of 2-6 kpc. We also suggest possible association of the pulsar with the nearby supernova remnant G54.4-0.3 and discuss its implications for the pulsar proper motion, age and distance.

Calvera is an unusual isolated neutron star with pure thermal X-ray spectrum typical for central compact objects in supernova remnants. On the other hand, its rotation period and spin-down rate are typical for ordinary rotation-powered pulsars. It was discovered and studied in X-rays and not yet detected in other spectral domains. We present deep optical imaging of the Calvera field obtained with the Gran Telescopio Canarias in $g'$ and $i'$ bands. Within $\approx 1^{\prime\prime}$ vicinity of Calvera, we detected two point-like objects invisible at previous shallow observations. However, accurate astrometry showed that none of them can be identified with the pulsar. We put new upper limits on its optical brightness of $g' > 27.87$ and $i' > 26.84$. We also reanalyzed all available archival X-ray data on Calvera. Comparison of the Calvera thermal emission parameters and upper limits on optical and non-thermal X-ray emission with respective data on rotation-powered pulsars shows that Calvera might belong to the class of ordinary middle-aged pulsars, if we assume that its distance is in the range of $1.5-5$ kpc.

We report the results of XMM-Newton observations of the Galactic mixed-morphology supernova remnant G350.0$-$2.0. Diffuse thermal X-ray emission fills the north-western part of the remnant surrounded by radio shell-like structures. We did not detect any X-ray counterpart of the latter structures, but found several bright blobs within the diffuse emission. The X-ray spectrum of the most part of the remnant can be described by a collisionally-ionized plasma model VAPEC with solar abundances and a temperature of $\approx 0.8$ keV. The solar abundances of plasma indicate that the X-ray emission comes from the shocked interstellar material. The overabundance of Fe was found in some of the bright blobs. We also analysed the brightest point-like X-ray source 1RXS J172653.4$-$382157 projected on the extended emission. Its spectrum is well described by the two-temperature optically thin thermal plasma model MEKAL typical for cataclysmic variable stars. The cataclysmic variable source nature is supported by the presence of a faint ($g\approx21$) optical source with non-stellar spectral energy distribution at the X-ray position of 1RXS J172653.4$-$382157. It was detected with the XMM-Newton optical/UV monitor in the $U$ filter and was also found in the archival H$\alpha$ and optical/near-infrared broadband sky survey images. On the other hand, the X-ray spectrum is also described by the power law plus thermal component model typical for a rotation powered pulsar. Therefore, the pulsar interpretation of the source cannot be excluded. For this source, we derived the upper limit for the pulsed fraction of 27 per cent.

PSR J1357$-$6429 is a young radio pulsar that was detected in X-rays and $\gamma$-rays. We present the high spatial resolution near-infrared imaging of the pulsar field in $J$, $H$ and $K_s$ bands obtained with the VLT/NaCo using the Adaptive Optic system. We found a faint source at the most precise pulsar radio position which we propose as the pulsar near-infrared counterpart candidate. It is confidently detected in the $J$ and $K_s$ bands, with $J$ = 23.51$\pm$0.24 and $K_s$ = 21.82$\pm$0.25. There is a hint of the source in the $H$ band with an upper limit $H$ $>$ 22.8. The dereddened source fluxes are compatible with the extrapolation of the pulsar X-ray spectrum towards the near-infrared. If the candidate is the true counterpart, by this property PSR J1357$-$6429 would be similar to the nearby middle-age pulsar PSR B0656+14. In this case, both pulsars demonstrate an unusually high near-infrared efficiency relative to the X-ray efficiency as compared to other pulsars detected in both ranges.

We analysed Chandra observations of the bright Fermi pulsar J0633+0632 and found evidence of an absorption feature in its spectrum at $804^{+42}_{-26}$ eV (the errors here and below are at 90% confidence) with equivalent width of $63^{+47}_{-36}$ eV. In addition, we analysed in detail the X-ray spectral continuum taking into account correlations between the interstellar absorption and the distance to the source. We confirm early findings by Ray et al. (2011) that the spectrum contains non-thermal and thermal components. The latter is equally well described by the blackbody and magnetised atmosphere models and can be attributed to the emission from the bulk of the stellar surface in both cases. The distance to the pulsar is constrained in a range of 1--4 kpc from the spectral fits. We infer the blackbody surface temperature of $108^{+22}_{-14}$ eV, while for the atmosphere model, the temperature, as seen by a distant observer, is $53^{+12}_{-7}$ eV. In the latter case J0633+0632 is one of the coldest middle-aged isolated neutron stars with measured temperatures. Finally, it powers an extended pulsar wind nebula whose shape suggests a high pulsar proper motion. Looking backwards the direction of the presumed proper motion we found a likely birthplace of the pulsar -- the Rosette nebula, a 50-Myr-old active star-forming region located at about 1$.\!\!^\circ$5 from the pulsar. If true, this constrains the distance to the pulsar in the range of 1.2--1.8 kpc.

PSR J1357$-$6429 is a young and energetic radio pulsar detected in X-rays and $\gamma$-rays. It powers a compact pulsar wind nebula with a jet visible in X-rays and a large scale plerion detected in X-ray and TeV ranges. Previous multiwavelength studies suggested that the pulsar has a significant proper motion of about 180 mas yr$^{-1}$ implying an extremely high transverse velocity of about 2000 km s$^{-1}$. In order to verify that, we performed radio-interferometric observations of PSR J1357$-$6429 with the the Australia Telescope Compact Array (ATCA) in the 2.1 GHz band. We detected the pulsar with a mean flux density of $212\pm5$ $\mu$Jy and obtained the most accurate pulsar position, RA = 13:57:02.525(14) and Dec = $-$64:29:29.89(15). Using the new and archival ATCA data, we did not find any proper motion and estimated its 90 per cent upper limit $\mu < 106$ mas yr$^{-1}$. The pulsar shows a highly polarised single pulse, as it was earlier observed at 1.4 GHz. Spectral analysis revealed a shallow spectral index $\alpha_{\nu}$ = $0.5 \pm 0.1$. Based on our new radio position of the pulsar, we disclaim its optical counterpart candidate reported before.

We present spectral and timing analyses of the X-ray emission from the pulsar wind nebula DA 495 and its central object, J1952.2+2925, suggested to be the pulsar, using archival Chandra and XMM-Newton data. J1952.2+2925 has a pure thermal spectrum which is equally well fitted either by the blackbody model with a temperature of $\approx 215$ eV and an emitting area radius of $\approx 0.6$ km or by magnetized neutron star atmosphere models with temperatures of 80-90 eV. In the latter case the thermal emission can come from the entire neutron star surface which temperature is consistent with standard neutron star cooling scenarios. We place also an upper limit on the J1952.2+2925 nonthermal flux. The derived spectral parameters are generally compatible with published ones based only on the Chandra data, but they are much more accurate due to the inclusion of XMM-Newton data. No pulsations were found and we placed an upper limit for the J1952.2+2925 pulsed emission fraction of 40 per cent. Utilizing the interstellar absorption-distance relation, we estimated the distance to DA 495, which can be as large as 5 kpc if J1952.2+2925 emission is described by the atmosphere models. We compiled possible multi-wavelength spectra of the nebula including radio data; they depend on the spectral model of the central object. Comparing the results with other pulsar plus wind nebula systems we set reasonable constraints on the J1952.2+2925 spin-down luminosity and age. We suggest that the Fermi source 3FGL J1951.6+2926 is the likely $\gamma$-ray counterpart of J1952.2+2925.

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.

We present results of the spectral analysis of the X-ray emission from the middle-aged Fermi pulsar J1741$-$2054 using all Chandra archival data collected in 2010 and 2013. We confirm early findings by Romani et al. 2010 that the pulsar spectrum contains a thermal emission component. The component is best described by the blackbody model with the temperature $\approx$ 60 eV and the emitting area radius $\approx$ 17 $D_{\rm kpc}$ km. The thermal emission likely originates from the entire surface of the cooling neutron star if the distance to the pulsar is $\approx$ 0.8 kpc. The latter is supported by a large absorbing column density inferred from the X-ray fit and empirical optical extinction--distance relations along the pulsar line of sight. The neutron star surface temperature and characteristic age make it similar to the well studied middle-aged pulsar B1055$-$52. Like the latter PSR J1741$-$2054 is hotter than the standard cooling scenario predicts.

A middle-aged radio-quiet pulsar J0357+3205 was discovered in gamma-rays with $Fermi$ and later in X-rays with $Chandra$ and $XMM$-$Newton$ observatories. It produces an unusual thermally-emitting pulsar wind nebula observed in X-rays. Deep optical observations were obtained to search for the pulsar optical counterpart and its nebula using the Gran Telescopio Canarias (GTC). The direct imaging mode in the Sloan $g'$ band was used. Archival X-ray data were reanalysed and compared with the optical data. No pulsar optical counterpart was detected down to $g'\geq~28_{\cdotp}^{\text{m}}1$. No pulsar nebula was either identified in the optical. We confirm early results that the X-ray spectrum of the pulsar consists of a nonthermal power-law component of the pulsar magnetospheric origin dominating at high energies and a soft thermal component from the neutron star surface. Using magnetised partially ionised hydrogen atmosphere models in X-ray spectral fits we found that the thermal component can come from entire surface of the cooling neutron star with a temperature of 36$^{+9}_{-6}$ eV, making it one of the coldest among cooling neutron stars known. The surface temperature agrees with the standard neutron star cooling scenario. The optical upper limit does not put any additional constraints on the thermal component, however it implies a strong spectral break for the nonthermal component between the optical and X-rays as is observed in other middle-aged pulsars. The thermal emission from the entire surface of the neutron star likely dominates over the nonthermal emission in the UV range. Observations of the PSR J0357+3205 in this range are promising to put more stringent constraints on its thermal properties.

We report the first high spatial resolution near-infrared imaging of the Vela pulsar in the $K_s$ band obtained with the new adaptive optics system recently mounted on the Gemini-South telescope. For the first time, we have firmly detected the pulsar in this band with $K_s$ $\approx$ 21\fm8, and have resolved in detail an extended feature barely detected previously in the immediate vicinity of the pulsar in the $J_sH$ bands. The pulsar $K_s$ flux is fully consistent with the extension of the flat optical spectrum of the pulsar towards the infrared and does not confirm the strong infrared flux excess in the pulsar emission suggested earlier by the low spatial resolution data. The extended feature is about two times brighter than the pulsar and is likely associated with its X-ray counter-jet. It extends $\sim$ 2\asec southwards of the pulsar along the X-ray ounter-jet and shows knot-like structures and a red spectrum.

The young radio pulsar J1124-5916 is associated with a Cas A like supernova remnant G292.0+1.8. It powers a compact torus-like pulsar wind nebula with a jet first detected in X-rays and then identified in the optical and mid-infrared. We carried out deep near-infrared observations of the pulsar field to identify the pulsar and its nebula in this range. The direct imaging mode of the NACO adaptive optics instrument at the ESO VLT in the H and Ks bands was used. In both bands we detected a faint, H=21.30(10) and Ks=20.45(10), extended elliptical object, whose center position is consistent with the X-ray position of the pulsar. The morphology of the object and the orientation of its major axis are in a good agreement with those observed for the pulsar torus-like nebula in the mid-infrared, optical, and X-rays. This suggests that it is the near-infrared counterpart of the nebula. The measured fluxes compiled with the data in other ranges show a complicated unabsorbed power law spectrum of the torus-like nebula with several steep breaks between the near-infrared and mid-infrared, the optical and X-rays, and possibly in the mid-infrared. This implies a multiple relativistic particle population responsible for the synchrotron emission of the nebula in different spectral ranges. We have not resolved the pulsar counterpart from its nebula and place only upper limits on its brightness, H >= 23.9$ and Ks >= 22.7. Based on that, its contribution to the total near-infrared flux of the pulsar+nebula system is <= 10%, which is comparable with the expected contribution in the optical.

We present and discuss new visual wavelength-range observations of the inner regions of the supernova remnant SNR 0540-69.3 that is located in the Large Magellanic Cloud (LMC). These observations provide us with more spatial and spectral information than were previously available for this object. We use these data to create a detailed three-dimensional model of the remnant, assuming linear expansion of the ejecta. With the observations and the model we study the general three-dimensional structure of the remnant, and the influence of an active region in the remnant - a "blob" - that we address in previous papers. We used the fibre-fed integral-field spectrograph VIMOS at the Very Large Telescope of the European Southern Observatory. The observations provide us with three-dimensional data in [OIII]5007 and [SII]6717,6731 at an 0.33x0.33 spatial sampling and a velocity resolution of about 35 km/s. We decomposed the two, partially overlapping, sulphur lines and used them to calculate electron densities across the remnant at high signal-to-noise ratio. Our analysis reveals a structure that stretches from the position of the "blob", and into the plane of the sky at a position angle of about 60 degrees. We speculate that the pulsar is positioned along this activity axis, where it has a velocity along the line of sight of a few hundred km/s. The "blob" is most likely a region of shock activity, as it is mainly bright in [SII]; future observations of [OII]3729 would be useful to test whether the S/O abundance ratio is higher than average for that location in the remnant. The striking resemblance in X-rays between the pulsar wind nebula (PWN) of SNR 0540-69.3 and the Crab, in combination with our findings in this paper, suggests that the symmetry axis is part of a torus in the PWN. (abridged)

(abridged) PSR J1048-5832 is a young Vela-like pulsar that was recently detected in gamma-rays with Fermi, and also in X-rays with Chandra and XMM-Newton. It powers a compact PWN visible in X-rays. We present deep optical observations with the ESO Very Large Telescope to search for optical counterparts of the pulsar and its nebula and to explore their multi-wavelength emission properties. The data were obtained in V and R bands and compared with archival data in other spectral domains. We do not detect the pulsar in the optical and derive informative upper limits of R > 28.1 mag and V > 28.4 mag for its brightness. Using a red-clump star method, we estimate an interstellar extinction towards the pulsar of 2 mag, which is consistent with the absorbing column density derived form X-rays. The respective distance agrees with the dispersion measure distance. We reanalyse the Chandra X-ray data and compare the dereddened upper limits with the unabsorbed X-ray spectrum of the pulsar. We find that regarding its optical-X-ray spectral properties this gamma-ray pulsar is not distinct from other pulsars detected in both ranges. However, like the Vela pulsar, it is very inefficient in the optical and X-rays. Among a dozen optical sources overlapping with the pulsar X-ray nebula we find one with V = 26.9 and R = 26.3 mag, whose colour is slightly bluer then that of the field stars and consistent with the peculiar colours typical for pulsar nebula features. It positionally coincides with a relatively bright feature of the pulsar X-ray nebula, resembling the Crab wisp and locating in 2 arcsec from the pulsar. We suggest this source as a counterpart candidate to the feature. Based on the substantial interstellar extinction towards the pulsar and its optical inefficiency, further optical studies should be carried out at longer wavelengths.

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.

The Vela and Geminga pulsars are rotation powered neutron stars, which have been identified in various spectral domains, from the near-infrared to hard $\gamma$-rays. In the near-infrared they exhibit tentative emission excesses, as compared to the optical range. To check whether these features are real, we analysed archival mid-infrared broadband images obtained with the Spitzer Space Telescope in the 3.6--160 $\mu$m range and compared them with the data in other spectral domains. In the 3.6 and 5.8 $\mu$m bands we detected at $\sim$ (4--5)$\sigma$ significance level a point-like object, that is likely to be the counterpart of the Vela pulsar. Its position coincides with the pulsar at < 0.4 arcsec 1$\sigma$-accuracy level. Combining the measured fluxes with the available multiwavelength spectrum of the pulsar shows a steep flux increase towards the infrared, confirming the reality of the near-infrared excess reported early, and, hence, the reality of the suggested mid-infrared pulsar identification. Geminga is also identified, but only at a marginal 2$\sigma$ detection level in one 3.6 $\mu$m band. This needs a farther confirmation by deeper observations, while the estimated flux is also compatible with the near-infrared Geminga excess. The detection of the infrared excess is in contrast to the Crab pulsar, where it is absent, but is similar to the two magnetars, 4U 0142+61 and 1E 2259+586, showing similar features. We discuss X-ray irradiated fall-back discs around the pulsars, unresolved pulsar nebula structures, and pulsar magnetospheres as possible origins of the excesses. We note also possible infrared signatures of an extended tail behind Geminga and of the Vela plerion radio lobes.

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 present high spatial resolution optical imaging and polarization observations of the PSR B0540-69.3 and its highly dynamical pulsar wind nebula (PWN) performed with HST, and compare them with X-ray data obtained with the Chandra X-ray Observatory. We have studied the bright region southwest of the pulsar where a bright "blob" is seen in 1999. We show that it may be a result of local energy deposition around 1999, and that the emission from this then faded away. Polarization data from 2007 show that the polarization properties show dramatic spatial variations at the 1999 blob position arguing for a local process. Several other positions along the pulsar-"blob" orientation show similar changes in polarization, indicating previous recent local energy depositions. In X-rays, the spectrum steepens away from the "blob" position, faster orthogonal to the pulsar-"blob" direction than along this axis of orientation. This could indicate that the pulsar-"blob" orientation is an axis along where energy in the PWN is mainly injected, and that this is then mediated to the filaments in the PWN by shocks. We highlight this by constructing an [S II]-to-[O III]-ratio map. We argue, through modeling, that the high [S II]/[O III] ratio is not due to time-dependent photoionization caused by possible rapid Xray emission variations in the "blob" region. We have also created a multiwavelength energy spectrum for the "blob" position showing that one can, to within 2sigma, connect the optical and X-ray emission by a single power law. We obtain best power-law fits for the X-ray spectrum if we include "extra" oxygen, in addition to the oxygen column density in the interstellar gas of the Large Magellanic Cloud and the Milky Way. This oxygen is most naturally explained by the oxygen-rich ejecta of the supernova remnant. The oxygen needed likely places the progenitor mass in the 20 - 25 Msun range.

G292.0+1.8 is a Cas A-like supernova remnant that contains the young pulsar PSR J1124-5916 powering a compact torus-like pulsar wind nebula visible in X-rays. A likely counterpart to the nebula has been detected in the optical VRI bands. To confirm the counterpart candidate nature, we examined archival mid-infrared data obtained with the Spitzer Space Telescope. Broad-band images taken at 4.5, 8, 24, and 70 microns were analyzed and compared with available optical and X-ray data. The extended counterpart candidate is firmly detected in the 4.5 and 8 micron bands. It is brighter and more extended in the bands than in the optical, and its position and morphology agree well with the coordinates and morphology of the torus-like pulsar wind nebula in X-rays. The source is not visible in 24 and 70 micron images, which are dominated by bright emission from the remnant shell and filaments. We compiled the infrared fluxes of the nebula, which probably contains a contribution from an unresolved pulsar in its center, with the optical and X-ray data. The resulting unabsorbed multiwavelength spectrum is described by power laws of significantly steeper slope in the infrared-optical than in X-rays, implying a double-knee spectral break between the optical and X-rays. The 24 and 70 microns flux upper limits suggest a second break and a flatter spectrum at the long wavelength limit. These features are common to two other pulsar wind nebulae associated with the remnants B0540-69.3 and 3C 58 and observed in all three ranges. The position, morphology, and spectral properties of the detected source allow us to comfirm that it is the infrared-optical counterpart to both the pulsar and its wind nebula system in the G292.0+1.8 supernova remnant.

G292.0+1.8 is the Cas A-like supernova remnant containing the young pulsar PSR J1124-5916, which powers a compact torus-like pulsar wind nebula with a jet visible in X-rays. We have performed deep optical observations of the pulsar field to detect the optical counterpart of the pulsar and its nebula. The observations were carried out using the direct imaging mode of FORS2 at the ESO VLT/UT1 telescope in the V, R, and I bands. We also analyzed archival images obtained with the Chandra/ACIS-I, ACIS-S, and HRC-S in X-rays. In all three optical bands we detect a faint elliptical nebulosity, whose brightness peak and center position are consistent at a sub-arcsecond level with the X-ray position of the pulsar. The field is densely packed with background stars, but after subtraction of these stars the morphology of the object and the orientation of its major axis appear to be in a good agreement with the brightest inner part of the pulsar nebula torus region seen almost edge on in X-rays. Within the nebulosity we do not resolve any point-like optical object that could be identified with the pulsar and estimate its contribution to the observed nebulosity flux as <~20%. Extracting the X-ray spectrum from the physical region equivalent to the optical source position and extent and combining that with the measured optical fluxes, we compile a tentative multi-wavelength spectrum of the inner part of the nebula. Within uncertainties of the interstellar extinction towards G292.0+1.8 it is reminiscent of either the Crab or PSR B540-69 and J0205+6449 pulsar wind nebula spectra. The position, morphology, and spectral properties of the detected nebulosity suggest that it is the likely optical counterpart of the pulsar plus its wind nebula system in G292.0+1.8.

We have performed a deep optical imaging of 3C 58 SNR with the NOT in the B and V bands to detect the optical counterpart of the associated pulsar J0295+6449 and its torus-like wind nebula visible in X-rays. We analyzed our data together with the archival data obtained with the Chandra in X-rays and with the Spitzer in the mid-IR. We detect a faint extended elliptical object with B=24.06 and V=23.11 whose peak brightness and center position are consistent at the sub-arcsecond level with the position of the pulsar. Its morphology and orientation are in excellent agreement with the torus-like pulsar nebula, seen almost edge on in X-rays although its extension is only about a half of that in X-rays. In the optical we likely see only the brightest central part of the torus with the pulsar. The object is identical to the counterpart of the torus recently detected in the mid-IR. The estimated pulsar contribution to the optical flux is less than 10%. Combinig the optical/mid-IR fluxes and X-ray power-law spectrum extracted from the spatial region constrained by the optical/IR source extent we compile a tentative multi-wavelength spectrum of the central part of the nebula. Within uncertainties of the interstellar extinction it is reminiscent of either the Crab or B0540-69 pulsar wind nebula spectra. The properties of the object strongly suggest it to be the optical counterpart of the 3C 58 pulsar + its wind nebula system, making 3C 58 the third member of such a class of the torus-like systems identified in the optical and mid-IR.

Aims: We performed deep optical observations of the field of an old, fast-moving radio pulsar PSR B1133+16 in an attempt to detect its optical counterpart and a bow shock nebula. Methods: The observations were carried out using the direct imaging mode of FORS1 at the ESO VLT/UT1 telescope in the B, R, and H_alpha bands. We also used archival images of the same field obtained with the VLT in the B band and with the Chandra/ACIS in X-rays. Results: In the B band we detected a faint (B=28.1+/-0.3) source that may be the optical counterpart of PSR B1133+16, as it is positionally consistent with the radio pulsar and with the X-ray counterpart candidate published earlier. Its upper limit in the R band implies a color index B-R <0.5, which is compatible with the index values for most pulsars identified in the optical range. The derived optical luminosity and its ratio to the X-ray luminosity of the candidate are consistent with expected values derived from a sample of pulsars detected in both spectral domains. No Balmer bow shock was detected, implying a low density of ambient matter around the pulsar. However, in the X-ray and H_alpha images we found the signature of a trail extending ~4"-5" behind the pulsar and coinciding with the direction of its proper motion. If confirmed by deeper studies, this is the first time such a trail has been seen in the optical and X-ray wavelengths. Conclusions: Further observations at later epochs are necessary to confirm the identification of the pulsar by the candidate's proper motion measurements.

We have obtained the spectrum of a middle-aged PSR B0656+14 in the 4300-9000 AA range with the ESO/VLT/FORS2. Preliminary results show that at 4600-7000 AA the spectrum is almost featureless and flat with a spectral index $\alpha_nu ~ -0.2 that undergoes a change to a positive value at longer wavelengths. Combining with available multiwavelength data suggests two wide, red and blue, flux depressions whose frequency ratio is about 2 and which could be the 1st and 2nd harmonics of electron/positron cyclotron absorption formed at magnetic fields ~10^8G in upper magnetosphere of the pulsar.

We carried out a deep subarcsecond BRI imaging of the two middle-aged pulsars to establish their properties in the optical range. Both pulsars are detected at >10 sigma level. Geminga is for the first time reliably detected in the I band with a magnitude of 25.10+/-0.14. We also reanalyze archival ESO/NTT and HST broadband data and find that some published fluxes for Geminga were estimated inaccurately. The resulting dereddened broadband spectra of both pulsars are remarkably similar to each other and show significant flux increases towards the far-UV and near-IR, and a wide flux excess in V-I bands. This suggests a multicomponent structure of the optical emission. The nonthermal power law component of the pulsar magnetospheric origin dominates in the most part of the optical range. For PSR B0656+14 it is compatible with a low energy extension of the power law tail seen in hard X-rays. For Geminga the respective extension overshoots by a factor of 100 the nonthermal optical flux, which has a less steep spectral slope than in X-rays. This implies a spectral break at a photon energy of about 1 keV. The flux increases towards the far-UV are compatible with contributions of the Rayleigh-Jeans parts of the blackbody components from whole surfaces of the neutron stars dominating in soft X-rays. The V-I excess, which is most significant for PSR B0656+14, suggests a third spectral component of still unidentified origin. Faint, a few arcseconds in size nebulae extended perpendicular to the proper motion directions of the pulsars, are seen around both objects in our deepest I band images. They can be optical counterparts of the bow-shock head of Geminga and of the tentative pulsar wind nebula of PSR B0656+14 observed in X-rays.

We have observed the middle-aged pulsar PSR B0656+14 with the prism and the NUV MAMA detector of the Space Telescope Imaging Spectrograph (STIS) to measure the pulsar spectrum and periodic pulsations in the near-ultraviolet (NUV). The pulsations are clearly detected, double-peaked and very similar to the optical pulse profile. The NUV pulsed fraction is 70 +/- 12%. The spectral slope of the dereddened phase-integrated spectrum in the 1800 - 3200 A range is 0.35 +/- 0.5 which together with the high pulse fraction indicates a non-thermal origin for the NUV emission. The total flux in the range 1700-3400 A is estimated to be 3.4 +/- 0.3e-15 erg/s/cm2 when corrected for E(B-V)=0.03. At a distance of 288 pc this corresponds to a luminosity 3.4e28 erg/s assuming isotropy of the emission. We compare the NUV pulse profile with observations from radio to gamma-rays. The first NUV sub-pulse is in phase with the gamma-ray pulse marginally detected with the EGRET, while the second NUV sub-pulse is similar both in shape and in phase with the non- thermal pulse in hard X-rays. This indicates a single origin of the non-thermal emission in the optical-NUV and in the X-rays. This is also supported by the observed NUV spectral slope, which is compatible with a blackbody plus power-law fit extended from the X-ray range, but dominated by the power-law component in most of the NUV range.

We have used the ESO NTT/EMMI and VLT/FORS1 instruments to examine the LMC supernova remnant 0540-69.3 as well as its pulsar (PSR B0540-69) and pulsar-powered nebula in the optical range.Spectroscopic observations of the remnant covering the range of 3600-7350 A centered on the pulsar produced results consistent with those of Kirshner et al. (1989), but also revealed many new emission lines. The most important are [Ne III] 3869, 3967 and Balmer lines of hydrogen. In both the central part of the remnant, as well as in nearby H II regions, the [O III] temperature is higher than about 20 000 K, but lower than previously estimated. For PSR B0540-69, previous optical data are mutually inconsistent: HST/FOS spectra indicate a significantly higher absolute flux and steeper spectral index than suggested by early time-resolved groundbased UBVRI photometry. We show that the HST and VLT spectroscopic data for the pulsar have more then about 50% nebular contamination, and that this is the reason for the previous difference. Using HST/WFPC2 archival images obtained in various bands from the red part of the optical to the NUV range we have performed an accurate photometric study of the pulsar, and find that the spectral energy distribution of the pulsar emission has a negative slope with 1.07 +/- 0.2. This is steeper than derived from previous UBVRI photometry, and also different from the almost flat spectrum of the Crab pulsar. We also estimate that the proper motion of the pulsar is 4.9 +/- 2.3 mas/yr, corresponding to a transverse velocity of 1190 +/- 560 km/s, projected along the southern jet of the pulsar nebula.

Using available multiwavelength data for the pulsars of different ages detected in the optical range, we analyze the efficiencies of the conversion of the pulsar spindown power \dot{E} into the observed non-thermal luminosity L in different spectral domains. We find that these pulsars show significantly non-monotonic dependencies of the optical and X-ray efficiencies (\eta = L/\dot{E}) versus pulsar age with a pronounced minimum at the beginning of the middle-age epoch and comparably higher efficiencies of younger and older pulsars. In addition, we find a strong correlation between the optical and 2--10 keV X-ray luminosities of these pulsars that implies the same origin of their nonthermal emission in both spectral domains.

The young PSR B0540-69 (B0540) in the LMC is the only pulsar (except the Crab pulsar) for which a near-UV spectrum has been obtained. However, the absolute flux and spectral index of previously published HST/FOS data are significantly higher than suggested by broadband groundbased UBVRI photometry. Using our ESO/VLT/FORS1 spectral observations and HST/WFPC2 archival images we show that the old HST and new VLT spectral data are >50% contaminated by the Pulsar Wind Nebula (PWN) and that this is the reason for the above mentioned difference. We find that the broadband HST spectrum for the range 3300-8000 A is clearly non- thermal and has a negative spectral index of 1.07(+0.20/-0.19). This is dif- ferent from the almost flat spectrum of the Crab pulsar. The PWN of B0540 shows a clear asymmetry of the surface brightness distribution along the major axis of the PWN torus-like structure with respect to the pulsar position, also seen in Chandra X-ray images. This can be linked to the asymmetry of the surrounding SN ejecta. We find no significant spectral index variation over the PWN. Using HST archival images we estimate the proper motion of B0540 to be 4.9+/-2.3 mas/yr, i.e. a transverse velocity of 1190+/-560 km/s along the southern jet of the PWN. This can make PSR B0540 the third pulsar with a proper motion aligned with the jet axis of its PWN, which poses constraints on pulsar kick models. We discuss the interstellar absorption toward B0540 including the contributions from the Milky Way, LMC and the supernova ejecta, and compare unabsorbed multi- wavelength spectra of B0540 and the Crab pulsar, and their PWNs. Compared with the Crab, B0540 and its PWN show a significant depression in the optical range.

We report on the first detection of the Vela pulsar in the near-infrared with the VLT/ISAAC in the Js and H bands. The pulsar magnitudes are Js=22.71 +/- 0.10 and H=22.04 +/- 0.16. We compare our results with the available multiwavelength data and show that the dereddened phase-averaged optical spectrum of the pulsar can be fitted with a power law F_nu propto nu^(-alpha_nu) with alpha_nu = 0.12 +/- 0.05, assuming the color excess E(B-V)=0.055 +/-0.005 based on recent spectral fits of the emission of the Vela pulsar and its supernova remnant in X-rays. The negative slope of the pulsar spectrum is different from the positive slope observed over a wide optical range in the young Crab pulsar spectrum. The near-infrared part of the Vela spectrum appears to have the same slope as the phase-averaged spectrum in the high energy X-ray tail, obtained in the 2-10 keV range with the RXTE. Both of these spectra can be fitted with a single power law suggesting their common origin. Because the phase-averaged RXTE spectrum in this range is dominated by the second X-ray peak of the pulsar light curve, coinciding with the second main peak of its optical pulse profile, we suggest that this optical peak can be redder than the first one. We also detect two faint extended structures in the 1.5''-3.1'' vicinity of the pulsar, projected on and aligned with the south-east jet and the inner arc of the pulsar wind nebula, detected in X-rays with Chandra. We discuss their possible association with the nebula.

We report on deep BVR-imaging of the field of the nearby millisecond pulsar PSR J0030+0451 obtained with the ESO/VLT/FORS2. We do not detect any optical counterpart down to B>27.3, V>27.0 and R>27.0 in the immediate vicinity of the radio pulsar position. The closest detected sources are offset by >3'', and they are excluded as counterpart candidates by our astrometry. Using our upper limits in the optical, and including recent XMM-Newton X-ray data we show that any nonthermal power-law spectral component of neutron star magnetospheric origin, as suggested by the interpretation of X-ray data, must be suppressed by at least a factor of ~500 in the optical range. This either rules out the nonthermal interpretation or suggests a dramatic spectral break in the 0.003-0.1 keV range of the power-law spectrum. Such a situation has never been observed in the optical/X-ray spectral region of ordinary pulsars, and the origin of such a break is unclear. An alternative interpretation with a purely thermal X-ray spectrum is consistent with our optical upper limits. In this case the X-ray emission is dominated by hot polar caps of the pulsar.

We report the B band optical observations of an old (17.5 Myr) radiopulsar PSR B0950+08 obtained with the Suprime-Cam at the Subaru telescope. We detected a faint object, B=27.07(16). Within our astrometrical accuracy it coincides with the radio position of the pulsar and with the object detected earlier by Pavlov et al. (1996) in UV with the HST/FOC/F130LP. The positional coincidence and spectral properties of the object suggest that it is the optical counterpart of PSR B0950+08. Its flux in the B band is two times higher than one would expect from the suggested earlier Rayleigh-Jeans interpretation of the only available HST observations in the adjacent F130LP band. Based on the B and F130LP photometry of the suggested counterpart and on the available X-ray data we argue in favour of nonthermal origin of the broad-band optical spectrum of PSR B0950+08, as it is observed for the optical emission of the younger, middle-aged pulsars PSR B0656+14 and Geminga. At the same time, the optical efficiency of PSR B0950+08, estimated from its spin-down power and the detected optical flux, is by several orders of magnitude higher than for these pulsars, and comparable with that for the much younger and more energetic Crab pulsar. We cannot exclude the presence of a compact, about 1'', faint pulsar nebula around PSR B0950+08, elongated perpendicular to the vector of its proper motion, unless it is not a projection of a faint extended object on the pulsar position.

We study the thermodynamic properties of a partially ionized hydrogen plasma in strong magnetic fields, B ~ 10^{12}-10^{13} G, typical of neutron stars. The properties of the plasma depend significantly on the quantum-mechanical sizes and binding energies of the atoms, which are strongly modified by thermal motion across the field. We use new fitting formulas for the atomic binding energies and sizes, based on accurate numerical calculations and valid for any state of motion of the atom. In particular, we take into account decentered atomic states, neglected in previous studies of thermodynamics of magnetized plasmas. We also employ analytic fits for the thermodynamic functions of nonideal fully ionized electron-ion Coulomb plasmas. This enables us to construct an analytic model of the free energy. An ionization equilibrium equation is derived, taking into account the strong magnetic field effects and the nonideality effects. This equation is solved by an iteration technique. Ionization degrees, occupancies, and the equation of state are calculated.

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.

Fields of Geminga, PSR B0950+08, PSR J1908+0734 and PSR J0108-1431 have been observed within the program of optical investigation of nearby neutron stars carried out with the 6m telescope. Our multicolour BVR_cI_c photometry has yielded the following magnitudes of Geminga: V=25.3+/-0.4; R_c=25.4+/-0.3. The estimates of the 3\sigma upper limits in the B and I_c bands are 25.4 and 24.9, respectively. This is in agreement with previous observations of Geminga. We have suggested candidate for optical counterpart of PSR B0950+08, with R_c=25.4+/-0.3. When compared to the HST data obtained in UV its flux shows that this old pulsar can be brighter in near-IR than in near-UV. We have obtained for the first time images of the fields of PSR J1908+0734 and PSR J0108-1431. Although no optical counterparts have been detected, the magnutude upper limits obtained suggest further observations of these pulsars.

We study the heat capacity and neutrino emission reactions (direct and modified Urca processes, nucleon-nucleon bremsstrahlung, Cooper pairing of nucleons) in matter of supranuclear density of the neutron star cores with superfluid neutrons and protons. Various superfluidity types are analysed (singlet-state pairing and two types of triplet-state pairing, without and with nodes of the gap at a nucleon Fermi surface). The results are used for cooling simulations of isolated neutron stars. Both, the standard cooling and the cooling enhanced by the direct Urca process, are strongly affected by nucleon superfluidity. Comparison of cooling theory of isolated neutron stars with observations of their thermal radiation may give stringent constraints on the critical temperatures of the neutron and proton superfluidities in the neutron star cores.

We study partially ionized hydrogen plasma at temperatures T = 10^{5.5}-10^{6.5} K and magnetic fields B = 10^{12}-10^{13} G, typical of photospheres of isolated neutron stars. We construct an analytic model of the plasma free energy and derive a generalized Saha equation which is used to obtain the opacities.

Cooling of neutron stars (NSs) with superfluid cores is simulated taking into account neutrino emission produced by Cooper pairing of nucleons. The critical temperatures of neutron and proton superfluidities, $T_{cn}$ and $T_{cp}$, are assumed to be constant over the NS core, and treated as free parameters. They are constrained using the surface temperatures $T_s$ of isolated NSs (RX J0822-43, PSR 1055-52, 1E 1207-52, Vela, Geminga, PSR 0656+14, RX J0002+62), obtained by interpretation of observed thermal radiation either with black body spectrum or with hydrogen atmosphere models.

We study hydrogen plasmas at magnetic fields B ~ 10^{12}-10^{13} Gauss, densities ~ 10^{-3}-10^3 g/cc and temperatures T ~ 10^{5.5}-10^{6.5} K, typical of photospheres of middle-aged cooling neutron stars. We construct an analytical free energy model of the partially ionized plasma, including into consideration the decentred atomic states, which arise due to the thermal motion across the strong field. We show that these states, neglected in previous studies, may contribute appreciably into thermodynamics of the outer atmospheric layers at density below 1 g/cc and typical B and T. We take into account Coulomb non-ideality of the ionized component of the plasma affected by intense magnetic field. Ionization degree, occupancies and equation of state are calculated, and their dependences on the temperature, density and magnetic field are studied.

Cooling of neutron stars with dipole magnetic fields is simulated using a realistic model of the anisotropic surface temperature distribution produced by magnetic fields. Suppression of the electron thermal conductivity of outer stellar layers across the field increases thermal isolation of these layers near the magnetic equator. Enhancement of the radiative and longitudinal electron thermal conductivities in quantizing magnetic fields reduces thermal isolation near the magnetic poles. The equatorial increase of the isolation is pronounced for $B \ga 10^{10}$ G, while the polar decrease -- for $B \ga 10^{12}$ G. The effects compensate partly each other, and the actual influence of the magnetic fields on the cooling is weaker than predicted by the traditional theories where the equatorial effects have been neglected.

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.