Coherent radio emission in pulsars is excited due to instabilities in a relativistically streaming non-stationary plasma flow, which is generated from sparking discharges in the inner acceleration region (IAR) near the stellar surface. A number of detailed works have shown the IAR to be a partially screened gap (PSG) dominated by non-dipolar magnetic fields with continuous outflow of ions from the surface. The phenomenon of subpulse drifting is expected to originate due to variable $\mathbf{E}\times\mathbf{B}$ drift of the sparks in PSG, where the sparks lag behind corotation velocity of the pulsar. Detailed observations show a wide variety of subpulse drifting behaviour where subpulses in different components of the profile have different phase trajectories. But the drifting periodicity is seen to be constant, within measurement errors, across all components of the profile. Using the concept of sparks lagging behind corotation speed in PSG as well as the different orientations of the surface non-dipolar magnetic fields we have simulated the expected single pulse behaviour in a representative sample of pulsars. Our results show that the different types of drifting phase behaviour can be reproduced using these simple assumptions of spark dynamics in a non-dipolar IAR.

The dynamic activity of the Sun, governed by its cycle of sunspots -- strongly magnetized regions that are observed on its surface -- modulate our solar system space environment creating space weather. Severe space weather leads to disruptions in satellite operations, telecommunications, electric power grids and air-traffic on polar routes. Forecasting the cycle of sunspots, however, has remained a challenging problem. We use reservoir computing -- a model-free, neural--network based machine-learning technique -- to forecast the upcoming solar cycle, sunspot cycle 25. The standard algorithm forecasts that solar cycle 25 is going to last about ten years, the maxima is going to appear in the year 2024 and the maximum number of sunspots is going to be 113 ($\pm15$). We also develop a novel variation of the standard algorithm whose forecasts for duration and peak timing matches that of the standard algorithm, but whose peak amplitude forecast is 124 ($\pm2$) -- within the upper bound of the standard reservoir computing algorithm. We conclude that sunspot cycle 25 is likely to be a weak, lower than average solar cycle, somewhat similar in strength to sunspot cycle 24.

The first stages of planet formation take place in protoplanetary disks that are largely made up of gas. Understanding how the gas affects planetesimals in the protoplanetary disk is therefore essential. In this paper, we discuss whether or not gas flow can erode planetesimals. We estimate how much shear stress is exerted onto the planetesimal surface by the gas as a function of disk and planetesimal properties. To determine whether erosion can take place, we compare this with previous measurements of the critical stress that a pebble-pile planetesimal can withstand before erosion begins. If erosion takes place, we estimate the erosion time of the affected planetesimals. We also illustrate our estimates with two-dimensional numerical simulations of flows around planetesimals using the lattice Boltzmann method. We find that the wall shear stress can overcome the critical stress of planetesimals in an eccentric orbit within the innermost regions of the disk. The high eccentricities needed to reach erosive stresses could be the result of shepherding by migrating planets. We also find that if a planetesimal erodes, it does so on short timescales. For planetesimals residing outside of $1 \ \rm{au}$, we find that they are mainly safe from erosion, even in the case of highly eccentric orbits.

The partially screened vacuum gap model (PSG) for the inner acceleration region in normal radio pulsars, a variant of the pure vacuum gap model, attempts to account for the observed thermal X-ray emission from polar caps and the subpulse drifting timescales. We have used this model to explain the presence of death lines, and extreme location of PSR J2144$-$3933 in the $P-\dot{P}$ diagram. This model requires maintaining the polar cap near a critical temperature and the presence of non-dipolar surface magnetic field to form the inner acceleration region. In the PSG model, thermostatic regulation is achieved by sparking discharges which are a feature of all vacuum gap models. We demonstrate that non-dipolar surface magnetic field reduces polar cap area in PSR J2144$-$3933 such that only one spark can be produced and is sufficient to sustain the critical temperature. This pulsar has a single component profile over a wide frequency range. Single-pulse polarimetric observations and the rotating vector model confirm that the observer's line-of-sight traverses the emission beam centrally. These observations are consistent with a single spark operating within framework of the PSG model leading to single-component emission. Additionally, single-pulse modulations of this pulsar, including lack of subpulse drifting, presence of single-period nulls and microstructure, are compatible with a single spark either in PSG or in general vacuum gap models.

Periodic modulations are seen in normal period pulsars ($P >$ 0.1 sec) over timescales ranging from a few seconds to several minutes. Such modulations have usually been associated with the phenomenon of subpulse drifting. A number of recent studies have shown subpulse drifting to exhibit very specific physical characteristics : i) drifting is seen only in conal components of the pulse profile and is absent in central core emission; ii) drifting pulsars are distributed over a narrow range of spin-down energy loss ($\dot{E}$), where pulsars with $\dot{E} < 2\times$10$^{32}$ erg s$^{-1}$ show this behaviour, iii) drifting periodicity ($P_3$) is anti-correlated with $\dot{E}$, such that pulsars with lower values of $\dot{E}$ tend to have longer $P_3$. These detailed characterisations of drifting behaviour on the other hand also revealed the presence of other distinct periodic modulations, which can be broadly categorised into two types, periodic nulling and periodic amplitude modulation. In contrast to drifting these periodic phenomena are seen across the entire profile in both the core and conal components simultaneously and are not restricted to any specific $\dot{E}$ range. In this work we have assembled an exhaustive list of around 70 pulsars which show such periodic modulations, 22 of which were newly detected using observations from the Giant Meterwave Radio Telescope and the remaining compiled from past publications. The presence of such a significant group in the pulsar population suggests that periodic nulling and periodic amplitude modulations are newly emergent phenomena in pulsars with their physical origin distinct from subpulse drifting.

Multi-wavelength observations of pulsar emission properties are powerful means to constrain their magnetospheric activity and magnetic topology. Usually a star centred magnetic dipole model is invoked to explain the main characteristics of this radiation. However in some particular pulsars where observational constraints exist, such simplified models are unable to predict salient features of their multi-wavelength emission. This paper aims to carefully model the radio and X-ray emission of PSR~J1136+1551 with an off-centred magnetic dipole to reconcile both wavelength measurements. We simultaneously fit the radio pulse profile with its polarization and the thermal X-ray emission from the polar cap hot spots of PSR~J1136+1551. We are able to pin down the parameters of the non-dipolar geometry (which we have assumed to be an offset dipole) and the viewing angle, meanwhile accounting for the time lag between X-rays and radio emission. Our model fits the data if the off-centred magnetic dipole lies about 20\% below the neutron star surface. We also expect very asymmetric polar cap shapes and sizes, implying non antipodal and non identical thermal emission from the hot spots. We conclude that a non-dipolar surface magnetic field is an essential feature to explain the multi-wavelength aspects of PSR~J1136+1551 and other similar pulsars.

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.

PSR J0108$-$1431 is an old pulsar where the X-ray emission is expected to have a thermal component from the polar cap and a non-thermal component from the magnetosphere. Although the phase-integrated spectra are fit best with a single non-thermal component modeled with a power-law (PL) of photon index $\Gamma=2.9$, the X-ray pulse profiles do show the presence of phase-separated thermal and non-thermal components. The spectrum extracted from half the rotational phase away from the X-ray peak fits well with either a single blackbody (BB) or a neutron star atmosphere (NA) model, whereas, the spectrum from the rest of the phase range is dominated by a PL. From Bayesian analysis, the estimated BB area is smaller than the expected polar cap area for a dipolar magnetic field with a probability of 86% whereas the area estimate from the NA model is larger with a probability of 80%. Due to the ambiguity in the thermal emission model, the polar cap area cannot be reliably estimated and hence cannot be used to understand the nature of the surface magnetic field. Instead, we can infer the presence of multipolar magnetic field from the misalignment between the pulsar's thermal X-ray peak and the radio emission peak. For J0108$-$1431, we estimated a phase-offset $\Delta\phi > 0.1$ between the thermal polar cap emission peak and the radio emission peak and argue that this is best explained by the presence of a multipolar surface magnetic field.

We report on a detailed analysis of the radio emission during the different modes of the pulsar J0826+2637 (B0823+26), observed using the Giant Meterwave Radio Telescope at 306-339 MHz observing frequencies. The pulsar profile has a postcursor and interpulse emission in addition to the main pulse. The single pulses showed the presence of nulling, periodic fluctuation in the emission as well as two prominent modes. In addition the pulsar also showed the presence of a null state where no emission was seen for roughly an hour which was immediately followed by a short duration ($\sim$5 minutes) bright state termed the Q-bright state. The nulling varied significantly in the two modes, from a few percent nulls in B-mode to more than 90 percent nulling during the Q-mode. Additionally, the pulsar showed the presence of low level emission in both the interpulse and postcursor components when the main pulse nulled in B-mode. We detected periodic fluctuations in both the main pulse and postcursor during B-mode which were most likely a form of periodic amplitude modulation unrelated to subpulse drifting. We have also detected the appearance of periodicity during the transitions from the null to the burst states in the Q-mode, which was longer than the B-mode modulations. Our analysis further revealed a significant increase in the main pulse and post-cursor intensity during the transition from the Q-mode to the short duration Q-bright mode. On the other hand no commensurate variation was visible in the interpulse intensity.

We report a detailed analysis of the emission behaviour of the five component, core-double cone, pulsar J2006$-$0807 (B2003$-$08). The single pulses revealed the presence of the three major phenomena of subpulse drifting, nulling and mode changing. The pulsar switched between four different emission modes, two of which showed systematic drifting with prominent drift bands, and were classified as modes A and B respectively. The drifting was seen primarily in the conal components and exhibited the rare phenomenon of bi-drifting, where the drift direction in the second component was opposite to the fourth component. This made PSR J2006$-$0807 the only known example where systematic drift bands were seen around a central core emission. The emission showed a gradual decrease in intensity during mode A which stabilised to a relatively constant level in the subsequent mode B. The presence of a low frequency, weak and wide structure in the fluctuation spectra was also seen primarily in the core component during modes A and B. The core component vanished during mode C and was most prominent during the fourth mode D. Both these modes were frequently interspersed with null pulses. No detectable drifting was seen during modes C and D, but the pulsar showed short duration periodic nulling in the core as well as the conal components. In addition to the four emission modes the pulsar also nulled for long durations lasting up to hundred rotation periods.

We solve for waves in an isothermal, stratified medium with a magnetic field that points along a direction perpendicular to that of gravity and varies exponentially in the direction of gravity. For waves propagating along the magnetic field, we calculate approximate dispersion relation as a function of the magnetic field strength. We also find exact solutions for two different cases: (a) waves propagating along the direction of the magnetic field and (b) waves propagating along the direction of gravity. In each of these cases, we find solutions in terms of either confluent hypergeometric functions or Gauss hypergeometric functions depending on whether the ratio of the scale height of the magnetic field over the density scale height is equal to two or not.

We consider the kinematic fluctuation dynamo problem in a flow that is random, white-in-time, with both solenoidal and potential components. This model is a generalization of the well-studied Kazantsev model. If both the solenoidal and potential parts have the same scaling exponent, then, as the compressibility of the flow increases, the growth rate decreases but remains positive. If the scaling exponents for the solenoidal and potential parts differ, in particular if they correspond to typical Kolmogorov and Burgers values, we again find that an increase in compressibility slows down the growth rate but does not turn it off. The slow down is, however, weaker and the critical magnetic Reynolds number is lower than when both the solenoidal and potential components display the Kolmogorov scaling. Intriguingly, we find that there exist cases, when the potential part is smoother than the solenoidal part, for which an increase in compressibility increases the growth rate. We also find that the critical value of the scaling exponent above which a dynamo is seen is unity irrespective of the compressibility. Finally, we realize that the dimension $d = 3$ is special, since for all other values of $d$ the critical exponent is higher and depends on the compressibility.

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.

In this study we propose a classification scheme for the phenomenon of subpulse drifting in pulsars. We have assembled an exhaustive list of pulsars which exhibit subpulse drifting from previously published results as well as recent observations using the Giant Meterwave Radio Telescope. We have estimated detailed phase variations corresponding to the drifting features. Based on phase behaviour the drifting population was classified into four groups : coherent phase-modulated drifting, switching phase-modulated drifting, diffuse phase-modulated drifting and low-mixed phase-modulated drifting. We have re-established the previous assertion that the subpulse drifting is primarily associated with the conal components in pulsar profile. The core components generally do not show the drifting phenomenon. However, in core emission of certain pulsars longer periodic fluctuations are seen, which are similar to periodic nulling, and likely arise due to a different physical phenomenon. In general the nature of the phase variations of the drifting features across the pulsar profile appears to be associated with specific pulsar profile classes, but we also find several examples that show departures from this trend. It has also been claimed in previous works that the spin-down energy loss is anti-correlated with the drifting periodicity. We have verified this dependence using a larger sample of drifting measurements.

A potential resolution for the generation of coherent radio emission in pulsar plasma is the existence of relativistic charge solitons, which are solutions of nonlinear Schr\"{o}dinger equation (NLSE). In an earlier study, Melikidze et al. (2000) investigated the nature of these charge solitons; however, their analysis ignored the effect of nonlinear Landau damping, which is inherent in the derivation of the NLSE in the pulsar pair plasma. In this paper we include the effect of nonlinear Landau damping and obtain solutions of the NLSE by applying a suitable numerical scheme. We find that for reasonable parameters of the cubic nonlinearity and nonlinear Landau damping, soliton-like intense pulses emerge from an initial disordered state of Langmuir waves and subsequently propagate stably over sufficiently long times, during which they are capable of exciting the coherent curvature radiation in pulsars. We emphasize that this emergence of {\em stable} intense solitons from a disordered state does not occur in a purely cubic NLSE; thus, it is {\em caused} by the nonlinear Landau damping.

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.

We report a detailed observational study of the single pulses from the pulsar J1822$-$2256. The pulsar shows the presence of subpulse drifting, nulling as well as multiple emission modes. During these observations the pulsar existed primarily in two modes; mode A with prominent drift bands and mode B which was more disorderly without any clear subpulse drifting. A third mode C was also seen for a short duration with a different drifting periodicity compared to mode A. The nulls were present throughout the observations but were more frequent during the disorderly B mode. The nulling also exhibited periodicity with a clear peak in the fluctuation spectra. Before the transition from mode A to nulling the pulsar switched to a third drifting state with periodicity different from both mode A and C. The diversity seen in the single pulse behaviour of the pulsar J1822$-$2256 provides an unique window into the emission physics.

We report a detailed study of subpulse drifting in four long period pulsars. These pulsars were observed in the Meterwavelength Single-pulse Polarimetric Emission Survey and the presence of phase modulated subpulse drifting was reported in each case. We have carried out longer duration and more sensitive observations lasting 7000-12000 periods, between frequency range of 306 and 339 MHz. The drifting features were characterised in great detail including the phase variations across the pulse window. In two pulsars J0820$-$1350 and J1720$-$2933 the phases changed steadily across the pulse window. The pulsar J1034$-$3224 has five components. The leading component was very weak and was barely detectable in our observations. The four trailing components showed the presence of subpulse drifting. The phase variations changed in alternate components with a reversal in the sign of the gradient. This phenomenon is known as bi-drifting. The pulsar J1555$-$3134 showed the presence of two distinct peak frequencies of comparable strengths in the fluctuation spectrum. The two peaks did not appear to be harmonically related and were most likely a result of different physical processes. Additionally, the long observations enabled us to explore the temporal variations of the drifting features. The subpulse drifting was largely constant with time but small fluctuations around a mean value was seen.

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.

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.

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.

The pulsar radio emission originates from regions below 10% of the light cylinder radius. This requires a mechanism where coherent emission is excited in relativistic pair plasma with frequency $\nu_{cr}$ which is below the plasma frequency $\nu_{\circ}$ i.e. $\nu_{cr} < \nu_{\circ}$. A possible model for the emission mechanism is charged bunches (charged solitons) moving relativistically along the curved open dipolar magnetic field lines capable of exciting coherent curvature radio emission. In this article we review the results from high quality observations in conjunction with theoretical models to unravel the nature of coherent curvature radio emission in pulsars.

A detailed analysis of nulling was conducted for the pulsars studied in the Meterwavelength Single-pulse Polarimetric Emission Survey. We characterized nulling in 36 pulsars including 17 pulsars where the phenomena were reported for the first time. The most dominant nulls lasted for short durations, less than five periods. The longer duration nulls extending to hundreds of periods were also seen in some cases. A careful analysis showed the presence of periodicities in the transition from the null to the burst states in 11 pulsars. In our earlier work fluctuation spectrum analysis showed multiple periodicities in 6 of these 11 pulsars. We demonstrate that the longer periodicity in each case was associated with nulling. The shorter periodicities usually originate due to subpulse drifting. The nulling periodicities were more aligned with the periodic amplitude modulation indicating a possible common origin for both. Most prevalent nulling lasts for a single period and can be potentially explained using random variations affecting the plasma processes in the pulsar magnetosphere. On the other hand, the longer duration nulls require changes in the pair production processes that need an external triggering mechanism for the change. The presence of periodic nulling puts an added constrain on the triggering mechanism which also needs to be periodic.

We study the role of turbulence in magnetic reconnection, within the framework of magneto-hydrodynamics, using three-dimensional direct numerical simulations. For small turbulent intensity we find that the reconnection rate obeys Sweet-Parker scaling. For large enough turbulent intensity reconnection rate departs significantly from Sweet-Parker behaviour, becomes almost a constant as a function of the Lundquist number. We further study energisation of test-particles in the same setup. We find that the speed of the energised particles obeys a Maxwellian distribution, whose variance also obeys Sweet-Parker scaling for small turbulent intensity but depends weakly on the Lundquist number for large turbulent intensity. Furthermore, the variance is found to increase with the strength of the reconnecting magnetic field.

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

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.

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.

Cosmic ray electrons (CREs) are a crucial part of the ISM and are observed via synchrotron emission. While much modelling has been carried out on the CRE distribution and propagation of the Milky Way, little has been done on normal external star-forming galaxies. Recent spectral data from a new generation of radio telescopes enable us to find more robust estimations of the CRE propagation. We model the synchrotron spectral index of M 51 using the time-dependent diffusion energy-loss equation and to compare the model results with the observed spectral index determined from recent low-frequency observations with LOFAR. This is the first time that this model for CRE propagation has been solved for a realistic distribution of CRE sources, which we derive from the observed star formation rate, in an external galaxy. The radial variation of the synchrotron spectral index and scale-length produced by the model are compared to recent LOFAR and older VLA observational data and also to new observations of M 51 at 325MHz obtained with the GMRT. We find that propagation of CREs by diffusion alone is sufficient to reproduce the observed spectral index distribution in M 51. An isotropic diffusion coefficient with a value of $6.6\pm0.2\,\times 10^{28}\mathrm{cm}^2\,\mathrm{s}^{-1}$ is found to fit best and is similar to what is seen in the Milky Way. We estimate an escape time of $11\,\mathrm{Myr}$ from the central galaxy to $88\,\mathrm{Myr}$ in the extended disk. It is found that an energy dependence of the diffusion coefficient is not important for CRE energies in the range $0.01\,\mathrm{GeV}$--$3\,\mathrm{GeV}$. We are able to reproduce the dependence of the observed synchrotron scale-lengths on frequency, with $l \propto \nu^{-1/4}$ in the outer disk and $l \propto \nu^{-1/8}$ in the inner disk.

We consider strongly stratified forced turbulence in a plane-parallel layer with helicity and corresponding large-scale dynamo action in the lower part and non-helical turbulence in the upper. The magnetic field is found to develop strongly concentrated bipolar structures near the surface. They form elongated bands with a sharp interface between opposite polarities. Unlike earlier experiments with imposed magnetic field, the inclusion of rotation does not strongly suppress the formation of these structures. We perform a systematic numerical study of this phenomenon by varying magnetic Reynolds number, scale separation ratio, and Coriolis number. We focus on the formation of a current sheet between bipolar regions where reconnection of oppositely oriented field lines occurs. We determine the reconnection rate by measuring either the inflow velocity in the vicinity of the current sheet or by measuring the electric field in the reconnection region. We demonstrate that for large Lundquist numbers, S>10^3, the reconnection rate is nearly independent of S in agreement with results of recent numerical simulations performed by other groups in simpler settings.

We review a physical model where the high brightness temperature of 10$^{25}-10^{30}$ K observed in pulsar radio emission is explained by coherent curvature radiation excited in the relativistic electron-positron plasma in the pulsar magnetosphere.

Recent work by Mitra et al. (2014) has shown that in strongly stratified forced two-layer turbulence with helicity and corresponding large-scale dynamo action in the lower layer, a magnetic field occurs in the upper layer in the form of sharply bounded bipolar magnetic spots. Here we extend this model to spherical wedge geometry covering the northern hemisphere up to 75{\deg} latitude and an azimuthal extent of 180{\deg}. The kinetic helicity and therefore also the large-scale magnetic field are strongest at low latitudes. For moderately strong stratification, several bipolar spots form that fill eventually the full longitudinal extent. At early times, the polarity of spots reflects the orientation of the underlying azimuthal field, as expected from {\Omega}-shaped flux loops. At late times their tilt changes such that there is a radial field of opposite orientation at different latitudes separated by about 10{\deg}. Our model demonstrates the spontaneous formation of spots of sizes much larger than the pressure scale height. Their tendency to produce filling factors close to unity is argued to be reminiscent of highly active stars. We confirm that strong stratification and strong scale separation are essential ingredients behind magnetic spot formation, which appears to be associated with downflows at larger depths.

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.

We present the measurements of scatter broadening time-scales ($\tau_{sc}$) for 124 pulsars at 327 MHz, using the upgraded Ooty Radio Telescope (ORT). These pulsars lie in the dispersion measure range of 37 $-$ 503 pc cm$^{-3}$ and declination ($\delta$) range of $-$57$^{\circ} < \delta< 60^{\circ}$. New $\tau_{sc}$ estimates for 58 pulsars are presented, increasing the sample of all such measurements by about 40% at 327 MHz. Using all available $\tau_{sc}$ measurements in the literature, we investigate the dependence of $\tau_{sc}$ on dispersion measure. Our measurements, together with previously reported values for $\tau_{sc}$, affirm that the ionized interstellar medium upto 3 kpc is consistent with Kolmogorov spectrum, while it deviates significantly beyond this distance.

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.

We perform direct numerical simulations of the equations of magnetohydrodynamics with external random forcing and in the presence of gravity. The domain is divided into two parts: a lower layer where the forcing is helical and an upper layer where the helicity of the forcing is zero with a smooth transition in between. At early times, a large-scale helical dynamo develops in the bottom layer. At later times the dynamo saturates, but the vertical magnetic field continues to develop and rises to form dynamic bipolar structures at the top, which later disappear and reappear. Some of the structures look similar to $\delta$ spots observed in the Sun. This is the first example of magnetic flux concentrations, owing to strong density stratification, from self-consistent dynamo simulations that generate bipolar, super-equipartition strength, magnetic structures whose energy density can exceeds the turbulent kinetic energy by even a factor of ten.

The pulsar radio emission is generated in the near magnetosphere of the neutron star and it has to propagate through the rest of it to emerge into the interstellar medium. An important issue is whether this propagation affects the planes of polarization of the generated radiation. Observationally, there is a sufficient evidence that the emerging radiation is polarized parallel or perpendicular to the magnetic field line planes that should be associated with the ordinary O and extraordinary X plasma modes respectively, excited by some radiative process. This strongly suggests that the excited X- and O-modes are not affected by the so-called Adiabatic Walking that causes a slow rotation of polarization vectors. In this paper we demonstrate that the conditions for the Adiabatic Walking are not fulfilled within the soliton model of pulsar radio emission, in which the coherent curvature radiation occurs at frequencies much lower than the characteristic plasma frequency. The X-mode propagates freely and observationally represents the primary polarization mode. The O-mode has difficulties in escaping from the pulsar plasma, nonetheless it is sporadically observed as a weaker secondary polarization mode. We discuss a possible scenario under which the O-mode can also escape from the plasma and reach an observer.

Accurate knowledge of time-variation in meridional flow-speed and profile is crucial for estimating a solar cycle's features, which are ultimately responsible for causing space climate variations. However, no consensus has been reached yet about the Sun's meridional circulation pattern observations and theories. By implementing an Ensemble Kalman Filter (EnKF) data assimilation in a Babcock-Leighton solar dynamo model using Data Assimilation Research Testbed (DART) framework, we find that the best reconstruction of time-variation in meridional flow-speed can be obtained when ten or more observations are used with an updating time of 15 days and a $\le 10\%$ observational error. Increasing ensemble-size from 16 to 160 improves reconstruction. Comparison of reconstructed flow-speed with "true-state" reveals that EnKF data assimilation is very powerful for reconstructing meridional flow-speeds and suggests that it can be implemented for reconstructing spatio-temporal patterns of meridional circulation.

We solve for the motion of charged particles in a helical time-periodic ABC (Arnold-Beltrami-Childress) magnetic field. The magnetic field lines of a stationary ABC field with coefficients $A=B=C=1$ are chaotic, and we show that the motion of a charged particle in such a field is also chaotic at late times with positive Lyapunov exponent. We further show that in time-periodic ABC fields, the kinetic energy of a charged particle can increase indefinitely with time. At late times the mean kinetic energy grows as a power law in time with an exponent that approaches unity. For an initial distribution of particles, whose kinetic energy is uniformly distributed within some interval, the PDF of kinetic energy is, at late times, close to a Gaussian but with steeper tails.

In the presence of strong density stratification, turbulence can lead to a large-scale instability of a horizontal magnetic field if its strength is in a suitable range (within a few percent of the turbulent equipartition value). This instability is related to a suppression of the turbulent pressure so that the turbulence contribution to the mean magnetic pressure becomes negative. This results in the excitation of a negative effective magnetic pressure instability (NEMPI). This instability has so far only been studied for an imposed magnetic field. We want to know how NEMPI works when the mean magnetic field is generated self-consistently by an $\alpha^2$ dynamo, whether it is affected by global spherical geometry, and whether it can influence the properties of the dynamo itself. We adopt the mean-field approach which has previously been shown to provide a realistic description of NEMPI in direct numerical simulations. We assume axisymmetry and solve the mean-field equations with the Pencil-Code for an adiabatic stratification at a total density contrast in the radial direction of approximately 4 orders of magnitude. NEMPI is found to work when the dynamo-generated field is about 4% of the equipartition value, which is achieved through strong $\alpha$ quenching. This instability is excited in the top 5% of the outer radius provided the density contrast across this top layer is at least 10. NEMPI is found to occur at lower latitudes when the mean magnetic field is stronger. For weaker fields, NEMPI can make the dynamo oscillatory with poleward migration. NEMPI is a viable mechanism for producing magnetic flux concentrations in a strongly stratified spherical shell in which a magnetic field is generated by a strongly quenched $\alpha$ effect dynamo.

In astrophysics, turbulent diffusion is often used in place of microphysical diffusion to avoid resolving the small scales. However, we expect this approach to break down when time and length scales of the turbulence become comparable with other relevant time and length scales in the system. Turbulent diffusion has previously been applied to the magneto-rotational instability (MRI), but no quantitative comparison of growth rates at different turbulent intensities has been performed. We investigate to what extent turbulent diffusion can be used to model the effects of small-scale turbulence on the kinematic growth rates of the MRI, and how this depends on angular velocity and magnetic field strength. We use direct numerical simulations in three-dimensional shearing boxes with periodic boundary conditions in the spanwise direction and additional random plane-wave volume forcing to drive a turbulent flow at a given length scale. We estimate the turbulent diffusivity using a mixing length formula and compare with results obtained with the test-field method. It turns out that the concept of turbulent diffusion is remarkably accurate in describing the effect of turbulence on the growth rate of the MRI. No noticeable breakdown of turbulent diffusion has been found, even when time and length scales of the turbulence become comparable with those imposed by the MRI itself. On the other hand, quenching of turbulent magnetic diffusivity by the magnetic field is found to be absent. Turbulence reduces the growth rate of the MRI in a way that is the same as microphysical magnetic diffusion.

IC443 is a relatively large Galactic (~45 arcmin) SNR with a high radio surface brightness. It has fine scale structure down to arcsec scales, and so is difficult to image on all angular scales with a single instrument. Here observations of IC443 at 151 MHz made with both the GMRT and the CLFST are combined to give a composite image of IC443 on all scales from >45 arcmin down to ~20 arcsec.

We consider mean-field dynamo models with fluctuating \alpha effect, both with and without shear. The \alpha effect is chosen to be Gaussian white noise with zero mean and given covariance. We show analytically that the mean magnetic field does not grow, but, in an infinitely large domain, the mean-squared magnetic field shows exponential growth of the fastest growing mode at a rate proportional to the shear rate, which agrees with earlier numerical results of Yousef et al (2008) and recent analytical treatment by Heinemann et al (2011) who use a method different from ours. In the absence of shear, an incoherent \alpha^2 dynamo may also be possible. We further show by explicit calculation of the growth rate of third and fourth order moments of the magnetic field that the probability density function of the mean magnetic field generated by this dynamo is non-Gaussian.

As a test bed for the growth of protoplanetary bodies in a turbulent circumstellar disk we examine the fate of a boulder using direct numerical simulations of particle seeded gas flowing around it. We provide an accurate description of the flow by imposing no-slip and non-penetrating boundary conditions on the boulder surface using the immersed boundary method pioneered by Peskin (2002). Advected by the turbulent disk flow, the dust grains collide with the boulder and we compute the probability density function (PDF) of the normal component of the collisional velocity. Through this examination of the statistics of collisional velocities we test the recently developed concept of collisional fusion which provides a physical basis for a range of collisional velocities exhibiting perfect sticking. A boulder can then grow sufficiently rapidly to settle into a Keplerian orbit on disk evolution time scales.

We have examined the cyclotron resonance instability developing in the relativistic out-flowing plasma in the pulsar magnetosphere. The instability condition leads to radio emission in the sub-GHz frequency regime which is likely to be seen as off-pulse emission. Recent studies have shown the presence of off-pulse emission in long period pulsars, and we demonstrate this plasma process to be an energetically viable mechanism.

Using direct numerical simulations, we verify that "flow IV" of Roberts (1972) exhibits dynamo action dominated by horizontally averaged large-scale magnetic field. With the test-field method we compute the turbulent magnetic diffusivity and find that it is negative and overcomes the molecular diffusivity, thus explaining quantitatively the large-scale dynamo for magnetic Reynolds numbers above $\approx8$. As expected for a dynamo of this type, but contrary to $\alpha$-effect dynamos, the two horizontal field components grow independently of each other and have arbitrary amplitude ratios and phase differences. Small length scales of the mean magnetic field are shown to be stabilized by the turbulent magnetic diffusivity becoming positive at larger wavenumbers. Oscillatory decaying or growing solutions have also been found in certain wavenumber intervals and sufficiently large values of the magnetic Reynolds number. For magnetic Reynolds numbers below $\approx0.5$ the turbulent magnetic diffusivity is confirmed to be positive, as expected for all incompressible flows. Earlier claims of a dynamo driven by a modified Taylor-Green flow through negative eddy diffusivity could not be confirmed.

Aims. To investigate the flux density modulation from pulsars and the existence of specific behaviour of modulation index versus frequency. Methods. Several pulsars have been observed with the Effelsberg radio telescope at 8.35 GHz. Their flux density time series have been corrected for interstellar scintillation effects. Results. We present the measurement of modulation indices for 8 pulsars. We confirm the presence of a critical frequency at ~1 GHz for these pulsars (including 3 new ones from this study). We derived intrinsic modulation indices for the resulting flux density time series. Our data analysis revealed strong single pulses detected from 5 pulsars.

Searching for the physical mechanism that can excite the coherent radio emission in pulsars is still an enigmatic problem. A wealth of high quality observations exist, which over the years have been instrumental in putting stringent constraints to pulsar emission models. In this article we will discuss the observational results that strongly suggests that pulsar radio emission is excited by coherent curvature radiation. We will also mention issues that remain to be resolved.

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 surface layers of the Sun are strongly stratified. In the presence of turbulence with a weak mean magnetic field, a large-scale instability resulting in the formation of non-uniform magnetic structures, can be excited over the scale of many turbulent eddies or convection cells. This instability is caused by a negative contribution of turbulence to the effective (mean-field) magnetic pressure and has previously been discussed in connection with the formation of active regions and perhaps sunspots. We want to understand the effects of rotation on this instability in both two and three dimensions. We use mean-field magnetohydrodynamics in a parameter regime in which the properties of the negative effective magnetic pressure instability have previously been found to be in agreement with those of direct numerical simulations. We find that the instability is suppressed already for relatively slow rotation with Coriolis numbers (i.e. inverse Rossby numbers) around 0.2. The suppression is strongest at the equator. In the nonlinear regime, we find traveling wave solutions with propagation in the prograde direction at the equator with additional poleward migration away from the equator. The prograde rotation of the magnetic pattern near the equator is argued to be a possible explanation for the faster rotation speed of magnetic tracers found on the Sun. In the bulk of the domain, kinetic and current helicities are negative in the northern hemisphere and positive in the southern.

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.

In Basu et al. 2011 we reported the detection of Off-pulse emission from two long period pulsars B0525+21 and B2045-16. The pulsars were observed at a single epoch using the 325 MHz frequency band of the Giant Meterwave Radio Telescope (GMRT). In this paper we report a detailed study of the Off-pulse emission from these two pulsars using multiple observations at two different frequencies, 325 MHz and 610 MHz bands of GMRT. We report detection of Off-pulse emission during each observation and based on the scintillation effects and spectral index of Off-pulse emission we conclude a magnetospheric origin. The magnetospheric origin of Off-pulse emission gives rise to various interesting possibilities about its emission mechanism and raises questions about the structure of the magnetosphere.

We present evidence for finite magnetic helicity density in the heliosphere and numerical models thereof, and relate it to the magnetic field properties of the dynamo in the solar convection zone. We use simulations and solar wind data to compute magnetic helicity either directly from the simulations, or indirectly using time series of the skew-symmetric components of the magnetic correlation tensor. We find that the solar dynamo produces negative magnetic helicity at small scales and positive at large scales. However, in the heliosphere these properties are reversed and the magnetic helicity is now positive at small scales and negative at large scales. We explain this by the fact that a negative diffusive magnetic helicity flux corresponds to a positive gradient of magnetic helicity, which leads to a change of sign from negative to positive values at some radius in the northern hemisphere.

We study spontaneous breakdown of chiral symmetry during the nonlinear evolution of the Tayler instability. We start with an initial steady state of zero helicity. Within linearized perturbation calculations, helical perturbations of this initial state have the same growth rate for either sign of helicity. Direct numerical simulations (DNS) of the fully nonlinear equations, however, show that an infinitesimal excess of one sign of helicity in the initial perturbation gives rise to a saturated helical state. We further show that this symmetry breaking can be described by weakly nonlinear finite--amplitude equations with undetermined coefficients which can be deduced solely from symmetry consideration. By fitting solutions of the amplitude equations to data from DNS we further determine the coefficients of the amplitude equations.

We study the radio--FIR correlation between the nonthermal (synchrotron) radio continuum emission at \lambda 90 cm (333 MHz) and the far infrared emission due to cool (~20 K) dust at \lambda 70\mu m in spatially resolved normal galaxies at scales of ~1 kpc. The slope of the radio--FIR correlation significantly differs between the arm and interarm regions. However, this change is not evident at a lower wavelength of \lambda 20 cm (1.4 GHz). We find the slope of the correlation in the arm to be 0.8 \pm 0.12 and we use this to determine the coupling between equipartition magnetic field (B_{eq}) and gas density (\rho_{gas}) as B_{eq} \propto \rho_{gas}^{0.51 \pm 0.12}. This is close to what is predicted by MHD simulations of turbulent ISM, provided the same region produces both the radio and far infrared emission. We argue that at 1 kpc scales this condition is satisfied for radio emission at 1.4 GHz and may not be satisfied at 333 MHz. Change of slope observed in the interarm region could be caused by propagation of low energy (~1.5 GeV) and long lived (~ 10^8 yr) cosmic ray electrons at 333 MHz.

In two-dimensional forced Navier-Stokes turbulence, energy cascades to the largest scales in the system to form a pair of coherent vortices known as the Bose condensate. We show, both numerically and analytically, that the energy condensation saturates and the system reaches a statistically stationary state. The time scale of saturation is inversely proportional to the viscosity and the saturation energy level is determined by both the viscosity and the force. We further show that, without sufficient resolution to resolve the small-scale enstrophy spectrum, numerical simulations can give a spurious result for the saturation energy level. We also find that the movement of the condensate is similar to the motion of an inertial particle with an effective drag force. Furthermore, we show that the profile of the saturated coherent vortices can be described by a Gaussian core with exponential wings.

The negative effective magnetic pressure instability operates on scales encompassing many turbulent eddies and is here discussed in connection with the formation of active regions near the surface layers of the Sun. This instability is related to the negative contribution of turbulence to the mean magnetic pressure that causes the formation of large-scale magnetic structures. For an isothermal layer, direct numerical simulations and mean-field simulations of this phenomenon are shown to agree in many details in that their onset occurs at the same depth. This depth increases with increasing field strength, such that the maximum growth rate of this instability is independent of the field strength, provided the magnetic structures are fully contained within the domain. A linear stability analysis is shown to support this finding. The instability also leads to a redistribution of turbulent intensity and gas pressure that could provide direct observational signatures.

The negative effective magnetic pressure instability discovered recently in direct numerical simulations (DNS) may play a crucial role in the formation of sunspots and active regions in the Sun and stars. This instability is caused by a negative contribution of turbulence to the effective mean Lorentz force (the sum of turbulent and non-turbulent contributions) and results in formation of large-scale inhomogeneous magnetic structures from initial uniform magnetic field. Earlier investigations of this instability in DNS of stably stratified, externally forced, isothermal hydromagnetic turbulence in the regime of large plasma beta are now extended into the regime of larger scale separation ratios where the number of turbulent eddies in the computational domain is about 30. Strong spontaneous formation of large-scale magnetic structures is seen even without performing any spatial averaging. These structures encompass many turbulent eddies. The characteristic time of the instability is comparable to the turbulent diffusion time, L^2/eta_t, where eta_t is the turbulent diffusivity and L is the scale of the domain. DNS are used to confirm that the effective magnetic pressure does indeed become negative for magnetic field strengths below the equipartition field. The dependence of the effective magnetic pressure on the field strength is characterized by fit parameters that seem to show convergence for larger values of the magnetic Reynolds number.

We report the detection of off-pulse emission from two long period pulsars PSR B0525+21 and PSR B2045-16 using the Giant Metrewave Radio Telescope (GMRT). We recorded high time resolution interferometric data at 325 MHz and gated the data offline to separate the on- and off-pulse sections of the pulsar. On imaging the two temporal sections separately, we detected a point source in the off-pulse images of both the pulsars coincident with the pulsar locations in the on-pulse images. The long periods (3.75 and 1.96 s respectively) and low spin-down energies (E ~ 10exp(31) erg s-1) of the two pulsars argue against a pulsar wind nebula (PWN) origin for the off-pulse emission. The derived estimates of the ISM particle density required to drive a PWN do not support such an interpretation. A magnetospheric origin for the off-pulse emission raises questions regarding the location of the emission region, which is an important input into pulsar emission models.

The chiral symmetry breaking properties of the Tayler instability are discussed. Effective amplitude equations are determined in one case. This model has three free parameters that are determined numerically. Comparison with chiral symmetry breaking in biochemistry is made.

We present the first demonstration of the negative effective magnetic pressure instability in direct numerical simulations of stably stratified, externally forced, isothermal hydromagnetic turbulence in the regime of large plasma beta. By the action of this instability, an initially uniform horizontal magnetic field forms flux concentrations whose scale is large compared to the turbulent scale. We further show that the magnetic energy of these large-scale structures is only weakly dependent on the magnetic Reynolds number, provided its value is large enough for the instability to be excited. Our results support earlier mean-field calculations and analytic work which identified this instability. Applications to the formation of active regions in the Sun are discussed.

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.

We report Giant Meterwave Radio Telescope (GMRT) continuum observations of six nearby normal galaxies at 333 MHz. The galaxies are observed with angular resolutions better than ~20" (corresponding to a linear scale of about 0.4 - 1 kpc). These observations are sensitive to all the angular scales of interest, since the resolution of the shortest baseline in GMRT is greater than the angular size of the galaxies. Further, for five of these galaxies we show that at 333 MHz, the mean thermal fraction is less than 5%. Using archival data at about 1 GHz, we estimate the mean thermal fraction to be about 10% at that frequency. We also find that the nonthermal spectral index is generally steeper in regions with low thermal fraction and/or located in the outer parts of the galaxy. In regions of high thermal fraction, the nonthermal spectral index is flatter, and has a narrow distribution peaking at ~ -0.78 with a spread of 0.16, putting stringent constraints on the physical conditions for generation, diffusion and energy losses of cosmic ray electrons at scales of ~ 1 kpc.

Evidence for a parity-breaking nature of the magnetic buoyancy instability in a stably stratified gas is reported. In the absence of rotation, no helicity is produced, but the non-helical state is found to be unstable to small helical perturbations during the development of the instability. The parity-breaking nature of an instability in magnetohydrodynamics appears to be the first of its kind and is similar to chiral symmetry breaking in biochemistry. Applications to the production of mean fields in galaxy clusters are being discussed.

We extend earlier models of turbulent dynamos with an upper, nearly force-free exterior to spherical geometry, and study how flux emerges from lower layers to the upper ones without being driven by magnetic buoyancy. We also study how this affects the possibility of plasmoid ejection. A spherical wedge is used that includes northern and southern hemispheres up to mid-latitudes and a certain range in longitude of the Sun. In radius, we cover both the region that corresponds to the convection zone in the Sun and the immediate exterior up to twice the radius of the Sun. Turbulence is driven with a helical forcing function in the interior, where the sign changes at the equator between the two hemispheres. An oscillatory large-scale dynamo with equatorward migration is found to operate in the turbulence zone. Plasmoid ejections occur in regular intervals, similar to what is seen in earlier Cartesian models. These plasmoid ejections are tentatively associated with coronal mass ejections. The magnetic helicity is found to change sign outside the turbulence zone, which is in agreement with recent findings for the solar wind.

We present a natural framework for studying the persistence problem in two-dimensional fluid turbulence by using the Okubo-Weiss parameter $\Lambda$ to distinguish between vortical and extensional regions. We then use a direct numerical simulation (DNS) of the two-dimensional, incompressible Navier--Stokes equation with Ekman friction to study probability distribution functions (PDFs) of the persistence times of vortical and extensional regions by employing both Eulerian and Lagrangian measurements. We find that, in the Eulerian case, the persistence-time PDFs have exponential tails; by contrast, this PDF for Lagrangian particles, in vortical regions, has a power-law tail with an exponent $\theta=2.9\pm0.2$.

We study the ability of magnetic helicity expulsion to alleviate catastrophic $\alpha$-quenching in mean field dynamos in two--dimensional spherical wedge domains. Motivated by the physical state of the outer regions of the Sun, we consider $\alpha^2\Omega$ mean field models with a dynamical $\alpha$ quenching. We include two mechanisms which have the potential to facilitate helicity expulsion, namely advection by a mean flow ("solar wind") and meridional circulation. We find that a wind alone can prevent catastrophic quenching, with the field saturating at finite amplitude. In certain parameter ranges, the presence of a large-scale meridional circulation can reinforce this alleviation. However, the saturated field strengths are typically below the equipartition field strength. We discuss possible mechanisms that might increase the saturated field.

Magnetic helicity fluxes are investigated in a family of gauges in which the contribution from ideal magnetohydrodynamics takes the form of a purely advective flux. Numerical simulations of magnetohydrodynamic turbulence in this advective gauge family exhibit instabilities triggered by the build-up of unphysical irrotational contributions to the magnetic vector potential. As a remedy, the vector potential is evolved in a numerically well behaved gauge, from which the advective vector potential is obtained by a gauge transformation. In the kinematic regime, the magnetic helicity density evolves similarly to a passive scalar when resistivity is small and turbulent mixing is mild, i.e. when the fluid Reynolds number is not too large. In the dynamical regime, resistive contributions to the magnetic helicity flux in the advective gauge are found to be significant owing to the development of small length scales in the irrotational part of the magnetic vector potential.

We present a unified three-dimensional model of the convection zone and upper atmosphere of the Sun in spherical geometry. In this model, magnetic fields, generated by a helically forced dynamo in the convection zone, emerge without the assistance of magnetic buoyancy. We use an isothermal equation of state with gravity and density stratification. Recurrent plasmoid ejections, which rise through the outer atmosphere, is observed. In addition, the current helicity of the small--scale field is transported outwards and form large structures like magnetic clouds.

A strong toroidal field can exist in form of a magnetic layer in the overshoot region below the solar convection zone. This motivates a more detailed study of the magnetic buoyancy instability with rotation. We calculate the alpha effect due to helical motions caused by a disintegrating magnetic layer in a rotating density-stratified system with angular velocity Omega making an angle theta with the vertical. We also study the dependence of the alpha effect on theta and the strength of the initial magnetic field. We carry out three-dimensional hydromagnetic simulations in Cartesian geometry. A turbulent EMF due to the correlations of the small scale velocity and magnetic field is generated. We use the test-field method to calculate the transport coefficients of the inhomogeneous turbulence produced by the layer. We show that the growth rate of the instability and the twist of the magnetic field vary monotonically with the ratio of thermal conductivity to magnetic diffusivity. The resulting alpha effect is inhomogeneous and increases with the strength of the initial magnetic field. It is thus an example of an "anti-quenched" alpha effect. The alpha effect is nonlocal, requiring around 8--16 Fourier modes to reconstruct the actual EMF based on the actual mean field.

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 direct numerical simulations of the equations of compressible magnetohydrodynamics in a wedge-shaped spherical shell, without shear, but with random helical forcing which has negative (positive) helicity in the northern (southern) hemisphere. We find a large-scale magnetic field that is nearly uniform in the azimuthal direction and approximately antisymmetric about the equator. Furthermore, the large-scale field in each hemisphere oscillates on nearly dynamical time scales with reversals of polarity and equatorward migration. Corresponding mean-field models also show similar migratory oscillations with a frequency that is nearly independent of the magnetic Reynolds number. This mechanism may be relevant for understanding equatorward migration seen in the solar dynamo.

Self-consistent convective dynamo simulations in wedge-shaped spherical shells are presented. Differential rotation is generated by the interaction of convection with rotation. Equatorward acceleration and dynamo action are obtained only for sufficiently rapid rotation. The angular velocity tends to be constant along cylinders. Oscillatory large-scale fields are found to migrate in the poleward direction. Comparison with earlier simulations in full spherical shells and Cartesian domains is made.

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 carry out systematic and high-resolution studies of dynamo action in a shell model for magnetohydrodynamic (MHD) turbulence over wide ranges of the magnetic Prandtl number $Pr_{\rm M}$ and the magnetic Reynolds number $Re_{\rm M}$. Our study suggests that it is natural to think of dynamo onset as a nonequilibrium, first-order phase transition between two different turbulent, but statistically steady, states. The ratio of the magnetic and kinetic energies is a convenient order parameter for this transition. By using this order parameter, we obtain the stability diagram (or nonequilibrium phase diagram) for dynamo formation in our MHD shell model in the $(Pr^{-1}_{\rm M}, Re_{\rm M})$ plane. The dynamo boundary, which separates dynamo and no-dynamo regions, appears to have a fractal character. We obtain hysteretic behavior of the order parameter across this boundary and suggestions of nucleation-type phenomena.

We use direct numerical simulations of forced MHD turbulence with a forcing function that produces two different signs of kinetic helicity in the upper and lower parts of the domain. We show that the mean flux of magnetic helicity from the small-scale field between the two parts of the domain can be described by a Fickian diffusion law with a diffusion coefficient that is approximately independent of the magnetic Reynolds number and about one third of the estimated turbulent magnetic diffusivity. The data suggest that the turbulent diffusive magnetic helicity flux can only be expected to alleviate catastrophic quenching at Reynolds numbers of more than several thousands. We further calculate the magnetic helicity density and its flux in the domain for three different gauges. We consider the Weyl gauge, in which the electrostatic potential vanishes, the pseudo-Lorenz gauge, where the speed of light is replaced by the sound speed, and the `resistive gauge' in which the Laplacian of the magnetic vector potential acts as resistive term. We find that, in the statistically steady state, the time-averaged magnetic helicity density and the magnetic helicity flux are the same in all three gauges.

We use three-dimensional direct numerical simulations of the helically forced magnetohydrodynamic equations in spherical shell segments in order to study the effects of changes in the geometrical shape and size of the domain on the growth and saturation of large-scale magnetic fields. We inject kinetic energy along with kinetic helicity in spherical domains via helical forcing using Chandrasekhar-Kendall functions. We take perfect conductor boundary conditions for the magnetic field to ensure that no magnetic helicity escapes the domain boundaries. We find dynamo action giving rise to magnetic fields at scales larger than the characteristic scale of the forcing. The magnetic energy exceeds the kinetic energy over dissipative time scales, similar to that seen earlier in Cartesian simulations in periodic boxes. As we increase the size of the domain in the azimuthal direction we find that the nonlinearly saturated magnetic field organizes itself in long-lived cellular structures with aspect ratios close to unity. These structures tile the domain along the azimuthal direction, thus resulting in very small longitudinally averaged magnetic fields for large domain sizes. The scales of these structures are determined by the smallest scales of the domain, which in our simulations is usually the radial scale. We also find that increasing the meridional extent of the domains produces little qualitative change, except a marginal increase in the large-scale field. We obtain qualitatively similar results in Cartesian domains with similar aspect ratios.

Forty years have passed since the discovery of pulsars, yet the physical mechanism of their coherent radio emission is a mystery. Recent observational and theoretical studies strongly suggest that the radiation outcoming from the pulsar magnetosphere consists mainly of extraordinary waves polarized perpendicular to the planes of pulsar dipolar magnetic field. However, the fundamental question whether these waves are excited by maser or coherent curvature radiation, remains open. High quality single pulse polarimetry is required to distinguish between these two possible mechanisms. Here we showcase such {\it decisive} strong single pulses from 10 pulsars observed with the GMRT, showing extremely high linear polarization with the position angle following locally the mean position angle traverse. These pulses, which are relatively free from depolarization, must consist of exclusively single polarization mode. We associate this mode with the extraordinary wave excited by the coherent curvature radiation. This crucial observational signature enables us to argue, for the first time, in favor of the coherent curvature emission mechanism, excluding the maser mechanism.

The magnetospheric locations of pulsar radio emission region are not well known. The actual form of the so--called radius--to--frequency mapping should be reflected in the aberration--retardation (A/R) effects that shift and/or delay the photons depending on the emission height in the magnetosphere. Recent studies suggest that in a handful of pulsars the A/R effect can be discerned w.r.t the peak of the central core emission region. To verify these effects in an ensemble of pulsars we launched a project analysing multi--frequency total intensity pulsar profiles obtained from the new observations from the Giant Meterwave Radio Telescope (GMRT), Arecibo Observatory (AO) and archival European Pulsar Network (EPN) data. For all these profiles we measure the shift of the outer cone components with respect to the core component which is necessary for establishing the A/R effect. Within our sample of 23 pulsars 7 show the A/R effects, 12 of them (doubtful cases) show a tendency towards this effect, while the remaining 4 are obvious counter examples. The counter--examples and doubtful cases may arise from uncertainties in determination of the location of the meridional plane and/or the core emission component. It hence appears that the A/R effects are likely to operate in most pulsars from our sample. We conclude that in cases where those effects are present the core emission has to originate below the conal emission region.

Simulations of stochastically forced shear-flow turbulence in a shearing-periodic domain are used to study the spontaneous generation of large-scale flow patterns in the direction perpendicular to the plane of the shear. Based on an analysis of the resulting large-scale velocity correlations it is argued that the mechanism behind this phenomenon could be the mean-vorticity dynamo effect pioneered by Elperin, Kleeorin, and Rogachevskii in 2003 (Phys. Rev. E 68, 016311). This effect is based on the anisotropy of the eddy viscosity tensor. One of its components may be able to replenish cross-stream mean flows by acting upon the streamwise component of the mean flow. Shear, in turn, closes the loop by acting upon the cross-stream mean flow to produce stronger streamwise mean flows. The diagonal component of the eddy viscosity is found to be of the order of the rms turbulent velocity divided by the wavenumber of the energy-carrying eddies.

We study the dependence of turbulent transport coefficients, such as the components of the $\alpha$ tensor ($\alpha_{ij}$) and the turbulent magnetic diffusivity tensor ($\eta_{ij}$), on shear and magnetic Reynolds number in the presence of helical forcing. We use three-dimensional direct numerical simulations with periodic boundary conditions and measure the turbulent transport coefficients using the kinematic test field method. In all cases the magnetic Prandtl number is taken as unity. We find that with increasing shear the diagonal components of $\alpha_{ij}$ quench, whereas those of $\eta_{ij}$ increase. The antisymmetric parts of both tensors increase with increasing shear. We also propose a simple expression for the turbulent pumping velocity (or $\gamma$ effect). This pumping velocity is proportional to the kinetic helicity of the turbulence and the vorticity of the mean flow. For negative helicity, i.e. for a positive trace of $\alpha_{ij}$, it points in the direction of the mean vorticity, i.e. perpendicular to the plane of the shear flow. Our simulations support this expression for low shear and magnetic Reynolds number. The transport coefficients depend on the wavenumber of the mean flow in a Lorentzian fashion, just as for non-shearing turbulence.

We have observed a total of 67 pulsars at five frequencies ranging from 243 to 3100 MHz. Observations at the lower frequencies were made at the Giant Metre Wave Telescope in India and those at higher frequencies at the Parkes Telescope in Australia. We present profiles from 34 of the sample with the best signal to noise ratio and the least scattering. The general `rules' of pulsar profiles are seen in the data; profiles get narrower, the polarization fraction declines and outer components become more prominent as the frequency increases. Many counterexamples to these rules are also observed, and pulsars with complex profiles are especially prone to rule breaking. We hypothesise that the location of pulsar emission within the magnetosphere evolves with time as the the pulsar spins down. In highly energetic pulsars, the emission comes from a confined range of high altitudes, in the middle range of spin down energies the emission occurs over a wide range of altitudes whereas in pulsars with low spin-down energies it is confined to low down in the magnetosphere.

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.

The motivation for considering distributed large scale dynamos in the solar context is reviewed in connection with the magnetic helicity constraint. Preliminary accounts of 3-dimensional direct numerical simulations (in spherical shell segments) and simulations of 2-dimensional mean field models (in spherical shells) are presented. Interesting similarities as well as some differences are noted.

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.

In an earlier paper (Ahuja, et al, 2005), based on simultaneous multi-frequency observations with the Giant Metrewave Radio Telescope (GMRT), we reported the variation of pulsar dispersion measures (DMs) with frequency. A few different explanations are possible for such frequency dependence, and a possible candidate is the effect of pulse shape evolution on the DM estimation technique. In this paper we describe extensive simulations we have done to investigate the effect of pulse profile evolution on pulsar DM estimates. We find that it is only for asymmetric pulse shapes that the DM estimate is significantly affected due to profile evolution with frequency. Using multi-frequency data sets from our earlier observations, we have carried out systematic analyses of PSR B0329+54 and PSR B1642-03. Both these pulsars have central core dominated emission which does not show significant asymmetric profile evolution with frequency. Even so, we find that the estimated DM shows significant variation with frequency for these pulsars. We also report results from new, simultaneous multi-frequency observations of PSR B1133+16 carried out using the GMRT in phased array mode. This pulsar has an asymmetric pulse profile with significant evolution with frequency. We show that in such a case, amplitude of the observed DM variations can be attributed to profile evolution with frequency. We suggest that genuine DM variations with frequency could arise due to propagation effects through the interstellar medium and/or the pulsar magnetosphere.

PSR B0031-07 is well known to exhibit three different modes of drifting sub-pulses (mode A, B and C). It has recently been shown that in a multifrequency observation, consisting of 2700 pulses, all driftmodes were visible at low frequencies, while at 4.85 GHz only mode-A drift or non-drifting emission was detected. This suggests that modes A and B are emitted in sub-beams, rotating at a fixed distance from the magnetic axis, with the mode-B sub-beams being closer to the magnetic axis than the mode-A sub-beams. Diffuse emission between the sub-beams can account for the non-drifting emission. Using the results of an analysis of simultaneous multifrequency observations of PSR B0031-07, we set out to construct a geometrical model that includes emission from both sub-beams and diffuse emission and describes the regions of the radio emission of PSR B0031-07 at each emission frequency for driftmodes A and B. Based on the vertical spacing between driftbands, we have determined the driftmode of each sequence of drift. To restrict the model, we calculated average polarisation and intensity characteristics for each driftmode and at each frequency. The model reproduces the observed polarisation and intensity characteristics, suggesting that diffuse emission plays an important role in the emission properties of PSR B0031-07. The model further suggests that the emission heights of this pulsar range from a few kilometers to a little over 10 kilometers above the pulsar surface. We also find that the relationships between height and frequency of emission that follow from curvature radiation and from plasma-frequency emission could not be used to reproduce the observed frequency dependence of the width of the average intensity profiles.

Here we present the results from an analysis of a multifrequency simultaneous observation of PSR B0031$-$07. We have constructed a geometrical model, based on an empirical relationship between height and frequency of emission, that reproduces many of the observed characteristics. The model suggests very low emission altitudes for this pulsar of only a few kilometers above the star's surface.

We report the discovery of a young pulsar associated with the supernova remnant G21.5-0.9, using the Giant Metrewave Radio Telescope (GMRT) located near Pune, India. Discovered at a frequency of 610 MHz, J1833-1034 has a period of 61.86 ms and a period derivative of $2.0 \times 10^{-13}$, making it similar to other known young pulsars. The characteristic age of the pulsar is $\approx 4900$ yr, somewhat higher than estimates for the age of the remnant, but not incompatible with it. The pulsar has a spin-down luminosity of $3.3 \times 10^{37}$ erg s$^{-1}$, which is the second highest amongst all the known Galactic pulsars.

The combination of dispersion measures of pulsars, distances from the model of Cordes and Lazio (2002) and emission measures from the WHAM survey enabled a statistical study of electron densities and filling factors of the diffuse ionized gas (DIG) in the Milky Way. The emission measures were corrected for absorption and contributions from beyond the pulsar distance. For a sample of 157 pulsars at |b| > 5 degrees, mainly located in interarm regions within about 3 kpc from the Sun, we find that: (1) The average volume filling factor along the line of sight is inversely proportional to the mean electron density in clouds. (2) The average volume filling factor increases towards larger distances from the Galactic plane. (3) The local volume filling factor may reach a maximum near |z| = 0.9 kpc, whereas the local electron density continues to decrease at higher |z|, thus causing the observed flattening in the distribution of dispersion measures perpendicular to the plane above this height. (4) The scale heights of the electron density, the volume filling factor and the emission measure are the same and in the range 250-500 pc.

The well known drifter PSR B0031-07 is known to exhibit drifting subpulses where the spacing between the drift bands (P_3) shows three distinct modes A, B and C corresponding to 12, 6 and 4 seconds respectively. We have investigated periodicities and polarisation properties of PSR B0031-07 for a sequence of 2700 single pulses taken simultaneously at 328 MHz and 4.85 GHz. We found that mode A occurs simultaneously at these frequencies, while modes B and C only occur at 328 MHz. However, when the pulsar is emitting in mode B at the lower frequency there is still emission at the higher frequency, hinting towards the presence of mode B emission at a weaker level. Further, we have established that modes A and B are associated with two orthogonal modes of polarisation, respectively. Based on these observations, we suggest a geometrical model where modes A and B at a given frequency are emitted in two concentric rings around the magnetic axis with mode B being nested inside mode A. Further, it is evident that this nested configuration is preserved across frequency with the higher frequency arising closer to the stellar surface compared to the lower one, consistent with the well known radius-to-frequency mapping operating in pulsars.

In this paper, we describe a novel experiment for the accurate estimation of pulsar dispersion measures using the Giant Metre-wave Radio Telescope. This experiment was carried out for a sample of twelve pulsars, over a period of more than one year (January 2001 to May 2002) with observations about once every fortnight. At each epoch, the pulsar DMs were obtained from simultaneous dual frequency observations, without requiring any absolute timing information. The DM estimates were obtained from both the single pulse data streams and from the average profiles. The accuracy of the DM estimates at each epoch is ~ 1 part in 10^4 or better, making the data set useful for many different kinds of studies. The time series of DM shows significant variations on time scales of weeks to months for most of the pulsars. A comparison of the mean DM values from these data show significant deviations from catalog values (as well as from other estimates in literature) for some of the pulsars, with PSR B1642-03 showing the most notable changes. From our analysis results it appears that constancy of pulsar DMs (at the level of 1 in 10^3 or better) can not be taken for granted. For PSR B2217+47, we see evidence for a large-scale DM gradient over a one year period, which is modeled as being due to a blob of enhanced electron density sampled by the line of sight. For some pulsars, including pulsars with fairly simple profiles like PSR B1642-03, we find evidence for small changes in DM values for different frequency pairs of measurement, a result that needs to be investigated in detail. Another interesting result is that we find significant differences in DM values obtained from average profiles and single pulse data.

We present results of timing measurements of the binary millisecond pulsar PSR J2145-0750. Combining timing data obtained with the Effelsberg and Lovell radio telescopes we measure a significant timing parallax of 2.0(6) mas placing the system at 500 pc distance to the solar system. The detected secular change of the projected semi-major axis of the orbit $\dot x=1.8(6)\times 10^{-14}$ lt-s s$^{-1}$, where $x=(a_{\rm p}\sin i)/c$, is caused by the proper motion of the system. With this measurement we can constrain the orbital inclination angle to $i<61\degr$, with a median likelihood value of $46\degr$ which is consistent with results from polarimetric studies of the pulsar magnetosphere. This constraint together with the non-detection of Shapiro delay rules out certain combinations of the companion mass, $m_2$, and the inclination, $i$. For typical neutron star masses and using optical observations of the carbon/oxygen-core white dwarf we derive a mass range for the companion of $0.7 M_\odot\leq m_2\leq 1.0 M_\odot$. We apply evolutionary white dwarf cooling models to revisit the cooling age of the companion. Our analysis reveals that the companion has an effective temperature of $T_{\rm eff}=5750\pm600$ K and a cooling age of $\tau_{\rm cool}=3.6(2)$ Gyr, which is roughly a factor of three lower than the pulsar's characteristic age of 10.4 Gyr. The cooling age implies an initial spin period of $P_0=13.0(5)$ ms, which is very close to the current period.

In this paper we report multi-frequency measurements of pulse broadening times (tau_d) for nine medium dispersion measure (DM $\approx 150-400$ pc cm$^{-3}$) pulsars observed over a wide frequency range. The low frequency data at 243, 325 and 610 MHz are new observations done with the Giant Metrewave Radio Telescope (GMRT). The frequency dependence of tau_d for all but one (PSR B1933+16) of our sources is consistent with the Kolmogorov spectrum of electron density fluctuations in a turbulent medium. PSR B1933+16, however, shows a very flat spectrum as previously observed for high DM pulsars. Our observations combined with earlier published results enable us to study the spectral index of tau_d over the whole observed DM range. While the spectral properties are generally consistent with a Kolmogorov spectrum, pulsars seen along line-of-sights towards the inner Galaxy or complex regions often show deviations from this expected behaviour.

The linear polarization position angle (PPA) of pulsar emission shows a wide variety of structures, which deviates from the classical rotating vector model. In this study we investigate the effect on the PPA due to emission arising from a range of heights. We demonstrate that by appropriately choosing the emission height we can reconstruct the PPA sweep across the integrated pulse profiles.

We use a set of carefully selected published average multifrequency polarimetric observations for six bright cone dominated pulsars and devise a method to combine the multifrequency polarization position angle (PPA) sweep traverses. We demonstrate that the PPA traverse is in excellent agreement with the rotating vector model over this broad frequency range confirming that radio emission emanates from perfectly dipolar field lines. For pulsars with central core emission in our sample, we find the peak of central core component to lag the steepest gradient of the PPA traverse at several frequencies. Also significant frequency evolution of the core width is observed over this frequency range. The above facts strongly suggest: (a) the peak core emission does not lie on the fiducial plane containing the dipole magnetic axis and the rotation axis, and (b) the core emission does not originate from the polar cap surface.

Combining dispersion measures, distances and emission measures for 157 pulsars lying above \mid b \mid > 5 degree and between 60 degree < l < 360 degree we find the mean volume filling factor (\bar{f_v}) of the diffused ionized gas in the Milky Way. This filling factor is inversely related to the mean electron density (\bar{n_c}) in the clouds, \bar{f_v} = (0.0184 +/- 0.0014) \bar{n_c}^{-1.07 +/- 0.03}, implying a nearly constant average electron density distribution within 3 kpc of the sun.

We present a study of single pulses from nine bright northern pulsars to investigate the behaviour of circular polarisation, V. The observations were conducted with the Effelsberg 100-m radio telescope at 1.41 GHz and 4.85 GHz and the Westerbork radio telescope at 352 MHz. For the first time, we present the average profile of the absolute circular polarisation |V| in the single pulses. We demonstrate that the average profile of |V| is the distinguishing feature between pulse components that exhibit low V in the single pulses and components that exhibit high V of either handedness, despite both cases resulting in a low mean. We also show that the |V| average profile remains virtually constant with frequency, which is not generally the case for V, leading us to the conclusion that |V| is a key quantity in the pulsar emission problem.

We have obtained new rotation measure for 11 pulsars observed with the Effelsberg 100-m radio telescope, in the direction of the Perseus arm. Using a combination of 34 published and the 11 newly measured pulsar rotation measures we study the magnetic field structure towards the Perseus arm. We find that two pulsars towards l$\sim$ 149$^{\circ}$ (Region 1) and four pulsars towards l$\sim113^{\circ}$ (Region 2) lie behind HII regions which seriously affects the pulsar rotation measures. The rotation measure of PSR J2337+6151 seem to be affected by its passage through the supernova remnant G114.3+0.3. For Region 1, we are able to constrain the random component of the magnetic field to $5.7\mu$G. For the large-scale component of the Galactic magnetic field we determine a field strength of $1.7\pm1.0\mu$G. This average field is constant on Galactic scales lying within the Galactic longitude range of $85^{\circ} <$ l $ < 240^{\circ}$ and we find no evidence for large scale field reversal upto 5-6 kpc. We highlight the great importance to include the effects of foreground emission in any systematic study.

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.

PSR B0943+10 is known to show remarkably stable drifting subpulses, which can be interpreted in terms of a circumferential motion of 20 sparks, each completing one circulation around the periphery of the polar cap in 37 pulsar periods. We use this observational constraint and argue that the vacuum gap model can adequately describe the observed drift patterns. Further we demonstrate that {\em only} the presence of strong non-dipolar surface magnetic field can favor such vacuum gap formation. Subsequently, for the first time we are able to constrain the parameters of the surface magnetic field, and model the expected magnetic structure on the polar cap of PSR B0943+10 considering the inverse Compton scattering photon dominated vacuum gap.

We propose a vacuum gap (VG) model which can be applied uniformly for normal and high magnetic field pulsars. The model requires strong and non-dipolar surface magnetic field near the pulsar polar cap. We assume that the actual surface magnetic field in pulsars results from a superposition of global dipole field and crust-anchored small scale magnetic anomaly. We provide a numerical formalism for modelling such structures of surface magnetic field and explore it within the framework of VG model, which requires strong surface fields more than 10^{13} G.

We report multifrequency measurements of scatter broadening times for nine highly dispersed pulsars over a wide frequency range (0.6 -- 4.9 GHz). We find the scatter broadening times to be larger than expected and to scale with frequency with an average power-law index of $3.44\pm 0.13$, i.e. significantly less than that expected from standard theories. Such possible discrepancies have been predicted very recently by Cordes & Lazio.

We have measured the scatter broadening of pulsars in the direction of the Gum nebula. For the first time, our observations show clear variations of scattering properties across the Gum nebula. The IRAS-Vela shell is shown to be a high scattering region. Our revised estimations of distances to these pulsars are consistently less by a factor of 2--3, which has very important consequences for the deduced values of radio luminosity and transverse velocity of pulsars.

In this paper we revisit the radio pulsar death line problem within the framework of curvature radiation and/or inverse compton scattering induced vacuum gap model above neutron star polar caps. Our special interest is in the recently detected pulsar PSR J2144-3933 with extremal period 8.5 seconds, which lies far beyond conventional death lines. We argue, that formation of vacuum gaps requires a complicated multipolar surface magnetic field, with a strenght $B_s$ much higher than the surface dipolar component $B_d$, and radii of curvature ${\cal R}$ much smaller than the neutron star radius $R=10^6$ cm. Such a multipolar surface field is also consistent with death lines including the extremal pulsar PSR J2144-3933. Since vacuum gap models produce sparks, our paper naturally supports the spark related models of subpulse drift phenomenon as well as to the spark associated models of coherent pulsar radio emission.

The evolution of the multipolar structure of the magnetic field of isolated neutron stars is studied assuming the currents to be confined to the crust. Lower orders ($\le 25$) of multipole are seen to evolve in a manner similar to the dipole suggesting little or no evolution of the expected pulse shape. We also study the multifrequency polarization position angle traverse of PSR B0329+54 and find a significant frequency dependence above 2.7 GHz. We interpret this as an evidence of strong multipolar magnetic field present in the radio emission region.

The evolution of the multipolar structure of the magnetic field of isolated neutron stars is studied assuming the currents to be confined to the crust. We find that except for multipoles of very high order ($l\gsim 25$) the evolution is similar to that of a dipole. Therefore no significant evolution is expected in pulse shape of isolated radio pulsars due to the evolution of the multipole structure of the magnetic field.

Characterizing the shape and evolution of pulsar radio emission beams is important for understanding the observed emission. The various attempts by earlier workers investigating beam shapes have resulted in widely different conclusions. Using a carefully selected subset of the recently published multifrequency polarimetry observations of 300 radio pulsars (Gould & Lyne, 1998), we attempt to model the shape of pulsar beams. Assuming that the beam shape is elliptical, in general, and that it may depend on the angle between the rotation and the magnetic axes, we seek a consistent model where we also solve for the dependence of the beam size on frequency. From the six-frequency data on conal triple and multiple component profiles, we show that a) the pulsar emission beams follow a nested cone structure with at least three distinct cones, although only one or more of the cones may be active in a given pulsar; b) each emission cone is illuminated in the form of an annular ring of width typically about 20% of the cone radius. Although some slight preference is evident for a model where the beam is circular for an aligned rotator & latitudinally compressed for an orthogonal rotator, the possibility that the beam shape is circular at all inclinations is found to be equally consistent with the data. While the overall size scales as $P^{-0.5}$ (where $P$ is the pulsar period) as expected from the notion of dipolar open field lines, we see no evidence in support of the beam shape evolution with pulsar period.