-
The ORT and the uGMRT Pulsar Monitoring Program : Pulsar Timing Irregularities & the Gaussian Process Realization
Authors:
Himanshu Grover,
Bhal Chandra Joshi,
Jaikhomba Singha,
Erbil Gügercinoğlu,
Paramasivan Arumugam,
Debades Bandyopadhyay,
James O. Chibueze,
Shantanu Desai,
Innocent O. Eya,
Anu Kundu,
Johnson O. Urama
Abstract:
The spin-down law of pulsars is generally perturbed by two types of timing irregularities: glitches and timing noise. Glitches are sudden changes in the rotational frequency of pulsars, while timing noise is a discernible stochastic wandering in the phase, period, or spin-down rate of a pulsar. We present the timing results of a sample of glitching pulsars observed using the Ooty Radio Telescope (…
▽ More
The spin-down law of pulsars is generally perturbed by two types of timing irregularities: glitches and timing noise. Glitches are sudden changes in the rotational frequency of pulsars, while timing noise is a discernible stochastic wandering in the phase, period, or spin-down rate of a pulsar. We present the timing results of a sample of glitching pulsars observed using the Ooty Radio Telescope (ORT) and the upgraded Giant Metrewave Radio Telescope (uGMRT). Our findings include timing noise analysis for 17 pulsars, with seven being reported for the first time. We detected five glitches in four pulsars and a glitch-like event in J1825-0935. The frequency evolution of glitch in pulsars, J0742-2822 and J1740-3015, is presented for the first time. Additionally, we report timing noise results for three glitching pulsars. The timing noise was analyzed separately in the pre-glitch region and post-glitch regions. We observed an increase in the red noise parameters in the post-glitch regions, where exponential recovery was considered in the noise analysis. Timing noise can introduce ambiguities in the correct evaluation of glitch observations. Hence, it is important to consider timing noise in glitch analysis. We propose an innovative glitch verification approach designed to discern between a glitch and strong timing noise. The novel glitch analysis technique is also demonstrated using the observed data.
△ Less
Submitted 23 May, 2024;
originally announced May 2024.
-
Low-frequency pulse-jitter measurement with the uGMRT I : PSR J0437$-$4715
Authors:
Tomonosuke Kikunaga,
Shinnosuke Hisano,
Neelam Dhanda Batra,
Shantanu Desai,
Bhal Chandra Joshi,
Manjari Bagchi,
T. Prabu,
Keitaro Takahashi,
Swetha Arumugam,
Adarsh Bathula,
Subhajit Dandapat,
Debabrata Deb,
Churchil Dwivedi,
Yashwant Gupta,
Shebin Jose Jacob,
Fazal Kareem,
Nobleson K,
Pragna Mamidipaka,
Avinash Kumar Paladi,
Arul Pandian B,
Prerna Rana,
Jaikhomba Singha,
Aman Srivastava,
Mayuresh Surnis,
Pratik Tarafdar
Abstract:
High-precision pulsar timing observations are limited in their accuracy by the jitter noise that appears in the arrival time of pulses. Therefore, it is important to systematically characterise the amplitude of the jitter noise and its variation with frequency. In this paper, we provide jitter measurements from low-frequency wideband observations of PSR J0437$-$4715 using data obtained as part of…
▽ More
High-precision pulsar timing observations are limited in their accuracy by the jitter noise that appears in the arrival time of pulses. Therefore, it is important to systematically characterise the amplitude of the jitter noise and its variation with frequency. In this paper, we provide jitter measurements from low-frequency wideband observations of PSR J0437$-$4715 using data obtained as part of the Indian Pulsar Timing Array experiment. We were able to detect jitter in both the 300 - 500 MHz and 1260 - 1460 MHz observations of the upgraded Giant Metrewave Radio Telescope (uGMRT). The former is the first jitter measurement for this pulsar below 700 MHz, and the latter is in good agreement with results from previous studies. In addition, at 300 - 500 MHz, we investigated the frequency dependence of the jitter by calculating the jitter for each sub-banded arrival time of pulses. We found that the jitter amplitude increases with frequency. This trend is opposite as compared to previous studies, indicating that there is a turnover at intermediate frequencies. It will be possible to investigate this in more detail with uGMRT observations at 550 - 750 MHz and future high sensitive wideband observations from next generation telescopes, such as the Square Kilometre Array. We also explored the effect of jitter on the high precision dispersion measure (DM) measurements derived from short duration observations. We find that even though the DM precision will be better at lower frequencies due to the smaller amplitude of jitter noise, it will limit the DM precision for high signal-to-noise observations, which are of short durations. This limitation can be overcome by integrating for a long enough duration optimised for a given pulsar.
△ Less
Submitted 18 April, 2024; v1 submitted 4 December, 2023;
originally announced December 2023.
-
Using low-frequency scatter-broadening measurements for precision estimates of dispersion measures
Authors:
Jaikhomba Singha,
Bhal Chandra Joshi,
M. A. Krishnakumar,
Fazal Kareem,
Adarsh Bathula,
Churchil Dwivedi,
Shebin Jose Jacob,
Shantanu Desai,
Pratik Tarafdar,
P. Arumugam,
Swetha Arumugam,
Manjari Bagchi,
Neelam Dhanda Batra,
Subhajit Dandapat,
Debabrata Deb,
Jyotijwal Debnath,
A Gopakumar,
Yashwant Gupta,
Shinnosuke Hisano,
Ryo Kato,
Tomonosuke Kikunaga,
Piyush Marmat,
K. Nobleson,
Avinash K. Paladi,
Arul Pandian B.
, et al. (6 additional authors not shown)
Abstract:
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…
▽ More
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, $α$, 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 $α$ 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 $α$ 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.
△ Less
Submitted 28 September, 2023;
originally announced September 2023.
-
Comparing recent PTA results on the nanohertz stochastic gravitational wave background
Authors:
The International Pulsar Timing Array Collaboration,
G. Agazie,
J. Antoniadis,
A. Anumarlapudi,
A. M. Archibald,
P. Arumugam,
S. Arumugam,
Z. Arzoumanian,
J. Askew,
S. Babak,
M. Bagchi,
M. Bailes,
A. -S. Bak Nielsen,
P. T. Baker,
C. G. Bassa,
A. Bathula,
B. Bécsy,
A. Berthereau,
N. D. R. Bhat,
L. Blecha,
M. Bonetti,
E. Bortolas,
A. Brazier,
P. R. Brook,
M. Burgay
, et al. (220 additional authors not shown)
Abstract:
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 PTA…
▽ More
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σ$. 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.
△ Less
Submitted 1 September, 2023;
originally announced September 2023.
-
The second data release from the European Pulsar Timing Array: IV. Implications for massive black holes, dark matter and the early Universe
Authors:
J. Antoniadis,
P. Arumugam,
S. Arumugam,
P. Auclair,
S. Babak,
M. Bagchi,
A. -S. Bak Nielsen,
E. Barausse,
C. G. Bassa,
A. Bathula,
A. Berthereau,
M. Bonetti,
E. Bortolas,
P. R. Brook,
M. Burgay,
R. N. Caballero,
C. Caprini,
A. Chalumeau,
D. J. Champion,
S. Chanlaridis,
S. Chen,
I. Cognard,
M. Crisostomi,
S. Dandapat,
D. Deb
, et al. (89 additional authors not shown)
Abstract:
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 sup…
▽ More
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.
△ Less
Submitted 15 May, 2024; v1 submitted 28 June, 2023;
originally announced June 2023.
-
The second data release from the European Pulsar Timing Array V. Search for continuous gravitational wave signals
Authors:
J. Antoniadis,
P. Arumugam,
S. Arumugam,
S. Babak,
M. Bagchi,
A. S. Bak Nielsen,
C. G. Bassa,
A. Bathula,
A. Berthereau,
M. Bonetti,
E. Bortolas,
P. R. Brook,
M. Burgay,
R. N. Caballero,
A. Chalumeau,
D. J. Champion,
S. Chanlaridis,
S. Chen,
I. Cognard,
S. Dandapat,
D. Deb,
S. Desai,
G. Desvignes,
N. Dhanda-Batra,
C. Dwivedi
, et al. (75 additional authors not shown)
Abstract:
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 o…
▽ More
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.
△ Less
Submitted 25 June, 2024; v1 submitted 28 June, 2023;
originally announced June 2023.
-
The second data release from the European Pulsar Timing Array II. Customised pulsar noise models for spatially correlated gravitational waves
Authors:
J. Antoniadis,
P. Arumugam,
S. Arumugam,
S. Babak,
M. Bagchi,
A. S. Bak Nielsen,
C. G. Bassa,
A. Bathula,
A. Berthereau,
M. Bonetti,
E. Bortolas,
P. R. Brook,
M. Burgay,
R. N. Caballero,
A. Chalumeau,
D. J. Champion,
S. Chanlaridis,
S. Chen,
I. Cognard,
S. Dandapat,
D. Deb,
S. Desai,
G. Desvignes,
N. Dhanda-Batra,
C. Dwivedi
, et al. (73 additional authors not shown)
Abstract:
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 su…
▽ More
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.
△ Less
Submitted 28 June, 2023;
originally announced June 2023.
-
The second data release from the European Pulsar Timing Array III. Search for gravitational wave signals
Authors:
J. Antoniadis,
P. Arumugam,
S. Arumugam,
S. Babak,
M. Bagchi,
A. -S. Bak Nielsen,
C. G. Bassa,
A. Bathula,
A. Berthereau,
M. Bonetti,
E. Bortolas,
P. R. Brook,
M. Burgay,
R. N. Caballero,
A. Chalumeau,
D. J. Champion,
S. Chanlaridis,
S. Chen,
I. Cognard,
S. Dandapat,
D. Deb,
S. Desai,
G. Desvignes,
N. Dhanda-Batra,
C. Dwivedi
, et al. (73 additional authors not shown)
Abstract:
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…
▽ More
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σ$ 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.
△ Less
Submitted 28 June, 2023;
originally announced June 2023.
-
Multi-band Extension of the Wideband Timing Technique
Authors:
Avinash Kumar Paladi,
Churchil Dwivedi,
Prerna Rana,
Nobleson K,
Abhimanyu Susobhanan,
Bhal Chandra Joshi,
Pratik Tarafdar,
Debabrata Deb,
Swetha Arumugam,
A Gopakumar,
M A Krishnakumar,
Neelam Dhanda Batra,
Jyotijwal Debnath,
Fazal Kareem,
Paramasivan Arumugam,
Manjari Bagchi,
Adarsh Bathula,
Subhajit Dandapat,
Shantanu Desai,
Yashwant Gupta,
Shinnosuke Hisano,
Divyansh Kharbanda,
Tomonosuke Kikunaga,
Neel Kolhe,
Yogesh Maan
, et al. (5 additional authors not shown)
Abstract:
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 sp…
▽ More
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.
△ Less
Submitted 8 November, 2023; v1 submitted 25 April, 2023;
originally announced April 2023.
-
Noise analysis of the Indian Pulsar Timing Array data release I
Authors:
Aman Srivastava,
Shantanu Desai,
Neel Kolhe,
Mayuresh Surnis,
Bhal Chandra Joshi,
Abhimanyu Susobhanan,
Aurélien Chalumeau,
Shinnosuke Hisano,
Nobleson K.,
Swetha Arumugam,
Divyansh Kharbanda,
Jaikhomba Singha,
Pratik Tarafdar,
P Arumugam,
Manjari Bagchi,
Adarsh Bathula,
Subhajit Dandapat,
Lankeswar Dey,
Churchil Dwivedi,
Raghav Girgaonkar,
A. Gopakumar,
Yashwant Gupta,
Tomonosuke Kikunaga,
M. A. Krishnakumar,
Kuo Liu
, et al. (6 additional authors not shown)
Abstract:
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 vari…
▽ More
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.
△ Less
Submitted 16 June, 2023; v1 submitted 21 March, 2023;
originally announced March 2023.
-
Nanohertz Gravitational Wave Astronomy during the SKA Era: An InPTA perspective
Authors:
Bhal Chandra Joshi,
Achamveedu Gopakumar,
Arul Pandian,
Thiagaraj Prabu,
Lankeswar Dey,
Manjari Bagchi,
Shantanu Desai,
Pratik Tarafdar,
Prerna Rana,
Yogesh Maan,
Neelam Dhanda Batra,
Raghav Girgaonkar,
Nikita Agarwal,
Paramasivan Arumugam,
Sarmistha Banik,
Avishek Basu,
Adarsh Bathula,
Subhajit Dandapat,
Yashwant Gupta,
Shinnosuke Hisano,
Ryo Kato,
Divyansh Kharbanda,
Tomonosuke Kikunaga,
Neel Kolhe,
M. A. Krishnakumar
, et al. (12 additional authors not shown)
Abstract:
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 pulsa…
▽ More
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.
△ Less
Submitted 13 July, 2022;
originally announced July 2022.
-
Pulsar timing irregularities and the Neutron Star interior in the era of SKA: An Indian Outlook
Authors:
Jaikhomba Singha,
Bhal Chandra Joshi,
Debades Bandyopadhyay,
Himanshu Grover,
Shantanu Desai,
P. Arumugam,
Sarmistha Banik
Abstract:
There are two types of timing irregularities seen in pulsars: glitches and timing noise. Both of these phenomena can help us probe the interior of such exotic objects. This article presents a brief overview of the observational and theoretical aspects of pulsar timing irregularities and the main results from the investigations of these phenomena in India. The relevance of such Indian programs for…
▽ More
There are two types of timing irregularities seen in pulsars: glitches and timing noise. Both of these phenomena can help us probe the interior of such exotic objects. This article presents a brief overview of the observational and theoretical aspects of pulsar timing irregularities and the main results from the investigations of these phenomena in India. The relevance of such Indian programs for monitoring of young pulsars with the Square Kilometer Array (SKA) is presented, highlighting possible contributions of the Indian neutron star community to the upcoming SKA endeavour.
△ Less
Submitted 13 July, 2022;
originally announced July 2022.
-
The Indian Pulsar Timing Array: First data release
Authors:
Pratik Tarafdar,
Nobleson K.,
Prerna Rana,
Jaikhomba Singha,
M. A. Krishnakumar,
Bhal Chandra Joshi,
Avinash Kumar Paladi,
Neel Kolhe,
Neelam Dhanda Batra,
Nikita Agarwal,
Adarsh Bathula,
Subhajit Dandapat,
Shantanu Desai,
Lankeswar Dey,
Shinnosuke Hisano,
Prathamesh Ingale,
Ryo Kato,
Divyansh Kharbanda,
Tomonosuke Kikunaga,
Piyush Marmat,
B. Arul Pandian,
T. Prabu,
Aman Srivastava,
Mayuresh Surnis,
Sai Chaitanya Susarla
, et al. (13 additional authors not shown)
Abstract:
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 collabor…
▽ More
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.
△ Less
Submitted 25 October, 2022; v1 submitted 18 June, 2022;
originally announced June 2022.
-
Constraining the parameterized neutron star equation of state with astronomical observations
Authors:
Jaikhomba Singha,
S. Mullai Vaneshwar,
Ankit Kumar
Abstract:
We utilise the phenomenologically parameterized piecewise polytropic equations of state to study various neutron star properties. We investigate the compliance of these equations of state with several astronomical observations. We also demonstrate that the theoretical estimates of the fractional moment of inertia cannot explain all the pulsar glitches observed. We model the crust as a solid sphero…
▽ More
We utilise the phenomenologically parameterized piecewise polytropic equations of state to study various neutron star properties. We investigate the compliance of these equations of state with several astronomical observations. We also demonstrate that the theoretical estimates of the fractional moment of inertia cannot explain all the pulsar glitches observed. We model the crust as a solid spheroidal shell to calculate the fractional moment of inertia of fast-spinning neutron stars. We also show that the braking index obtained in a simple magnetic dipole radiation model with a varying moment of inertia deviates significantly from the observed data. Future developments in both theory and observations may allow us to use the fractional moment of inertia and braking index as observational constraints for neutron star equation of state.
△ Less
Submitted 21 February, 2022;
originally announced March 2022.
-
Low-frequency wideband timing of InPTA pulsars observed with the uGMRT
Authors:
K Nobleson,
Nikita Agarwal,
Raghav Girgaonkar,
Arul Pandian,
Bhal Chandra Joshi,
M A Krishnakumar,
Abhimanyu Susobhanan,
Shantanu Desai,
T Prabu,
Adarsh Bathula,
Timothy T Pennucci,
Sarmistha Banik,
Manjari Bagchi,
Neelam Dhanda Batra,
Arpita Choudhary,
Subhajit Dandapat,
Lankeswar Dey,
Yashwant Gupta,
Shinnosuke Hisano,
Ryo Kato,
Divyansh Kharbanda,
Tomonosuke Kikunaga,
Neel Kolhe,
Yogesh Maan,
Piyush Marmat
, et al. (7 additional authors not shown)
Abstract:
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…
▽ More
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-μs precision and the root-mean-square of the post-fit ToA residuals are in the sub-μ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.
△ Less
Submitted 23 February, 2022; v1 submitted 13 December, 2021;
originally announced December 2021.
-
Evidence for profile changes in PSR J1713+0747 using the uGMRT
Authors:
Jaikhomba Singha,
Mayuresh P Surnis,
Bhal Chandra Joshi,
Pratik Tarafdar,
Prerna Rana,
Abhimanyu Susobhanan,
Raghav Girgaonkar,
Neel Kolhe,
Nikita Agarwal,
Shantanu Desai,
T Prabu,
Adarsh Bathula,
Subhajit Dandapat,
Lankeswar Dey,
Shinnosuke Hisano,
Ryo Kato,
Divyansh Kharbanda,
Tomonosuke Kikunaga,
Piyush Marmat,
Sai Chaitanya Susarla,
Manjari Bagchi,
Neelam Dhanda Batra,
Arpita Choudhury,
A Gopakumar,
Yashwant Gupta
, et al. (7 additional authors not shown)
Abstract:
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…
▽ More
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ν^{+1.34}$.
△ Less
Submitted 16 August, 2021; v1 submitted 9 July, 2021;
originally announced July 2021.
-
A real-time Automated Glitch Detection Pipeline at Ooty Radio Telescope
Authors:
Jaikhomba Singha,
Avishek Basu,
M. A. Krishnakumar,
Bhal Chandra Joshi,
P. Arumugam
Abstract:
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 Absolut…
▽ More
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.
△ Less
Submitted 8 July, 2021; v1 submitted 3 June, 2021;
originally announced June 2021.
-
Phase-locking in $k$-partite networks of delay-coupled oscillators
Authors:
Joydeep Singha,
Ramakrishna Ramaswamy
Abstract:
We examine the dynamics of an ensemble of phase oscillators that are divided in $k$ sets, with time-delayed coupling interactions {\em only} between oscillators in different sets or partitions. The network of interactions thus form a $k-$partite graph. We observe a variety of phase-locked states that include, in addition to the in-phase fully synchronized solution, a variety of splay cluster solut…
▽ More
We examine the dynamics of an ensemble of phase oscillators that are divided in $k$ sets, with time-delayed coupling interactions {\em only} between oscillators in different sets or partitions. The network of interactions thus form a $k-$partite graph. We observe a variety of phase-locked states that include, in addition to the in-phase fully synchronized solution, a variety of splay cluster solutions; all oscillators within a partition are synchronised and the phase differences between oscillators in different partitions are multiples of $2π/k$. Such solutions exist independent of the delay and we determine the generalised stability criteria for the existence of these phase-locked solutions for the $k-$partite system. With increase in time-delay, there is an increase in multistability: the above generic solutions coexist with a number of other partially synchronized solutions. We apply the Ott-Antonsen ansatz for the special case of a symmetric $k-$partite graph to obtain a single time-delayed differential equation for the attracting synchronization manifold. The agreement with numerical results for the specific case of oscillators on a tripartite lattice (the $k=3$ case) is excellent.
△ Less
Submitted 25 October, 2021; v1 submitted 19 May, 2021;
originally announced May 2021.
-
High Precision Measurements of Interstellar Dispersion Measure with the upgraded GMRT
Authors:
M. A. Krishnakumar,
P. K. Manoharan,
Bhal Chandra Joshi,
Raghav Girgaonkar,
Shantanu Desai,
Manjari Bagchi,
K. Nobleson,
Lankeswar Dey,
Abhimanyu Susobhanan,
Sai Chaitanya Susarla,
Mayuresh P. Surnis,
Yogesh Maan,
A. Gopakumar,
Avishek Basu,
Neelam Dhanda Batra,
Arpita Choudhary,
Kishalay De,
Yashwant Gupta,
Arun Kumar Naidu,
Dhruv Pathak,
Jaikhomba Singha,
T. Prabu
Abstract:
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 wo…
▽ More
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.
△ Less
Submitted 1 July, 2021; v1 submitted 13 January, 2021;
originally announced January 2021.
-
pinta: The uGMRT Data Processing Pipeline for the Indian Pulsar Timing Array
Authors:
Abhimanyu Susobhanan,
Yogesh Maan,
Bhal Chandra Joshi,
T. Prabu,
Shantanu Desai,
K. Nobleson,
Sai Chaitanya Susarla,
Raghav Girgaonkar,
Lankeswar Dey,
Neelam Dhanda Batra,
Yashwant Gupta,
A. Gopakumar,
Manjari Bagchi,
Avishek Basu,
Suryarao Bethapudi,
Arpita Choudhary,
Kishalay De,
M. A. Krishnakumar,
P. K. Manoharan,
Arun Kumar Naidu,
Dhruv Pathak,
Jaikhomba Singha,
Mayuresh P. Surnis
Abstract:
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…
▽ More
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.
△ Less
Submitted 16 April, 2021; v1 submitted 6 July, 2020;
originally announced July 2020.
-
Chimera states for a globally coupled sine circle map lattice: spatiotemporal intermittency and hyperchaos
Authors:
Joydeep Singha,
Neelima Gupte
Abstract:
We study the existence of chimera states, i.e. mixed states, in a globally coupled sine circle map lattice, with different strengths of inter-group and intra-group coupling. We find that at specific values of the parameters of the CML, a completely random initial condition evolves to chimera states, having a phase synchronised and a phase desynchronised group, where the space time variation of the…
▽ More
We study the existence of chimera states, i.e. mixed states, in a globally coupled sine circle map lattice, with different strengths of inter-group and intra-group coupling. We find that at specific values of the parameters of the CML, a completely random initial condition evolves to chimera states, having a phase synchronised and a phase desynchronised group, where the space time variation of the phases of the maps in the desynchronised group shows structures similar to spatiotemporally intermittent regions. Using the complex order parameter we obtain a phase diagram that identifies the region in the parameter space which supports chimera states of this type, as well as other types of phase configurations such as globally phase synchronised states, two phase clustered states and fully phase desynchronised states. We estimate the volume of the basin of attraction of each kind of solution. The STI chimera region is studied in further detail via numerical and analytic stability analysis, and the Lyapunov spectrum is calculated. This state is identified to be hyperchaotic as the two largest Lyapunov exponents are found to be positive. The distributions of laminar and burst lengths in the incoherent region of the chimera show exponential behaviour. The average fraction of laminar/burst sites is identified to be the important quantity which governs the dynamics of the chimera. After an initial transient, these settle to steady values which can be used to reproduce the phase diagram in the chimera regime.
△ Less
Submitted 6 July, 2019;
originally announced July 2019.
-
Chimera states in coupled map lattices: Spatiotemporally intermittent behavior and an equivalent cellular automaton
Authors:
Joydeep Singha,
Neelima Gupte
Abstract:
We study the existence of chimera states, i.e. mixed states, in a globally coupled sine circle map lattice, with different strengths of inter-group and intra-group coupling. We find that at specific values of the parameters of the CML, a completely random initial condition evolves to chimera states, having a phase synchronised and a phase desynchronised group, where the space time variation of the…
▽ More
We study the existence of chimera states, i.e. mixed states, in a globally coupled sine circle map lattice, with different strengths of inter-group and intra-group coupling. We find that at specific values of the parameters of the CML, a completely random initial condition evolves to chimera states, having a phase synchronised and a phase desynchronised group, where the space time variation of the phases of the maps in the desynchronised group shows structures similar to spatiotemporally intermittent regions. Using the complex order parameter we obtain a phase diagram that identifies the region in the parameter space which supports chimera states of this type, as well as other types of phase configurations. An equivalent cellular automaton is obtained which shows space time behaviour similar to the CML. We also derive a mean field equation for this cellular automaton and compare its solutions with the corresponding phase configurations obtained in the parameter space, and show that it matches with the behavior seen for the CML.
△ Less
Submitted 29 November, 2018; v1 submitted 9 October, 2018;
originally announced October 2018.
-
Hand Gesture Recognition Based on Karhunen-Loeve Transform
Authors:
Joyeeta Singha,
Karen Das
Abstract:
In this paper, we have proposed a system based on K-L Transform to recognize different hand gestures. The system consists of five steps: skin filtering, palm cropping, edge detection, feature extraction, and classification. Firstly the hand is detected using skin filtering and palm cropping was performed to extract out only the palm portion of the hand. The extracted image was then processed using…
▽ More
In this paper, we have proposed a system based on K-L Transform to recognize different hand gestures. The system consists of five steps: skin filtering, palm cropping, edge detection, feature extraction, and classification. Firstly the hand is detected using skin filtering and palm cropping was performed to extract out only the palm portion of the hand. The extracted image was then processed using the Canny Edge Detection technique to extract the outline images of palm. After palm extraction, the features of hand were extracted using K-L Transform technique and finally the input gesture was recognized using proper classifier. In our system, we have tested for 10 different hand gestures, and recognizing rate obtained was 96%. Hence we propose an easy approach to recognize different hand gestures.
△ Less
Submitted 11 June, 2013;
originally announced June 2013.
-
Recognition of Indian Sign Language in Live Video
Authors:
Joyeeta Singha,
Karen Das
Abstract:
Sign Language Recognition has emerged as one of the important area of research in Computer Vision. The difficulty faced by the researchers is that the instances of signs vary with both motion and appearance. Thus, in this paper a novel approach for recognizing various alphabets of Indian Sign Language is proposed where continuous video sequences of the signs have been considered. The proposed syst…
▽ More
Sign Language Recognition has emerged as one of the important area of research in Computer Vision. The difficulty faced by the researchers is that the instances of signs vary with both motion and appearance. Thus, in this paper a novel approach for recognizing various alphabets of Indian Sign Language is proposed where continuous video sequences of the signs have been considered. The proposed system comprises of three stages: Preprocessing stage, Feature Extraction and Classification. Preprocessing stage includes skin filtering, histogram matching. Eigen values and Eigen Vectors were considered for feature extraction stage and finally Eigen value weighted Euclidean distance is used to recognize the sign. It deals with bare hands, thus allowing the user to interact with the system in natural way. We have considered 24 different alphabets in the video sequences and attained a success rate of 96.25%.
△ Less
Submitted 6 June, 2013;
originally announced June 2013.
-
Indian Sign Language Recognition Using Eigen Value Weighted Euclidean Distance Based Classification Technique
Authors:
Joyeeta Singha,
Karen Das
Abstract:
Sign Language Recognition is one of the most growing fields of research today. Many new techniques have been developed recently in these fields. Here in this paper, we have proposed a system using Eigen value weighted Euclidean distance as a classification technique for recognition of various Sign Languages of India. The system comprises of four parts: Skin Filtering, Hand Cropping, Feature Extrac…
▽ More
Sign Language Recognition is one of the most growing fields of research today. Many new techniques have been developed recently in these fields. Here in this paper, we have proposed a system using Eigen value weighted Euclidean distance as a classification technique for recognition of various Sign Languages of India. The system comprises of four parts: Skin Filtering, Hand Cropping, Feature Extraction and Classification. Twenty four signs were considered in this paper, each having ten samples, thus a total of two hundred forty images was considered for which recognition rate obtained was 97 percent.
△ Less
Submitted 4 March, 2013;
originally announced March 2013.
