The core component width in normal pulsars, with periods ($P$) $>$ 0.1 seconds, measured at the half-power point at 1 GHz has a lower boundary line (LBL) which closely follows the $P^{-0.5}$ scaling relation. This result is of fundamental importance for understanding the emission process and requires extended studies over a wider frequency range. In this paper we have carried out a detailed study of the profile component widths of 123 normal pulsars observed in the Meterwavelength Single-pulse Polarimetric Emission Survey at 333 and 618 MHz. The components in the pulse profile were separated into core and conal classes. We found that at both frequencies the core as well as the conal component widths versus period had a LBL which followed the $P^{-0.5}$ relation with a similar lower boundary. The radio emission in normal pulsars have been observationally shown to arise from a narrow range of heights around a few hundred kilometers above the stellar surface. In the past the $P^{-0.5}$ relation has been considered as evidence for emission arising from last open dipolar magnetic field lines. We show that the $P^{-0.5}$ dependence only holds if the trailing and leading half-power points of the component are associated with the last open field line. In such a scenario we do not find any physical motivation which can explain the $P^{-0.5}$ dependence for both core and conal components as evidence for dipolar geometry in normal pulsars. We believe the period dependence is a result of an yet unexplained physical phenomenon.

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.

A "pulsar timing array" (PTA), in which observations of a large sample of pulsars spread across the celestial sphere are combined, allows investigation of "global" phenomena such as a background of gravitational waves or instabilities in atomic timescales that produce correlated timing residuals in the pulsars of the array. The Parkes Pulsar Timing Array (PPTA) is an implementation of the PTA concept based on observations with the Parkes 64-m radio telescope. A sample of 20 millisecond pulsars is being observed at three radio-frequency bands, 50cm (~700 MHz), 20cm (~1400 MHz) and 10cm (~3100 MHz), with observations at intervals of 2 - 3 weeks. Regular observations commenced in early 2005. This paper describes the systems used for the PPTA observations and data processing, including calibration and timing analysis. The strategy behind the choice of pulsars, observing parameters and analysis methods is discussed. Results are presented for PPTA data in the three bands taken between 2005 March and 2011 March. For ten of the 20 pulsars, rms timing residuals are less than 1 microsec for the best band after fitting for pulse frequency and its first time derivative. Significant "red" timing noise is detected in about half of the sample. We discuss the implications of these results on future projects including the International Pulsar Timing Array (IPTA) and a PTA based on the Square Kilometre Array. We also present an "extended PPTA" data set that combines PPTA data with earlier Parkes timing data for these pulsars.

This work is a continuation of two previous papers of a series, in which we examined the pulse-width statistics of normal radio pulsars. In the first paper we compiled the largest ever database of pulsars with interpulses in their mean profiles. In the second one we confirmed the existence of the lower boundary in the scatter plot of core component pulse-widths versus pulsar period W50 sim 2.5 P^{-0.5}[deg], first discovered by Rankin using much smaller number of interpulse cases. In this paper we show that the same lower boundary also exists for conal profile components. Rankin proposed a very simple method of estimation of pulsar inclination angle based on comparing the width W50 of its core component with the period dependent value of the lower boundary. We claim that this method can be extended to conal components as well. To explain an existence of the lower boundary Rankin proposed that the core emission originates at or near the polar cap surface. We demonstrated clearly that no coherent pulsar radio emission can originate at altitudes lower than 10 stellar radii, irrespective of the actual mechanism of coherence. We argue that the lower boundary reflects the narrowest angular structures that can be distinguished in the average pulsar beam. These structures represent the core and the conal components in mean pulsar profiles. The P^{-0.5} dependence follows from the dipolar nature of magnetic field lines in the radio emission region, while the numerical factor of about 2.5 deg reflects the curvature radius of a non-dipolar surface magnetic field in the partially screened gap above the polar cap, where dense electron-positron plasma is created. Both core and conal emission should originate at altitudes of about 50 stellar radii in a typical pulsar, with a possibility that the core beam is emitted at a slightly lower heights than the conal ones.

We performed a statistical analysis of half-power pulse-widths of the core components in average pulsar profiles. We confirmed an existence of the lower bound of the distribution of half-power pulse-width versus the pulsar period W50~2.45deg P^(-0.5) found by Rankin (1990). Using our much larger database we found W50= (2.51 +/- 0.08)deg P^(-0.50 +/-0.02) for 21 pulsars with double-pole interpulses for which measurement of the core component width was possible. On the other hand, all single-pole interpulse cases lie in the swarm of pulsars above the boundary line. Using the Monte Carlo simulations based on exact geometrical calculations we found that the Rankin's method of estimation of the inclination angle alpha ~ asin(2.45deg P^(-0.5)/W50) in pulsars with core components is quite good an approximation, except for very small angles alpha in almost aligned rotators.

We performed Monte Carlo simulations of different properties of pulsar radio emission, such as: pulsar periods, pulse-widths, inclination angles and rates of occurrence of interpulse emission (IP). We used recently available large data sets of the pulsar periods P, the pulse profile widths W and the magnetic inclination angle alpha. We also compiled the largest ever database of pulsars with interpulse emission, divided into the double-pole (DP-IP) and the single-pole (SP-IP) cases. Their distribution on the P - Pdot diagram strongly suggests a secular alignment of the magnetic axis from the originally random orientation. We derived possible parent distribution functions of important pulsar parameters by means of the Kolmogorov-Smirnov significance test using the available data sets (P, W, alpha and IP), different models of pulsar radio beam rho = rho(P) as well as different trial distribution functions of pulsar period and the inclination angles. The best suited parent period distribution function is the log-normal distribution, although the gamma function distribution cannot be excluded. The strongest constraint on derived model distribution functions was the requirement that the numbers of interpulses were exactly (within 1sigma errors) at the observed level of occurrences. We found that a suitable model distribution function for the inclination angle is the complicated trigonometric function which has two local maxima, one near 0 deg and the other near 90 deg. The former and the latter implies the right rates of IP occurrence. It is very unlikely that the pulsar beam deviates significantly from the circular cross-section. We found that the upper limit for the average beaming factor fb describing a fraction of the full sphere (called also beaming fraction) covered by a pulsar beam is about 10%. This implies that the number of the neutron stars in the Galaxy might be underestimated.

The Monte Carlo simulations of pulsar periods, pulse-widths and magnetic inclination angles are performed. Using the available observational data sets we study a possible trial parent distribution functions by means of the Kolmogorov-Smirnov significance tests. We also use an additional condition that the numbers of generated interpulses, whether from both magnetic poles or from single pole, are at the observed levels. We conclude that the parent distribution function of magnetic inclination angles is neither flat nor cosine but it is a more complicated function with a local maximum near alpha=25deg and another weaker one near alpha=90deg. The plausible distribution function of pulsar periods is represented by the gamma function. The beaming fraction describing the fraction of observable radio pulsars is about 0.12.