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 describe system verification tests and early science results from the pulsar processor (PTUSE) developed for the newly-commissioned 64-dish SARAO MeerKAT radio telescope in South Africa. MeerKAT is a high-gain (~2.8 K/Jy) low-system temperature (~18 K at 20cm) radio array that currently operates from 580-1670 MHz and can produce tied-array beams suitable for pulsar observations. This paper presents results from the MeerTime Large Survey Project and commissioning tests with PTUSE. Highlights include observations of the double pulsar J0737-3039A, pulse profiles from 34 millisecond pulsars from a single 2.5h observation of the Globular cluster Terzan 5, the rotation measure of Ter5O, a 420-sigma giant pulse from the Large Magellanic Cloud pulsar PSR J0540-6919, and nulling identified in the slow pulsar PSR J0633-2015. One of the key design specifications for MeerKAT was absolute timing errors of less than 5 ns using their novel precise time system. Our timing of two bright millisecond pulsars confirm that MeerKAT delivers exceptional timing. PSR J2241-5236 exhibits a jitter limit of <4 ns per hour whilst timing of PSR J1909-3744 over almost 11 months yields an rms residual of 66 ns with only 4 min integrations. Our results confirm that the MeerKAT is an exceptional pulsar telescope. The array can be split into four separate sub-arrays to time over 1000 pulsars per day and the future deployment of S-band (1750-3500 MHz) receivers will further enhance its capabilities.

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

We report the detection of a bright fast radio burst, FRB\,191108, with Apertif on the Westerbork Synthesis Radio Telescope (WSRT). The interferometer allows us to localise the FRB to a narrow $5\arcsec\times7\arcmin$ ellipse by employing both multibeam information within the Apertif phased-array feed (PAF) beam pattern, and across different tied-array beams. The resulting sight line passes close to Local Group galaxy M33, with an impact parameter of only 18\,kpc with respect to the core. It also traverses the much larger circumgalactic medium of M31, the Andromeda Galaxy. We find that the shared plasma of the Local Group galaxies could contribute $\sim$10\% of its dispersion measure of 588\,pc\,cm$^{-3}$, but not detectable scintillation, temporal scattering, or significant Faraday rotation in this case. FRB\,191108 has a Faraday rotation measure of +474\,$\pm\,3$\,rad\,m$^{-2}$, which is too large to be explained by either the Milky Way or the intergalactic medium. This indicates a dense local magneto-ionic environment in the source host galaxy---as has been argued for other FRBs. We found no accompanying persistent radio sources in the Apertif imaging survey data.

Repeating fast radio bursts (FRBs) present excellent opportunities to identify FRB progenitors and host environments, as well as decipher the underlying emission mechanism. Detailed studies of repeating FRBs might also hold clues to the origin of FRBs as a population. We aim to detect the first two repeating FRBs: FRB 121102 (R1) and FRB 180814.J0422+73 (R2), and characterise their repeat statistics. We also want to significantly improve the sky localisation of R2. We use the Westerbork Synthesis Radio Telescope to conduct extensive follow-up of these two repeating FRBs. The new phased-array feed system, Apertif, allows covering the entire sky position uncertainty of R2 with fine spatial resolution in one pointing. We characterise the energy distribution and the clustering of detected R1 bursts. We detected 30 bursts from R1. Our measurements indicate a dispersion measure of 563.5(2) pc cm$^{-3}$, suggesting a significant increase in DM over the past few years. We place an upper limit of 8% on the linear polarisation fraction of the brightest burst. We did not detect any bursts from R2. A single power-law might not fit the R1 burst energy distribution across the full energy range or widely separated detections. Our observations provide improved constraints on the clustering of R1 bursts. Our stringent upper limits on the linear polarisation fraction imply a significant depolarisation, either intrinsic to the emission mechanism or caused by the intervening medium, at 1400 MHz that is not observed at higher frequencies. The non-detection of any bursts from R2 implies either a highly clustered nature of the bursts, a steep spectral index, or a combination of both. Alternatively, R2 has turned off completely, either permanently or for an extended period of time.

Binary millisecond pulsars (MSPs) provide several opportunities for research of fundamental physics. However, finding them can be challenging. Several subdwarf B (sdB) binary systems with possible neutron star companions have been identified, allowing us to perform a targeted search for MSPs within these systems. Six sdBs with companions in the neutron star mass range, as determined from their optical light curves, were observed with the Green Bank and Westerbork radio telescopes. The data were searched for periodic signals as well as single pulses. No radio pulsations from sdB systems were detected, down to an average sensitivity limit of 0.11 mJy. We did, however, discover a pulsar in the field of sdB HE0532-4503. Follow-up observations with the Giant Metrewave Radio Telescope showed that this pulsar, J0533-4524, is not spatially coincident with the sdB system. The pulsar has a relatively low magnetic field but still emits giant pulses. We place an upper limit of three to the number of radio pulsars in the six sdB systems. The non-detections may be explained by a combination of the MSP beaming fraction, luminosity, and a recycling fraction <0.5. Alternatively, the assumption of co-rotation between the MSP and sdB may break down, which implies the systems are more edge-on than previously thought. This would shift the predicted companion masses into the white dwarf range. It would also explain the relative lack of edge-on sdB systems with massive companions.

XTE J1810-197 (PSR J1809-1943) was the first ever magnetar which was found to emit transient radio emission. It has recently undergone another radio and high-energy outburst. This is only the second radio outburst that has been observed from this source. We observed J1810-197 soon after its recent radio outburst at low radio frequencies using the Giant Metrewave Radio Telescope. We present the 650 MHz flux density evolution of the source in the early phases of the outburst, and its radio spectrum down to frequencies as low as 300 MHz. The magnetar also exhibits radio emission in the form of strong, narrow bursts. We show that the bursts have a characteristic intrinsic width of the order of 0.5-0.7 ms, and discuss their properties in the context of giant pulses and giant micropulses from other pulsars. We also show that the bursts exhibit spectral structures which cannot be explained by interstellar propagation effects. These structures might indicate a phenomenological link with the repeating fast radio bursts which also show interesting, more detailed frequency structures. While the spectral structures are particularly noticeable in the early phases of the outburst, these seem to be less prominent as well as less frequent in the later phases, suggesting an evolution of the underlying cause of these spectral structures.

Pulsars typically exhibit radio emission in the form of narrow pulses originated from confined regions of their magnetospheres. A potential presence of magnetospherically originated emission outside this region, the so-called off-pulse emission, would challenge the existing theories. Detection of significant off-pulse emission has been reported so far from only two pulsars, B0525+21 and B2045-16, at 325 and 610 MHz. However, the nature of this newly uncovered off-pulse emission remains unclear. To probe its origin we conducted very high resolution radio observations of B0525+21 and B2045-16 with the European VLBI Network (EVN) at 1.39 GHz. Whereas the pulsed emission is detected at a level consistent with previous observations, we report absence of any off-pulse emission above $42$ and $96\ \mathrm{\mu Jy\ beam^{-1}}$ (three times the rms noise levels) for B0525+21 and B2045-16, respectively. Our stringent upper limits imply the off-pulse emission to be less than $0.4$ and $0.3\%$ of the period-averaged pulsed flux density, i.e., much fainter than the previously suggested values of $1$-$10\%$. Since the EVN data are most sensitive to extremely compact angular scales, our results suggest a non-magnetospheric origin for the previously reported off-pulse emission. Presence of extended emission that is resolved out on these milliarcsecond scales still remains plausible. In this case, we constrain the emission to arise from structures with sizes of $\sim (0.61$-$19) \times 10^3\ \mathrm{au}$ for B0525+21 and $\sim (0.48$-$8.3) \times 10^3\ \mathrm{au}$ for B2045-16. These constraints might indicate that the two pulsars are accompanied by compact bow-shock pulsar wind nebulae. Future observations probing intermediate angular scales ($\sim 0.1$-$5\ \mathrm{arcsec}$) will help in clarifying the actual origin of the off-pulse emission.

We present multi-frequency scatter broadening evolution of 29 pulsars observed with the LOw Frequency ARray (LOFAR) and Long Wavelength Array (LWA). We conducted new observations using LOFAR Low Band Antennae (LBA) as well as utilized the archival data from LOFAR and LWA. This study has increased the total of all multi-frequency or wide-band scattering measurements up to a dispersion measure (DM) of 150~pc\,cm$^{-3}$ by 60\%. The scatter broadening timescale ($\tau_{sc}$) measurements at different frequencies are often combined by scaling them to a common reference frequency of 1\,GHz. Using our data, we show that the $\tau_{sc}$--DM variations are best fitted for reference frequencies close to 200--300\,MHz, and scaling to higher or lower frequencies results in significantly more scatter in data. We suggest that this effect might indicate a frequency dependence of the scatter broadening scaling index ($\alpha$). However, a selection bias due to our chosen observing frequencies can not be ruled out with the current data set. Our data did not favour any particular model of the DM -- $\tau_{sc}$ relations, and we do not see a statistically significant break at the low DM range in this relation. The turbulence spectral index ($\beta$) is found to be steeper than that is expected from a Kolmogorov spectrum. This indicates that the local ISM turbulence may have a low wave-number cutoff or presence of large scale inhomogeneities in the line of sight to some of the reported pulsars.

Subpulse modulation has been regarded as one of the most insightful and intriguing aspects of pulsar radio emission. This phenomenon is generally explained by the presence of a carousel of sparks in the polar acceleration gap region that rotates around the magnetic axis due to the E$\times$B drift. While there have been extensive single pulse studies, geometric signatures of the underlying carousel, or lack thereof, in simultaneous multi-frequency observations have remained largely unexplored. This work presents a theoretical account of such expected signatures, particularly that of a geometry induced phase-offset in subpulse modulation, including various formulae that can be readily applied to observations. A notable result is a method to resolve aliasing in the measured subpulse modulation period without relying on knowledge of the viewing geometry parameters. It is also shown in detail that the geometry induced phase-offset enables critical tests of various observed phenomena as well as proposed hypotheses, e.g., multi-altitude emission, magnetic field twisting, pseudo-nulls, etc., in addition to that of the carousel model itself. Finally, a detailed analysis of a 327 MHz pulse-sequence of PSR B1237+25 is presented as a case study to test the single-frequency multi-altitude emission scenario. The analysis provides a firm evidence of inner and outer conal components of this pulsar to have originated from the same carousel of subbeams and emitted at different heights.

We report on a time-domain search for pulsars in 44 steep spectrum radio sources originally identified from recent imaging surveys. The time-domain search was conducted at 327 MHz using the Ooty radio telescope, and utilized a semi-coherent dedispersion scheme retaining the sensitivity even for sub-millisecond periods up to reasonably high dispersion measures. No new pulsars were found. We discuss the nature of these steep spectrum sources and argue that majority of the sources in our sample should either be pulsars or a new category of Galactic sources. Several possibilities that could hinder detection of these sources as pulsars, including anomalously high scattering or alignment of the rotation and magnetic axes, are discussed in detail, and we suggest unconventional search methods to further probe these possibilities.

Sub-pulse drifting has been regarded as one of the most insightful aspects of the pulsar radio emission. The phenomenon is generally explained with a system of emission sub-beams rotating around the magnetic axis, originating from a carousel of sparks near the pulsar surface (the carousel model). Since the observed radio emission at different frequencies is generated at different altitudes in the pulsar magnetosphere, corresponding sampling of the carousel on the polar cap differs slightly in magnetic latitude. When this aspect is considered, it is shown here that the carousel model predicts important observable effects in multi-frequency or wide-band observations. Also presented here are brief mentions of how this aspect can be exploited to probe the electrodynamics in the polar cap by estimating various physical quantities, and correctly interpret various carousel related phenomena, in addition to test the carousel model itself.

With the installation of a new phased array system called Apertif, the instantaneous field of view of the Westerbork Synthesis Radio Telescope (WSRT) has increased to 8.7$\,$deg$^2$. This system has turned the WSRT in to an highly effective telescope to conduct Fast Radio Burst (FRB) and pulsar surveys. To exploit this advantage, an advanced and real-time backend, called the Apertif Radio Transient System (ARTS), is being developed and commissioned at the WSRT. In addition to the real-time detection of FRBs, ARTS will localize the events to about 1/2600 of the field of view --- essential information for identifying the nature of FRBs. ARTS will also trigger real-time follow up with LOFAR of newly detected FRBs, to achieve localization at arcsecond precision. We review the upcoming time-domain surveys with Apertif, and present the current status of the ongoing commissioning of the time domain capabilities of Apertif.

Although originally discovered as a radio-quiet gamma-ray pulsar, J1732-3131 has exhibited intriguing detections at decameter wavelengths. We report an extensive follow-up of the pulsar at 327 MHz with the Ooty radio telescope. Using the previously observed radio characteristics, and with an effective integration time of 60 hrs, we present a detection of the pulsar at a confidence level of 99.82%. The 327 MHz mean flux density is estimated to be 0.5-0.8 mJy, which establishes the pulsar to be a steep spectrum source and one of the least-luminous pulsars known to date. We also phase-aligned the radio and gamma-ray profiles of the pulsar, and measured the phase-offset between the main peaks in the two profiles to be 0.24$\pm$0.06. We discuss the observed phase-offset in the context of various trends exhibited by the radio-loud gamma-ray pulsar population, and suggest that the gamma-ray emission from J1732-3131 is best explained by outer magnetosphere models. Details of our analysis leading to the pulsar detection, and measurements of various parameters and their implications relevant to the pulsar's emission mechanism are presented.

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.

We present the results obtained from imaging observations, and search for radio pulsations towards the magnetar SGR J1935+2154 made using the Giant Metrewave Radio Telescope, and the Ooty Radio Telescope. We present the high resolution radio image of the supernova remnant (SNR) G57.2+0.8, which is positionally associated with SGR J1935+2154. We did not detect significant periodic radio pulsations from the magnetar, with 8$\sigma$ upper limits on its flux density of 0.4, and 0.2 mJy at 326.5, and 610 MHz, respectively, for an assumed duty cycle of 10\%. The corresponding 6$\sigma$ upper limits at the two frequencies for any burst emission with an assumed width of 10 ms are 0.5 Jy, and 63 mJy, respectively. No continuum radio point source was detected at the position of SGR J1935+2154 with a 3$\sigma$ upper limit of 1.2 mJy. We also did not detect significant diffuse radio emission in a radius of 70 arc seconds in coincidence with the diffuse X-ray emission reported recently, with a 3$\sigma$ upper limit of 4.5 mJy. Using the archival HI spectra, we estimate the distance of SNR G57.2+0.8 to be 11.7 $\pm$ 2.8 kpc. Based on measured HI column density (N$_H$) along this line of sight, we argue that the magnetar could be physically associated with SNR G57.2+0.8. Based on present data, we can not rule out either a pulsar wind nebula or a dust scattering halo origin for the diffuse X-ray emission seen around the magnetar.

We report discovery of several energetic radio bursts at 34 MHz, using the Gauribidanur radio telescope. The radio bursts exhibit two important properties associated with the propagation of astronomical signals through the interstellar medium: (i) frequency dependent dispersive delays across the observing bandwidth, and (ii) Faraday rotation of the plane of linear polarization. These bursts sample a range of dispersion measures (DM; 1.4--3.6$~{\rm pc}~{\rm cm}^{-3}$), and show DM-variation at timescales of the order of a minute. Using groups of bursts having a consistent DM, we show that the bursts have originated from the radio-quiet gamma-ray pulsar Geminga. Detection of these bursts supports the existence of occasional radio emission from Geminga. The rare occurrence of these bursts, and the short timescale variation in their DM (if really caused by the intervening medium or the pulsar magnetosphere), might provide clues as to why the pulsar has not been detected in earlier sensitive searches. We present details of the observations and search procedure used to discover these bursts, a detailed analysis of their properties, and evidences of these bursts being associated with Geminga pulsar, and discuss briefly the possible emission mechanism of these bursts.

We report the results of (a) extensive follow-up observations of the gamma-ray pulsar J1732-3131 that has been recently detected at decameter wavelengths, and (b) deep searches for counterparts of 9 other radio-quiet gamma-ray pulsars at 34 MHz, using the Gauribidanur radio telescope. No periodic signal from J1732-3131 could be detected above a detection threshold of $8\sigma$, even with an effective integration time of more than 40 hours. However, the average profile obtained by combining data from several epochs, at a dispersion measure of 15.44 pc/cc, is found to be consistent with that from the earlier detection of this pulsar at a confidence level of 99.2 %. We present this consistency between the two profiles as an evidence that J1732-3131 is a faint radio pulsar with an average flux density of 200--400 mJy at 34 MHz. Detection sensitivity of our deep searches, despite the extremely bright sky background at such low frequencies, is generally comparable to that of higher frequency searches for these pulsars, when scaled using reasonable assumptions about the underlying pulsar spectrum. We provide details of our deep searches, and put stringent upper limits on the decameter wavelength flux densities of several radio-quiet gamma-ray pulsars.

The origin of two distinct pairs of conal emission components in pulsars, associated with the "outer" and the "inner" emission cones, as well as the marked difference in their observed spectral properties, is poorly understood. The sub-pulse modulation in the corresponding conal components, if mapped back to the underlying system of sub-beams rotating around the magnetic axis in the polar cap, as envisioned by Ruderman & Sutherland (1975), provides a potential way to investigate the emission morphologies in the two conal regions, and more importantly, any inter-relationship between them. The bright pulsar B1237+25 with its special viewing geometry where the sightline traverses almost through the magnetic axis, along with a rich variety in pulse-to-pulse fluctuations, provides an excellent, but challenging opportunity to map the underlying emission patterns across the full transverse slice of its polar emission region. We present here our analysis on a number of pulse-sequences from this star to map and study any relationship between the underlying patterns responsible for emission in the two pairs of presumed conal-components and a core-component of this pulsar. The results from our correlation analysis of the two conal emission patterns strongly support the view that the two cones of this pulsar (the outer and the inner cone) originate from a common system of sub-beams. We also see evidence for a twist in the emission columns, most likely associated with a corresponding twist in the magnetic field structure. We discuss these results, and their implications, including a possibility that the core component of this pulsar shares its origin partly with the conal counterparts.

We report the design and development of a self-contained multi-band receiver (MBR) system, intended for use with a single large aperture to facilitate sensitive & high time-resolution observations simultaneously in 10 discrete frequency bands sampling a wide spectral span (100-1500 MHz) in a nearly log-periodic fashion. The development of this system was primarily motivated by need for tomographic studies of pulsar polar emission regions. Although the system design is optimized for the primary goal, it is also suited for several other interesting astronomical investigations. The system consists of a dual-polarization multi-band feed (with discrete responses corresponding to the 10 bands pre-selected as relatively RFI-free), a common wide-band RF front-end, and independent back-end receiver chains for the 10 individual sub-bands. The raw voltage time-sequences corresponding to 16 MHz bandwidth each for the two linear polarization channels and the 10 bands, are recorded at the Nyquist rate simultaneously. We present the preliminary results from the tests and pulsar observations carried out with the Green Bank Telescope using this receiver. The system performance implied by these results, and possible improvements are also briefly discussed.

Broadband antenna feeds are of particular interest to existing and future radio telescopes for multi-frequency studies of astronomical sources. Although a 1:15 range in frequency is difficult to achieve, the well-known Eleven feed design offers a relatively uniform response over such a range, and reasonably well-matched responses in E & H planes. However, given the severe Radio Frequency Interference in several bands over such wide spectral range, one desires to selectively reject the corresponding bands. With this view, we have explored the possibilities of having a multi-band feed antenna spanning a wide frequency range, but which would have good response only in a number of pre-selected (relatively) RFI-free windows (for a particular telescope-site). The designs we have investigated use the basic configuration of pairs of dipoles as in the Eleven feed, but use simple wire dipoles instead of folded dipoles used in the latter. From our study of the two designs we have investigated, we find that the design with feed-lines constructed using co-axial lines shows good rejection in the unwanted parts of the spectrum and control over the locations of resonant bands.

We report our search for and a possible detection of periodic radio pulses at 34.5 MHz from the Fermi-LAT pulsar J1732-3131. The candidate detection has been possible in only one of the many sessions of observations made with the low-frequency array at Gauribidanur, India, when the otherwise radio weak pulsar may have apparently brightened many folds. The candidate dispersion measure along the sight-line, based on the broad periodic profiles from about 20 minutes of data, is estimated to be 15.44 +/-0.32 pc/cc. We present the details of our periodic & single-pulse search, and discuss the results and their implications relevant to both, the pulsar and the intervening medium.

Although several existing and upcoming telescopes have imaging as their primary mode, they also have a sensitive phased-array mode with a multiple-beam forming capability enabling high time resolution studies of several types of objects, including pulsars. For example, the potentially wide coverage in frequency, combined with its collecting area, makes the MWA-LFD a unique instrument for low-frequency detection and studies of pulsars and transients. A software data-processing pipeline is being developed by the Raman Research Institute for this purpose. We describe the various issues relevant to the detection strategies, illustrated with real data at low radio frequencies.