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 study the bi-drifting pulsar B1839-04, where the observed subpulse drift direction in the two leading pulse components is opposite from that in the two trailing components. Such diametrically opposed apparent motions challenge our understanding of an underlying structure. We find that for the geometry spanned by the observer and the pulsar magnetic and rotation axes, the observed bi-drifting in B1839-04 can be reproduced assuming a non-dipolar configuration of the surface magnetic field. Acceptable solutions are found to either have relatively weak $(\sim 10^{12} \,{\rm G})$ or strong $(\sim 10^{14} \,{\rm G})$ surface magnetic fields. Our single pulse modeling shows that a global electric potential variation at the polar cap that leads to a solid-body-like rotation of spark forming regions is favorable in reproducing the observed drift characteristics. This variation of the potential additionally ensures that the variability is identical in all pulse components resulting in the observed phase locking of subpulses. Thorough and more general studies of pulsar geometry show that a low ratio of impact factor to opening angle $(\beta / \rho)$ increases the likelihood of bi-drifting to be observed. We thus conclude that bi-drifting is visible when our line of sight crosses close to the magnetic pole.

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.

The unique and highly unusual drift feature reported for PSR J0815+0939, wherein one component's subpulses drift in the direction opposite of the general trend, is a veritable challenge to pulsar theory. In this paper, we observationally quantify the drift direction throughout its profile, and find that the second component is the only one that exhibits "bi-drifting", meaning that only second component moves in the direction opposite of the others. We here present a model that shows that the observed bi-drifting phenomenon follows from the insight that the discharging regions, i.e. sparks, do not rotate around the magnetic axis per se, but rather around the point of electric potential extremum at the polar cap: minimum in the pulsar case ($\mathbf{\Omega} \cdot \mathbf{B} < 0$) and maximum in the antipulsar case ($\mathbf{\Omega} \cdot \mathbf{B} > 0$). We show that a purely dipolar surface magnetic field cannot exhibit bi-drifting behaviour, though certain non-dipolar configurations can. We can distinguish two types of solutions, with relatively low ($\sim 10^{12} \,{\rm G}$) and high ($\sim 10^{14} \,{\rm G}$) surface magnetic fields. Depending on the strength of the surface magnetic field, the radius of the curvature of magnetic field lines ranges from $10^{5} \,{\rm cm}$ to $10^{7}\,{\rm cm}$. Pulsar \J~allows us to gain an understanding of the polar-cap conditions essential for plasma generation processes in the inner acceleration region, by linking the observed subpulse shift to the underlying spark motion.

PSR J0922+0638 (B0919+06) shows unexplained anomalous variations in the on-pulse phase, where the pulse appears to episodically move to an earlier longitude for a few tens of rotations before reverting to the usual phase for approximately several hundred to more than a thousand rotations. These events, where the pulse moves in phase by up to 5$^{\circ}$, have been previously detected in observations from $\sim$300 to 2000 MHz. We present simultaneous observations from the Effelsberg 100-m radio telescope at 1350 MHz and the Bornim (Potsdam) station of the LOw Frequency ARray at 150 MHz. Our observations present the first evidence for an absence of the anomalous phase-shifting behaviour at 150 MHz. Instead, the observed intensity at the usual pulse-phase typically decreases, often showing a pseudo-nulling feature corresponding to the times when phase shifts are observed at 1350 MHz. The presence of weak emission at the usual pulse-phase supports the theory that these shifts may result from processes similar to the 'profile-absorption' expected to operate for PSR J0814+7429 (B0809+74). A possible mechanism for this could be intrinsic variations of the emission within the pulsar's beam combined with absorption by expanding shells of electrons in the line of sight.

We constrain the X-ray properties of the nearby $(360\,{\rm pc})$, old ($5\,{\rm Myr}$) pulsar B1133+16 with $\sim 100\,{\rm ks}$ effective exposure time by {\it XMM-Newton}. The observed pulsar flux in the 0.2-3 keV energy range is $\sim 10^{-14} \, {\rm erg \, cm}^{-2} \, {\rm s}^{-1}$, which results in the recording of $\sim 600$ source counts with the EPIC pn and MOS detectors. The X-ray radiation is dominated by nonthermal radiation and is well described by both a single power-law model (PL) and a sum of blackbody and power-law emission (BB+PL). The BB+PL model results in a spectral photon index $\Gamma=2.4^{+0.4}_{-0.3}$ and a nonthermal flux in the 0.2-3 keV energy range of $(7\pm 2) \times 10^{-15}\, {\rm erg \, cm}^{-2} \, {\rm s}^{-1}$. The thermal emission is consistent with the blackbody emission from a small hot spot with a radius of $R_{\rm pc} \approx 14^{+7}_{-5} \, {\rm m}$ and a temperature of $T_{\rm s} = 2.9^{+0.6}_{-0.4} \, {\rm MK}$. Assuming that the hot spot corresponds to the polar cap of the pulsar, we can use the magnetic flux conservation law to estimate the magnetic field at the surface $B_{\rm s} \approx 3.9 \times 10^{14} \, {\rm G}$. The observations are in good agreement with the predictions of the partially screened gap model, which assumes the existence of small-scale surface magnetic field structures in the polar cap region.

We report on simultaneous X-ray and radio observations of the radio-mode-switching pulsar PSR B1822-09 with ESA's XMM-Newton and the WSRT, GMRT and Lovell radio telescopes. PSR B1822-09 switches between a radio-bright and radio-quiet mode, and we discovered a relationship between the durations of its modes and a known underlying radio-modulation timescale within the modes. We discovered X-ray (energies 0.2-1.4 keV) pulsations with a broad sinusoidal pulse, slightly lagging the radio main pulse in phase by 0.094 +/- 0.017, with an energy-dependent pulsed fraction varying from ~0.15 at 0.3 keV to ~0.6 at 1 keV. No evidence is found for simultaneous X-ray and radio mode switching. The total X-ray spectrum consists of a cool component (T ~ 0.96 x 10^6 K, hot-spot radius R ~ 2.0 km) and a hot component (T ~ 2.2 x 10^6 K, R ~ 100 m). The hot component can be ascribed to the pulsed emission and the cool component to the unpulsed emission. The high-energy characteristics of PSR B1822-09 resemble those of middle-aged pulsars such as PSR B0656+14, PSR B1055-52 and Geminga, including an indication for pulsed high-energy gamma-ray emission in Fermi LAT data. Explanations for the high pulsed fraction seem to require different temperatures at the two poles of this orthogonal rotator, or magnetic anisotropic beaming effects in its strong magnetic field. In the X-ray skymap we found a harder source at only (5.1+/- 0.5 )arcsec from PSR B1822-09, which might be a pulsar wind nebula.

A large sample of pulsars was observed as part of the Meterwavelength Single-pulse Polarimetric Emission Survey. We carried out a detailed fluctuation spectral analysis which revealed periodic features in 46% pulsars including 22 pulsars where drifting characteristics were reported for the first time. The pulsar population can be categorized into three distinct groups, pulsars which show systematic drift motion within the pulse window, the pulsars showing no systematic drift but periodic amplitude fluctuation and pulsars with no periodic variations. We discovered the dependence of the drifting phenomenon on the spin down energy loss ($\dot{E}$), with the three categories occupying distinctly different regions along the $\dot{E}$ axis. The estimation of the drift periodicity ($P_3$) from the peak frequency in the fluctuation spectra is ambiguous due to the aliasing effect. However, using basic physical arguments we were able to determine $P_3$ in pulsars showing systematic drift motion. The estimated $P_3$ values in these pulsars were anti-correlated with $\dot{E}$ which favoured the Partially Screened Gap model of Inner Acceleration Region in pulsars.

We have conducted the Meterwavelength Single-pulse Polarimetric Emission Survey to study the radio emission properties of normal pulsars. A total of 123 pulsars with periods between 0.1 seconds and 8.5 seconds were observed in the survey at two different frequencies, 105 profiles at 333 MHz, 118 profiles at 618 MHz and 100 pulsars at both. In this work we concentrate primarily on the time-averaged properties of the pulsar emission. The measured widths of the pulsar profiles in our sample usually exhibit the radius to frequency mapping. We validate the existence of lower bounds for the distribution of profile widths with pulsar period ($P$), which is seen for multiple definitions of the width, viz. a lower boundary line (LBL) at $2.7^{\circ} P^{-0.5}$ with width measured at 50% level of profile peak, a LBL at $5.7^{\circ} P^{-0.5}$ for 10% level of peak and LBL at $6.3^{\circ} P^{-0.5}$ for width defined as 5$\sigma$ above the baseline level. In addition we have measured the degree of linear polarization in the average profile of pulsars and confirmed their dependence on pulsar spindown energy loss ($\dot{E}$). The single pulse polarization data show interesting trends with the polarization position angle (PPA) distribution exhibiting the simple rotating vector model for high $\dot{E}$ pulsars while the PPA becomes more complex for medium and low $\dot{E}$ pulsars. The single pulse total intensity data is useful for studying a number of emission properties from pulsars like subpulse drifting, nulling and mode changing which is being explored in separate works.

Magnetars are the most magnetized objects in the known universe. Powered by the magnetic energy, and not by the rotational energy as in the case of radio pulsars, they have long been regarded as a completely different class of neutron stars. The discovery of pulsed radio emission from a few magnetars weakened the idea of a clean separation between magnetars and normal pulsars. We use the partially screened gap (PSG) model to explain radio emission of magnetars. The PSG model requires that the temperature of the polar cap is equal to the so-called critical value, i.e., the temperature at which the thermal ions outflowing from the stellar surface screen the acceleration gap. We show that a magnetar has to fulfill the temperature, power, and visibility conditions in order to emit radio waves. First, in order to form PSG, the residual temperature of the surface has to be lower than the critical value. Second, since the radio emission is powered by the rotational energy, it has to be high enough to enable heating of the polar cap by backstreaming particles to the critical temperature. Finally, the structure of the magnetic field has to be altered by magnetospheric currents in order to widen a radio beam and increase the probability of detection. Our approach allows us to predict whether a magnetar can emit radio waves using only its rotational period, period derivative, and surface temperature in the quiescent mode.

The analysis of X-ray observations suggest an ultrastrong ($B\gtrsim 10^{14} \,{\rm G}$) surface magnetic field at the polar cap of pulsars (Szary, 2013). On the other hand, the temperature of the polar caps is about a few millions Kelvin. Based on these two facts we use the Partially Screened Gap (PSG) model to describe the Inner Acceleration Region (IAR). The PSG model assumes that the temperature of the actual polar cap is equal to the so-called critical value, i.e. the temperature at which the outflow of thermal ions from the surface screens the gap completely. We have found that, depending on the conditions above the polar cap, the generation of high energetic photons in IAR can be caused either by Curvature Radiation (CR) or by Inverse Compton Scattering (ICS). Completely different properties of both processes result in two different scenarios of breaking the acceleration gap: the so-called PSG-off mode for the gap dominated by CR and the PSG-on mode for the gap dominated by ICS. The existence of two different mechanisms of gap breakdown naturally explains the mode-changing phenomenon. Different characteristics of plasma generated in the acceleration region for both modes also explain the pulse nulling phenomenon.

We investigate radio emission efficiency $\xi$ of pulsars and report a near linear inverse correlation between $\xi$ and the spindown power $\dot E$, as well as a near linear correlation between $\xi$ and pulsar age $\tau$. This is a consequence of very weak, if any, dependences of radio luminosity $L$ on pulsar period $P$ and period derivative $\dot{P}$, in contrast to X-ray or $\gamma$-ray emission luminosities. The analysis of radio fluxes suggests that these correlations are not due to a selection effect, but are intrinsic to the pulsar radio emission physics. We have found that, although with a large variance, the radio luminosity of pulsars is $\left<L\right>\approx 10^{29} \,{\rm erg/s}$, regardless of the position in the $P-\dot P$ diagram. Within such a picture, a model-independent statement can be made that the death line of radio pulsars corresponds to an upper limit in the efficiency of radio emission. If we introduce the maximum value for a radio efficiency into Monte Carlo-based population syntheses we can reproduce the observed sample using the random luminosity model. The Kolmogorov-Smirnov test on a synthetic flux distribution shows high probability of reproducing the observed distribution. Our results suggests that the plasma responsible for generating radio emission is produced under similar conditions regardless of pulsar age, dipolar magnetic field strength, and spin-down rate. The magnetic fields near the pulsar surface are likely dominated by crust-anchored magnetic anomalies, which do not significantly differ among pulsars, leading to similar conditions for generating electron-positron pairs necessary to power radio emission.

Despite the fact that pulsars have been observed for almost half a century, many questions have remained unanswered. We use the analysis of X-ray observations in order to study the polar cap region of radio pulsars. The size of the hot spots implies that the magnetic field configuration just above the stellar surface differs significantly from a purely dipole one. We can estimate the surface magnetic field as of the order of $10^{14}\,{\rm G}$. On the other hand, the temperature of the hot spots is about a few million Kelvins. Based on these two facts the Partially Screened Gap (PSG) model was proposed to describe the Inner Acceleration Region (IAR). The PSG model assumes that the temperature of the actual polar cap is equal to the so-called critical value, i.e. the temperature at which the outflow of thermal ions from the surface screens the gap completely. We have found that, depending on the conditions above the polar cap, the generation of high energetic photons in IAR can be caused either by Curvature Radiation (CR) or by Inverse Compton Scattering (ICS). This results in two different scenarios of breaking the acceleration gap: the so-called PSG-off mode for CR-dominated gaps and the PSG-on mode ICS-dominated gaps. The existence of two different mechanisms of gap breakdown naturally explains the mode-changing and the pulse nulling. Furthermore, the mode changes of the IAR may explain the anti-correlation of radio and X-ray emission in very recent observations of PSR B0943+10 (Hermsen et al., 2013). Simultaneous analysis of X-ray and radio properties have allowed to develop a model which explains the drifting subpulse phenomenon. According to this model the drift takes place when the charge density in IAR differs from the Goldreich-Julian co-rotational density. The proposed model allows to verify both the radio drift parameters and X-ray efficiency of the observed pulsars.

Observations of the thermal X-ray emission from old radio pulsars implicate that the size of hot spots is much smaller then the size of the polar cap that follows from the purely dipolar geometry of pulsar magnetic field. Plausible explanation of this phenomena is an assumption that the magnetic field at the stellar surface differs essentially from the purely dipolar field. Using the conservation of the magnetic flux through the area bounded by open magnetic field lines we can estimate the surface magnetic field as of the order of $10^{14}$G. Based on observations that the hot spot temperature is about a few million Kelvins the Partially Screened Gap (PSG) model was proposed which assumes that the temperature of the actual polar cap equals to the so called critical temperature. We discuss correlation between the temperature and corresponding area of the thermal X-ray emission for a number of pulsars. We have found that depending on the conditions in a polar cap region the gap breakdown can be caused either by the Curvature Radiation (CR) or by the Inverse Compton Scattering (ICS). When the gap is dominated by ICS the density of secondary plasma with Lorentz factors $10^{2}-10^{3}$ is at least an order of magnitude higher then in a CR scenario. We believe that two different gap breakdown scenarios can explain the mode-changing phenomenon and in particular the pulse nulling. Measurements of the characteristic spacing between sub-pulses ($P_{2}$) and the period at which a pattern of pulses crosses the pulse window ($P_{3}$) allowed us to determine more strict conditions for avalanche pair production in the PSG.

Observations of the thermal X-ray emission from radio pulsars implicate that the size of hot spots is much smaller then the size of the polar cap that follows from the purely dipolar geometry of pulsar magnetic field. Most plausible explanation of this phenomena is an assumption that the magnetic field at the stellar surface differs essentially from the purely dipolar field. We can determine magnetic field at the surface by the conservation of the magnetic flux through the area bounded by open magnetic field lines. Then the value of the surface magnetic field can be estimated as of the order of $10^{14}$ G. On the other hand observations show that the temperature of the hot spot is about a few million Kelvins. Based on these observations the Partially Screened Gap (PSG) model was proposed which assumes that the temperature of the actual polar cap (hot spot) equals to the so called critical temperature. We discuss correlation between the temperature and corresponding area of the thermal X-ray emission for a number of pulsars. The results of our analysis show that the PSG model is suitable to explain both cases: when the hot spot is smaller and larger then conventional polar cap. We argue that in the second case structure and curvature of field lines allow pair creation in the closed field lines region thus the secondary particles can heat the stellar surface outside the actual polar cap. We have found that the Curvature Radiation (CR) plays dominant role in avalanche pair production in the PSG. We studied dependence of the PSG parameters on the pulsar period, the magnetic field strength and the curvature of field lines.

We propose a model, which naturally explains an unusual thermal X-ray emission, observed from PSR J1119-6127. The model is based on the assumption that the pulsar magnetic field near the stellar surface differs significantly from the pure dipole field. We show that the structure and curvature of the field lines can be of the kind that allows the pair creation in the closed field line region. The created pairs propagate along the closed field lines and heat the stellar surface beyond the local poles. It is demonstrated that such a configuration can be easily realized.