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

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