Search

High-precision pulsar timing is highly dependent on precise and accurate modeling of any effects that impact the data. It was shown that commonly used Solar Wind models do not accurately account for variability in the amplitude of the Solar wind on both short and long time scales. In this study, we test and validate a new, cutting-edge Solar wind modeling method included in the \texttt{enterprise} software suite through extended simulations, and we apply it to investigate temporal variability in LOFAR data. Our model testing scheme in itself provides an invaluable asset for pulsar timing array (PTA) experiments. As improperly accounting for the solar wind signature in pulsar data can induce false-positive signals, it is of fundamental importance to include in any such investigations. We employ a Bayesian approach utilizing a continuously varying Gaussian process to model the solar wind referred to as Solar Wind Gaussian Process (SWGP). We conduct noise analysis on eight pulsars from the LOFAR dataset with most pulsars having a timespan of $\sim 11$ years encompassing one full solar activity cycle. Our analysis reveals a strong correlation between the electron density at 1 AU and the ecliptic latitude (ELAT) of the pulsar. Pulsars with $|ELAT|< 3^{\circ}$ exhibit significantly higher average electron densities. We observe distinct temporal patterns in electron densities in different pulsars. In particular, pulsars within $|ELAT|< 3^{\circ}$ exhibit similar temporal variations, while the electron densities of those outside this range correlate with the solar activity cycle. The continuous variability in electron density offered in this model represents a substantial improvement over previous models, which assume a single value for piece-wise bins of time. This advancement holds promise for solar wind modeling in future International Pulsar Timing Array data combinations., Comment: Accepted in Section 9. Sun and the Heliosphere, Astronomy and Astrophysics

A pulsar's pulse profile gets broadened at low frequencies due to dispersion along the line of sight or due to multi-path propagation. The dynamic nature of the interstellar medium makes both of these effects time-dependent and introduces slowly varying time delays in the measured times-of-arrival similar to those introduced by passing gravitational waves. In this article, we present a new method to correct for such delays by obtaining unbiased dispersion measure (DM) measurements by using low-frequency estimates of the scattering parameters. We evaluate this method by comparing the obtained DM estimates with those, where scatter-broadening is ignored using simulated data. A bias is seen in the estimated DMs for simulated data with pulse-broadening with a larger variability for a data set with a variable frequency scaling index, $\alpha$, as compared to that assuming a Kolmogorov turbulence. Application of the proposed method removes this bias robustly for data with band averaged signal-to-noise ratio larger than 100. We report, for the first time, the measurements of the scatter-broadening time and $\alpha$ from analysis of PSR J1643$-$1224, observed with upgraded Giant Metrewave Radio Telescope as part of the Indian Pulsar Timing Array experiment. These scattering parameters were found to vary with epoch and $\alpha$ was different from that expected for Kolmogorov turbulence. Finally, we present the DM time-series after application of the new technique to PSR J1643$-$1224., Comment: 9 pages, 8 figures, Submitted to MNRAS

The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTAs that constitute the International Pulsar Timing Array (IPTA). We show that despite making different modeling choices, there is no significant difference in the GWB parameters that are measured by the different PTAs, agreeing within $1\sigma$. The pulsar noise parameters are also consistent between different PTAs for the majority of the pulsars included in these analyses. We bridge the differences in modeling choices by adopting a standardized noise model for all pulsars and PTAs, finding that under this model there is a reduction in the tension in the pulsar noise parameters. As part of this reanalysis, we "extended" each PTA's data set by adding extra pulsars that were not timed by that PTA. Under these extensions, we find better constraints on the GWB amplitude and a higher signal-to-noise ratio for the Hellings and Downs correlations. These extensions serve as a prelude to the benefits offered by a full combination of data across all pulsars in the IPTA, i.e., the IPTA's Data Release 3, which will involve not just adding in additional pulsars, but also including data from all three PTAs where any given pulsar is timed by more than as single PTA., Comment: 21 pages, 9 figures, submitted to ApJ

Annual variations of interstellar scintillation can be modelled to constrain parameters of the ionized interstellar medium. If a pulsar is in a binary system, then investigating the orbital parameters is possible through analysis of the orbital variation of scintillation. In observations carried out from 2011 January to 2020 August by the European Pulsar Timing Array radio telescopes, PSRs~J0613$-$0200 and J0636+5128 show strong annual variations in their scintillation velocity, while the former additionally exhibits an orbital fluctuation. Bayesian theory and Markov-chain-Monte-Carlo methods are used to interpret these periodic variations. We assume a thin and anisotropic scattering screen model, and discuss the mildly and extremely anisotropic scattering cases. PSR~J0613$-$0200 is best described by mildly anisotropic scattering, while PSR~J0636+5128 exhibits extremely anisotropic scattering. We measure the distance, velocity and degree of anisotropy of the scattering screen for our two pulsars, finding that scattering screen distances from Earth for PSRs~J0613$-$0200 and J0636+5128 are 316$^{+28}_{-20}$\,pc and 262$^{+96}_{-38}$\,pc, respectively. The positions of these scattering screens are coincident with the shell of the Local Bubble towards both pulsars. These associations add to the growing evidence of the Local Bubble shell as a dominant region of scattering along many sightlines., Comment: Accepted by SCIENCE CHINA Physics, Mechanics & Astronomy ( SCPMA)

Pulsar Timing Array experiments probe the presence of possible scalar or pseudoscalar ultralight dark matter particles through decade-long timing of an ensemble of galactic millisecond radio pulsars. With the second data release of the European Pulsar Timing Array, we focus on the most robust scenario, in which dark matter interacts only gravitationally with ordinary baryonic matter. Our results show that ultralight particles with masses $10^{-24.0}~\text{eV} \lesssim m \lesssim 10^{-23.3}~\text{eV}$ cannot constitute $100\%$ of the measured local dark matter density, but can have at most local density $\rho\lesssim 0.3$ GeV/cm$^3$., Comment: 5 pages + acknowledgements + refs, 2 figures. Text and figures match the version published in PRL

The European Pulsar Timing Array (EPTA) and Indian Pulsar Timing Array (InPTA) collaborations have measured a low-frequency common signal in the combination of their second and first data releases respectively, with the correlation properties of a gravitational wave background (GWB). Such signal may have its origin in a number of physical processes including a cosmic population of inspiralling supermassive black hole binaries (SMBHBs); inflation, phase transitions, cosmic strings and tensor mode generation by non-linear evolution of scalar perturbations in the early Universe; oscillations of the Galactic potential in the presence of ultra-light dark matter (ULDM). At the current stage of emerging evidence, it is impossible to discriminate among the different origins. Therefore, in this paper, we consider each process separately, and investigate the implications of the signal under the hypothesis that it is generated by that specific process. We find that the signal is consistent with a cosmic population of inspiralling SMBHBs, and its relatively high amplitude can be used to place constraints on binary merger timescales and the SMBH-host galaxy scaling relations. If this origin is confirmed, this is the first direct evidence that SMBHBs merge in nature, adding an important observational piece to the puzzle of structure formation and galaxy evolution. As for early Universe processes, the measurement would place tight constraints on the cosmic string tension and on the level of turbulence developed by first-order phase transitions. Other processes would require non-standard scenarios, such as a blue-tilted inflationary spectrum or an excess in the primordial spectrum of scalar perturbations at large wavenumbers. Finally, a ULDM origin of the detected signal is disfavoured, which leads to direct constraints on the abundance of ULDM in our Galaxy., Comment: 30 pages, 23 figures, replaced to match the version published in Astronomy & Astrophysics, note the change in the numbering order in the series (now paper IV)

We present the results of a search for continuous gravitational wave signals (CGWs) in the second data release (DR2) of the European Pulsar Timing Array (EPTA) collaboration. The most significant candidate event from this search has a gravitational wave frequency of 4-5 nHz. Such a signal could be generated by a supermassive black hole binary (SMBHB) in the local Universe. We present the results of a follow-up analysis of this candidate using both Bayesian and frequentist methods. The Bayesian analysis gives a Bayes factor of 4 in favor of the presence of the CGW over a common uncorrelated noise process, while the frequentist analysis estimates the p-value of the candidate to be 1%, also assuming the presence of common uncorrelated red noise. However, comparing a model that includes both a CGW and a gravitational wave background (GWB) to a GWB only, the Bayes factor in favour of the CGW model is only 0.7. Therefore, we cannot conclusively determine the origin of the observed feature, but we cannot rule it out as a CGW source. We present results of simulations that demonstrate that data containing a weak gravitational wave background can be misinterpreted as data including a CGW and vice versa, providing two plausible explanations of the EPTA DR2 data. Further investigations combining data from all PTA collaborations will be needed to reveal the true origin of this feature., Comment: 13 figures, 15 pages, accepted

The nanohertz gravitational wave background (GWB) is expected to be an aggregate signal of an ensemble of gravitational waves emitted predominantly by a large population of coalescing supermassive black hole binaries in the centres of merging galaxies. Pulsar timing arrays, ensembles of extremely stable pulsars, are the most precise experiments capable of detecting this background. However, the subtle imprints that the GWB induces on pulsar timing data are obscured by many sources of noise. These must be carefully characterized to increase the sensitivity to the GWB. In this paper, we present a novel technique to estimate the optimal number of frequency coefficients for modelling achromatic and chromatic noise and perform model selection. We also incorporate a new model to fit for scattering variations in the pulsar timing package temponest and created realistic simulations of the European Pulsar Timing Array (EPTA) datasets that allowed us to test the efficacy of our noise modelling algorithms. We present an in-depth analysis of the noise properties of 25 millisecond pulsars (MSPs) that form the second data release (DR2) of the EPTA and investigate the effect of incorporating low-frequency data from the Indian PTA collaboration. We use enterprise and temponest packages to compare noise models with those reported with the EPTA DR1. We find that, while in some pulsars we can successfully disentangle chromatic from achromatic noise owing to the wider frequency coverage in DR2, in others the noise models evolve in a more complicated way. We also find evidence of long-term scattering variations in PSR J1600$-$3053. Through our simulations, we identify intrinsic biases in our current noise analysis techniques and discuss their effect on GWB searches. The results presented here directly help improve sensitivity to the GWB and are already being used as part of global PTA efforts., Comment: 20 pages, 6 figures, 9 tables

Pulsar timing arrays offer a probe of the low-frequency gravitational wave spectrum (1 - 100 nanohertz), which is intimately connected to a number of markers that can uniquely trace the formation and evolution of the Universe. We present the dataset and the results of the timing analysis from the second data release of the European Pulsar Timing Array (EPTA). The dataset contains high-precision pulsar timing data from 25 millisecond pulsars collected with the five largest radio telescopes in Europe, as well as the Large European Array for Pulsars. The dataset forms the foundation for the search for gravitational waves by the EPTA, presented in associated papers. We describe the dataset and present the results of the frequentist and Bayesian pulsar timing analysis for individual millisecond pulsars that have been observed over the last ~25 years. We discuss the improvements to the individual pulsar parameter estimates, as well as new measurements of the physical properties of these pulsars and their companions. This data release extends the dataset from EPTA Data Release 1 up to the beginning of 2021, with individual pulsar datasets with timespans ranging from 14 to 25 years. These lead to improved constraints on annual parallaxes, secular variation of the orbital period, and Shapiro delay for a number of sources. Based on these results, we derived astrophysical parameters that include distances, transverse velocities, binary pulsar masses, and annual orbital parallaxes., Comment: 29 pages, 9 figures, 13 tables, Astronomy & Astrophysics in press

We present the results of the search for an isotropic stochastic gravitational wave background (GWB) at nanohertz frequencies using the second data release of the European Pulsar Timing Array (EPTA) for 25 millisecond pulsars and a combination with the first data release of the Indian Pulsar Timing Array (InPTA). We analysed (i) the full 24.7-year EPTA data set, (ii) its 10.3-year subset based on modern observing systems, (iii) the combination of the full data set with the first data release of the InPTA for ten commonly timed millisecond pulsars, and (iv) the combination of the 10.3-year subset with the InPTA data. These combinations allowed us to probe the contributions of instrumental noise and interstellar propagation effects. With the full data set, we find marginal evidence for a GWB, with a Bayes factor of four and a false alarm probability of $4\%$. With the 10.3-year subset, we report evidence for a GWB, with a Bayes factor of $60$ and a false alarm probability of about $0.1\%$ ($\gtrsim 3\sigma$ significance). The addition of the InPTA data yields results that are broadly consistent with the EPTA-only data sets, with the benefit of better noise modelling. Analyses were performed with different data processing pipelines to test the consistency of the results from independent software packages. The inferred spectrum from the latest EPTA data from new generation observing systems is rather uncertain and in mild tension with the common signal measured in the full data set. However, if the spectral index is fixed at 13/3, the two data sets give a similar amplitude of ($2.5\pm0.7)\times10^{-15}$ at a reference frequency of $1\,{\rm yr}^{-1}$. By continuing our detection efforts as part of the International Pulsar Timing Array (IPTA), we expect to be able to improve the measurement of spatial correlations and better characterise this signal in the coming years., Comment: 21 pages, 14 figures, 4 appendix figures, accepted for publication in A&A

We search for a stochastic gravitational wave background (SGWB) generated by a network of cosmic strings using six millisecond pulsars from Data Release 2 (DR2) of the European Pulsar Timing Array (EPTA). We perform a Bayesian analysis considering two models for the network of cosmic string loops, and compare it to a simple power-law model which is expected from the population of supermassive black hole binaries. Our main strong assumption is that the previously reported common red noise process is a SGWB. We find that the one-parameter cosmic string model is slightly favored over a power-law model thanks to its simplicity. If we assume a two-component stochastic signal in the data (supermassive black hole binary population and the signal from cosmic strings), we get a $95\%$ upper limit on the string tension of $\log_{10}(G\mu) < -9.9$ ($-10.5$) for the two cosmic string models we consider. In extended two-parameter string models, we were unable to constrain the number of kinks. We test two approximate and fast Bayesian data analysis methods against the most rigorous analysis and find consistent results. These two fast and efficient methods are applicable to all SGWBs, independent of their source, and will be crucial for analysis of extended data sets., Comment: 14 pages, 6 figures; typo corrected in (5)

The wideband timing technique enables the high-precision simultaneous estimation of pulsar Times of Arrival (ToAs) and Dispersion Measures (DMs) while effectively modeling frequency-dependent profile evolution. We present two novel independent methods that extend the standard wideband technique to handle simultaneous multi-band pulsar data incorporating profile evolution over a larger frequency span to estimate DMs and ToAs with enhanced precision. We implement the wideband likelihood using the libstempo python interface to perform wideband timing in the tempo2 framework. We present the application of these techniques to the dataset of fourteen millisecond pulsars observed simultaneously in Band 3 (300 - 500 MHz) and Band 5 (1260 - 1460 MHz) of the upgraded Giant Metrewave Radio Telescope (uGMRT) with a large band gap of 760 MHz as a part of the Indian Pulsar Timing Array (InPTA) campaign. We achieve increased ToA and DM precision and sub-microsecond root mean square post-fit timing residuals by combining simultaneous multi-band pulsar observations done in non-contiguous bands for the first time using our novel techniques., Comment: Published in MNRAS

The Indian Pulsar Timing Array (InPTA) collaboration has recently made its first official data release (DR1) for a sample of 14 pulsars using 3.5 years of uGMRT observations. We present the results of single-pulsar noise analysis for each of these 14 pulsars using the InPTA DR1. For this purpose, we consider white noise, achromatic red noise, dispersion measure (DM) variations, and scattering variations in our analysis. We apply Bayesian model selection to obtain the preferred noise models among these for each pulsar. For PSR J1600$-$3053, we find no evidence of DM and scattering variations, while for PSR J1909$-$3744, we find no significant scattering variations. Properties vary dramatically among pulsars. For example, we find a strong chromatic noise with chromatic index $\sim$ 2.9 for PSR J1939+2134, indicating the possibility of a scattering index that doesn't agree with that expected for a Kolmogorov scattering medium consistent with similar results for millisecond pulsars in past studies. Despite the relatively short time baseline, the noise models broadly agree with the other PTAs and provide, at the same time, well-constrained DM and scattering variations., Comment: Accepted for publication in PRD, 30 pages, 17 figures, 4 tables

Decades long monitoring of millisecond pulsars, which exhibit highly stable rotational periods, in pulsar timing array experiments is on the threshold of discovering nanohertz stochastic gravitational wave background. This paper describes the Indian Pulsar timing array (InPTA) experiment, which employs the upgraded Giant Metrewave Radio Telescope (uGMRT) for timing an ensemble of millisecond pulsars for this purpose. We highlight InPTA's observation strategies and analysis methods, which are relevant for a future PTA experiment with the more sensitive Square Kilometer Array (SKA) telescope. We show that the unique multi-sub-array multi-band wide-bandwidth frequency coverage of the InPTA provides Dispersion Measure estimates with unprecedented precision for PTA pulsars, e.g., ~ 2 x 10{-5} pc-cm{-3} for PSR J1909-3744. Configuring the SKA-low and SKA-mid as two and four sub-arrays respectively, it is shown that comparable precision is achievable, using observation strategies similar to those pursued by the InPTA, for a larger sample of 62 pulsars requiring about 26 and 7 hours per epoch for the SKA-mid and the SKA-low telescopes respectively. We also review the ongoing efforts to develop PTA-relevant general relativistic constructs that will be required to search for nanohertz gravitational waves from isolated super-massive black hole binary systems like blazar OJ 287. These efforts should be relevant to pursue persistent multi-messenger gravitational wave astronomy during the forthcoming era of the SKA telescope, the Thirty Meter Telescope, and the next-generation Event Horizon Telescope., Comment: Accepted for publication in Journal of Astronomy and Astrophysics for Special Issue on Indian Participation in the SKA (Editors : Abhirup Datta, Nirupam Roy, Preeti Kharb and Tirthankar Roy Choudhury)

We present the pulse arrival times and high-precision dispersion measure estimates for 14 millisecond pulsars observed simultaneously in the 300-500 MHz and 1260-1460 MHz frequency bands using the upgraded Giant Metrewave Radio Telescope (uGMRT). The data spans over a baseline of 3.5 years (2018-2021), and is the first official data release made available by the Indian Pulsar Timing Array collaboration. This data release presents a unique opportunity for investigating the interstellar medium effects at low radio frequencies and their impact on the timing precision of pulsar timing array experiments. In addition to the dispersion measure time series and pulse arrival times obtained using both narrowband and wideband timing techniques, we also present the dispersion measure structure function analysis for selected pulsars. Our ongoing investigations regarding the frequency dependence of dispersion measures have been discussed. Based on the preliminary analysis for five millisecond pulsars, we do not find any conclusive evidence of chromaticity in dispersion measures. Data from regular simultaneous two-frequency observations are presented for the first time in this work. This distinctive feature leads us to the highest precision dispersion measure estimates obtained so far for a subset of our sample. Simultaneous multi-band uGMRT observations in Band 3 and Band 5 are crucial for high-precision dispersion measure estimation and for the prospect of expanding the overall frequency coverage upon the combination of data from the various Pulsar Timing Array consortia in the near future. Parts of the data presented in this work are expected to be incorporated into the upcoming third data release of the International Pulsar Timing Array., Comment: 23 pages, 21 figures, 3 tables. Published in PASA

High-precision measurements of the pulsar dispersion measure (DM) are possible using telescopes with low-frequency wideband receivers. We present an initial study of the application of the wideband timing technique, which can simultaneously measure the pulsar times of arrival (ToAs) and DMs, for a set of five pulsars observed with the upgraded Giant Metrewave Radio Telescope (uGMRT) as part of the Indian Pulsar Timing Array (InPTA) campaign. We have used the observations with the 300-500 MHz band of the uGMRT for this purpose. We obtain high precision in DM measurements with precisions of the order 10^{-6}cm^{-3}pc. The ToAs obtained have sub-{\mu}s precision and the root-mean-square of the post-fit ToA residuals are in the sub-{\mu}s range. We find that the uncertainties in the DMs and ToAs obtained with this wideband technique, applied to low-frequency data, are consistent with the results obtained with traditional pulsar timing techniques and comparable to high-frequency results from other PTAs. This work opens up an interesting possibility of using low-frequency wideband observations for precision pulsar timing and gravitational wave detection with similar precision as high-frequency observations used conventionally., Comment: Accepted for publication in MNRAS

PSR J1713+0747 is one of the most precisely timed pulsars in the international pulsar timing array experiment. This pulsar showed an abrupt profile shape change between April 16, 2021 (MJD 59320) and April 17, 2021 (MJD 59321). In this paper, we report the results from multi-frequency observations of this pulsar carried out with the upgraded Giant Metrewave Radio Telescope (uGMRT) before and after the event. We demonstrate the profile change seen in Band 5 (1260 MHz - 1460 MHz) and Band 3 (300 MHz - 500 MHz). The timing analysis of this pulsar shows a disturbance accompanying this profile change followed by a recovery with a timescale of $\sim 159$ days. Our data suggest that a model with chromatic index as a free parameter is preferred over models with combinations of achromaticity with DM bump or scattering bump. We determine the frequency dependence to be $\sim\nu^{+1.34}$., Comment: Accepted for publication in MNRAS-Letters

Glitches are the observational manifestations of superfluidity inside neutron stars. The aim of this paper is to describe an automated glitch detection pipeline, which can alert the observers on possible real-time detection of rotational glitches in pulsars. Post alert, the pulsars can be monitored at a higher cadence to measure the post-glitch recovery phase. Two algorithms namely, Median Absolute Deviation (MAD) and polynomial regression have been explored to detect glitches in real time. The pipeline has been optimized with the help of simulated timing residuals for both the algorithms. Based on the simulations, we conclude that the polynomial regression algorithm is significantly more effective for real time glitch detection. The pipeline has been tested on a few published glitches. This pipeline is presently implemented at the Ooty Radio Telescope. In the era of upcoming large telescopes like SKA, several hundreds of pulsars will be observed regularly and such a tool will be useful for both real-time detection as well as optimal utilization of observation time for such glitching pulsars., Comment: 10 pages, 6 figures, Accepted for publication in MNRAS

Pulsar radio emission undergoes dispersion due to the presence of free electrons in the interstellar medium (ISM). The dispersive delay in the arrival time of pulsar signal changes over time due to the varying ISM electron column density along the line of sight. Correcting for this delay accurately is crucial for the detection of nanohertz gravitational waves using Pulsar Timing Arrays. In this work, we present in-band and inter-band DM estimates of four pulsars observed with uGMRT over the timescale of a year using two different template alignment methods. The DMs obtained using both these methods show only subtle differences for PSR 1713+0747 and J1909$-$3744. A considerable offset is seen in the DM of PSR J1939+2134 and J2145$-$0750 between the two methods. This could be due to the presence of scattering in the former and profile evolution in the latter. We find that both methods are useful but could have a systematic offset between the DMs obtained. Irrespective of the template alignment methods followed, the precision on the DMs obtained is about $10^{-3}$ pc cm$^{-3}$ using only BAND3 and $10^{-4}$ pc cm$^{-3}$ after combining data from BAND3 and BAND5 of the uGMRT. In a particular result, we have detected a DM excess of about $5\times10^{-3}$ pc cm$^{-3}$ on 24 February 2019 for PSR J2145$-$0750. This excess appears to be due to the interaction region created by fast solar wind from a coronal hole and a coronal mass ejection (CME) observed from the Sun on that epoch. A detailed analysis of this interesting event is presented., Comment: 11 pages, 6 figures, 2 tables. Accepted by A&A

We introduce pinta, a pipeline for reducing the upgraded Giant Metre-wave Radio Telescope (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. We also discuss a novel and independent determination of the relative time offsets between the different back-end modes of uGMRT and the interpretation of the uGMRT observation frequency settings, and their agreement with results obtained from engineering tests. Further, we demonstrate the capability of pinta to generate data products which can produce high-precision TOAs using PSR J1909-3744 as an example. These results are crucial for performing precision pulsar timing with the uGMRT.

The abrupt change in the pulse period of a pulsar is called a pulsar glitch. In this paper, we present eleven pulsar glitches detected using the Ooty Radio Telescope (ORT) and the upgraded Giant Metrewave Radio Telescope (uGMRT) in high cadence timing observations of 8 pulsars. The measured relative amplitude of glitches ($\Delta \nu/\nu$) from our data ranges from $10^{-6}$ to $10^{-9}$. Among these glitches, three are new discoveries, being reported for the first time. We also reanalyze the largest pulsar glitch in the Crab pulsar (PSR J0534+2200) by fitting the ORT data to a new phenomenological model including the slow rise in the post-glitch evolution. We measure an exponential recovery of 30 days after the Vela glitch detected on MJD 57734 with a healing factor $Q=5.8\times 10^{-3}$. Further, we report the largest glitch ($\Delta \nu/\nu = 3147.9 \times 10^{-9}$) so far in PSR J1731$-$4744., Comment: 11 pages, 12 figures, 3 tables; Accepted for publication in MNRAS

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., Comment: Accepted for publication in ApJ

We report on the results of multi-frequency follow-up observations of three pulsars (PSRs J0026+6320, J2208+5500 and J2217+5733) discovered with the Giant Metrewave Radio Telescope (GMRT). These observations were carried out with the GMRT and the Ooty Radio Telescope (ORT). We report improved timing solutions for all three pulsars. For PSR J2208+5500, we estimate the nulling fraction to be 53(3)\%. The steep spectrum of this pulsar, its single component profile, and narrow pulse width suggest its single component to be a core component. If so, this significant cessation of emission in a core component is inconsistent with a geometric origin of nulls, such as those due to `empty' sightline traverses, and more likely due to intrinsic changes in the pulsar magnetosphere. We have measured scatter-broadening timescales at 325 and 610 MHz for PSRs J0026+6320 and J2217+5733. The implied scatter-broadening frequency scaling index of $-$2.9 for both pulsars is different from that expected assuming Kolmogorov turbulence in the interstellar medium. We also report spectral indices, obtained from imaging observations, for all three pulsars for the first time. The spectra for two of these pulsars indicate a possible spectral turnover between 100$-$300 MHz. Multi-frequency timing analyses carried out for these pulsars have enabled us to determine dispersion measures (DMs) with accuracies of 0.01 pc cm\sups{--3}. This demonstrates the usefulness of quasi-simultaneous multi-frequency multi-epoch timing observations with the GMRT and the ORT for studying variations in DM for millisecond pulsars., Comment: Accepted for publication in ApJ

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., Comment: Accepted for publication in ApJ

The radio as well as the high energy emission mechanism in pulsars is yet not understood properly. A multi-wavelength study is likely to help in better understanding of such processes. The first Indian space-based observatory, ASTROSAT, has five instruments aboard, which cover the electromagnetic spectrum from infra-red (1300 $\AA$) to hard X-ray (380 KeV). Cadmium Zinc Telluride Imager (CZTI), one of the five instruments is a hard X-ray telescope functional over an energy range of 20-380 KeV. We aim to estimate the timing offset introduced in the data acquisition pipeline of the instrument, which will help in time alignment of high energy time series with those from two other ground-based observatories, viz. the Giant Meterwave Radio Telescope (GMRT) and the Ooty Radio Telescope (ORT). PSR B0531+21 is a well-studied pulsar with nearly aligned radio and hard X-ray pulse profiles. We use simultaneous observations of this pulsar with the ASTROSAT, the ORT and the GMRT. The pulsar was especially observed using the ORT with almost daily cadence to obtain good timing solutions. We also supplement the ORT data with archival FERMI data for estimation of timing noise. The timing offset of ASTROSAT instruments was estimated from fits to arrival time data at the ASTROSAT and the radio observatories. We estimate the offset between the GMRT and the ASTROSAT-CZTI to be -4716 $\pm$ 50 $\mu s$. The corresponding offset with the ORT was -29639 $\pm$ 50 $\mu s$. The offsets between the GMRT and Fermi-LAT -5368 $\pm$ 56 $\mu s$. (Abridged), Comment: 6 pages, 5 figures, 2 tables, Revised and Updated, accepted for publication in A&A

We report the results of a blind pulsar survey carried out with the Giant Metrewave Radio Telescope (GMRT) at 325 MHz. The survey covered about 10\% of the region between Galactic longitude 45$^{\circ} < l <$ 135$^{\circ}$ and Galactic latitude 1$^{\circ} < |b| <$ 10$^{\circ}$ with a dwell time of 1800 s, resulting in the detection of 28 pulsars. One of these, PSR J1838+1523, was previously unknown and has a period of 549 ms and a dispersion measure of 68 pc cm\sups{$-$3}. We also present the timing solution of this pulsar obtained from multi-frequency timing observations carried out with the GMRT and the Ooty Radio Telescope. The measured flux density of this pulsar is 4.3$\pm$1.8 and 1.2$\pm$0.7 mJy at 325 and 610 MHz, respectively. This implies a spectral index of $-$2$\pm$0.8, thus making the expected flux density at 1.4 GHz to be about 0.2 mJy, which would be just detectable in the high frequency pulsar surveys like the Northern High Time Resolution Universe pulsar survey. This discovery underlines the importance of low frequency pulsar surveys in detecting steep spectrum pulsars, thus providing complementary coverage of the pulsar population., Comment: Accepted for publication in MNRAS

Single pulse observations, characterizing in detail, the nulling behaviour of PSR B1706$-$16 are being reported for the first time in this paper. Our regular long duration monitoring of this pulsar reveals long nulls of 2 to 5 hours with an overall nulling fraction of 31$\pm$2\%. The pulsar shows two distinct phases of emission. It is usually in an active phase, characterized by pulsations interspersed with shorter nulls, with a nulling fraction of about 15 \%, but it also rarely switches to an inactive phase, consisting of long nulls. The nulls in this pulsar are concurrent between 326.5 and 610 MHz. Profile mode changes accompanied by changes in fluctuation properties are seen in this pulsar, which switches from mode A before a null to mode B after the null. The distribution of null durations in this pulsar is bimodal. With its occasional long nulls, PSR B1706$-$16 joins the small group of intermediate nullers, which lie between the classical nullers and the intermittent pulsars. Similar to other intermediate nullers, PSR B1706$-$16 shows high timing noise, which could be due to its rare long nulls if one assumes that the slowdown rate during such nulls is different from that during the bursts., Comment: Accepted for publication in MNRAS

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., Comment: 8 pages, 6 figures; Accepted for publication in MNRAS

We present multi-wavelength scatter broadening observations of 47 pulsars, made with the Giant Metre-wave Radio Telescope (GMRT), Ooty Radio Telescope (ORT) and Long Wavelength Array (LWA). The GMRT observations have been made in the phased array mode at 148, 234, and 610 MHz and the ORT observations at 327 MHz. The LWA data sets have been obtained from the LWA pulsar data archive. The broadening of each pulsar as a function of observing frequency provides the frequency scaling index, $\alpha$. The estimations of $\alpha$ have been obtained for 39 pulsars, which include entirely new estimates for 31 pulsars. This study increases the total sample of pulsars available with $\alpha$ estimates by $\sim$50\%. The overall distribution of $\alpha$ with the dispersion measure (DM) of pulsar shows interesting variations, which are consistent with the earlier studies. However, for a given value of DM a range of $\alpha$ values are observed, indicating the characteristic turbulence along each line of sight. For each pulsar, the estimated level of turbulence, $C^{2}_{n_e}$, has also been compared with $\alpha$ and DM. Additionally, we compare the distribution of $\alpha$ with the theoretically predicated model to infer the general characteristics of the ionized interstellar medium (ISM). Nearly 65\% of the pulsars show a flatter index (i.e., $\alpha < 4.4$) than that is expected from the Kolmogorov turbulence model. Moreover, the group of pulsars having flatter index is typically associated with an enhanced value of $C^{2}_{n_e}$ than those with steeper index., Comment: 13 pages, 4 figures, 3 tables. Accepted for publication in ApJ

The use of beta-solenoid proteins as functionalizable, nanoscale, self-assembling molecular building blocks may have many applications, including templating the growth of wires or higher-dimensional structures. By understanding their mechanical strengths, we can efficiently design the proteins for specific functions. We present a study of the mechanical properties of seven beta-solenoid proteins using GROMACS molecular dynamics software to produce force/torque-displacement data, implement umbrella sampling of bending/twisting trajectories, produce Potentials of Mean Force (PMFs), extract effective spring constants, and calculate rigidities for two bending and two twisting directions for each protein. We examine the differences between computing the strength values from force/torque-displacement data alone and PMF data, and show how higher precision estimates can be obtained from the former. In addition to the analysis of the methods, we report estimates for the bend/twist persistence lengths for each protein, which range from 0.5-3.4 $\mu$m. We note that beta-solenoid proteins with internal disulfide bridges do not enjoy enhanced bending or twisting strength, and that the strongest correlate with bend/twist rigidity is the number of hydrogen bonds per turn. In addition, we compute estimates of the Young's modulus ($Y$) for each protein, which range from $Y$ = 3.5 to 7.2 GPa.

We performed simultaneous observations at 326.5 MHz with the Ooty Radio Telescope and at 326, 610 and 1308 MHz with the Giant Meterwave Radio Telescope for a sample of 12 pulsars, where frequency dependent single pulse behaviour was reported. The single pulse sequences were analysed with fluctuation analysis, sensitive to both the average fluctuation properties (using longitude resolved fluctuation spectrum and two-dimensional fluctuation spectrum) as well as temporal changes in these (using sliding two-dimensional fluctuation spectrum ) to establish concurrent changes in subpulse drifting over the multiple frequencies employed. We report subpulse drifting in PSR J0934$-$5249 for the first time. We also report pulse nulling measurements in PSRs J0934$-$5249, B1508+55, J1822$-$2256, B1845$-$19 and J1901$-$0906 for the first time. Our measurements of subpulse drifting and pulse nulling for the rest of the pulsars are consistent with previously reported values. Contrary to previous belief, we find no evidence for a frequency dependent drift pattern in PSR B2016+28 implied by non-simultaneous observations by Oster et al. (1977). In PSRs B1237+25, J1822$-$2256, J1901$-$0906 and B2045$-$16, our longer and more sensitive observations reveal multiple drift rates with distinct P3. We increase the sample of pulsars showing concurrent nulling across multiple frequencies by more than 100 percent, adding 4 more pulsars to this sample. Our results confirm and further strengthen the understanding that the subpulse drifting and pulse nulling are broadband consistent with previous studies (Gajjar et al. 2014a; Rankin 1986; Weltevrede et al. 2007) and are closely tied to physics of polar gap., Comment: 22 pages, 44 figures, Single pulse studies of pulsars, accepted by A&A

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., Comment: 20 pages, 4 figures. Accepted for publication in ApJ

This paper describes a new real-time versatile backend, the Pulsar Ooty Radio Telescope New Digital Efficient Receiver (PONDER), which has been designed to operate along with the legacy analog system of the Ooty Radio Telescope (ORT). PONDER makes use of the current state of the art computing hardware, a Graphical Processing Unit (GPU) and sufficiently large disk storage to support high time resolution real-time data of pulsar observations, obtained by coherent dedispersion over a bandpass of 16 MHz. Four different modes for pulsar observations are implemented in PONDER to provide standard reduced data products, such as time-stamped integrated profiles and dedispersed time series, allowing faster avenues to scientific results for a variety of pulsar studies. Additionally, PONDER also supports general modes of interplanetary scintillation (IPS) measurements and very long baseline interferometry data recording. The IPS mode yields a single polarisation correlated time series of solar wind scintillation over a bandwidth of about four times larger (16 MHz) than that of the legacy system as well as its fluctuation spectrum with high temporal and frequency resolutions. The key point is that all the above modes operate in real time. This paper presents the design aspects of PONDER and outlines the design methodology for future similar backends. It also explains the principal operations of PONDER, illustrates its capabilities for a variety of pulsar and IPS observations and demonstrates its usefulness for a variety of astrophysical studies using the high sensitivity of the ORT., Comment: 25 pages, 14 figures, Accepted by Experimental Astronomy

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