We describe system verification tests and early science results from the pulsar processor (PTUSE) developed for the newly-commissioned 64-dish SARAO MeerKAT radio telescope in South Africa. MeerKAT is a high-gain (~2.8 K/Jy) low-system temperature (~18 K at 20cm) radio array that currently operates from 580-1670 MHz and can produce tied-array beams suitable for pulsar observations. This paper presents results from the MeerTime Large Survey Project and commissioning tests with PTUSE. Highlights include observations of the double pulsar J0737-3039A, pulse profiles from 34 millisecond pulsars from a single 2.5h observation of the Globular cluster Terzan 5, the rotation measure of Ter5O, a 420-sigma giant pulse from the Large Magellanic Cloud pulsar PSR J0540-6919, and nulling identified in the slow pulsar PSR J0633-2015. One of the key design specifications for MeerKAT was absolute timing errors of less than 5 ns using their novel precise time system. Our timing of two bright millisecond pulsars confirm that MeerKAT delivers exceptional timing. PSR J2241-5236 exhibits a jitter limit of <4 ns per hour whilst timing of PSR J1909-3744 over almost 11 months yields an rms residual of 66 ns with only 4 min integrations. Our results confirm that the MeerKAT is an exceptional pulsar telescope. The array can be split into four separate sub-arrays to time over 1000 pulsars per day and the future deployment of S-band (1750-3500 MHz) receivers will further enhance its capabilities.

We report here on initial results from the Thousand Pulsar Array (TPA) programme, part of the Large Survey Project "MeerTime" on the MeerKAT telescope. The interferometer is used in tied-array mode in the band from 856 to 1712~MHz, and the wide band coupled with the large collecting area and low receiver temperature make it an excellent telescope for the study of radio pulsars. The TPA is a 5 year project which aims to observe (a) more than 1000 pulsars to obtain high-fidelity pulse profiles, (b) some 500 of these pulsars over multiple epochs, (c) long sequences of single-pulse trains from several hundred pulsars. The scientific outcomes from the programme will include determination of pulsar geometries, the location of the radio emission within the pulsar magnetosphere, the connection between the magnetosphere and the crust and core of the star, tighter constraints on the nature of the radio emission itself as well as interstellar medium studies. First results presented here include updated dispersion measures, 26 pulsars with Faraday rotation measures derived for the first time and a description of interesting emission phenomena observed thus far.

We report on the energy-resolved timing and phase-resolved spectral analysis of X-ray emission from PSR J0659+1414 observed with XMM-Newton and NuSTAR. We find that the new data rule out the previously suggested model of the phase-dependent spectrum as a three-component (2 blackbodies + power-law) continuum, which shows large residuals between $0.3-0.7$ keV. Fitting neutron star atmosphere models or several blackbodies to the spectrum does not provide a better description of the spectrum, and requires spectral model components with unrealistically large emission region sizes. The fits improve significantly if we add a phase-dependent absorption feature with central energy $0.5-0.6$ keV and equivalent width up to $\approx 50$ eV. We detected the feature for about half of the pulse cycle. Energy-resolved pulse profiles support the description of the spectrum with a three-component continuum and an absorption component. The absorption feature could be interpreted as an electron cyclotron line originating in the pulsar magnetosphere and broadened by the non-uniformity of the magnetic field along the line of sight. The significant phase-variability in the thermal emission from the entire stellar surface may indicate multi-polar magnetic fields and a non-uniform temperature distribution. The strongly pulsed non-thermal spectral component detected with NuSTAR in the $3-20$ keV range is well fit by a power-law model with a photon index $\Gamma=1.5\pm0.2$.

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

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

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

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

Empirical calibrations of the stellar age-rotation-activity relation (ARAR) rely on observations of the co-eval populations of stars in open clusters. We used the Chandra X-ray Observatory to study M37, a 500-Myr-old open cluster that has been extensively surveyed for rotation periods ($P_{\rm rot}$). M37 was observed almost continuously for five days, for a total of 440.5 ksec, to measure stellar X-ray luminosities ($L_{\mathrm{X}}$), a proxy for coronal activity, across a wide range of masses. The cluster's membership catalog was revisited to calculate updated membership probabilities from photometric data and each star's distance to the cluster center. The result is a comprehensive sample of 1699 M37 members: 426 with $P_{\rm rot}$, 278 with X-ray detections, and 76 with both. We calculate Rossby numbers, $R_o = P_{\rm rot}/\tau$, where $\tau$ is the convective turnover time, and ratios of the X-ray-to-bolometric luminosity, $L_{\rm X}/L_{\rm bol}$, to minimize mass dependencies in our characterization of the rotation-coronal activity relation at 500 Myr. We find that fast rotators, for which $R_o<0.09\pm0.01$, show saturated levels of activity, with log($L_{\rm X}/L_{\rm bol}$)$=-3.06\pm0.04$. For $R_o\geq0.09\pm0.01$, activity is unsaturated and follows a power law of the form $R_o^\beta$, where $\beta$=$-2.03_{-0.14}^{+0.17}$. This is the largest sample available for analyzing the dependence of coronal emission on rotation for a single-aged population, covering stellar masses in the range 0.4$-$1.3 $M_{\odot}$, $P_{\rm rot}$ in the range 0.4$-$12.8 d, and $L_{\rm X}$ in the range 10$^{28.4-30.5}$ erg s$^{-1}$. Our results make M37 a new benchmark open cluster for calibrating the ARAR at ages of $\approx$500 Myr.

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

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

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

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

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

Supernova fallback disks around neutron stars have been discussed to influence the evolution of the diverse neutron star populations. Slowly rotating neutron stars are most promising to find such disks. Searching for the cold and warm debris of old fallback disks, we carried out Herschel PACS (70 $\mu$m, 160 $\mu$m) and Spitzer IRAC (3.6 $\mu$m, 4.5 $\mu$m) observations of eight slowly rotating ($P\approx 3 - 11$ s) nearby ($<1$ kpc) isolated neutron stars. Herschel detected 160 $\mu$m emission ($>5\sigma$) at locations consistent with the positions of the neutron stars RX J0806.4-4123 and RX J2143.0+0654. No other significant infrared emission was detected from the eight neutron stars. We estimate probabilities of 63%, 33% and 3% that, respectively, none, one, or both Herschel PACS 160 $\mu$m detections are unrelated excess sources due to background source confusion or an interstellar cirrus. If the 160 $\mu$m emission is indeed related to cold (10 K to 22 K) dust around the neutron stars, this dust is absorbing and re-emitting $\sim 10$% to $\sim 20$% of the neutron stars' X-rays. Such high efficiencies would be at least three orders of magnitude larger than the efficiencies of debris disks around nondegenerate stars. While thin dusty disks around the neutron stars can be excluded as counterparts of the 160 $\mu$m emission, dusty asteroid belts constitute a viable option.

We present the Chasing the Identification of ASCA Galactic Objects (ChIcAGO) survey, which is designed to identify the unknown X-ray sources discovered during the ASCA Galactic Plane Survey (AGPS). Little is known about most of the AGPS sources, especially those that emit primarily in hard X-rays (2-10 keV) within the F_x ~ 10^-13 to 10^-11 erg cm^-2 s^-1 X-ray flux range. In ChIcAGO, the subarcsecond localization capabilities of Chandra have been combined with a detailed multi-wavelength follow-up program, with the ultimate goal of classifying the >100 unidentified sources in the AGPS. Overall to date, 93 unidentified AGPS sources have been observed with Chandra as part of the ChIcAGO survey. A total of 253 X-ray point sources have been detected in these Chandra observations within 3' of the original ASCA positions. We have identified infrared and optical counterparts to the majority of these sources, using both new observations and catalogs from existing Galactic plane surveys. X-ray and infrared population statistics for the X-ray point sources detected in the Chandra observations reveal that the primary populations of Galactic plane X-ray sources that emit in the F_x ~ 10^-13 to 10^-11 erg cm^-2 s^-1 flux range are active stellar coronae, massive stars with strong stellar winds that are possibly in colliding-wind binaries, X-ray binaries, and magnetars. There is also a fifth population that is still unidentified but, based on its X-ray and infrared properties, likely comprise partly of Galactic sources and partly of active galactic nuclei.

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

We report on an X-ray observation of the 166 Myr old radio pulsar J0108-1431 with XMM-Newton. The X-ray spectrum can be described by a power-law model with a relatively steep photon index Gamma~3 or by a combination of thermal and non-thermal components, e.g., a power-law component with fixed photon index Gamma~2 plus a blackbody component with a temperature of kT=0.11 keV. The two-component model appears more reasonable considering different estimates for the hydrogen column density. The non-thermal X-ray efficiency in the single power-law model is eta^PL (1-10 keV) = L^PL (1-10 keV) / Edot ~ 0.003, higher than in most other X-ray detected pulsars. In the case of the combined model, the non-thermal and thermal X-ray efficiencies are even higher, eta^PL (1-10 keV) ~ eta^bb ~ 0.006. We detected X-ray pulsations at the radio period of P=0.808s with significance of 7sigma. The pulse shape in the folded X-ray lightcurve (0.15-2 keV) is asymmetric, with statistically significant contributions from up to 5 leading harmonics. Pulse profiles at two different energy ranges differ slightly: the profile is asymmetric at low energies, 0.15-1 keV, while at higher energies, 1-2 keV, it has a nearly sinusodial shape. The radio pulse peak leads the 0.15-2 keV X-ray pulse peak by delta phi = 0.06 +/- 0.03.

We present multi-wavelength observations of the radio magnetar PSR J1622-4950 and its environment. Observations of PSR J1622-4950 with Chandra (in 2007 and 2009) and XMM (in 2011) show that the X-ray flux of PSR J1622-4950 has decreased by a factor of ~50 over 3.7 years, decaying exponentially with a characteristic time of 360 +/- 11 days. This behavior identifies PSR J1622-4950 as a possible addition to the small class of transient magnetars. The X-ray decay likely indicates that PSR J1622-4950 is recovering from an X-ray outburst that occurred earlier in 2007, before the 2007 Chandra observations. Observations with the Australia Telescope Compact Array show strong radio variability, including a possible radio flaring event at least one and a half years after the 2007 X-ray outburst that may be a direct result of this X-ray event. Radio observations with the Molonglo Observatory Synthesis Telescope reveal that PSR J1622-4950 is 8' southeast of a diffuse radio arc, G333.9+0.0, which appears non-thermal in nature and which could possibly be a previously undiscovered supernova remnant. If G333.9+0.0 is a supernova remnant then the estimates of its size and age, combined with the close proximity and reasonable implied velocity of PSR J1622-4950, suggests that these two objects could be physically associated.

We present 35 ks Chandra ACIS observations of the 42 Myr old radio pulsar PSR B1451-68. A point source is detected 0.32" +/- 0.73" from the expected radio pulsar position. It has ~200 counts in the 0.3-8 keV energy range. We identify this point source as the X-ray counterpart of the radio pulsar. PSR B1451-68 is located close to a 2MASS point source, for which we derive 7% as the upper limit on the flux contribution to the measured pulsar X-ray flux. The pulsar spectrum can be described by either a power-law model with photon index Gamma=2.4 (+0.4/-0.3) and a unrealistically high absorbing column density N(H)= (2.5 (+1.2/-1.3)) * 10^(21) cm^-2, or by a combination of a kT=0.35 (+0.12/-0.07) keV blackbody and a Gamma = 1.4 +/- 0.5 power-law component for N(H)[DM]= 2.6 * 10^(20) cm^-2, estimated from the pulsar dispersion measure. At the parallactic, Lutz-Kelker bias corrected distance of 480 pc, the non-thermal X-ray luminosities in the 0.3-8 keV energy band are either Lx(nonth)= (11.3 +/- 1.7) * 10^(29) erg/s or Lx(nonth)= (5.9 (+4.9/-5.0)) * 10^(29) erg/s, respectively. This corresponds to non-thermal X-ray efficiencies of either eta(nonth)= Lx(nonth) / (dE/dt) ~ 0.005 or 0.003, respectively.

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

We report on follow-up observations of candidate X-ray bright, radio-quiet isolated neutron stars (INSs) identified from correlations of the ROSAT All-Sky Survey (RASS) and the Sloan Digital Sky Survey (SDSS) Data Release 4 in Ag\"ueros et al. (2006). We obtained Chandra X-ray Telescope exposures for 13 candidates in order to pinpoint the source of X-ray emission in optically blank RASS error circles. These observations eliminated 12 targets as good INS candidates. We discuss subsequent observations of the remaining candidate with the XMM-Newton X-ray Observatory, the Gemini North Observatory, and the Apache Point Observatory. We identify this object as a likely extragalactic source with an unusually high log(fX/fopt) ~ 2.4. We also use an updated version of the population synthesis models of Popov et al. (2010) to estimate the number of RASS-detected INSs in the SDSS Data Release 7 footprint. We find that these models predict ~3-4 INSs in the 11,000 square deg imaged by SDSS, which is consistent with the number of known INSs that fall within the survey footprint. In addition, our analysis of the four new INS candidates identified by Turner et al. (2010) in the SDSS footprint implies that they are unlikely to be confirmed as INSs; together, these results suggest that new INSs are not likely to be found from further correlations of the RASS and SDSS.

We present X-ray, infrared, optical and radio observations of four previously unidentified Galactic plane X-ray sources, AX J163252-4746, AX J184738-0156, AX J144701-5919 and AX J144547-5931. Detection of each source with the Chandra X-ray Observatory has provided sub-arcsecond localizations, which we use to identify bright infrared counterparts to all four objects. Infrared and optical spectroscopy of these counterparts demonstrate that all four X-ray sources are extremely massive stars, with spectral classifications Ofpe/WN9 (AX J163252-4746), WN7 (AX J184738-0156 = WR121a), WN7-8h (AX J144701-5919) and OIf+ (AX J144547-5931). AX J163252-4746 and AX J184738-0156 are both luminous, hard, X-ray emitters with strong Fe XXV emission lines in their X-ray spectra at ~6.7 keV. The multi-wavelength properties of AX J163252-4746 and AX J184738-0156 are not consistent with isolated massive stars or accretion onto a compact companion; we conclude that their X-ray emission is most likely generated in a colliding-wind binary system. For both AX J144701-5919 and AX J144547-5931, the X-ray emission is an order of magnitude less luminous and with a softer spectrum. These properties are consistent with a colliding-wind binary interpretation for these two sources also, but other mechanisms for the generation of X-rays cannot be excluded. There are many other as yet unidentified X-ray sources in the Galactic plane, with X-ray properties similar to those seen for AX J163252-4746, AX J184738-0156, AX J144701-5919 and AX J144547-5931. This may indicate a substantial population of X-ray-emitting massive stars and colliding-wind binaries in the Milky Way.

We perform population synthesis studies of different types of neutron stars taking into account the magnetic field decay. For the first time, we confront our results with observations using {\it simultaneously} the Log N -- Log S distribution for nearby isolated neutron stars, the Log N -- Log L distribution for magnetars, and the distribution of radio pulsars in the $P$ -- $\dot P$ diagram. We find that our theoretical model is consistent with all sets of data if the initial magnetic field distribution function follows a log-normal law with $<\log (B_0/[G]) > \sim 13.25$ and $\sigma_{\log B_0}\sim 0.6$. The typical scenario includes about 10% of neutron stars born as magnetars, significant magnetic field decay during the first million years of a NS life. Evolutionary links between different subclasses may exist, although robust conclusions are not yet possible. We apply the obtained field distribution and the model of decay to study long-term evolution of neuton stars till the stage of accretion from the interstellar medium. It is shown that though the subsonic propeller stage can be relatively long, initially highly magnetized neutron stars ($B_0 > \sim 10^{13}$ G) reach the accretion regime within the Galactic lifetime if their kick velocities are not too large. The fact that in previous studies made $>$10 years ago, such objects were not considered results in a slight increase of the Accretor fraction in comparison with earlier conclusions. Most of the neutron stars similar to the Magnificent seven are expected to become accreting from the interstellar medium after few billion years of their evolution. They are the main predecestors of accreting isolated neutron stars.

We have previously analyzed sensitive mid-infrared observations to establish that the Pipe Nebula has a very low star-formation efficiency. That study focused on YSOs with excess infrared emission (i.e, protostars and pre-main sequence stars with disks), however, and could have missed a population of more evolved pre-main sequence stars or Class III objects (i.e., young stars with dissipated disks that no longer show excess infrared emission). Evolved pre-main sequence stars are X-ray bright, so we have used ROSAT All-Sky Survey data to search for diskless pre-main sequence stars throughout the Pipe Nebula. We have also analyzed archival XMM-Newton observations of three prominent areas within the Pipe: Barnard 59, containing a known cluster of young stellar objects; Barnard 68, a dense core that has yet to form stars; and the Pipe molecular ring, a high-extinction region in the bowl of the Pipe. We additionally characterize the X-ray properties of YSOs in Barnard 59. The ROSAT and XMM-Newton data provide no indication of a significant population of more evolved pre-main sequence stars within the Pipe, reinforcing our previous measurement of the Pipe's very low star formation efficiency.

Context: Isolated cooling neutron stars with thermal X-ray emission remain rarely detected objects despite many searches investigating the ROSAT data. Aims: We simulate the population of close-by young cooling neutron stars to explain the current observational results. Given the inhomogeneity of the neutron star distribution on the sky it is particularly interesting to identify promising sky regions with respect to on-going and future searches. Methods: Applying a population synthesis model the inhomogeneity of the progenitor distribution and the inhomogeneity of the X-ray absorbing interstellar medium are considered for the first time. The total number of observable neutron stars is derived with respect to ROSAT count rates. In addition, we present sky maps of neutron star locations and discuss age and distance distributions of the simulated neutron stars. Implications for future searches are discussed. Results: With our advanced model we can successfully explain the observed logN - logS distribution of close-by neutron stars. Cooling neutron stars will be most abundant in the directions of rich OB associations. New candidates are expected to be identified behind the Gould Belt, in particular in the Cygnus-Cepheus region. They are expected to be on average younger and then hotter than the known population of isolated cooling neutron stars. In addition, we propose to use data on runaway stars to search for more radio-quiet cooling neutron stars.

We report on submillimetre bolometer observations of the isolated neutron star RX J1856.5--3754 using the LABOCA bolometer array on the Atacama Pathfinder Experiment (APEX) Telescope. No cold dust continuum emission peak at the position of RX J1856.5--3754 was detected. The 3 sigma flux density upper limit of 5 mJy translates into a cold dust mass limit of a few earth masses. We use the new submillimetre limit, together with a previously obtained H-band limit, to constrain the presence of a gaseous, circumpulsar disc. Adopting a simple irradiated-disc model, we obtain a mass accretion limit of dM/dt less than 10^{14} g/s, and a maximum outer disc radius of around 10^{14} cm. By examining the projected proper motion of RX J1856.5--3754, we speculate about a possible encounter of the neutron star with a dense fragment of the CrA molecular cloud a few thousand years ago.

We perform population synthesis studies of different types of neutron stars (thermally emitting isolated neutron stars, normal radio pulsars, magnetars) taking into account the magnetic field decay and using results from the most recent advances in neutron star cooling theory. For the first time, we confront our results with observations using {\it simultaneously} the Log N -- Log S distribution for nearby isolated neutron stars, the Log N -- Log L distribution for magnetars, and the distribution of radio pulsars in the $P$ -- $\dot P$ diagram. For this purpose, we fix a baseline neutron star model (all microphysics input), and other relevant parameters to standard values (velocity distribution, mass spectrum, birth rates ...), allowing to vary the initial magnetic field strength. We find that our theoretical model is consistent with all sets of data if the initial magnetic field distribution function follows a log-normal law with $<\log (B_0/[G])>\sim 13.25$ and $\sigma_{\log B_0}\sim 0.6$. The typical scenario includes about 10% of neutron stars born as magnetars, significant magnetic field decay during the first million years of a NS life (only about a factor of 2 for low field neutron stars but more than an order of magnitude for magnetars), and a mass distribution function dominated by low mass objects. This model explains satisfactorily all known populations. Evolutionary links between different subclasses may exist, although robust conclusions are not yet possible.

The ROSAT All-Sky Survey (RASS) was the first imaging X-ray survey of the entire sky. While X-ray source counterparts are known to range from distant quasars to nearby M dwarfs, the RASS data alone are often insufficient to determine the nature of an X-ray source. As a result, large-scale follow-up programs are required to construct samples of known X-ray emitters. We use optical data produced by the Sloan Digital Sky Survey (SDSS) to identify 709 stellar X-ray emitters cataloged in the RASS and falling within the SDSS Data Release 1 footprint. Most of these are bright stars with coronal X-ray emission unsuitable for SDSS spectroscopy, which is designed for fainter objects (g > 15 mag). Instead, we use SDSS photometry, correlations with the Two Micron All Sky Survey and other catalogs, and spectroscopy from the Apache Point Observatory 3.5 m telescope to identify these stellar X-ray counterparts. Our sample of 707 X-ray-emitting F, G, K, and M stars is one of the largest X-ray-selected samples of such stars. We derive distances to these stars using photometric parallax relations appropriate for dwarfs on the main sequence, and use these distances to calculate LX. We also identify a previously unknown cataclysmic variable (CV) as a RASS counterpart. Separately, we use correlations of the RASS and the SDSS spectroscopic catalogs of CVs and white dwarfs (WDs) to study the properties of these rarer X-ray-emitting stars. We examine the relationship between (fX/fg) and the equivalent width of the Hbeta emission line for 46 X-ray-emitting CVs and discuss tentative classifications for a subset based on these quantities. We identify 17 new X-ray-emitting DA (hydrogen) WDs, of which three are newly identified WDs. (abridged)

(abridged) We study a previously discovered protostellar source that is deeply embedded and drives an energetic molecular outflow. The source, UYSO1, is located close to IRAS 07029-1215 at a distance of ~1 kpc. The multi-wavelength observations resulted in the detection of a double intermediate-mass protostar at the location of UYSO1. In addition to the associated molecular outflow, with a projected size of 0.25 pc, two intersecting near-infrared jets with projected sizes of 0.4 pc and 0.2 pc were found. However, no infrared counterparts to the driving sources could be detected in sensitive near- to far-infrared observations. In interferometric millimeter observations, UYSO1 was resolved into two continuum sources with high column densities and gas masses of 3.5 Mo and 1.2 Mo, with a linear separation of 4200 AU. We report the discovery of a H2O maser toward one of the two sources. The total luminosity is roughly estimated to be ~50 Lo, shared by the two components, one of which is driving the molecular outflow that has a dynamical timescale of less than a few thousand years. The jets of the two individual components are not aligned. Submillimeter observations show that the region lacks typical hot-core chemistry. We thus find two protostellar objects, whose associated circumstellar and parent core masses are high enough to suggest that they may evolve into intermediate-mass stars. This is corroborated by their association with a very massive and energetic CO outflow, suggesting high protostellar accretion rates. The short dynamical timescale of the outflow, the pristine chemical composition of the cloud core and absence of hot core tracers, the absence of detectable radio continuum emission, and the very low protostellar luminosity argue for an extremely early evolutionary stage.

Context: Only two planetary systems around old ms-pulsars are currently known. Young radio pulsars and radio-quiet neutron stars cannot be covered by the usually-applied radio pulse timing technique. However, finding substellar companions around these neutron stars would be of great interest -- not only because of the companion's possible exotic formation but also due to the potential access to neutron star physics. Aims: We investigate the closest young neutron stars to search for substellar companions around them. Methods: Young, thus warm substellar companions are visible in the Near Infrared while the neutron star itself is much fainter. Four young neutron stars are moving fast enough to enable a common proper motion search for substellar companions within few years. Results. For Geminga, RX J0720.4-3125, RX J1856.6-3754, and PSR J1932+1059 we did not find any co-moving companion down to 12, 15, 11, 42 Jupiter masses for assumed ages of 1, 1, 1, 3.1 Myrs and distances of 250, 361, 167, 361 pc, respectively. Near Infrared limits are presented for these four as well as five other neutron stars for which we currently have only observations at one epoch. Conclusions: We conclude that young isolated neutron stars rarely have brown dwarf companions.

We present our new advanced model for population synthesis of close-by cooling NSs. Detailed treatment of the initial spatial distribution of NS progenitors and a detailed ISM structure up to 3 kpc give us an opportunity to discuss the strategy to look for new isolated cooling NSs. Our main results in this respect are the following: new candidates are expected to be identified behind the Gould Belt, in directions to rich OB associations, in particular in the Cygnus-Cepheus region; new candidates, on average, are expected to be hotter than the known population of cooling NS. Besides the usual approach (looking for soft X-ray sources), the search in 'empty' $\gamma$-ray error boxes or among run-away OB stars may yield new X-ray thermally emitting NS candidates.

Multi-wavelength (X-ray to radio) monitoring of Young Stellar Objects (YSOs) can provide important information about physical processes at the stellar surface, in the stellar corona, and/or in the inner circumstellar disk regions. While coronal processes should mainly cause variations in the X-ray and radio bands, accretion processes may be traced by time-correlated variability in the X-ray and optical/infrared bands. Several multi-wavelength studies have been successfully performed for field stars and approx. 1-10 Myr old T Tauri stars, but so far no such study succeeded in detecting simultaneous X-ray to radio variability in extremely young objects like class I and class 0 protostars. Here we present the first simultaneous X-ray, radio, near-infrared, and optical monitoring of YSOs, targeting the Coronet cluster in the Corona Australis star-forming region, which harbors at least one class 0 protostar, several class I objects, numerous T Tauri stars, and a few Herbig AeBe stars. [...] Seven objects are detected simultaneously in the X-ray, radio, and optical/infrared bands; they constitute our core sample. While most of these sources exhibit clear variability in the X-ray regime and several also display optical/infrared variability, none of them shows significant radio variability on the timescales probed. We also do not find any case of clearly time-correlated optical/infrared and X-ray variability. [...] The absence of time-correlated multi-wavelength variability suggests that there is no direct link between the X-ray and optical/infrared emission and supports the notion that accretion is not an important source for the X-ray emission of these YSOs. No significant radio variability was found on timescales of days.

We discuss a set of tests which confront observations of cooling compact objects and theories of their thermal evolution. As an example we apply the recently developed $\mathrm{Log N}$-$\mathrm{Log S}$ test of compact star cooling theories to hybrid stars with a color superconducting quark matter core, we also apply and discuss other existing tests. While there is not yet a microscopically founded superconducting quark matter phase which would fulfill constraints from cooling phenomenology, we explore the hypothetical 2SC+X phase and show that the magnitude and density-dependence of the X-gap can be chosen to satisfy a set of tests: the temperature -- age ($\mathrm{T}$-$\mathrm{t}$) test, the $\mathrm{Log N}$-$\mathrm{Log S}$ test, the brightness constraint, and the mass spectrum constraint. Some recent modifications of the population synthesis model used to obtain the $\mathrm{Log N}$-$\mathrm{Log S}$ distribution are briefly discussed. In addition, we propose to use the age-distance diagram as a new tool to study the local population of young isolated neutron stars.

The Magnificent Seven have all been discovered by their exceptional soft X-ray spectra and high ratios of X-ray to optical flux. They all are considered to be nearby sources. Searching for similar neutron stars with larger distances, one expects larger interstellar absorption resulting in harder X-ray counterparts. Current interstellar absorption treatment depends on chosen abundances and scattering cross-sections of the elements as well as on the 3D distribution of the interstellar medium. After a discussion of these factors we use the comprehensive 3D measurements of the Local Bubble by Lallement et al. 2003 to construct two simple models of the 3D distribution of the hydrogen column density. We test these models by using a set of soft X-ray sources with known distances. Finally, we discuss possible applications for distance estimations and population synthesis studies.

Only seven radio-quiet isolated neutron stars (INSs) emitting thermal X rays are known, a sample that has yet to definitively address such fundamental issues as the equation of state of degenerate neutron matter. We describe a selection algorithm based on a cross-correlation of the ROSAT All-Sky Survey (RASS) and the Sloan Digital Sky Survey (SDSS) that identifies X-ray error circles devoid of plausible optical counterparts to the SDSS g~22 magnitudes limit. We quantitatively characterize these error circles as optically blank; they may host INSs or other similarly exotic X-ray sources such as radio-quiet BL Lacs, obscured AGN, etc. Our search is an order of magnitude more selective than previous searches for optically blank RASS error circles, and excludes the 99.9% of error circles that contain more common X-ray-emitting subclasses. We find 11 candidates, nine of which are new. While our search is designed to find the best INS candidates and not to produce a complete list of INSs in the RASS, it is reassuring that our number of candidates is consistent with predictions from INS population models. Further X-ray observations will obtain pinpoint positions and determine whether these sources are entirely optically blank at g~22, supporting the presence of likely isolated neutron stars and perhaps enabling detailed follow-up studies of neutron star physics.

Our search for the earliest stages of massive star formation turned up twelve massive pre-protocluster candidates plus a few protoclusters. For this search, we selected 47 FIR-bright IRAS sources in the outer Galaxy. We mapped regions of several square arcminutes around the IRAS source in the millimeter continuum in order to find massive cold cloud cores possibly being in a very early stage of massive star formation. Masses and densities are derived for the 128 molecular cloud cores found in the obtained maps. We present these maps together with near-infrared, mid-infrared, and radio data collected from the 2MASS, MSX, and NVSS catalogs. Further data from the literature on detections of high-density tracers, outflows, and masers are added. The multi-wavelength datasets are used to characterize each observed region. The massive cloud cores (M>100 M_sun) are placed in a tentative evolutionary sequence depending on their emission at the investigated wavelengths. Candidates for the youngest stages of massive star formation are identified by the lack of detections in the above-mentioned near-infrared, mid-infrared, and radio surveys. Twelve massive cores prominent in the millimeter continuum fulfill this requirement. Since neither FIR nor radio emission have been detected from these cloud cores massive protostars must be very deeply embedded in these cores. Some of these objects may actually Pre-Proto-cluster cores: an up to now rare object class, where the initial conditions of massive star formation can be studied.

In the course of a comprehensive mm/submm survey of massive star-forming regions, a particularly interesting object has been found in the surroundings of the bright FIR source IRAS 07029-1215, in a distance of 1 kpc. The object -- named UYSO 1 (Unidentified Young Stellar Object 1) -- is cold (T=40 K), it has a massive envelope, and it is associated with an energetic molecular outflow. No infrared point source has been detected at its position for wavelengths below 20 micron. Therefore, it is a very good candidate for a member of the long searched for group of massive protostars.

In the course of a comprehensive mm/submm survey of massive star-forming regions, the vicinities of a sample of 47 luminous IRAS sources were closely investigated with SCUBA and IRAM bolometers in order to search for massive protostellar candidates. A particularly interesting object has been found in the surroundings of the bright FIR source IRAS 07029-1215. Follow-up line observations show that the object is cold, that it has a massive envelope, and that it is associated with an energetic molecular outflow. No infrared point source has been detected at its position. Therefore, it is a very good candidate for a member of the long searched-for group of massive protostars.

We present feasibility studies to directly image stellar surface features, which are caused by magnetic activity, with the Very Large Telescope Interferometer (VLTI). We concentrate on late type magnetically active stars, for which the distribution of starspots on the surface has been inferred from photometric and spectroscopic imaging analysis. The study of the surface spot evolution during consecutive rotation cycles will allow first direct measurements (apart from the Sun) of differential rotation which is the central ingredient of magnetic dynamo processes. The VLTI will provide baselines of up to 200 m, and two scientific instruments for interferometric studies at near- and mid-infrared wavelengths. Imaging capabilities will be made possible by closure-phase techniques. We conclude that a realistically modeled cool surface spot can be detected on stars with angular diameters exceeding ~2 mas using the VLTI with the first generation instrument AMBER. The spot parameters can then be derived with reasonable accuracy. We discuss that the lack of knowledge of magnetically active stars of the required angular size, especially in the southern hemisphere, is a current limitation for VLTI observations of these surface features.