We introduce pinta, a pipeline for reducing the uGMRT raw pulsar timing data developed for the Indian Pulsar Timing Array experiment. We provide a detailed description of the workflow and usage of pinta, as well as its computational performance and RFI mitigation characteristics. Furthermore, the results of a calibration experiment carried out to determine the relative time offsets between different back-end modes and the correct interpretation of the observatory frequency settings at the uGMRT, which are crucial for performing precision pulsar timing, are also discussed.

We present an experiment performed to understand the polarization convention adopted for band 4 (550--900 MHz) of the upgraded Giant Metrewave Radio Telescope (uGMRT). For that we observed the pulsar B1702--19 in this band, both in interferometry and pulsar modes, and compare the results with its already known Stokes $I$, $Q$, $U$, $V$ profiles obtained with the Lovell telescope. We find that the results obtained from interferometry and pulsar modes of the uGMRT agree with each other. However, although the Stokes $U$ profile obtained with the uGMRT match with that obtained by the Lovell telescope, Stokes $Q$ and $V$ do not. This can be explained if the $X$ and $Y$ dipoles in this band, from which $R$ and $L$ are derived, are swapped w.r.t. the IAU convention. The swapping makes $RR^*$ and $LL^*$ of the uGMRT band 4 to be $LL^*$ and $RR^*$ respectively according to the IEEE convention. This implies that if we need to compare polarization measurements obtained in band 4 of the uGMRT with telescopes like Lovell, Parkes, Very Large Array etc. (all follow IEEE convention for defining right and left hand circular polarization), we must interchange $RR^*$ and $LL^*$, and change the sign of Stokes $Q$ for the uGMRT data. Note that this is the current convention for uGMRT band 4, and is likely to change in future once the swapping of the dipoles is taken care of. Once it is done, it will be notified in another technical report.

(abridged) In this work we present the results of one year of upgraded Giant Metrewave Radio Telescope timing measurements of PSR~J0514$-$4002A, a 4.99-ms pulsar in a 18.8-day, eccentric ($e \, =\, 0.89$) orbit with a massive companion located in the globular cluster NGC~1851. Combining these data with earlier Green Bank Telescope data, we greatly improve the precision of the rate of advance of periastron, $\dot{\omega} \, = \, 0.0129592(16)\, \deg \, \rm yr^{-1}$ which, assuming the validity of general relativity, results in a much refined measurement of the total mass of the binary, $M_{\rm tot} \, = \, 2.4730(6) \, M_{\odot}$. Additionally, we measure the Einstein delay parameter, $\gamma\,, = \, 0.0216(9) \, \rm s$. Furthermore, we measure the proper motion of the system ($\mu_{\alpha} \, = \, 5.19(22)$ and $\mu_{\delta} = -0.56(25)\rm~mas ~ yr^{-1}$), which is not only important for analyzing its motion in the cluster, but is also essential for a proper interpretation of $\gamma$, given the latter parameter's correlation with the variation of the projected semi-major axis. The measurements of $\gamma$ and the proper motion enable a separation of the system component masses: we obtain a pulsar mass of $M_{p} \, = \, 1.25^{+0.05}_{-0.06} \, M_{\odot}$ and a companion mass of $M_{c} \, = \, 1.22^{+0.06}_{-0.05} \, M_{\odot}$. This raises the possibility that the companion is also a neutron star. Searches for radio pulsations from the companion have thus far been unsuccessful, hence we cannot confirm the latter hypothesis. The low mass of this millisecond pulsar - one of the lowest ever measured for such objects - clearly indicates that the recycling process can be achieved with a relatively small amount of mass transfer.

With 30 antennas and a maximum baseline length of 25 km, the Giant Metrewave Radio Telescope (GMRT) is the premier low-frequency radio interferometer today. We have carried out a study of possible expansions of the GMRT, via adding new antennas and installing focal plane arrays (FPAs), to improve its point-source sensitivity, surface brightness sensitivity, angular resolution, field of view, and U-V coverage. We have carried out array configuration studies, aimed at minimizing the number of new GMRT antennas required to obtain a well-behaved synthesized beam over a wide range of angular resolutions for full-synthesis observations. This was done via two approaches, tomographic projection and random sampling, to identify the optimal locations for the new antennas. We report results for the optimal locations of the antennas of an expanded array (the "EGMRT"), consisting of the existing 30 GMRT antennas, 30 new antennas at short distances, $\leq 2.5$ km from the array centre, and 26 new antennas at long distances, $\approx 5-25$ km from the array centre. The collecting area and the field of view of the proposed EGMRT array would be larger by factors of, respectively, $\approx 3$ and $\approx 30$, than those of the GMRT. Indeed, the EGMRT continuum sensitivity and survey speed with 550-850 MHz FPAs installed on the 45 antennas within a distance of $\approx 2.5$ km of the array centre would be far better than those of any existing interferometer, and comparable to the sensitivity and survey speed of Phase-1 of the Square Kilometre Array.

We report the discovery of a diffuse radio emission around PSR J0855--4644 using an upgraded GMRT (uGMRT) observation at 1.35 GHz. The radio emission is spatially coincident with the diffuse X-ray pulsar wind nebula (PWN) seen with XMM but is much larger in extent compared to the compact axisymmetric PWN seen with Chandra. The morphology of the emission, with a bright partial ring-like structure and two faint tail-like features strongly resembles a bow shock nebula, and indicates a velocity of 100 km/s through the ambient medium. We conclude that the emission is most likely to be associated with the radio PWN of PSR J0855-4644. From the integrated flux density, we estimate the energetics of the PWN.

The Giant Metrewave Radio Telescope (GMRT) is being upgraded to increase the receiver sensitivity. This makes the receiver more susceptible to man-made Radio Frequency Interference (RFI). To improve the receiver performance in presence of RFI, real-time RFI excision (filtering) is incorporated in the GMRT wideband backend (GWB). The RFI filtering system is implemented on FPGA and CPU-GPU platforms to detect and remove broadband and narrowband RFI. The RFI is detected using a threshold-based technique where the threshold is computed using Median Absolute Deviation (MAD) estimator. The filtering is carried out by replacing the RFI samples by either noise samples or constant value or threshold. This paper describes the status of the real-time broadband RFI excision system in the wideband receiver chain of the upgraded GMRT (uGMRT). The test methodology for carrying out various tests to demonstrate the performance of broadband RFI excision at the system level and on radio astronomical imaging experiments are also described.

Microstructure emission, involving short time scale, often quasi-periodic, intensity fluctuations in subpulse emission, is well known in normal period pulsars. In this letter, we present the first detections of quasi-periodic microstructure emission from millisecond pulsars (MSPs), from Giant Metrewave Radio Telescope (GMRT) observations of two MSPs at 325 and 610 MHz. Similar to the characteristics of microstructure observed in normal period pulsars, we find that these features are often highly polarized, and exhibit quasi-periodic behavior on top of broader subpulse emission, with periods of the order of a few $\mu$s. By measuring their widths and periodicities from single pulse intensity profiles and their autocorrelation functions, we extend the microstructure timescale - rotation period relationship by more than an order of magnitude down to rotation periods $\sim$ 5 ms, and find it to be consistent with the relationship derived earlier for normal pulsars. The similarity of behavior is remarkable, given the significantly different physical properties of MSPs and normal period pulsars, and rules out several previous speculations about the possible different characteristics of microstructure in MSP radio emission. We discuss the possible reasons for the non-detection of these features in previous high time resolution MSP studies along with the physical implications of our results, both in terms of a geometric beam sweeping model and temporal modulation model for micropulse production.

It is an exceptionally opportune time for Astrophysics when a number of next-generation mega-instruments are poised to observe the universe across the entire electromagnetic spectrum with unprecedented data quality. The Square Kilometre Array (SKA) is undoubtedly one of the major components of this scenario. In particular, the SKA is expected to discover tens of thousands of new neutron stars giving a major fillip to a wide range of scientific investigations. India has a sizeable community of scientists working on different aspects of neutron star physics with immediate access to both the uGMRT (an SKA pathfinder) and the recently launched X-ray observatory Astrosat. The current interests of the community largely centre around studies of - a) the generation of neutron stars and the SNe connection}, b) the neutron star population and evolutionary pathways}, c) the evolution of neutron stars in binaries and the magnetic fields}, d) the neutron star equation of state}, e) the radio pulsar emission mechanism}, and, f) the radio pulsars as probes of gravitational physics}. Most of these studies are the main goals of the SKA first phase, which is likely to be operational in the next four years. This article summarises the science goals of the Indian neutron star community in the SKA era, with significant focus on coordinated efforts among the SKA and other existing/upcoming instruments.

The Square Kilometre Array Telescope Manager (SKA TM) will be responsible for assisting the SKA Operations and Observation Management, carrying out System diagnosis and collecting Monitoring & Control data from the SKA sub-systems and components. To provide adequate compute resources, scalability, operation continuity and high availability, as well as strict Quality of Service, the TM cyber-infrastructure (embodied in the Local Infrastructure - LINFRA) consists of COTS hardware and infrastructural software (for example: server monitoring software, host operating system, virtualization software, device firmware), providing a specially tailored Infrastructure as a Service (IaaS) and Platform as a Service (PaaS) solution. The TM infrastructure provides services in the form of computational power, software defined networking, power, storage abstractions, and high level, state of the art IaaS and PaaS management interfaces. This cyber platform will be tailored to each of the two SKA Phase 1 telescopes (SKA_MID in South Africa and SKA_LOW in Australia) instances, each presenting different computational and storage infrastructures and conditioned by location. This cyber platform will provide a compute model enabling TM to manage the deployment and execution of its multiple components (observation scheduler, proposal submission tools, M&C components, Forensic tools and several Databases, etc). In this sense, the TM LINFRA is primarily focused towards the provision of isolated instances, mostly resorting to virtualization technologies, while defaulting to bare hardware if specifically required due to performance, security, availability, or other requirement.

A fully real-time coherent dedispersion system has been developed for the pulsar back-end at the Giant Metrewave Radio Telescope (GMRT). The dedispersion pipeline uses the single phased array voltage beam produced by the existing GMRT software back-end (GSB) to produce coherently dedispersed intensity output in real time, for the currently operational bandwidths of 16 MHz and 32 MHz. Provision has also been made to coherently dedisperse voltage beam data from observations recorded on disk. We discuss the design and implementation of the real-time coherent dedispersion system, describing the steps carried out to optimise the performance of the pipeline. Presently functioning on an Intel Xeon X5550 CPU equipped with a NVIDIA Tesla C2075 GPU, the pipeline allows dispersion free, high time resolution data to be obtained in real-time. We illustrate the significant improvements over the existing incoherent dedispersion system at the GMRT, and present some preliminary results obtained from studies of pulsars using this system, demonstrating its potential as a useful tool for low frequency pulsar observations. We describe the salient features of our implementation, comparing it with other recently developed real-time coherent dedispersion systems. This implementation of a real-time coherent dedispersion pipeline for a large, low frequency array instrument like the GMRT, will enable long-term observing programs using coherent dedispersion to be carried out routinely at the observatory. We also outline the possible improvements for such a pipeline, including prospects for the upgraded GMRT which will have bandwidths about ten times larger than at present.

In this paper we present results from flux density measurements for 21 pulsars over a wide frequency range, using the Giant Metrewave Radio Telescope (GMRT) and the Effelsberg telescope. Our sample was a set of mostly newly discovered pulsars from the selection of candidates for gigahertz-peaked spectra (GPS) pulsars. Using the results of our observations along with previously published data, we identify two new GPS pulsars. One of them, PSR J1740+1000, with dispersion measure of 24 pc cm$^{-3}$, is the first GPS pulsar with such a low DM value.We also selected several strong candidates for objects with high frequency turnover in their spectra which require further investigation.We also revisit our source selection criteria for future searches for GPS pulsars.

Using the Giant Metrewave Radio Telescope (GMRT) we performed deep observations to search for radio pulsations in the directions of unidentified Fermi Large Area Telescope (LAT) gamma-ray sources. We report the discovery of an eclipsing black-widow millisecond pulsar, PSR J1544+4937, identified with the un-cataloged gamma-ray source Fermi J1544.2+4941. This 2.16 ms pulsar is in a 2.9 hours compact circular orbit with a very low-mass companion (Mc > 0.017 Msun). At 322 MHz this pulsar is found to be eclipsing for 13% of its orbit, whereas at 607 MHz the pulsar is detected throughout the low-frequency eclipse phase. Variations in the eclipse ingress phase are observed, indicating a clumpy and variable eclipsing medium. Moreover, additional short-duration absorption events are observed around the eclipse boundaries. Using the radio timing ephemeris we were able to detect gamma-ray pulsations from this pulsar, confirming it as the source powering the gamma-ray emission.

The GMRT Epoch of Reionization (EoR) experiment is an ongoing effort to measure the power spectrum from neutral hydrogen at high redshift. We have previously reported an upper limit of (70 mK)^2 at wavenumbers of k=0.65 h/Mpc using a basic piecewise-linear foreground subtraction. In this paper we explore the use of a singular value decomposition to remove foregrounds with fewer assumptions about the foreground structure. Using this method we also quantify, for the first time, the signal loss due to the foreground filter and present new power spectra adjusted for this loss, providing a revised measurement of a 2-sigma upper limit at (248 mK)^2 for k=0.50 h/Mpc. While this revised limit is larger than previously reported, we believe it to be more robust and still represents the best current constraints on reionization at z=8.6.

Next-generation radio arrays, including the Square Kilometre Array (SKA) and its pathfinders, will open up new avenues for exciting transient science at radio wavelengths. Their innovative designs, comprising a large number of small elements, pose several challenges in digital processing and optimal observing strategies. The Giant Metrewave Radio Telescope (GMRT) presents an excellent test-bed for developing and validating suitable observing modes and strategies for transient experiments with future arrays. Here we describe the first phase of the ongoing development of a transient detection system for GMRT that is planned to eventually function in a commensal mode with other observing programs. It capitalizes on the GMRT's interferometric and sub-array capabilities, and the versatility of a new software backend. We outline considerations in the plan and design of transient exploration programs with interferometric arrays, and describe a pilot survey that was undertaken to aid in the development of algorithms and associated analysis software. This survey was conducted at 325 and 610 MHz, and covered 360 deg$^2$ of the sky with short dwell times. It provides large volumes of real data that can be used to test the efficacies of various algorithms and observing strategies applicable for transient detection. We present examples that illustrate the methodologies of detecting short-duration transients, including the use of sub-arrays for higher resilience to spurious events of terrestrial origin, localization of candidate events via imaging and the use of a phased array for improved signal detection and confirmation. In addition to demonstrating applications of interferometric arrays for fast transient exploration, our efforts mark important steps in the roadmap toward SKA-era science.

We present the results from nearly three years of monitoring of the variations in dispersion measure (DM) along the line-of-sight to 11 millisecond pulsars using the Giant Metrewave Radio Telescope (GMRT). These results demonstrate accuracies of single epoch DM estimates of the order of 5x10^(-4) cm^(-3) pc. A preliminary comparison with the Parkes Pulsar Timing Array (PPTA) data shows that the measured DM fluctuations are comparable. We show effects of DM variations due to the solar wind and solar corona and compare with the existing models.

We have developed a multi-pixel beamformer technique, which can be used for enhancing the capabilities for studying pulsars using an interferometric array. Using the Giant Metrewave Radio Telescope (GMRT), we illustrate the application of this efficient technique, which combines the enhanced sensitivity of a coherent array beamformer with the wide field-of-view seen by an incoherent array beamformer. Multi-pixel beamformer algorithm is implemented using the recorded base-band data. With the optimisations in multi-pixelisation described in this paper, it is now possible to form 16 directed beams in real-time. We discuss a special application of this technique, where we use continuum imaging followed by the multi-pixel beamformer to obtain the precise locations of newly discovered millisecond pulsars with the GMRT. Accurate positions measured with single observations enable highly sensitive follow-up studies using coherent array beamformer and rapid follow up at higher radio frequencies and other wavelengths. Normally, such accurate positions can only be obtained from a long-term pulsar timing program. The multi-pixel beamformer technique can also be used for highly sensitive targeted pulsar searches in extended supernova remnants. In addition this method can provide optimal performance for the large scale pulsar surveys using multi-element arrays.

We present the results from timing observations with the GMRT of the young pulsar J1833-1034, in the galactic supernova remnant G21.5-0.9. We detect the presence of 4 glitches in this pulsar over a period of 5.5 years, making it one of a set of pulsars that show fairly frequent glitches. The glitch amplitudes, characterized by the fractional change of the rotational frequency, range from 1 \times 10^-9 to 7 \times 10^-9, with no evidence for any appreciable relaxation of the rotational frequency after the glitches. The fractional changes observed in the frequency derivative are of the order of 10-5 . We show conclusively that, in spite of having significant timing noise, the sudden irregularities like glitches detected in this pulsar can not be modeled as smooth timing noise. Our timing solution also provides a stable estimate of the second derivative of the pulsar spin-down model, and a plausible value for the braking index of 1.857, which, like the value for other such young pulsars, is much less than the canonical value of 3.0. PSR J1833-1034 appears to belong to a class of pulsars exhibiting fairly frequent occurrence of low amplitude glitches. This is further supported by an estimate of the glitch activity parameter, Ag = 1.53 \times 10^-15 s^-2, which is found to be significantly lower than the trend of glitch activity versus characteristic age (or spin frequency derivative) that a majority of the glitching pulsars follow. We present evidence for a class of such young pulsars, including the Crab, where higher internal temperature of the neutron star could be responsible for the nature of the observed glitch activity.

We present a summary of the Fermi Pulsar Search Consortium (PSC), an international collaboration of radio astronomers and members of the Large Area Telescope (LAT) collaboration, whose goal is to organize radio follow-up observations of Fermi pulsars and pulsar candidates among the LAT gamma-ray source population. The PSC includes pulsar observers with expertise using the world's largest radio telescopes that together cover the full sky. We have performed very deep observations of all 35 pulsars discovered in blind frequency searches of the LAT data, resulting in the discovery of radio pulsations from four of them. We have also searched over 300 LAT gamma-ray sources that do not have strong associations with known gamma-ray emitting source classes and have pulsar-like spectra and variability characteristics. These searches have led to the discovery of a total of 43 new radio millisecond pulsars (MSPs) and four normal pulsars. These discoveries greatly increase the known population of MSPs in the Galactic disk, more than double the known population of so-called `black widow' pulsars, and contain many promising candidates for inclusion in pulsar timing arrays.

We present the results of flux density monitoring of PSR B0329+54 at the frequency of 4.8 GHz using the 32-meter TCfA radiotelescope. The observations were conducted between 2002 and 2005. The main goal of the project was to find interstellar scintillation (ISS) parameters for the pulsar at the frequency at which it was never studied in detail. To achieve this the 20 observing sessions consisted of 3-minute integrations which on average lasted 24 hours. Flux density time series obtained for each session were analysed using structure functions. For some of the individual sessions as well as for the general average structure function we were able to identify two distinctive timescales present, the timescales of diffractive and refractive scintillations. To the best of our knowledge, this is the first case when both scintillation timescales, t_DISS=42.7 minutes and t_RISS=305 minutes, were observed simultaneously in a uniform data set and estimated using the same method. The obtained values of the ISS parameters combined with the data found in the literature allowed us to study the frequency dependence of these parameters over a wide range of observing frequencies, which is crucial for understanding the ISM turbulence. We found that the Kolmogorov spectrum is not best suited for describing the density fluctuations of the ISM, and a power-law spectrum with beta =4 seems to fit better with our results. We were also able to estimate the transition frequency (transition from strong to weak scintillation regimes) as 10.1 GHz, much higher than was previously predicted. We were also able to estimate the strength of scattering parameter u=2.67$ and the Fresnel scale as 6.7x10^8 meters.

We describe the serendipitous discovery of a radio point source in a 618 MHz image of the supernova remnant(SNR) G76.9+1.0. The SNR has a bipolar structure and the point source is located near a faint bridge of emission joining the two lobes of emission. The point source was also detected in follow-up higher frequency(1170 MHz) observations. The spectral index for the point source obtained from the GMRT observations is alpha = -2.1. The steep spectrum, as well as the location of the point source near the centre of the SNR establish the fact that it is indeed the pulsar J2022+3842 associated with this SNR. Consistent with this, subsequent analysis of archival Chandra X-ray data shows a point source coincident with the radio point source, as well as diffuse extended X-ray emission surrounding the unresolved source. However, no pulsed emission was detected despite deep searches at both 610 MHz and 1160 MHz although pulsed emission has been seen at 2 GHz with the GBT. It appears that the most likely reason for not detecting the pulsed signal at the GMRT is temporal broadening: for the estimated DM towards this SNR, the pulse broadening time could be as large as tens of milliseconds. The diffuse X-ray emission is elongated along the same direction as the bipolar structure seen in the radio. We interpret the radio lobes as having been formed from an equatorial wind. Although direct detection of pulsed signal has not been possible, we show convincingly that sensitive, high-resolution, radio imaging at multiple frequencies is a useful method to search for pulsar candidates.

We investigate a high frequency turn-over effect in radio spectra for pulsars with positive or flat spectral index. Using GMRT and Effelsberg observatory, we estimated the flux density to reconstruct pulsar spectra. We find objects that have a maximum flux in their spectrum above 1 GHz and whose spectral indices are positive at lower frequencies. Some pulsars with a turn-over in their spectrum at high frequencies are found to exist in very interesting environments. We call these objects gigahertz-peaked spectra pulsars.

We present a new upper limit to the 21cm power spectrum during the Epoch of Reionization (EoR) which constrains reionization models with an unheated IGM. The GMRT-EoR experiment is an ongoing effort to make a statistical detection of the power spectrum of 21cm neutral hydrogen emission at redshift z~9. Data from this redshift constrain models of the EoR, the end of the Dark Ages arising from the formation of the first bright UV sources, probably stars or mini-quasars. We present results from approximately 50 hours of observations at the Giant Metrewave Radio Telescope in India from December 2007. We describe radio frequency interference (RFI) localisation schemes which allow bright sources on the ground to be identified and physically removed in addition to automated flagging. Singular-value decomposition is used to remove remaining broadband RFI by identifying ground sources with large eigenvalues. Foregrounds are modelled using a piecewise linear filter and the power spectrum is measured using cross-correlations of foreground subtracted images.

The new era of software signal processing has a large impact on radio astronomy instrumentation. Our design and implementation of a 32 antennae, 33 MHz, dual polarization, fully real-time software backend for the GMRT, using only off-the-shelf components, is an example of this. We have built a correlator and a beamformer, using PCI-based ADC cards and a Linux cluster of 48 nodes with dual gigabit inter-node connectivity for real-time data transfer requirements. The highly optimized compute pipeline uses cache efficient, multi-threaded parallel code, with the aid of vectorized processing. This backend allows flexibility in final time and frequency resolutions, and the ability to implement algorithms for radio frequency interference rejection. Our approach has allowed relatively rapid development of a fairly sophisticated and flexible backend receiver system for the GMRT, which will greatly enhance the productivity of the telescope. In this paper we describe some of the first lights using this software processing pipeline. We believe this is the first instance of such a real-time observatory backend for an intermediate sized array like the GMRT.

We report on the unique nulling properties of PSR B0818-41, using the GMRT at 325 and 610 MHz. We find following interesting behaviour just before and after the nulls: (i) The pulsar's intensity does not switch off abruptly at the null, but fades gradually, taking ~ 10 P1. Just after nulls intensity rises to a maximum over a short (less than one period) time scale. (ii) While the last active pulses before nulls are dimmer, the first few active pulses just after nulls outshine normal ones. This effect is very clear for inner region of pulsar profile, where mean intensity of last few active pulses just after nulls is ~ 2.8 times more than that for last active pulses just before nulls. (iii) There is a significant evolution of shape of the pulsar's profile, around nulls, especially at beginning of bursts: an enhanced bump of intensity in inner region, a change in ratio of strengths of the leading and trailing peaks towards a more symmetric profile, an increase in profile width of about 10%, and a shift of profile centre towards later longitudes. (iv) Just before nulls, the apparent drift rate becomes slower, transitioning to an almost phase stationary drift pattern. Further, when the pulsar comes out of null, the increased intensity is very often accompanied by what looks like a disturbed drift rate behaviour, which settles down to the regular drift pattern as the pulsar intensity returns to normal. Thus, we find some very specific and well correlated changes in the radio emission properties of PSR B0818-41 when the emission restarts after a null. These could imply that the phenomenon of nulling is associated with some kind of a "reset" of the pulsar radio emission engine. We also present plausible explanations for some of the observed behaviour, using the Partially Screened Gap model of the inner pulsar accelerator.

We have developed a method to compute the possible distribution of radio emission regions in a typical pulsar magnetosphere, taking into account the viewing geometry and rotational effects of the neutron star. Our method can estimate the emission altitude and the radius of curvature of particle trajectory as a function of rotation phase for a given inclination angle, impact angle, spin-period, Lorentz factor, field line constant and the observation frequency. Further, using curvature radiation as the basic emission mechanism, we simulate the radio intensity profiles that would be observed from a given distribution of emission regions, for different values of radio frequency and Lorentz factor. We show clearly that rotation effects can introduce significant asymmetries into the observed radio profiles. We investigate the dependency of profile features on various pulsar parameters. We find that the radiation from a given ring of field lines can be seen over a large range of pulse longitudes, originating at different altitudes, with varying spectral intensity. Preferred heights of emission along discrete sets of field lines are required to reproduce realistic pulsar profiles, and we illustrate this for a known pulsar. Finally, we show how our model provides feasible explanations for the origin of core emission, and also for one-sided cones which have been observed in some pulsars.

We present the results from a study of wide profile pulsars using high sensitivity multifrequency observations with the GMRT. Since the line of sight samples a large region of the polar cap in case of the wide profile pulsars, presence of simultaneous multiple drift regions is quite probable (as seen in PSR B0826-34 and PSR B0818-41). We solve the aliasing problem of PSR B0818-41 using the observed phase relationship of the drift regions, and determine its pattern rotation period P4 to be ~ 10s, which makes it the fastest known carousel. We find that, for all the pulsars showing drifting in multiple rings of emission, the drift pattern from the rings are phase locked. This can constraint the theoretical models of pulsar emission as it favors a pan magnetospeheric radiation mechanism.

The GMRT reionization effort aims to map out the large scale structure of the Universe during the epoch of reionization (EoR). Removal of polarized Galactic emission is a difficult part of any 21 cm EoR program, and we present new upper limits to diffuse polarized foregrounds at 150 MHz. We find no high significance evidence of polarized emission in our observed field at mid galactic latitude (J2000 08h26m+26). We find an upper limit on the 2-dimensional angular power spectrum of diffuse polarized foregrounds of [l^2 C_l/(2 PI)]^{1/2}< 3K in frequency bins of width 1 MHz at 300<l<1000. The 3-dimensional power spectrum of polarized emission, which is most directly relevant to EoR observations, is [k^3 P_p(k)/(2 PI^2)]^{1/2}< 2K at k_perp > 0.03 h/Mpc, k < 0.1 h/Mpc. This can be compared to the expected EoR signal in total intensity of [k^3 P(k)/ (2 PI^2) ]^{1/2} ~ 10 mK. We find polarized structure is substantially weaker than suggested by extrapolation from higher frequency observations, so the new low upper limits reported here reduce the anticipated impact of these foregrounds on EoR experiments. We discuss Faraday beam and depth depolarization models and compare predictions of these models to our data. We report on a new technique for polarization calibration using pulsars, as well as a new technique to remove broadband radio frequency interference. Our data indicate that, on the edges of the main beam at GMRT, polarization squint creates ~ 3% leakage of unpolarized power into polarized maps at zero rotation measure. Ionospheric rotation was largely stable during these solar minimum night time observations.

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.

PSR B0818-41 is one of the few pulsars which show multiple drift regions having well defined phase relationship. In this paper we report new results from the multifrequency GMRT observations of this pulsar. Significant linear polarization is observed with depolarization at the edge of the profile, due to orthogonal polarization. Circular polarization changes sign near the middle of the pulse profile at 1060 MHz, but not observed at 325 and 610 MHz. Remarkable frequency evolution of polarization angle (PA) is observed. Based on the frequency evolution of average profile, PA swing and results from subpulse drifting, we propose two possible emission geometries, alpha~11 deg, beta~-5.4 deg and alpha~175.4 deg, beta~-6.9 deg. Simulation of the pulsar radiation pattern with both these geometries reproduces the observed features in the drift pattern quite well. In addition to the remarkable subpulse drifting observed at 325 MHz, we report subpulse drifting at 244 and 610 MHz. We observe changes of drift rates, transitions from negative to stationary and stationary to negative drift rates, with some connection with nulls. Though P3m is the same for all the drift regions, P2m and Delta(Phi) values are different for the inner and outer drift regions. The peak emission from the leading and the trailing outer regions are offset by ~ 9 P1. Utilising this information, we solve the aliasing problem and predict P4~ 10 s, which makes it the fastest known carousel. The drift pattern from the two rings are phase locked for PSR B0818-41. The same is found to be true for, all pulsars showing drifting in multiple rings of emission and puts constraints on the theoretical models of pulsar emission mechanism and, favors a pan magnetospeheric radiation mechanism.

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.

We present an overview for the reionization search at GMRT. The forecast sensitivities are promising for an early detection. RFI mitigation has been successful. Several hundred hours of telescope time have already been invested in this ongoing effort, and analysis of the data is in progress.

PSR J0514-4002A is a 5-ms pulsar is located in the globular cluster NGC 1851; it belongs to a highly eccentric (e = 0.888) binary system. It is one of the earliest known examples of a numerous and fast-growing class of eccentric binary MSPs recently discovered in globular clusters. Using the GBT, we have obtained a phase-coherent timing solution for the pulsar, which includes a measurement of the rate of advance of periastron: 0.01289(4) degrees per year, which if due completely to general relativity, implies a total system mass of 2.453(14) solar masses. We also derive m_p < 1.5 solar masses and m_c > 0.96 solar masses. The companion is likely to be a massive white dwarf star.

We report new results obtained from multi-frequency observations of PSR B0826-34 with the Giant Metrewave Radio Telescope (GMRT). (1) We find no evidence of weak emission during the typical long null state of this pulsar, simultaneously at 303 and 610 MHz, as well as individually at 157, 325, 610 and 1060 MHz at separate epochs. Our limit of non-detection is at ~ 1% or better of the peak of the active state profile, and corresponds to ~ 2 mJy at 610 MHz. (2) Significant correlation in the total intensity of the individual pulses between 303 and 610 MHz is reported from the simultaneous dual frequency observations, which is indicative of the broadband nature of the emission. We also report correlation between total energy in the main pulse and inter-pulse region from the high sensitivity single frequency observations at 610 and 1060 MHz. (3) Though we find the drift pattern to be very similar in the simultaneous 303 and 610 MHz data, we observe that the drift band separation (P2) evolves significantly between these two frequencies, and in a manner opposite to the average profile evolution. In addition, we confirm the dependence of P2 on pulse longitude at 303 MHz and find indications for the same at 610 MHz. We also present results for subpulse width at different frequencies, and as well as a function of pulse longitude. (4) As a natural out-come of the simultaneous dual frequency observations, we obtain an accurate DM value, equal to 52.2(6) pc/cc, for this pulsar.

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.

We have used the Green Bank Telescope to observe the millisecond pulsar PSR J0514-4002A on 43 occasions spread over 2 years. This 5-ms pulsar is located in the globular cluster NGC 1851; it belongs to a binary system and has a highly eccentric (e = 0.888) orbit. We have obtained a phase-coherent timing solution for this object, including very precise position, spin and orbital parameters. The pulsar is located 4.6" (about 1.3 core radii) from the center of the cluster, and is likely to lie on its more distant half. The non-detection of eclipses at superior conjunction can be used, given the peculiar geometry of this system, to rule out the possibility of an extended companion. We have measured the rate of advance of periastron for this binary system to be $\dot{\omega}$ = 0.01289(4) degrees per year, which if due completely to general relativity, implies a total system mass of 2.453(14) solar masses. Given the known mass function, the pulsar mass has to be < 1.5 solar masses, and the mass of the companion has to be > 0.96 solar masses, implying that it is a heavy white dwarf. The 350-MHz flux density of this pulsar varies between 0.2 and 1.4 mJy; the origin of these variations is not known.

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.

We report the discovery of a remarkable subpulse drift pattern in the relatively less studied wide profile pulsar, B0818-41, using high sensitivity GMRT observations. We find simultaneous occurrence of three drift regions with two different drift rates: an inner region with steeper apparent drift rate flanked on each side by a region of slower apparent drift rate. Furthermore, these closely spaced drift bands always maintain a constant phase relationship. Though these drift regions have significantly different values for the measured P2, the measured P3 value is the same and equal to 18.3 P1. We interpret the unique drift pattern of this pulsar as being created by the intersection of our line of sight (LOS) with two conal rings on the polar cap of a fairly aligned rotator (inclination angle alpha ~ 11 deg), with an ``inner'' LOS geometry (impact angle beta ~ -5.4 deg). We argue that both the rings have the same values for the carousel rotation periodicity P4 and the number of sparks Nsp. We find that Nsp is 19-21 and show that it is very likely that, P4 is the same as the measured P3, making it a truly unique pulsar. We present results from simulations of the radiation pattern using the inferred parameters, that support our interpretations and reproduce the average profile as well as the observed features in the drift pattern quite well.

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.

We present the first direct evidence for turn-over in pulsar radio spectra at high frequencies. Two pulsars are now shown to have a turn-over frequency > 1GHz. We also find some evidence that the peak frequency of turn-over in pulsar spectra appears to depend on dispersion measure and pulsar age.

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.

The main aim is to investigate the possibility of a high frequency turn-over in the radio spectrum of pulsars. Using the GMRT, multi-frequency flux density measurements of several candidate pulsars have been carried out and their spectra have been extended to lower frequencies. We present the first direct evidence for turn-over in pulsar radio spectra at high frequencies. A total of 3 pulsars (including 2 new ones from this study) are now shown to have a turn-over frequency > 1 GHz, and one is shown to have a turn-over at ~600 MHz.

We revisit the phenomenon of pulse nulling using high-quality single-pulse data of PSR B1133+16 from simultaneous multifrequency observations. Observations were made at 325, 610, 1400 and 4850 MHz as part of a joint program between the European Pulsar Network and the Giant Metrewave Radio Telescope. The pulse energy time series are analysed to derive improved statistics of nulling pulses as well as to investigate the frequency dependence of the phenomenon. The pulsar is observed to be in null state for approximately 15% of the time; however, we find that nulling does not always occur simultaneously at all four frequencies of observation. We characterise this "selective nulling'' as a function of frequency, separation in frequency, and combination of frequencies. The most remarkable case is a significantly large number of nulls ($\approx$6%) at lower frequencies, that are marked by the presence of a fairly narrow emission feature at the highest frequency of 4850 MHz. We refer to these as "low frequency (LF) nulls." Our analysis shows that this high frequency emission tends to occur preferentially over a narrow range in longitude and with pulse widths typically of the order of a few milliseconds. We discuss the implications of our results for the pulsar emission mechanism in general and for the broadbandness of nulling phenomenon in particular. Our results signify the presence of an additional process of emission which does not turn off when the pulsar nulls at low frequencies, and becomes more prominent at higher frequencies. Our analysis also hints at a possible outer gap origin for this new population of pulses, and thus a likely connection to some high-energy emission processes that occur in the outer parts of the pulsar magnetosphere.

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.

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.

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.

We present new results from high sensitivity GMRT observations of PSR B0826-34. We provide a model to explain the observed subpulse drift properties of this pulsar, including the apparent reversals of the drift direction. In this model, PSR B0826-34 is close to being an aligned rotator. We solve for the emission geometry of this pulsar and show that the angle between the rotation and the magnetic axes is less than 5 deg. We see evidence for as many as 6 to 7 drifting bands in the main pulse at 318 MHz, which are part of a circulating system of about 15 spark-associated subpulse emission beams. We provide quantitative treatments of the aliasing problem and various effects of geometry. The observed drift rate is an aliased version of the true drift rate, such that a subpulse drifts to the location of the adjacent subpulse (or a multiple thereof) in about one pulsar period. We show that small variations, of the order of 3-8%, in the mean drift rate are then enough to explain the apparent reversals of drift direction. We find the mean circulation time of the drift pattern to be significantly longer than the predictions of the original RS75 model and propose an explanation for this, based on modified models with temperature regulated partial ion flow in the polar vacuum gap. From the variation of the mean subpulse separation across the main pulse window, we show that the spark pattern is not centred around the dipole axis, but around a point much closer (within a degree or so) to the rotation axis -- we discuss the implication of this.

Using the Giant Metrewave Radio Telescope (GMRT), we have discovered PSR J0514-4002A, a binary millisecond pulsar in the globular cluster NGC 1851. This pulsar has a rotational period of 4.99 ms and the most eccentric pulsar orbit yet found: e = 0.89. The orbital period is 18.8 days, and companion has a minimum mass of 0.9 M_sun; its nature is presently unclear. After accreting matter from a low-mass stellar companion, this pulsar exchanged it for its more massive present companion. This system presents the strongest evidence to date of such a process.

We report the discovery of the binary millisecond pulsar J0514-4002A, which is the first known pulsar in the globular cluster NGC 1851 and the first pulsar discovered using the Giant Metrewave Radio Telescope (GMRT). The pulsar has a rotational period of 4.99 ms, an orbital period of 18.8 days, and the most eccentric pulsar orbit yet measured (e = 0.89). The companion has a minimum mass of 0.9 M_sun and its nature is presently unclear. After accreting matter from a low-mass companion star which spun it up to a (few) millisecond spin period, the pulsar eventually exchanged the low-mass star for its more massive present companion. This is exactly the same process that could form a system containing a millisecond pulsar and a black hole; the discovery of NGC 1851A demonstrates that such systems might exist in the Universe, provided that stellar mass black holes exist in globular clusters.

In this paper we demonstrate that a large, unexplored reservoir of information about pulsar emission exists, that is directly linked to the radiating particles and their radiation process: We present a study of flux density measurements of individual pulses simultaneously observed at four different frequencies. Correcting for effects caused by the interstellar medium, we derive intrinsic flux density spectra of individual radio pulses observed at several frequencies for the first time. Pulsar B0329+54 was observed at 238, 626, 1412 and 4850 MHz, while observations of PSR B1133+16 were made at 341, 626, 1412 and 4850 MHz. We derive intrinsic pulse-to-pulse modulation indices which show a minimum around 1 GHz. Correlations between the flux densities of different frequency pairs worsen as the frequency separation widens and also tend to be worse for outer profile components. The single pulse spectra of PSR B0329+54 resemble the spectra of the integrated profile. However, the spectral index distributions for the single pulses of PSR B1133+16 show significant deviations from a Gaussian. This asymmetry is caused by very strong pulses with flux densities exceeding the mean value by more than a factor of ten. These strong pulses occur preferentially at the trailing edge of the leading component and appear to be broadband in most cases. Their properties are similar to those of so-called giant pulses, suggesting that these phenomena are related.

We have conducted a detailed analysis of the emission geometry of a handful of radio pulsars that have prominent, multiple-component profiles at meter wavelengths. From careful determination of the total number of emission components and their locations in pulse longitude, we find that all of the six pulsars show clear evidence for retardation and aberration effects in the conal emission beams. Using this information, coupled with a dipolar field geometry, we obtain estimates of the height and transverse location in the magnetosphere, for each of the emitting cones in these pulsars. These results support our earlier conclusions for PSR B0329+54 in that we find successive outer cones (in cases of multi-cone pulsars) being emitted at higher altitudes in the magnetosphere. The range of inferred heights is from ~200 to ~2200 km. The set of ``active'' field lines from which the conal emissions originate are located in the region from ~0.22 to ~0.74 of the polar cap radius. At the neutron star surface, these conal rings map to radii of a few to several tens of meters and the separation between successive rings is about 10 to 20 meters. We discuss the implications of these findings for the understanding of the pulsar emission geometry and for current theories and models of the emission mechanism.

In this paper, we present a study of orthogonal polarization modes in the radio emission of PSR B1133+16, conducted within the frame of simultaneous, multi-frequency, single-pulse observations. Simultaneously observing at two frequencies (1.41 GHz and 4.85 GHz) provides the means to study the bandwidth of polarization features such as the polarization position angle. We find two main results. First, that there is a high degree of correlation between the polarization modes at the two frequencies. Secondly, the modes occur more equally and the fractional linear polarization decreases towards higher frequencies. We discuss this frequency evolution and propose propagation effects in the pulsar magnetosphere as its origin.

Recent pulsar scintillation measurements from Ooty, in conjunction with those from Parkes and other large radio telescopes, are used for a systematic investigation of the Local Interstellar Medium (LISM) towards the general direction of the Loop I Bubble. For several pulsars, clear evidence is found for an enhanced level of scattering which is over and above what can be accounted for by the enhanced scattering model for the Local Bubble. These results are interpreted in terms of enhanced scattering due to turbulent plasma associated with the Loop I shell. Useful constraints are obtained for the scattering properties of the shell. The inferred value for the scattering measure for the Loop I shell is found to be ~0.3 pc m^{-20/3}. Assuming a shell thickness ~5-10 pc, this implies an average strength of scattering in the shell that is ~100-200 times larger than that in the ambient ISM. An alternative explanation, where the enhanced level of scattering is due to a possible "interaction zone" between the Local Bubble and the Loop I Bubble, is also considered; it is found to be somewhat less satisfactory in explaining the observations. The best fit value for the scattering measure for such an interaction zone region is estimated to be ~1.1 pc m^{-20/3}. Further, several pulsars beyond ~1 kpc are found to show enhanced levels of scattering over and above that expected from this "two-bubble model." For some of the low-latitude pulsars, this is found to be due to enhanced scattering from plasma inside the intervening Sagittarius spiral arm. We discuss the implications of our results for the interpretation of scintillation data and for the general understanding of the LISM.

We report on new results from simultaneous, dual frequency, single pulse observation of PSR B0329+54 using the Giant Metrewave Radio Telescope. We find that the longitude separation of subpulses at two different frequencies (238 and 612 MHz) is less than that for the corresponding components in the average profile. A similar behaviour has been noticed before in a number of pulsars. We argue that subpulses are emitted within narrow flux tubes of the dipolar field lines and that the mean pulsar beam has a conal structure. In such a model the longitudes of profile components are determined by the intersection of the line of sight trajectory with subpulse-associated emission beams. Thus, we show that the difference in the frequency dependence of subpulse and profile component longitudes is a natural property of the conal model of pulsar emission beam. We support our conclusions by numerical modelling of pulsar emission, using the known parameters for this pulsar, which produce results that agree very well with our dual frequency observations.

We have analyzed high-quality single pulse data of PSR B0329+54 at 325 MHz and 606 MHz to study the structure of the emission beam. Using the window-threshold technique, which is suitable for detecting weak emission components, we have detected 4 additional emission components in the pulse window. Three of these are new components and the fourth is a confirmation of a recently proposed component. Hence PSR B0329+54 is now known to have 9 emission components - the highest among all known pulsars. The distribution of the pulse components around the central core component indicates that the emission beam consists of four nested cones. The asymmetry in the location of the conal components in the leading versus trailing parts of the profile is interpreted as being due to aberration and retardation in the pulsar magnetosphere. These measurements allow us to determine the precise location of the 4 conal rings of emission. We find that the successive outer cones are emitted at higher altitudes in the magnetosphere. Further, for any given cone, the emission height at the lower frequency is found to be more than that at the higher frequency. The inferred heights range from ~160 km to ~1150 km. The set of ``active'' field lines, from which most of the conal radiation appears to originate, are found to be confined to a region located within ~0.5 to ~0.6 of the polar cap radius. We discuss the implications of our new findings on our understanding of the pulsar emission geometry and its impact on the emission mechanisms.

Results from new observations of pulsars using the Ooty Radio Telescope (ORT) are used for investigating the structure of the Local Interstellar Medium (LISM). The observations show anomalous scintillation effects towards several nearby pulsars, and these are modeled in terms of large-scale spatial inhomogeneities in the distribution of plasma density fluctuations in the LISM. A 3-component model, where the Solar neighbourhood is surrounded by a shell of enhanced plasma turbulence, is proposed for the LISM. The inferred scattering structure is strikingly similar to the Local Bubble. Further, analysis based on recent scintillation measurements show evidence for enhanced scattering towards pulsars located in the general direction of the Loop I Superbubble. The model for the LISM has been further extended by incorporating the scattering due to turbulent plasma associated with Loop I.

Refractive Interstellar Scintillation (RISS) effects on pulsar signals are powerful techniques for discriminating between different models that have been proposed for the power spectrum of plasma density fluctuations in the Interstellar Medium. The nature of the spectrum is considered to be a major input for understanding the underlying mechanism of interstellar plasma turbulence. Data from our long-term pulsar scintillation observations using the Ooty Radio Telescope (ORT) at 327 MHz are used to investigate the nature of the spectrum in the Local Interstellar Medium. The emerging picture is that the spectrum is Kolmogorov-like ($ \beta \approx 11/3 $) in the spatial scale range $\sim 10^6$ m to $\sim 10^{11}$ m, and there is evidence for excess power at larger spatial scales ($ \sim 10^12 - 10^13 $ m).

To study the structure of emission beam, we have analysed the single pulse data of PSR B0329+-54 at 325 and 606 MHz. In order to unambiguously detect the weak emission components, we have developed a new data analysis technique, which we term ``window-thresholding''. By applying this technique to the data, we have detected three new emission components, and also confirmed the presence of a component which was proposed earlier. Hence our analysis indicates that PSR B0329+-54 has nine components, which is among the highest of all the known pulsars. The symmetric distribution of pulse components about the pulse centre, indicates that the emission beam is conal.

Refraction of pulsar radiation by electron density irregularities in the interstellar medium sometimes produces multiple imaging of pulsars which can lead to periodic oscillations of intensity in pulsar dynamic spectra records. Such events can be used as tools to resolve the emission regions in pulsar magnetospheres. Here we describe results from the recent observation of a double imaging event for PSR B1133+16, which place fairly tight constraints on the location of the emission regions. Our analysis constrains the location of the scattering object to the shell of the Local Bubble. The phase of the oscillations shows significant variations across the pulse. The minimum value for the transverse separation of the emitting regions at the two edges of the pulse is inferred to be 3 x 10^5 m. This translates to a minimum emission altitude of 2.6 x 10^6 m. The non-monotonic variations of the fringe phase with pulse longitude are interpreted as variations of the altitude of the emission regions for the orthogonal polarization modes of this pulsar. This is in agreement with theories where propagation effects, such as refraction, are responsible for the orthogonal modes.

Refractive interstellar scintillation (RISS) is thought to be the cause behind a variety of phenomena seen at radio wavelengths in pulsars and compact radio sources. Though there is substantial observational data to support several consequences of it, the quantitative predictions from theories have not been thoroughly tested. In this paper, data from our long-term scintillation study of 18 pulsars are used to test the predictions. The fluctuations of decorrelation bandwidth ($\nu_d$), scintillation time scale ($\tau_d$) and flux density (F) are examined for their cross-correlations and compared with the predictions. The theory predicts a strong correlation between $\nu_d$ and $\tau_d$, and strong anti-correlations between $\nu_d$ and F, and $\tau_d$ and F. For 5 pulsars, we see a reasonable agreement. There is considerable difficulty in reconciling the results for the rest of the pulsars. Our analysis shows the underlying noise sources can sometimes reduce the correlation, but cannot cause an absence of correlation. It is also unlikely that the poor flux correlations arise from a hitherto unrecognized intrinsic flux variations. For PSR B0834+06, which shows anomalous behaviour of persistent drift slopes, positive correlation is found between $\tau_d$ and the drift-corrected $\nu_d$. Many pulsars show an anti-correlation between $\nu_d$ and the drift slope, and this is in accordance with the simple models of RISS. The detections of correlated variations of observables and a reasonable agreement between the predicted and measured correlations for some pulsars confirm RISS as the primary cause of the observed fluctuations. However, the complexity seen with the detailed results suggests the necessity of more comprehensive theoretical treatments for describing refractive fluctuations and their correlations.

Refractive scintillation effects are powerful techniques for discriminating between different models proposed for the electron density fluctuation spectrum in the ISM. Data from our long-term scintillation study of 18 pulsars in the DM range 3-35 pc cm^{-3} (Paper I) are used to investigate two important observable effects of refractive scintillation: (i) modulations of DISS observables and flux, and (ii) drifting bands. Our data provide simultaneous measurements of decorrelation bandwidth, scintillation time scale, flux density and drift slope of patterns. The observed modulations of the first 3 are compared with the theoretical predictions, and constraints are placed on the density spectral form. The measured modulation indices are larger than predicted by a Kolmogorov form. Measurements of drift slopes and DISS parameters are used to independently estimate the slope of the density spectrum, which is found to be consistent with a Kolmogorov form for several pulsars.The contradictory results from these 2 independent methods are not reconcilable with the simple theoretical models. Our observations show anomalous behaviour like persistent drifting bands for some pulsars.This can be interpreted as an excess power in the low wavenumber range (~10^{-12} - 10^{-13} m^{-1}) compared to the Kolmogorov expectations or the existence of localized density structures.The results from our observations and those from the earlier studies are examined for understanding the overall nature of the density spectrum.The emerging picture is a Kolmogorov-like form in the wavenumber range ~10^{-6} m^{-1} to ~10^{-11} m^{-1}, which either steepens or has a bump near ~ 10^{-12} - 10^{-13} m^{-1}. We also discuss possible implications of our results for the current theoretical models.

We report long-term scintillation observations of 18 pulsars in the DM range 3-35 pc cm^{-3} carried out during 1993-1995 using the Ooty Radio Telescope at 327 MHz. These observations were made with the aim of studying refractive effects in pulsar scintillation, and obtaining reliable estimates of diffractive and refractive scintillation properties. Dynamic scintillation spectra of pulsars were regularly monitored at 10-90 epochs spanning 100-1000 days. Significant changes are observed in the dynamic spectra over time scales as short as a few days. Large-amplitude fluctuations are observed in quantities such as decorrelation bandwidth, scintillation time scale, drift rate, and flux density. Several pulsars show organized features such as drifting bands in a highly pronounced manner. For some pulsars, gradual and systematic variations are seen in the drift rate of patterns which undergo several sign reversals during the observing time spans. Anomalous behaviour such as persistent drifts lasting over many months are seen for PSRs B0834+06 & B1919+21. Four pulsars were studied for 2-4 well separated observing sessions, each lasting over ~ 100 days. In some cases, significant variations are seen in the average scintillation properties and/or flux densities between successive observing sessions. From our data, we have been able to obtain more accurate and reliable estimates of scintillation properties and flux densities than those from the earlier observations, by averaging out the fluctuations due to refractive scintillation effects. These measurements are used to derive parameters such as the strength of scattering and scintillation speeds. The present measurements are compared with earlier measurements and the long-term stability of scintillation properties and flux densities is discussed.

We present here the results from an extensive scintillation study of twenty pulsars in the dispersion measure (DM) range 3 - 35 pc cm^-3 caried out using the Ooty Radio Telescope (ORT) at 327 MHz, to investigate the distribution of ionized material in the local interstellar medium. Observations were made during the period January 1993 to August 1995, in which the dynamic scintillation spectra of these pulsars were regularly monitored over 10 - 90 epochs spanning 100 days. Reliable and accurate estimates of strengths of scattering have been deduced from the scintillation parameters averaged out for their long-term fluctuations arising from refractive scintillation (RISS) effects. Our analysis reveals several anomalies in the scattering strength, which suggest tht the distribution of scattering material in the Solar neighborhood is not uniform. We have modelled these anomalous scattering effects in terms of inhomogeneities in the distribution of electron dnsity fluctuations in the local interstellar medium (LISM). Our model suggests the presence of a low density bubble surrounded by a shell of much higher density fluctuations. We are able to put constraints on geometrical and scattering properties of such a structure, and find it to be morphologically similar to the Local Bubble known from other studies.

We present here the results from an extensive scintillation study of twenty pulsars in the dispersion measure (DM) range 3 - 35 pc cm^-3 carried out using the Ooty Radio Telescope, to investigate the distribution of ionized material in the local interstellar medium (LISM). Our analysis reveals several anomalies in the scattering strength, which suggest that the distribution of scattering material in the Solar neighborhood is not uniform. Our model suggests the presence of a low density bubble surrounded by a shell of much higher density fluctuations. We are able to put some constraints on geometrical and scattering properties of such a structure, and find it to be morphologically similar to the Local Bubble known from other studies.