The Parker Solar Probe was launched on 2018 August 12 and completed its second orbit on 2019 June 19 with perihelion of 35.7 solar radii. During this time, the Energetic particle Instrument-Hi (EPI-Hi, one of the two energetic particle instruments comprising the Integrated Science Investigation of the Sun, ISOIS) measured seven proton intensity increases associated with stream interaction regions (SIRs), two of which appear to be occurring in the same region corotating with the Sun. The events are relatively weak, with observed proton spectra extending to only a few MeV and lasting for a few days. The proton spectra are best characterized by power laws with indices ranging from -4.3 to -6.5, generally softer than events associated with SIRs observed at 1 au and beyond. Helium spectra were also obtained with similar indices, allowing He/H abundance ratios to be calculated for each event. We find values of 0.016-0.031, which are consistent with ratios obtained previously for corotating interaction region events with fast solar wind < 600 km s-1. Using the observed solar wind data combined with solar wind simulations, we study the solar wind structures associated with these events and identify additional spacecraft near 1 au appropriately positioned to observe the same structures after some corotation. Examination of the energetic particle observations from these spacecraft yields two events that may correspond to the energetic particle increases seen by EPI-Hi earlier.

A solar energetic particle event was detected by the Integrated Science Investigation of the Sun (ISOIS) instrument suite on Parker Solar Probe (PSP) on 2019 April 4 when the spacecraft was inside of 0.17 au and less than 1 day before its second perihelion, providing an opportunity to study solar particle acceleration and transport unprecedentedly close to the source. The event was very small, with peak 1 MeV proton intensities of ~0.3 particles (cm^2 sr s MeV)^-1, and was undetectable above background levels at energies above 10 MeV or in particle detectors at 1 au. It was strongly anisotropic, with intensities flowing outward from the Sun up to 30 times greater than those flowing inward persisting throughout the event. Temporal association between particle increases and small brightness surges in the extreme-ultraviolet observed by the Solar TErrestrial RElations Observatory, which were also accompanied by type III radio emission seen by the Electromagnetic Fields Investigation on PSP, indicates that the source of this event was an active region nearly 80 degrees east of the nominal PSP magnetic footpoint. This suggests that the field lines expanded over a wide longitudinal range between the active region in the photosphere and the corona.

A series of solar energetic particle (SEP) events were observed at Parker Solar Probe (PSP) by the Integrated Science Investigation of the Sun (\ISOIS) during the period from April 18, 2019 through April 24, 2019. The PSP spacecraft was located near 0.48 au from the Sun on Parker spiral field lines that projected out to 1 au within $\sim 25^\circ$ of near Earth spacecraft. These SEP events, though small compared to historically large SEP events, were amongst the largest observed thus far in the PSP mission and provide critical information about the space environment inside 1 au during SEP events. During this period the Sun released multiple coronal mass ejections (CMEs). One of these CMEs observed was initiated on April 20, 2019 at 01:25 UTC, and the interplanetary CME (ICME) propagated out and passed over the PSP spacecraft. Observations by the Electromagnetic Fields Investigation (FIELDS) show that the magnetic field structure was mostly radial throughout the passage of the compression region and the plasma that followed, indicating that PSP did not directly observe a flux rope internal to the ICME, consistent with the location of PSP on the ICME flank. Analysis using relativistic electrons observed near Earth by the Electron, Proton and Alpha Monitor (EPAM) on the Advanced Composition Explorer (ACE) demonstrates the presence of electron seed populations (40--300 keV) during the events observed. The energy spectrum of the \ISOIS~ observed proton seed population below 1 MeV is close to the limit of possible stationary state plasma distributions out of equilibrium. \ISOIS~ observations reveal the \revise{enhancement} of seed populations during the passage of the ICME, which \revise{likely indicates a key part} of the pre-acceleration process that occurs close to the Sun.

The Integrated Science Investigations of the Sun (IS$\odot$IS) instrument suite includes two Energetic Particle instruments: EPI-Hi, designed to measure ions from ~1-200 MeV/nuc, and EPI-Lo, designed to measure ions from ~20 keV/nuc to ~15 MeV/nuc. We present an analysis of eight energetic proton events observed across the energy range of both instruments during PSP's first two orbits in order to examine their combined energy spectra. Background corrections are applied to help resolve spectral breaks between the two instruments and are shown to be effective. In doing so we demonstrate that, even in the early stages of calibration, IS$\odot$IS is capable of producing reliable spectral observations across broad energy ranges. In addition to making groundbreaking measurements very near the Sun, IS$\odot$IS also characterizes energetic particle populations over a range of heliocentric distances inside 1 au. During the first two orbits, IS$\odot$IS observed energetic particle events from a single corotating interaction region (CIR) at three different distances from the Sun. The events are separated by two Carrington rotations and just 0.11 au in distance, however the relationship shown between proton intensities and proximity of the spacecraft to the source region shows evidence of the importance of transport effects on observations of energetic particles from CIRs. Future IS$\odot$IS observations of similar events over larger distances will help disentangle the effects of CIR-related acceleration and transport. We apply similar spectral analyses to the remaining five events, including four that are likely related to stream interaction regions (SIRs) and one solar energetic particle (SEP) event.

In mid-2012, a GMIR observed by Voyager 2 crossed through the heliosheath and collided with the heliopause, generating a pressure pulse that propagated into the very local interstellar medium. The effects of the transmitted wave were seen by Voyager 1 just 93 days after its own heliopause crossing. The passage of the transient was accompanied by long-lasting decreases in galactic cosmic ray intensities that occurred from ~2012.55 to ~2013.35 and ~2012.91 to ~2013.70 at Voyager 2 and Voyager 1, respectively. Omnidirectional (>20 MeV) proton-dominated measurements from each spacecraft's Cosmic Ray Subsystem reveal a remarkable similarity between these causally-related events, with a correlation coefficient of 91.2% and a time-lag of 130 days. Knowing the locations of the two spacecraft, we use the observed time-delay to calculate the GMIR's average speed through the heliosheath (inside the heliopause) as a function of temperature in the very local interstellar medium. This, combined with particle, field, and plasma observations enables us to infer previously unmeasured properties of the heliosheath, including a range of sound speeds and total effective pressures. For a nominal temperature of ~20,000 K just outside the heliopause, we find a sound speed of 314 (+/-) 32 km/s and total effective pressure of 267 (+/-) 55 fPa inside the heliopause. We compare these results with the Interstellar Boundary Explorer's data-driven models of heliosheath pressures derived from energetic neutral atom fluxes (the globally distributed flux) and present them as additional evidence that the heliosheath's dynamics are driven by suprathermal energetic processes.

We report on an Arecibo 4.5-GHz polarimetric single-pulse survey of the brightest pulsars at high frequency within its sky. The high frequency profiles are accompanied by a collection of both previously published and unpublished high quality 1.4- and 0.33-GHz observations. Here our analyses and discussion primarily involve the average and statistical properties of the 46 pulsar's polarimetric pulse sequences, profile classification and frequency evolution, and polarimetric profiles and peak-occurrence histograms. In most cases both the fractional linear polarization and profile widths decrease with frequency as expected, but there are some exceptions. Similarly, we were able to review and/or extend the profile classifications for this population of pulsars and work out their beaming characteristics quantitatively showing that almost all show properties compatible with the core/double-cone emission beam model. The entirety of these observation's average profiles are accessible for download.

Since crossing the heliopause on August 25, 2012, Voyager 1 observed reductions in galactic cosmic ray count rates caused by a time-varying depletion of particles with pitch angles near 90-deg, while intensities of particles with other pitch angles remain unchanged. Between late 2012 and mid-2017, three large-scale events occurred, lasting from ~100 to ~630 days. Omnidirectional and directional high-energy data from Voyager 1's Cosmic Ray Subsystem are used to report cosmic ray intensity variations. Omnidirectional (greater than ~20 MeV) proton-dominated measurements show up to a 3.8% intensity reduction. Bi-directional (greater than ~70 MeV) proton-dominated measurements taken from various spacecraft orientations provide insight about the depletion region's spatial properties. We characterize the anisotropy as a "notch" in an otherwise uniform pitch-angle distribution of varying depth and width centered about 90 degrees in pitch angle space. The notch averages 22-deg wide and 15% deep - signifying a depletion region that is broad and shallow. There are indications that the anisotropy is formed by a combination of magnetic trapping and cooling downstream of solar-induced transient disturbances in a region that is also likely influenced by the highly compressed fields near the heliopause.

The phenomenon of profile mode switching in pulsars, where the stable average pulse profile changes to another stable state on the timescale of a pulsar's period, remains poorly understood. We sought to understand how pulsars undergo profile mode switching through a comparative analysis of the polarization and geometry of the two different profile modes of PSR B0329+54. The polarization behavior and fitted parameters of the rotation-vector model remain constant between modes, and the emission height remains constant as well. These similarities lend support to a model of pair production in the surface plasma that would cause a change in the available electrons and therefore the differential emission intensity.

Simultaneous observations of PSR B0823+26 with ESA's XMM-Newton, the Giant Metrewave Radio Telescope and international stations of the Low Frequency Array revealed synchronous X-ray/radio switching between a radio-bright (B) mode and a radio-quiet (Q) mode. During the B mode we detected PSR B0823+26 in 0.2$-$2 keV X-rays and discovered pulsed emission with a broad sinusoidal pulse, lagging the radio main pulse by 0.208 $\pm$ 0.012 in phase, with high pulsed fraction of 70$-$80%. During the Q mode PSR B0823+26 was not detected in X-rays (2 $\sigma$ upper limit a factor ~9 below the B-mode flux). The total X-ray spectrum, pulse profile and pulsed fraction can globally be reproduced with a magnetized partially ionized hydrogen atmosphere model with three emission components: a primary small hot spot ($T$$\sim$3.6$\times10^6$ K, $R$$\sim$17 m), a larger cooler concentric ring ($T$$\sim$1.1$\times10^6$ K, $R$$\sim$280 m) and an antipodal hot spot ($T$$\sim$1.1$\times10^6 $ K, $R$$\sim$100 m), for the angle between the rotation axis and line of sight direction $\sim66^\circ$. The latter is in conflict with the radio derived value of $(84\pm0.7)^\circ$. The average X-ray flux within hours-long B-mode intervals varied by a factor $\pm$20%, possibly correlated with variations in the frequency and lengths of short radio nulls or short durations of weak emission. The correlated X-ray/radio moding of PSR B0823+26 is compared with the anti-correlated moding of PSR B0943+10, and the lack of X-ray moding of PSR B1822-09. We speculate that the X-ray/radio switches of PSR B0823+26 are due to variations in the rate of accretion of material from the interstellar medium through which it is passing.

New simultaneous X-ray and radio observations of the archetypal mode-switching pulsar PSR B0943+10 have been carried out with XMM-Newton and the LOFAR, LWA and Arecibo radio telescopes in November 2014. They allowed us to better constrain the X-ray spectral and variability properties of this pulsar and to detect, for the first time, the X-ray pulsations also during the X-ray-fainter mode. The combined timing and spectral analysis indicates that unpulsed non-thermal emission, likely of magnetospheric origin, and pulsed thermal emission from a small polar cap are present during both radio modes and vary in a correlated way.

In this paper we study a set of twelve pulsars that previously had not been characterized. Our timing shows that eleven of them are "normal" isolated pulsars, with rotation periods between 0.22 and 2.65 s, characteristic ages between 0.25 Myr and 0.63 Gyr, and estimated magnetic fields ranging from 0.05 to 3.8x 10^{12} G. The youngest pulsar in our sample, PSR~J0627+0706, is located near the Monoceros supernova remnant (SNR G205.5+0.5), but it is not the pulsar most likely to be associated with it. We also confirmed the existence of a candidate from an early Arecibo survey, PSR~J2053+1718, its subsequent timing and polarimetry are also presented here. It is an isolated pulsar with a spin period of 119 ms, a relatively small magnetic field of 5.8x10^9 G and a characteristic age of 6.7 Gyr; this suggests the pulsar was mildly recycled by accretion from a companion star which became unbound when that companion became a supernova. We report the results of single-pulse and average Arecibo polarimetry at both 327 and 1400 MHz aimed at understanding the basic emission properties and beaming geometry of these pulsars. Three of them (PSRs~J0943+2253, J1935+1159 and J2050+1259) have strong nulls and sporadic radio emission, several others exhibit interpulses (PSRs J0627+0706 and J0927+2345) and one shows regular drifting subpulses (J1404+1159).

The five-component profile of the 2.7-ms pulsar J0337+1715 appears to exhibit the best example to date of a core/double-cone emission-beam structure in a millisecond pulsar (MSP). Moreover, three other MSPs, the Binary Pulsar B1913+16, B1953+29 and J1022+1001, seem to exhibit core/single-cone profiles. These configurations are remarkable and important because it has not been clear whether MSPs and slow pulsars exhibit similar emission-beam configurations, given that they have considerably smaller magnetospheric sizes and magnetic field strengths. MSPs thus provide an extreme context for studying pulsar radio emission. Particle currents along the magnetic polar flux tube connect processes just above the polar cap through the radio-emission region to the light-cylinder and the external environment. In slow pulsars radio-emission heights are typically about 500 km around where the magnetic field is nearly dipolar, and estimates of the physical conditions there point to radiation below the plasma frequency and emission from charged solitons by the curvature process. We are able to estimate emission heights for the four MSPs and carry out a similar estimation of physical conditions in their much lower emission regions. We find strong evidence that MSPs also radiate by curvature emission from charged solitons.

Radio pulsar PSR B1946+35 is a classical example of a core/cone triple pulsar where the observer's line-of-sight cuts the emission beam centrally. In this paper we perform a detailed single-pulse polarimetric analysis of B1946+35 using sensitive Arecibo archival and new observations at 1.4 and 4.6 GHz to re-establish the pulsar's classification wherein a pair of inner conal "outriders" surround a central core component. The new 1.4 GHz observation consisted of a long single pulse sequence of 6678 pulses, and its fluctuation spectral analysis revealed that the pulsar shows a time-varying amplitude modulation, where for a thousand periods or so the spectra have a broad low frequency "red" excess and then at intervals they suddenly exhibit highly periodic longitude-stationary modulation of both the core and conal components for several hundred periods. The fluctuations of the leading conal and the core components are in phase, while those in the trailing conal component in counterphase. These fluctuation properties are consistant with shorter pulse sequence analyses reported in an earlier study by Weltevrede et al. (2006, 2007) as well as in our shorter pulse sequence data sets. We argue that this dual modulation of core and conal emission cannot be understood by a model where subpulse modulation is associated with the plasma {\bf E}$\times${\bf B} drift phenomenon. Rather the effect appears to represent a kind of periodic emission-pattern change over timescales of $\sim$18 s (or 25 pulsar periods), which has not been reported previously for any other pulsar.

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.

In pulsar timing, timing residuals are the differences between the observed times of arrival and the predictions from the timing model. A comprehensive timing model will produce featureless residuals, which are presumably composed of dominating noise and weak physical effects excluded from the timing model (e.g. gravitational waves). In order to apply the optimal statistical methods for detecting the weak gravitational wave signals, we need to know the statistical properties of the noise components in the residuals. In this paper we utilize a variety of non-parametric statistical tests to analyze the whiteness and Gaussianity of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 5-year timing data which are obtained from the Arecibo Observatory and the Green Bank Telescope from 2005 to 2010 (Demorest et al. 2013). We find that most of the data are consistent with white noise; Many data deviate from Gaussianity at different levels, nevertheless, removing outliers in some pulsars will mitigate the deviations.

We report on simultaneous X-ray and radio observations of the mode-switching pulsar PSR B0943+10 obtained with the XMM-Newton satellite and the LOFAR, LWA and Arecibo radio telescopes in November 2014. We confirm the synchronous X-ray/radio switching between a radio-bright (B) and a radio-quiet (Q) mode, in which the X-ray flux is a factor ~2.4 higher than in the B-mode. We discovered X-ray pulsations, with pulsed fraction of 38+/-5% (0.5-2 keV), during the B-mode, and confirm their presence in Q-mode, where the pulsed fraction increases with energy from ~20% up to ~65% at 2 keV. We found marginal evidence for an increase in the X-ray pulsed fraction during B-mode on a timescale of hours. The Q-mode X-ray spectrum requires a fit with a two-component model (either a power-law plus blackbody or the sum of two blackbodies), while the B-mode spectrum is well fit by a single blackbody (a single power-law is rejected). With a maximum likelihood analysis, we found that in Q-mode the pulsed emission has a thermal blackbody spectrum with temperature ~3.4x10^6 K and the unpulsed emission is a power-law with photon index ~2.5, while during B-mode both the pulsed and unpulsed emission can be fit by either a blackbody or a power law with similar values of temperature and photon index. A Chandra image shows no evidence for diffuse X-ray emission. These results support a scenario in which both unpulsed non-thermal emission, likely of magnetospheric origin, and pulsed thermal emission from a small polar cap (~1500 m^2) with a strong non-dipolar magnetic field (~10^{14} G), are present during both radio modes and vary in intensity in a correlated way. This is broadly consistent with the predictions of the partially screened gap model and does not necessarily imply global magnetospheric rearrangements to explain the mode switching.

A quasi-periodicity has been identified in the strange emission shifts in pulsar B1859+07 and possibly B0919+06. These events, first investigated by Rankin, Rodriguez & Wright in 2006, originally appeared disordered or random, but further mapping as well as Fourier analysis has revealed that they occur on a fairly regular basis of approximately 150 rotation periods in B1859+07 and perhaps some 700 in B0919+06. The events-which we now refer to as "swooshes"-are not the result of any known type of mode-changing, but rather we find that they are a uniquely different effect, produced by some mechanism other than any known pulse-modulation phenomenon. Given that we have yet to find another explanation for the swooshes, we have appealed to a last resort for periodicities in astrophysics: orbital dynamics in a binary system. Such putative "companions" would then have semi-major axes comparable to the light cylinder radius for both pulsars. However, in order to resist tidal disruption their densities must be at least some 10$^5$ grams/cm$^3$-therefore white-dwarf cores or something even denser might be indicated.

Radio pulsar B1933+16 is brightest core-radiation dominated pulsar in the Arecibo sky, and here we carry out a comprehensive high resolution polarimetric study of its radiation at both 1.5 and 4.6 GHz. At 1.5 GHz, the polarization is largely compatible with a rotating-vector model with $\alpha$ and $\beta$ values of 125 and --1.2$^{\circ}$, such that the core and conal regions can be identified with the primary and secondary polarization modes and plausibly with the extraordinary and ordinary propagation modes. Polarization modal segregation of profiles shows that the core is comprised of two parts which we associate with later X-mode and earlier O-mode emission. Analysis of the broad microstructures under the core shows that they have similar timescales to those of the largely conal radiation of other pulsars studied earlier. Aberration/retardation analysis was here possible for both the conal and core radiation and showed average physical emission heights of about 200 km for each. Comparison with other core-cone pulsars suggests that the core and conal emission arises from similar heights. Assuming the inner vacuum gap model, we note that at these emission altitudes the frequency of the observed radiation $\nu_{obs}$ is less than the plasma frequency $\nu_p$. We then conclude that the radio emission properties are consistent with the theory of coherent curvature radiation by charged solitons where the condition $\nu_{obs} < \nu_{p}$ is satisfied. However, the differences that exist between core and conal emission in their geometric locations within a pulse, polarization and modulation properties are yet to be understood.

We study the spin-down changes of PSR B1859$+$07 over a period of more than 28 years of radio observation. We identify that the time derivative of the rotational frequency ($\nu$) varies quasi-periodically with a period of $\sim$350 days, switching mainly between two spin-down states. The profile shape of the pulsar is correlated with the $\dot \nu$ variation, producing two slightly different profile shapes corresponding to high- and low-$\dot \nu$ states. In addition to these two normal emission states, we confirm the occasional flare-state of the pulsar, in which the emission appears early in spin phase compared to that of the common normal emission. The profile of the flare-state is significantly different from that of the two normal emission states. The correlation analysis further shows that the flare-state is not directly linked with the $\dot \nu$ changes. With a simple emission beam model, we estimate the emission altitude of the normal emission to be 240~km, and explain the origin of the flare-state as an emission height variation from the leading edge of the beam. We also argue that the emission of these states can be explained with a partially active beam model. In this scenario, the trailing portion of the radio beam is usually active and the normal emission is produced. The flare-state occurs when the leading edge of the beam becomes active while the trailing part is being blocked. This model estimates a fixed emission altitude of 360~km. However, the cause of the flare-state (i.e. the emission height variation, or the time-dependent activity across the radio beam) is not easily explained.

Bright single pulses of many radio pulsars show rapid intensity fluctuations (called microstructure) when observed with time resolutions of tens of microseconds. Here, we report an analysis of Arecibo 59.5 $\mu$sec-resolution polarimetric observations of 11 P-band and 32 L-band pulsars with periods ranging from 150 msec to 3.7 sec. These higher frequency observations forms the most reliable basis for detailed microstructure studies. Close inspection of individual pulses reveals that most pulses exhibit quasiperiodicities with a well-defined periodicity timescale ($P_{\mu}$). While we find some pulses with deeply modulating microstructure, most pulses show low-amplitude modulations on top of broad smooth subpulses features, thereby making it difficult to infer periodicities. We have developed a method for such low-amplitude fluctuations wherein a smooth subpulse envelope is subtracted from each de-noised subpulse; the fluctuating portion of each subpulse is then used to estimate $P_{\mu}$ via autocorrelation analysis. We find that the microstructure timescale $P_{\mu}$ is common across all Stokes parameters of polarized pulsar signals. Moreover, no clear signature of curvature radiation in vacuum in highly resolved microstructures was found. Our analysis further shows strong correlation between $P_\mu$ and the pulsar period $P$. We discuss implications of this result in terms of a coherent radiation model wherein radio emission arises due to formation and acceleration of electron-positron pairs in an inner vacuum gap over magnetic polar cap, and a subpulse corresponds to a series of non-stationary sparking discharges. We argue that in this model, $P_{\mu}$ reflects the temporal modulation of non-stationary plasma flow.

Two entwined problems have remained unresolved since pulsars were discovered nearly 50 years ago: the orientation of their polarized emission relative to the emitting magnetic field and the direction of putative supernova ``kicks' relative to their rotation axes. The rotational orientation of most pulsars can be inferred only from the (``fiducial') polarization angle of their radiation, when their beam points directly at the Earth and the emitting polar fluxtube field is $\parallel$ to the rotation axis. Earlier studies have been unrevealing owing to the admixture of different types of radiation (core and conal, two polarization modes), producing both $\parallel$ or $\perp$ alignments. In this paper we analyze the some 50 pulsars having three characteristics: core radiation beams, reliable absolute polarimetry, and accurate proper motions. The ``fiducial' polarization angle of the core emission, we then find, is usually oriented $\perp$ to the proper-motion direction on the sky. As the primary core emission is polarized $\perp$ to the projected magnetic field in Vela and other pulsars where X-ray imaging reveals the orientation, this shows that the proper motions usually lie $\parallel$ to the rotation axes on the sky. Two key physical consequences then follow: first, to the extent that supernova ``kicks' are responsible for pulsar proper motions, they are mostly $\parallel$ to the rotation axis; and second that most pulsar radiation is heavily processed by the magnetospheric plasma such that the lowest altitude ``parent' core emission is polarized $\perp$ to the emitting field, propagating as the extraordinary (X) mode.

The radio millisecond pulsar J1713+0747 is regarded as one of the highest-precision clocks in the sky, and is regularly timed for the purpose of detecting gravitational waves. The International Pulsar Timing Array collaboration undertook a 24-hour global observation of PSR J1713+0747 in an effort to better quantify sources of timing noise in this pulsar, particularly on intermediate (1 - 24 hr) timescales. We observed the pulsar continuously over 24 hr with the Arecibo, Effelsberg, GMRT, Green Bank, LOFAR, Lovell, Nancay, Parkes, and WSRT radio telescopes. The combined pulse times-of-arrival presented here provide an estimate of what sources of timing noise, excluding DM variations, would be present as compared to an idealized root-N improvement in timing precision, where N is the number of pulses analyzed. In the case of this particular pulsar, we find that intrinsic pulse phase jitter dominates arrival time precision when the S/N of single pulses exceeds unity, as measured using the eight telescopes that observed at L-band/1.4 GHz. We present first results of specific phenomena probed on the unusually long timescale (for a single continuous observing session) of tens of hours, in particular interstellar scintillation, and discuss the degree to which scintillation and profile evolution affect precision timing. This paper presents the data set as a basis for future, deeper studies.

Precursors and postcursors (PPCs) are rare emission components detected in a handful of pulsars that appear beyond the main pulse emission, in some cases far away from it. In this paper we attempt to characterize the PPC emission in relation to the pulsar main pulse geometry. In our analysis we find that PPC components have properties very different from that of outer conal emission. The separation of the PPC components from the main pulse center remains constant with frequency. In addition the beam opening angles corresponding to the separation of PPC components from the pulsar center are much larger than the largest encountered in conal emission. Pulsar radio emission is believed to originate within the magnetic polar flux tubes due to the growth of instabilities in the outflowing relativistic plasma. Observationally, there is strong evidence that the main pulse emission originates at altitudes of about 50 neutron star radii for a canonical pulsar. Currently, the most plausible radio emission model that can explain main pulse emission is the coherent curvature radiation mechanism, wherein relativistic charged solitons are formed in a non-stationary electron-positron-pair plasma. The wider beam opening angles of PPC require the emission to emanate from larger altitudes as compared to the main pulse, if both these components originate by the same emission mechanism. We explore this possibility and find that this emission mechanism is probably inapplicable at the height of the PPC emission. We propose that the PPC emission represents a new type of radiation from pulsars with a mechanism different from that of the main pulse.

Arecibo observations of the conal triple pulsar B1918+19 at 0.327- and 1.4-GHz are used to analyse its subpulse behaviour in detail. We confirm the presence of three distinct drift modes (A,B,C) plus a disordered mode (N) and show that they follow one another in specific cycles. Interpreting the pulsar's profile as resulting from a sightline traverse which cuts across an outer cone and tangentially grazes an inner cone, we demonstrate that the phase modulation of the inner cone is locked to the amplitude modulation of the outer cone in all the drift modes. The 9% nulls are found to be largely confined to the dominant B and N modes, and, in the N mode, create alternating bunches of nulls and emission in a quasi-periodic manner with an averaged fluctuation rate of about 12 rotation periods ($P_1$). We explore the assumption that the apparent drift is the first alias of a faster drift of subbeams equally spaced around the cones. This is shown to imply that all modes A, B and C have a common circulation time of 12 $P_1$ and differ only in the number of subbeams. This timescale is on the same order as predicted by the classic {\bf E}$\times${\bf B} drift and also coincides with the N-mode modulation. We therefore arrive at a picture where the circulation speed remains roughly invariant while the subbeams progressively diminish in number from modes A to B to C, and are then re-established during the N mode. We suggest that aliasing combined with subbeam loss may be responsible for apparently dramatic changes in drift rates in other pulsars.

Pulsars emit low-frequency radio waves through to high-energy gamma-rays that are generated anywhere from the surface out to the edges of the magnetosphere. Detecting correlated mode changes in the multi-wavelength emission is therefore key to understanding the physical relationship between these emission sites. Through simultaneous observations, we have detected synchronous switching in the radio and X-ray emission properties of PSR B0943+10. When the pulsar is in a sustained radio 'bright' mode, the X-rays show only an un-pulsed, non-thermal component. Conversely, when the pulsar is in a radio 'quiet' mode, the X-ray luminosity more than doubles and a 100%-pulsed thermal component is observed along with the non-thermal component. This indicates rapid, global changes to the conditions in the magnetosphere, which challenge all proposed pulsar emission theories.

The research presented here examines an 8-hour observation of pulsar B1822-09,taken by the Giant Metrewave Radio Telescope. B1822-09 has been known to exhibit two stable emission modes, the B-mode, where the precursor (PC) `turns-on', and the Q-mode, which is defined by interpulse (IP) emission. The results of our analysis, of this extremely long observation, have shown that B1822-09 exhibits at least three other emission behaviors that have not been seen before in other similar pulsars or in other observations of B1822-09. These three behaviors can be described as: Q-mode emission with PC emission, B-mode emission with IP emission, and instances where both the PC and IP are `on' when transitioning from one mode to the other. The pulse structure has been found to be more complex than previously thought. The MP has an inner cone/core triple (T) configuration together with a central sightline traverse. The IP is a 15/degr-wide region, that along with the MP originate from an open dipolar field. The PC emission comes from a still unknown source. We argue that the PC emission arises within the same region as the MP, but likely comes from higher in the magnetosphere. Overall, our analyses strongly suggest that mode changes allow information transfer between the two magnetic polar regions and contribute to global magnetospheric changes.

We report analysis of an 8 hr observation of PSR B0943+10 at 325 MHz performed at the Giant Metrewave Radio Telescope (GMRT) in India. B0943+10 is well known for displaying regular sub-pulse drifting and two emission modes. We investigate the modal behavior of B0943+10. By reconstructing an entire 'B' mode from two consecutive 'B' modes, we estimate that the pulsar spends roughly 7.5 hrs in the 'B' mode and about 2.2 hrs in the 'Q' mode, on average. Although the pulsar can switch modes within one pulse, the sub-pulse drift rate changes with a characteristic time of 1.2 hrs over the course of a 'B' mode. Under the subbeam carousel model we find the drift-rate changes are produced by a 10% increase in the average number of subbeams and a 16% increase in the carousel circulation time. We speculate that under the partially screened gap model the increase in circulation time should be related to a small increase in the neutron star surface temperature.

Lyne & Manchester (1988) identified a group of some 50 pulsars they called "partial cones" which they found difficult to classify and interpret. They were notable for their asymmetric average profiles and asymmetric polarization position-angle (PPA) traverses, wherein the steepest gradient (SG) point fell toward one edge of the total intensity profile. Over the last two decades, this population of pulsars has raised cautions regarding the core/cone model of the radio pulsar-emission beam which implies a high degree of order, symmetry and geometric regularity. In this paper we reinvestigate this population "partial cone" pulsars on the basis of new single pulse polarimetric observations of 39 of them, observed with the Giant Meterwave Radio Telescope in India and the Arecibo Observatory in Puerto Rico. These highly sensitive observations help us to establish that most of these "partial cones" exhibit a core/cone structure just as did the "normal" pulsars studied in the earlier papers of this series. In short, we find that many of these "partial cones" are partial in the sense that the emission above different areas of their polar caps can be (highly) asymmetric. However, when studied closely we find that their emission geometries are overall identical to core/double cone structure encountered earlier-that is, with specific conal dimensions scaling as the polar cap size.

We present a single pulse study of pulsar B1944+17, whose non-random nulls dominate nearly 70% of its pulses and usually occur at mode boundaries. When not in the null state, this pulsar displays four bright modes of emission, three of which exhibit drifting subpulses. B1944+17 displays a weak interpulse whose position relative to the main pulse we find to be frequency independent. Its emission is nearly 100% polarized, its polarization-angle traverse is very shallow and opposite in direction to that of the main pulse, and it nulls approximately two-thirds of the time. Geometric modeling indicates that this pulsar is a nearly aligned rotator whose alpha value is hardly 2 degrees--i.e., its magnetic axis is so closely aligned with its rotation axis that its sightline orbit remains within its conal beam. The star's nulls appear to be of two distinct types: those with lengths less than about 8 rotation periods appear to be pseudonulls--that is, produced by "empty" sightline traverses through the conal beam system; whereas the longer nulls appear to represent actual cessations of the pulsar's emission engine.

This paper reports new observations of pulsars B0943+10 and B1822--09 carried out with the Arecibo Observatory (AO) and the Giant Metrewave Radio Telescope (GMRT), respectively. Both stars exhibit two stable emission modes. We report the discovery in B0943+10 of a highly linearly polarized precursor component that occurs primarily in only one mode. This emission feature closely resembles B1822-09's precursor which also occurs brightly in only one mode. B0943+10's other mode is well known for its highly regular drifting subpulses that are apparently produced by a rotating carousel system of 20 beamlets. Similary, B1822-09 exhibits subpulse-modulation behavior only in the mode where its precursor is absent. We survey our 18 hours of B0943+10 observations and find that the sideband-modulation features, from which the carousel-rotation time can be directly determined, occur rarely--less than 5% of the time--but always indicating 20 beamlets. We present an analysis of B1822-09's modal modulation characteristics at 325-MHz and compare them in detail with B0943+10. The pulsar never seems to null, and we find a 43-rotation-period feature in the star's Q mode that modulates the interpulse as well as the conal features in the main pulse. We conclude that B1822-09 must have a nearly orthogonal geometry and that its carousel circulation time is long compared to the modal sub-sequences available in our observations, and the mainpulse/interpulse separation is almost exactly 180 degrees. We conclude the precursors for both stars are incompatible with core-cone emission. We assess the interesting suggestion by Dyks et al. that downward-going radiation produces B1822-09's precursor emission.

We present Arecibo time-aligned, total intensity profiles for 46 pulsars over an unusually wide range of radio frequencies and multi-frequency, polarization-angle density diagrams and/or polarization profiles for 58 pulsars pulsars at some or all of the frequencies 50, 111/130, 430 and 1400 MHz. The frequency-dependent dispersion delay has been removed In order to align the profiles for study of their spectral evolution and wherever possible the profiles of each pulsar are displayed on the same longitude scale. Most of the pulsars within Arecibo's declination range that are sufficiently bright for such spectral or single pulse analysis are included in this survey--and the calibrated pulse sequences are available by web download for further study.

Sensitive Arecibo observations provide an unprecedented ability to detect nulls for an accurate pulse-modulation quelling (PMQ) analysis. We demonstrate that a number of conal pulsars show "periodic nulling" similar to the phenomenon found earlier in pulsar B1133+16.

We present a single-pulse study of the four-component pulsar J1819+1305, whose ``null'' pulses bunch at periodic intervals of around 57 times the rotation period. The emission bursts between the null bunches exhibit characteristic modulations at two shorter periodicities of approximately 6.2 and 3 times the rotation period, the former found largely in the two outer components, and the latter only in the first component. Many bursts commence with bright emission in second component, exhibit positive six-period drift across the full profile width, and end with 3-period modulation in the leading component. The 57-period cycle can be modelled geometrically as a sparsely filled subbeam carousel with nulls appearing whenever our line of sight intersects a circulating empty region. This interpretation is compatible with other recent evidence for periodic, carousel-related nulling and appears to support the physics of a polar-gap emission model for ``drifting'' subpulses, but the subtle structure of the emission bursts defies an easy explanation.

New GMRT observations of the five-component pulsar B1857--26 provide detailed insight into its pulse-sequence modulation phenomena for the first time. The outer conal components exhibit a 7.4-rotation-period, longitude-stationary modulation. Several lines of evidence indicate a carousel circulation time $\P3hat$ of about 147 stellar rotations, characteristic of a pattern with 20 beamlets. The pulsar nulls some 20% of the time, usually for only a single pulse, and these nulls show no discernible order or periodicity. Finally, the pulsar's polarization-angle traverse raises interesting issues: if most of its emission is comprised of a single polarization mode, the full traverse exceeds 180\degr; or if both polarization modes are present, then the leading and the trailing halves of the profiles exhibit two different modes. In either case the rotating vector model fails to fit the polarization-angle traverse of the core component.

In this paper we report multifrequency single pulse polarization observations of the PSR B0329+54 normal mode using the Giant Meterwave Radio Telescope at 325 and 610 MHz and the Effelsberg Observatory at 2695 MHz. Our observations show that towards the central part of the polarization position-angle traverse there is an unusual ``arc''-like structure, which is comprised of a broadband ``kink'' and a frequency-dependent ``spot.'' The features are intimately connected with the intensity dependence of the core component: the stronger emission arrives earlier and its linear polarization is displaced farther along the ``kink''. Moreover, at high intensities, the circular polarization is --/+ antisymmetric; the nearly complete positive circular is characteristic of the weaker, later core subpulses. We find that the ``kink'' emission is associated with the extraordinary (X) propagation mode, and hence propagation effects do not appear capable of producing the core component's broadband, intensity-dependent emission. Rather, the overall evidence points to a largely geometric interpretation in which the ``kink'' provides a rare glimpse of the accelerating cascade or height-dependent amplifier responsible for the core radiation.

Averaged pulse profiles of three nearby pulsars: B1929+10, J0437-4715 and B0950+08 exhibit unusual `double notches'. These W-like looking features consist of two adjacent V-shaped dips that approach each other at increasing observation frequency nuobs roughly at a rate sep \propto nuobs^{-1/2}, where sep is the separation between the notches' minima. We show that basic properties of the notches, namely their W-like look and the rate of their converging can be understood within a narrow class of models of coherent radio emission from pulsars: the free electron maser models based on coherent inverse Compton scattering of parallel oscillations of ambient electric field. The observed properties of the pulsars imply that the Fourier spectrum of the wiggler-like oscillations is narrow and that the broad-band character of the radio emission reflects the width of the electron energy distribution. Such a model provides a natural explanation for the frequency-independent separation between the main pulse and interpulse of B0950+08 as well as for the lack of radius to frequency mapping in the conal-like emission of J0437-4715. The frequency behaviour of the main pulse in the profile of the first radio magnetar XTE J1810-197 can also be explained within this model.

We present a detailed study of the single radio pulses of PSR B0656+14. The emission can be characterized by two separate populations of pulses: bright pulses have a narrow ``spiky'' appearance in contrast to the underlying weaker broad pulses. The shape of the pulse profile requires an unusually long timescale to achieve stability (over 25,000 pulses at 327 MHz) caused by spiky emission. The extreme peak-fluxes of the brightest of these pulses indicates that PSR B0656+14, were it not so near, could only have been discovered as an RRAT source. The strongest bursts represent pulses from the bright end of an extended smooth pulse-energy distribution, which is unlike giant pulses, giant micropulses or the pulses of normal pulsars. Longer observations of the RRATs may reveal that they, like PSR B0656+14, emit weaker emission in addition to the bursts.

We present a detailed study of the single radio pulses of PSR B0656+14, a pulsar also known to be a strong pulsed source of high-energy emission. Despite the extensive studies at high-energy wavelengths, there is little or no published work on its single-pulse behaviour in the radio band. In this report we rectify this omission. We have found that the shape of the pulse profile of PSR B0656+14 requires an unusually long timescale to achieve stability (over 25,000 pulses at 327 MHz). This instability is caused by very bright and narrow pulses with widths and luminosities comparable to those observed for the RRATs. Many pulses are bright enough to qualify as "giant pulses", but are broader than those usually meant by this term. At 327 MHz the brightest pulse was about 116 times brighter than the average pulse. Although the most powerful pulses peak near the centre of the profile, occasional sudden strong pulses are also found on the extreme leading edge of the profile. One of them has a peak flux of about 2000 times the average flux at that pulse longitude. No "break" in the pulse-energy distributions is observed, but nevertheless there is evidence of two separate populations of pulses: bright pulses have a narrow "spiky" appearance consisting of short quasi-periodic bursts of emission with microstructure, in contrast to the underlying weaker broad pulses. Furthermore, the spiky pulses tend to appear in clusters which arise and dissipate over about 10 periods. We demonstrate that the spiky emission builds a narrow and peaked profile, whereas the weak emission produces a broad hump, which is largely responsible for the shoulders in the total emission profiles at both high and low frequencies.

The recently discovered RRAT sources are characterized by very bright radio bursts which, while being periodically related, occur infrequently. We find bursts with the same characteristics for the known pulsar B0656+14. These bursts represent pulses from the bright end of an extended smooth pulse-energy distribution and are shown to be unlike giant pulses, giant micropulses or the pulses of normal pulsars. The extreme peak-fluxes of the brightest of these pulses indicates that PSR B0656+14, were it not so near, could only have been discovered as an RRAT source. Longer observations of the RRATs may reveal that they, like PSR B0656+14, emit weaker emission in addition to the bursts.

Pulsar B0943+10 is well known for its `B' (burst) mode, characterized by accurately drifting subpulses, in contrast to its chaotic `Q' (quiet) mode. Six new Arecibo observations at 327 MHz with durations of 2+ hours each have shed considerable light on the modal dynamics of this pulsar. Of these, three were found to be exclusively `B' mode, and three were discovered to exhibit transitions from the `Q' to the `B' mode. One of these observations has permitted us to determine the circulation time of the subbeam carousel in the `Q' mode for the first time, at some 36.4$\pm$0.9 stellar rotation periods. The onset of the `B' mode is then observed to commence similarly in all three observations. The initial circulation time is about 36 periods and relaxes to nearly 38 periods in a roughly exponential fashion with a characteristic time of some 1.2 hours. This is the longest characteristic time ever found in a mode-switching pulsar. Moreover, just after the `B'-mode onset the pulsar exhibits a symmetrical resolved-double profile form with a somewhat stronger trailing component, but this second component slowly dies away leaving the usual single `B'-mode profile with the longitude of the magnetic axis falling at about its trailing half power point. Thus it would appear that Q-to-B- and B-to-Q-transitions have different characteristic times. Some speculations are given on the nature of this slow modal alternation.

A new single pulse behaviour has been identified in two pulsars, B0919+06 and B1859+07. Normally both starts emit bright subpulses in a region near the trailing edge of heir profile. However, occasionally these stars exhibit "events" wherein the emission longitude gradually decreases by about their profile width, remains in this position for typically tens of pulses and then gradually returns over a few pulses to the usual longitude. The effect bears some resemblance to a profile "mode change", but here the effect is gradual and episodic. When the separated profiles of the normal and "event" emission, they reveal a broad and complex profile structure in each pulsar - but one which can probably be understood correctly in terms of a single conical beam. Possibly the effect entails an extreme example of "absorption"-induced profile asymmetry, as suspected in other pulsars. Alternatively, shifting sources of illumination within the pulsar beam may be responsible.

The paper provides a new analysis of this famous five-component ({\bf M}) pulsar. In addition to the star's core-active ``abnormal'' mode, we find two distinct behaviors within its ``normal'' mode, a ``quiet-normal'' mode with regular 2.8-period subpulse modulation and little or no core activity, and a ``flare-normal'' mode, where the core is regularly bright and a nearly 4-period modulation is maintained. The ``flare-normal'' mode appears to be an intermediate state between the ``quiet normal'' and ``abnormal'' behaviors. Short 5--15-pulse ``flare-normal''-mode ``bursts'' and ``quiet normal'' intervals alternate with each other quasi-periodically, making a cycle some 60--80 pulses in duration. ``Abnormal''-mode intervals are interspersed within this overall cycle, usually persisting for only a few pulses, but occasionally lasting for scores or even many hundreds of pulses. Within subsequences where the core is exceptionally quiet, the pulsar provides a nearly ``textbook'' example of a central PA sightline traverse--with a PA rate measured to be at least some 180$\deg/\deg$. On this basis it is shown that the sightline impact angle $\beta$ must be about 0.25\degr or some 5% of the outer conal beam radius. The star's core component is found to be incomplete, despite the fact that the core's full antisymmetric circularly polarized signature is present. We find that the star's PA traverse is disrupted through the action of orthogonally polarized linear power in the longitude range of the core component. This behavior is modeled to show its effect. Finally, the star provides two well defined fiducial points from which emission-height estimates can be computed, the respective centers of both the linear PA traverse and the zero-crossing point of the circular polarization signature.

Analyses of multiple pulse sequences of B2303+30 reveal that this pulsar has two distinct emission modes. One mode (B) follows a steady even-odd pattern and is more intense. The second mode (Q) is characteristically weak, but has intermittent driftbands with a periodicity of approximately 3 P/cycles, and nulls much more frequently than the B mode. Both modes occur with roughly equal frequency and their profiles have a similar single-humped form with a slight asymmetry. Our observations and analyses strongly suggest that the subpulse driftrates in both modes are linked in a series of cycles, which can be modelled as relaxing oscillations in the underlying circulation rate.

The longitude separation between adjacent drifting subpulses, $P_2$, is roughly constant for many pulsars. It was then perplexing when pulsar B0809+74 was found to exhibit substantial variations in this measure, both with wavelength and with longitude position within the pulse window. We analyze these variations between 40 and 1400 MHz, and we show that they stem primarily from the incoherent superposition of the two orthogonal modes of polarization.

We report here the discovery of a significant source of systematic error in the rotation measure determinations of pulsars. Conventional analysis of high sensitivity polarimetric observations of PSR B2016+28 display variation of the rotation measure of $\pm$15 rad m$^{-2}$ (around the mean value of -34.6 rad m$^{-2}$) across the pulse profile. Analysis of single pulse data shows that this variation is an artifact of the incoherent superposition of quasi-orthogonal polarisation modes along with the frequency dependence of relative strength and/or quasi-orthogonality of the modes. Quasi-orthogonal polarization is common among pulsars, and therefore this effect needs to be taken into account in the interpretation of pulsar rotation measures.

We discuss the discovery of `notch-like' features in the mean pulse profile of the nearby, bright pulsar B0950+08. We compare these low-level features with those previously seen in the pulse profiles of pulsars J0437-4715 and B1929+10. While J0437-4715 is a binary millisecond pulsar and B0950+08 and B1929+10 are isolated, slow pulsars, all three pulsars are nearby and very bright. Furthermore, all three have detectable emission over an unusually wide range of pulse phase. We describe the similar properties of the notch features seen in all three pulsars and discuss possible interpretations.

The purpose of this paper is to suggest how detailed single-pulse observations of ``slow'' radio pulsars may be utilized to construct an empirical model for their emission. It links the observational synthesis developed in a series of papers by Rankin starting in the 1980s to the more recent empirical feedback model of Wright (2003a) by regarding the entire pulsar magnetosphere as a non-steady, non-linear interactive system with a natural built-in delay. It is argued that the enhanced role of the outer gap in such a system indicates an evolutionary link to younger pulsars, in which this region is thought to be highly active, and that pulsar magnetospheres should no longer be seen as being ``driven'' by events on the neutron star's polar cap, but as having more in common with planetary magnetospheres and auroral phenomena.

Both nulling and subpulse drifting are poorly understood phenomena. We probe their mechanisms by investigating how they interact in PSR B0809+74. We find that the subpulse drift is not aliased but directly reflects the actual motion of the subbeams. The carousel-rotation time must then be over 200 seconds, which is much longer than theoretically predicted. The drift pattern after nulls differs from the normal one, and using the absence of aliasing we determine the underlying changes in the subbeam-carousel geometry. We show that after nulls, the subbeam carousel is smaller, suggesting that we look deeper in the pulsar magnetosphere than we do normally. The many striking similarities with emission at higher frequencies, thought to be emitted lower too, confirm this. The emission-height change as well as the striking increase in carousel-rotation time can be explained by a post-null decrease in the polar gap height. This offers a glimpse of the circumstances needed to make the pulsar turn off so dramatically.

The average polarization properties of conal single and double profiles directly reflect the polarization-modal structure of the emission beams which produce them. These average properties require that the circulating subbeam systems which produce conal beams entail paired PPM and SPM emission elements which are offset from each other in both magnetic azimuth and magnetic colatitude. A clear delineation of the modal polarization topology of the conal beam promises to address fundamental questions about the nature and origin of this modal emission--and the modal parity at the outer beam edges is a fact of considerable significance. The different angular dependences of the modal ``beamlets'' suggests that the polarization modes are generated via propagation effects. This argument may prove much stronger if the modal emission is fundamentally only partially polarized. Several theories now promise quantitative comparison with the observations.

In this paper we return to the old problem of conal component-pair widths and profile dimensions. Observationally, we consider a set of 10 pulsars with prominent conal component pairs, for which well measured profiles exist over the largest frequency range now possible. Apart from some tendency to narrow at high frequency, the conal components exhibit almost constant widths. We use all three profile measures, the component separation as well as the outside half-power and 10% widths, to determine conal beam radii, which are the focus of our subsequent analysis. These radii at different frequencies are well fitted by a relationship introduced by Thorsett (1991), but the resulting parameters are highly correlated. Three different types of behavior are found: one group of stars exhibits a continuous variation of beam radius which can be extrapolated down to the stellar surface along the ``last open field lines''; a second group exhibits beam radii which asymptotically approach a minimum high frequency value that is 3--5 times larger; and a third set shows almost no spectral change in beam radius at all. The first two behaviors are associated with outer-cone component pairs; whereas the constant separation appears to reflect inner-cone emission.

We have found that there are two distinct emission modes in PSR B0809+74. Beside its normal and most common mode, the pulsar emits in a significantly different quasi-stable mode after most or possibly all nulls, occasionally for over 100 pulses. In this mode the pulsar is brighter, the subpulse separation is less, the subpulses drift more slowly and the pulse window is shifted towards earlier longitudes. We can now account for several previously unexplained phenomena associated with the nulling-drifting interaction: the unexpected brightness of the first active pulse and the low post-null driftrate. We put forward a new interpretation of the subpulse-position jump over the null, indicating that the speedup time scale of the post-null drifting is much shorter than previously thought. The speedup time scale we find is no longer discrepant with the time scales found for the subpulse-drift slowdown and the emission decay around the null.

Polarization-angle density displays are given for pulsars B0809+74 and B2303+30, which exhibit no evidence of the systematic polarization-angle rotation within individual subpulses previously reported for these two stars. The ``drifting'' subpulses of both pulsars exhibit strikingly linear and circular polarization which appears to reflect the characteristics of two nearly orthogonally polarized emission ``modes''--along which the severe average-profile depolarization that is characteristic of their admixture at comparable overall intensities.

The ``drifting'' subpulses exhibited by some radio pulsars have fascinated both observers and theorists for 30 years, and have been widely regarded as one of the most critical and potentially insightful aspects of their emission. Here, we report on detailed studies of pulsar B0943+10, whose nearly coherent sequences of ``drifting'' subpulses have permitted us to identify their origin as a system of subbeams that appear to circulate around the star's magnetic axis. We introduce several new techniques of analysis, and we find that both the primary and secondary features in the star's fluctuation spectra are aliases of their actual values. We have also developed a method of tracing the underlying pattern responsible for the observed sequences, using a ``cartographic'' transform and its inverse, permitting us to study the characteristics of the polar-cap emission ``map'' and to confirm that such a ``map'' in turn represents the observed sequence. We apply these techniques to the study of three different Arecibo observations. The ``B''-mode sequences are consistent in revealing that the emission pattern consists of 20 subbeams, which rotate around the magnetic axis in about 37 periods or 41 seconds. Even in the ``Q'' mode sequence, we find evidence of a compatible circulation time. The similarity of the subbeam patterns at different radio frequencies strongly suggests that the radiation is produced within a set of columns, which extend from close to the stellar surface up though the emission region and reflect some manner of a ``seeding''phenomenon at their base. The subbeam emission is then tied neither to the stellar surface nor to the field.

The beautiful sequences of ``drifting'' subpulses observed in some radio pulsars have been regarded as among the most salient and potentially instructive characteristics of their emission, not least because they have appeared to represent a system of subbeams in motion within the emission zone of the star. Numerous studies of these ``drift'' sequences have been published, and a model of their generation and motion articulated long ago by Ruderman & Sutherland (1975); but efforts thus far have failed to establish an illuminating connection between the drift phemomenon and the actual sites of radio emission. Through a detailed analysis of a nearly coherent sequence of ``drifting'' pulses from pulsar B0943+10, we have in fact identified a system of subbeams circulating around the magnetic axis of the star. A mapping technique, involving a ``cartographic'' transform and its inverse, permits us to study the character of the polar-cap emission ``map'' and then to confirm that it, in turn, represents the observed pulse sequence. On this basis, we have been able to trace the physical origin of the ``drifting-subpulse'' emission to a stably rotating and remarkably organized configuration of emission columns, in turn traceable possibly to the magnetic polar-cap ``gap'' region envisioned by some theories.