Showing 1–50 of 69 results for author: Sanidas, S

A study of two FRBs with low polarization fractions localized with the MeerTRAP transient buffer system

Authors: K. M. Rajwade, L. N. Driessen, E. D. Barr, I. Pastor-Marazuela, M. Berezina, F. Jankowski, A. Muller, L. Kahinga, B. W. Stappers, M. C. Bezuidenhout, M. Caleb, A. Deller, W. Fong, A. Gordon, M. Kramer, M. Malenta, V. Morello, J. X. Prochaska, S. Sanidas, M. Surnis, N. Tejos, S. Wagner

Abstract: Localisation of fast radio bursts (FRBs) to arcsecond and sub-arcsecond precision maximizes their potential as cosmological probes. To that end, FRB detection instruments are deploying triggered complex-voltage capture systems to localize FRBs, identify their host galaxy and measure a redshift. Here, we report the discovery and localisation of two FRBs (20220717A and 20220905A) that were captured… ▽ More Localisation of fast radio bursts (FRBs) to arcsecond and sub-arcsecond precision maximizes their potential as cosmological probes. To that end, FRB detection instruments are deploying triggered complex-voltage capture systems to localize FRBs, identify their host galaxy and measure a redshift. Here, we report the discovery and localisation of two FRBs (20220717A and 20220905A) that were captured by the transient buffer system deployed by the MeerTRAP instrument at the MeerKAT telescope in South Africa. We were able to localize the FRBs to a precision of $\sim$1 arc-second that allowed us to unambiguously identify the host galaxy for FRB 20220717A (posterior probability$\sim$0.97). FRB 20220905A lies in a crowded region of the sky with a tentative identification of a host galaxy but the faintness and the difficulty in obtaining an optical spectrum preclude a conclusive association. The bursts show low linear polarization fractions (10--17$\%$) that conform to the large diversity in the polarization fraction observed in apparently non-repeating FRBs akin to single pulses from neutron stars. We also show that the host galaxy of FRB 20220717A contributes roughly 15$\%$ of the total dispersion measure (DM), indicating that it is located in a plasma-rich part of the host galaxy which can explain the large rotation measure. The scattering in FRB 20220717A can be mostly attributed to the host galaxy and the intervening medium and is consistent with what is seen in the wider FRB population. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: 13 pages, 8 figures, 1 Appendix, revised version after addressing the reviewer's comments

Modeling non stationary noise in pulsar timing array data analysis

Authors: Mikel Falxa, J. Antoniadis, D. J. Champion, I. Cognard, G. Desvignes, L. Guillemot, H. Hu, G. Janssen, J. Jawor, R. Karuppusamy, M. J. Keith, M. Kramer, K. Lackeos, K. Liu, J. W. McKee, D. Perrodin, S. A. Sanidas, G. M. Shaifullah, G. Theureau

Abstract: Pulsar Timing Array (PTA) collaborations recently reported evidence for the presence of a gravitational wave background (GWB) in their datasets. The main candidate that is expected to produce such a GWB is the population of supermassive black hole binaries (SMBHB). Some analyses showed that the recovered signal may exhibit time-dependent properties, i.e. non-stationarity. In this paper, we propose… ▽ More Pulsar Timing Array (PTA) collaborations recently reported evidence for the presence of a gravitational wave background (GWB) in their datasets. The main candidate that is expected to produce such a GWB is the population of supermassive black hole binaries (SMBHB). Some analyses showed that the recovered signal may exhibit time-dependent properties, i.e. non-stationarity. In this paper, we propose an approximated non-stationary Gaussian process (GP) model obtained from the perturbation of stationary processes. The presented method is applied to the second data release of the European pulsar timing array to search for non-stationary features in the GWB. We analyzed the data in different time slices and showed that the inferred properties of the GWB evolve with time. We find no evidence for such non-stationary behavior and the Bayes factor in favor of the latter is $\mathcal{B}^{NS}_{S} = 1.5$. We argue that the evolution of the GWB properties most likely comes from the \mf{improvement of the observation cadence} with time and \mf{better} characterization of the noise of individual pulsars. Such non-stationary GWB could also be produced by the leakage of non-stationary features in the noise of individual pulsars or by the presence of an eccentric single source. △ Less

Submitted 6 May, 2024; originally announced May 2024.

Comments: 14 pages, 14 figures

Constraints on conformal ultralight dark matter couplings from the European Pulsar Timing Array

Authors: Clemente Smarra, Adrien Kuntz, Enrico Barausse, Boris Goncharov, Diana López Nacir, Diego Blas, Lijing Shao, J. Antoniadis, D. J. Champion, I. Cognard, L. Guillemot, H. Hu, M. Keith, M. Kramer, K. Liu, D. Perrodin, S. A. Sanidas, G. Theureau

Abstract: Millisecond pulsars are extremely precise celestial clocks: as they rotate, the beamed radio waves emitted along the axis of their magnetic field can be detected with radio telescopes, which allows for tracking subtle changes in the pulsars' rotation periods. A possible effect on the period of a pulsar is given by a potential coupling to dark matter, in cases where it is modeled with an "ultraligh… ▽ More Millisecond pulsars are extremely precise celestial clocks: as they rotate, the beamed radio waves emitted along the axis of their magnetic field can be detected with radio telescopes, which allows for tracking subtle changes in the pulsars' rotation periods. A possible effect on the period of a pulsar is given by a potential coupling to dark matter, in cases where it is modeled with an "ultralight" scalar field. In this paper, we consider a universal conformal coupling of the dark matter scalar to gravity, which in turn mediates an effective coupling between pulsars and dark matter. If the dark matter scalar field is changing in time, as expected in the Milky Way, this effective coupling produces a periodic modulation of the pulsar rotational frequency. By studying the time series of observed radio pulses collected by the European Pulsar Timing Array experiment, we present constraints on the coupling of dark matter, improving on existing bounds. These bounds can also be regarded as constraints on the parameters of scalar-tensor theories of the Fierz-Jordan-Brans-Dicke and Damour-Esposito-Farèse types in the presence of a (light) mass potential term. △ Less

Submitted 2 May, 2024; originally announced May 2024.

Comments: 17 pages, 5 figures

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.

Comments: 21 pages, 9 figures, submitted to ApJ

The second data release from the European Pulsar Timing Array: VI. Challenging the ultralight dark matter paradigm

Authors: Clemente Smarra, Boris Goncharov, Enrico Barausse, J. Antoniadis, S. Babak, A. -S. Bak Nielsen, C. G. Bassa, 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, G. Desvignes, M. Falxa, R. D. Ferdman, A. Franchini, J. R. Gair, E. Graikou, J. -M. Grie , et al. (46 additional authors not shown)

Abstract: 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 s… ▽ More 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 $ρ\lesssim 0.3$ GeV/cm$^3$. △ Less

Submitted 25 October, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 5 pages + acknowledgements + refs, 2 figures. Text and figures match the version published in PRL

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.

Comments: 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)

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.

Comments: 13 figures, 15 pages, accepted

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.

Comments: 20 pages, 6 figures, 9 tables

Journal ref: A&A 678, A49 (2023)

The second data release from the European Pulsar Timing Array I. The dataset and timing analysis

Authors: J. Antoniadis, S. Babak, A. -S. Bak Nielsen, C. G. Bassa, 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, G. Desvignes, M. Falxa, R. D. Ferdman, A. Franchini, J. R. Gair, B. Goncharov, E. Graikou, J. -M. Grießmeier, L. Guillemot, Y. J. Guo , et al. (44 additional authors not shown)

Abstract: 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 pu… ▽ More 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. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 29 pages, 9 figures, 13 tables, Astronomy & Astrophysics in press

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.

Comments: 21 pages, 14 figures, 4 appendix figures, accepted for publication in A&A

Journal ref: A&A 678, A50 (2023)

Variable Scintillation Arcs of Millisecond Pulsars observed with the Large European Array for Pulsars

Authors: R. A. Main, J. Antoniadis, S. Chen, I. Cognard, H. Hu, J. Jang, R. Karuppusamy, M. Kramer, K. Liu, Y. Liu, G. Mall, J. W. McKee, M. B. Mickaliger, D. Perrodin, S. A. Sanidas, B. W. Stappers, T. Sprenger, O. Wucknitz, C. G. Bassa, M. Burgay, R. Concu, M. Gaikwad, G. H. Janssen, K. J. Lee, A. Melis , et al. (4 additional authors not shown)

Abstract: We present the first large sample of scintillation arcs in millisecond pulsars, analysing 12 sources observed with the Large European Array for Pulsars (LEAP), and the Effelsberg 100\,m telescope. We estimate the delays from multipath propagation, measuring significant correlated changes in scattering timescales over a 10-year timespan. Many sources show compact concentrations of power in the seco… ▽ More We present the first large sample of scintillation arcs in millisecond pulsars, analysing 12 sources observed with the Large European Array for Pulsars (LEAP), and the Effelsberg 100\,m telescope. We estimate the delays from multipath propagation, measuring significant correlated changes in scattering timescales over a 10-year timespan. Many sources show compact concentrations of power in the secondary spectrum, which in PSRs J0613$-$0200 and J1600$-$3053 can be tracked between observations, and are consistent with compact scattering at fixed angular positions. Other sources such as PSRs J1643$-$1224 and J0621+1002 show diffuse, asymmetric arcs which are likely related to phase-gradients across the scattering screen. PSR B1937+21 shows at least three distinct screens which dominate at different times and evidence of varying screen axes or multi-screen interactions. We model annual and orbital arc curvature variations in PSR J0613$-$0200, providing a measurement of the longitude of ascending node, resolving the sense of the orbital inclination, where our best fit model is of a screen with variable axis of anisotropy over time, corresponding to changes in the scattering of the source. Unmodeled variations of the screen's axis of anisotropy are likely to be a limiting factor in determining orbital parameters with scintillation, requiring careful consideration of variable screen properties, or independent VLBI measurements. Long-term scintillation studies such as this serve as a complementary tool to pulsar timing, to measure a source of correlated noise for pulsar timing arrays, solve pulsar orbits, and to understand the astrophysical origin of scattering screens. △ Less

Submitted 23 June, 2023; originally announced June 2023.

Comments: Submitted to MNRAS, first revision. 18 pages, 16 figures

Practical approaches to analyzing PTA data: Cosmic strings with six pulsars

Authors: Hippolyte Quelquejay Leclere, Pierre Auclair, Stanislav Babak, Aurélien Chalumeau, Danièle A. Steer, J. Antoniadis, A. -S. Bak Nielsen, C. G. Bassa, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, D. J. Champion, S. Chanlaridis, S. Chen, I. Cognard, G. Desvignes, M. Falxa, R. D. Ferdman, A. Franchini, J. R. Gair, B. Goncharov, E. Graikou , et al. (47 additional authors not shown)

Abstract: 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 blac… ▽ More 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μ) < -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. △ Less

Submitted 3 May, 2024; v1 submitted 21 June, 2023; originally announced June 2023.

Comments: 14 pages, 6 figures; typo corrected in (5)

Journal ref: Phys. Rev. D 108 (2023), 123527

Discovery of an Extremely Intermittent Periodic Radio Source

Authors: M. P. Surnis, K. M. Rajwade, B. W. Stappers, G. Younes, M. C. Bezuidenhout, M. Caleb, L. N. Driessen, F. Jankowski, M. Malenta, V. Morello, S. Sanidas, E. Barr, M. Kramer, R. Fender, P. Woudt

Abstract: We report the serendipitous discovery of an extremely intermittent radio pulsar, PSR J1710-3452, with a relatively long spin period of 10.4 s. The object was discovered through the detection of 97 bright radio pulses in only one out of 66 epochs of observations spanning almost three years. The bright pulses have allowed the source to be localised to a precision of 0.5" through radio imaging. We ob… ▽ More We report the serendipitous discovery of an extremely intermittent radio pulsar, PSR J1710-3452, with a relatively long spin period of 10.4 s. The object was discovered through the detection of 97 bright radio pulses in only one out of 66 epochs of observations spanning almost three years. The bright pulses have allowed the source to be localised to a precision of 0.5" through radio imaging. We observed the source location with the Swift X-ray telescope but did not detect any significant X-ray emission. We did not identify any high-energy bursts or multi-frequency counterparts for this object. The solitary epoch of detection hinders the calculation of the surface magnetic field strength, but the long period and the microstructure in the single-pulses resembles the emission of radio-loud magnetars. If this is indeed a magnetar, it is located at a relatively high Galactic latitude (2.9 degree), making it potentially one of the oldest and the most intermittent magnetars known in the Galaxy. The very short activity window of this object is unique and may point towards a yet undetected population of long period, highly transient radio emitting neutron stars. △ Less

Submitted 26 June, 2023; v1 submitted 19 June, 2023; originally announced June 2023.

Comments: Accepted for publication in MNRAS Letters. Minor correction in the current version

A sample of Fast Radio Bursts discovered and localised with MeerTRAP at the MeerKAT telescope

Authors: F. Jankowski, M. C. Bezuidenhout, M. Caleb, L. N. Driessen, M. Malenta, V. Morello, K. M. Rajwade, S. Sanidas, B. W. Stappers, M. P. Surnis, E. D. Barr, W. Chen, M. Kramer, J. Wu, S. Buchner, M. Serylak, J. X. Prochaska

Abstract: We present a sample of well-localised Fast Radio Bursts (FRBs) discovered by the MeerTRAP project at the MeerKAT telescope in South Africa. We discovered the three FRBs in single coherent tied-array beams and localised them to an area of ~1 arcmin$^2$. We investigate their burst properties, scattering, repetition rates, and localisations in a multi-wavelength context. FRB 20201211A shows hints of… ▽ More We present a sample of well-localised Fast Radio Bursts (FRBs) discovered by the MeerTRAP project at the MeerKAT telescope in South Africa. We discovered the three FRBs in single coherent tied-array beams and localised them to an area of ~1 arcmin$^2$. We investigate their burst properties, scattering, repetition rates, and localisations in a multi-wavelength context. FRB 20201211A shows hints of scatter broadening but is otherwise consistent with instrumental dispersion smearing. For FRB 20210202D, we discovered a faint post-cursor burst separated by ~200 ms, suggesting a distinct burst component or a repeat pulse. We attempt to associate the FRBs with host galaxy candidates. For FRB 20210408H, we tentatively (0.35 - 0.53 probability) identify a compatible host at a redshift ~0.5. Additionally, we analyse the MeerTRAP survey properties, such as the survey coverage, fluence completeness, and their implications for the FRB population. Based on the entire sample of 11 MeerTRAP FRBs discovered by the end of 2021, we estimate the FRB all-sky rates and their scaling with the fluence threshold. The inferred FRB all-sky rates at 1.28 GHz are $8.2_{-4.6}^{+8.0}$ and $2.1_{-1.1}^{+1.8} \times 10^3$ sky$^{-1}$ d$^{-1}$ above 0.66 and 3.44 Jy ms for the coherent and incoherent surveys, respectively. The scaling between the MeerTRAP rates is flatter than at higher fluences at the 1.4-$σ$ level. There seems to be a deficit of low-fluence FRBs, suggesting a break or turn-over in the rate versus fluence relation below 2 Jy ms. We speculate on cosmological or progenitor-intrinsic origins. The cumulative source counts within our surveys appear consistent with the Euclidean scaling. △ Less

Submitted 5 July, 2023; v1 submitted 20 February, 2023; originally announced February 2023.

Comments: 23 pages, 10 figures, 3 tables. Accepted for publication in MNRAS. Supplementary data are available on Zenodo (https://doi.org/10.5281/zenodo.6047539) and via the Transient Name Server

FRB 20210405I: a nearby Fast Radio Burst localised to sub-arcsecond precision with MeerKAT

Authors: Laura Nicole Driessen, Ewan Barr, David Buckley, Manisha Caleb, Hao Chen, Weiwei Chen, Mariusz Gromadzki, Fabian Jankowski, Renee Kraan-Korteweg, Michael Kramer, Jesse Palmerio, Kaustubh Rajwade, Ben Stappers, Evangelia Tremou, Susanna Vergani, Patrick Woudt, Mechiel Christiaan Bezuidenhout, Mateusz Malenta, Vincent Morello, Sotiris Sanidas, Mayuresh Surnis, Rob Fender

Abstract: We present the first sub-arcsecond localised Fast Radio Burst (FRB) detected using MeerKAT. FRB 20210405I was detected in the incoherent beam using the MeerTRAP pipeline on 2021 April 05 with a signal to noise ratio of 140.8 and a dispersion measure of 565.17 pc cm$^{-3}$. It was detected while MeerTRAP was observing commensally with the ThunderKAT large survey project, and was sufficiently bright… ▽ More We present the first sub-arcsecond localised Fast Radio Burst (FRB) detected using MeerKAT. FRB 20210405I was detected in the incoherent beam using the MeerTRAP pipeline on 2021 April 05 with a signal to noise ratio of 140.8 and a dispersion measure of 565.17 pc cm$^{-3}$. It was detected while MeerTRAP was observing commensally with the ThunderKAT large survey project, and was sufficiently bright that we could use the ThunderKAT 8s images to localise the FRB. Two different models of the dispersion measure in the Milky Way and halo suggest that the source is either right at the edge of the Galaxy, or outside. This highlights the uncertainty in the Milky Way dispersion measure models, particularly in the Galactic Plane, and the uncertainty of Milky Way halo models. Further investigation and modelling of these uncertainties will be facilitated by future detections and localisations of nearby FRBs. We use the combined localisation, dispersion measure, scattering, specific luminosity and chance coincidence probability information to find that the origin is most likely extra-galactic and identify the likely host galaxy of the FRB: 2MASS J1701249$-$4932475. Using SALT spectroscopy and archival observations of the field, we find that the host is a disk/spiral galaxy at a redshift of $z=0.066$. △ Less

Submitted 29 October, 2023; v1 submitted 20 February, 2023; originally announced February 2023.

Comments: 15 pages, 4 tables, 10 figures. Accepted to MNRAS

A sub-arcsec localised fast radio burst with a significant host galaxy dispersion measure contribution

Authors: M. Caleb, L. N. Driessen, A. C. Gordon, N. Tejos, L. Bernales, H. Qiu, J. O. Chibueze, B. W. Stappers, K. M. Rajwade, F. Cavallaro, Y. Wang, P. Kumar, W. A. Majid, R. S. Wharton, C. J. Naudet, M. C. Bezuidenhout, F. Jankowski, M. Malenta, V. Morello, S. Sanidas, M. P. Surnis, E. D. Barr, W. Chen, M. Kramer, W. Fong , et al. (7 additional authors not shown)

Abstract: We present the discovery of FRB 20210410D, with the MeerKAT radio interferometer in South Africa, as part of the MeerTRAP commensal project. FRB 20210410D has a dispersion measure DM = 578.78 +/- 2 pc cm-3, and was localised to sub-arcsec precision in the 2s images made from the correlation data products. The localisation enabled the association of the FRB with an optical galaxy at z = 0.1415, whi… ▽ More We present the discovery of FRB 20210410D, with the MeerKAT radio interferometer in South Africa, as part of the MeerTRAP commensal project. FRB 20210410D has a dispersion measure DM = 578.78 +/- 2 pc cm-3, and was localised to sub-arcsec precision in the 2s images made from the correlation data products. The localisation enabled the association of the FRB with an optical galaxy at z = 0.1415, which when combined with the DM places it above the 3sigma scatter of the Macquart relation. We attribute the excess DM to the host galaxy after accounting for contributions from the Milky Way's interstellar medium and halo, and the combined effects of the intergalactic medium and intervening galaxies. This is the first FRB that is not associated with a dwarf galaxy, to exhibit a likely large host galaxy DM contribution. We do not detect any continuum radio emission at the FRB position or from the host galaxy down to a 3sigma RMS of 14.4 uJy/beam. The FRB has a scattering delay of 29.4 ms at 1 GHz, and exhibits candidate subpulses in the spectrum, which hint at the possibility of it being a repeating FRB. Although not constraining, we note that this FRB has not been seen to repeat in 7.28h at 1.3 GHz with MeerKAT, 3h at 2.4 GHz with Murriyang and 5.7h at simultaneous 2.3 GHz and 8.4 GHz observations with the Deep Space Network. We encourage further follow-up to establish a possible repeating nature. △ Less

Submitted 13 June, 2023; v1 submitted 19 February, 2023; originally announced February 2023.

Comments: Accepted for publication in MNRAS. 14 pages, 10 figures

Detection of quasi-periodic micro-structure in three millisecond pulsars with the Large European Array for Pulsars

Authors: K. Liu, J. Antoniadis, C. G. Bassa, S. Chen, I. Cognard, M. Gaikwad, H. Hu, J. Jang, G. H. Janssen, R. Karuppusamy, M. Kramer, K. J. Lee, R. A. Main, G. Mall, J. W. McKee, M. B. Mickaliger, D. Perrodin, S. A. Sanidas, B. W. Stappers, L. Wang, W. W. Zhu, M. Burgay, R. Concu, A. Corongiu, A. Melis , et al. (2 additional authors not shown)

Abstract: We report on the detection of quasi-periodic micro-structure in three millisecond pulsars (MSPs), PSRs J1022+1001, J2145-0750 and J1744-1134, using high time resolution data acquired with the Large European Array for Pulsars at a radio frequency of 1.4 GHz. The occurrence rate of quasi-periodic micro-structure is consistent among pulses with different peak flux densities. Using an auto-correlation… ▽ More We report on the detection of quasi-periodic micro-structure in three millisecond pulsars (MSPs), PSRs J1022+1001, J2145-0750 and J1744-1134, using high time resolution data acquired with the Large European Array for Pulsars at a radio frequency of 1.4 GHz. The occurrence rate of quasi-periodic micro-structure is consistent among pulses with different peak flux densities. Using an auto-correlation analysis, we measure the periodicity and width of the micro-structure in these three pulsars. The detected micro-structure from PSRs J1022+1001 and J1744-1134 is often highly linearly polarised. In PSR J1022+1001, the linear polarisation position angles of micro-structure pulses are in general flat with a small degree of variation. Using these results, we further examine the frequency and rotational period dependency of micro-structure properties established in previous work, along with the angular beaming and temporal modulation models that explains the appearance of micro-structure. We also discuss a possible link of micro-structure to the properties of some of the recently discovered fast radio bursts which exhibit a very similar emission morphology. △ Less

Submitted 20 June, 2022; originally announced June 2022.

Comments: 8 pages, 6 figures, 1 table, published in MNRAS

Discovery of a radio emitting neutron star with an ultra-long spin period of 76 seconds

Authors: Manisha Caleb, Ian Heywood, Kaustubh Rajwade, Mateusz Malenta, Benjamin Stappers, Ewan Barr, Weiwei Chen, Vincent Morello, Sotiris Sanidas, Jakob van den Eijnden, Michael Kramer, David Buckley, Jaco Brink, Sara Elisa Motta, Patrick Woudt, Patrick Weltevrede, Fabian Jankowski, Mayuresh Surnis, Sarah Buchner, Mechiel Christiaan Bezuidenhout, Laura Nicole Driessen, Rob Fender

Abstract: The radio-emitting neutron star population encompasses objects with spin periods ranging from milliseconds to tens of seconds. As they age and spin more slowly, their radio emission is expected to cease. We present the discovery of an ultra-long period radio-emitting neutron star, J0901-4046, with spin properties distinct from the known spin and magnetic-decay powered neutron stars. With a spin-pe… ▽ More The radio-emitting neutron star population encompasses objects with spin periods ranging from milliseconds to tens of seconds. As they age and spin more slowly, their radio emission is expected to cease. We present the discovery of an ultra-long period radio-emitting neutron star, J0901-4046, with spin properties distinct from the known spin and magnetic-decay powered neutron stars. With a spin-period of 75.88 s, a characteristic age of 5.3 Myr, and a narrow pulse duty-cycle, it is uncertain how radio emission is generated and challenges our current understanding of how these systems evolve. The radio emission has unique spectro-temporal properties such as quasi-periodicity and partial nulling that provide important clues to the emission mechanism. Detecting similar sources is observationally challenging, which implies a larger undetected population. Our discovery establishes the existence of ultra-long period neutron stars, suggesting a possible connection to the evolution of highly magnetized neutron stars, ultra-long period magnetars, and fast radio bursts △ Less

Submitted 2 June, 2022; originally announced June 2022.

Comments: Published in Nature Astronomy - https://www.nature.com/articles/s41550-022-01688-x

First discoveries and localisations of Fast Radio Bursts with MeerTRAP: a real-time, commensal MeerKAT survey

Authors: K. M. Rajwade, M. C. Bezuidenhout, M. Caleb, L. N. Driessen, F. Jankowski, M. Malenta, V. Morello, S. Sanidas, B. W. Stappers, M. P. Surnis, E. D. Barr, W. Chen, M. Kramer, J. Wu, S. Buchner, M. Serylak, F. Combes, W. Fong, N. Gupta, P. Jagannathan, C. D. Kilpatrick, J. -K. Krogager, P. Noterdaeme, C. Núnez, J. Xavier Prochaska , et al. (2 additional authors not shown)

Abstract: We report on the discovery and localization of fast radio bursts (FRBs) from the MeerTRAP project, a commensal fast radio transient-detection programme at MeerKAT in South Africa. Our hybrid approach combines a coherent search with an average field-of-view of 0.4 $\rm deg^{2}$ with an incoherent search utilizing a field-of-view of $\sim$1.27 $\rm deg^{2}$ (both at 1284~MHz). Here, we present resul… ▽ More We report on the discovery and localization of fast radio bursts (FRBs) from the MeerTRAP project, a commensal fast radio transient-detection programme at MeerKAT in South Africa. Our hybrid approach combines a coherent search with an average field-of-view of 0.4 $\rm deg^{2}$ with an incoherent search utilizing a field-of-view of $\sim$1.27 $\rm deg^{2}$ (both at 1284~MHz). Here, we present results on the first three FRBs: FRB 20200413A (DM=1990.05 pc cm$^{-3}$), FRB 20200915A (DM=740.65 pc cm$^{-3}$), and FRB 20201123A (DM=433.55 pc cm$^{-3}$). FRB 20200413A was discovered only in the incoherent beam. FRB 20200915A (also discovered only in the incoherent beam) shows speckled emission in the dynamic spectrum which cannot be explained by interstellar scintillation in our Galaxy or plasma lensing, and might be intrinsic to the source. FRB 20201123A shows a faint post-cursor burst about 200 ms after the main burst and warrants further follow-up to confirm whether it is a repeating FRB. FRB 20201123A also exhibits significant temporal broadening consistent with scattering by a turbulent medium. The broadening exceeds that predicted for medium along the sightline through our Galaxy. We associate this scattering with the turbulent medium in the environment of the FRB in the host galaxy. Within the approximately $1'$ localization region of FRB 20201123A, we identify one luminous galaxy ($r \approx 15.67$; J173438.35$-$504550.4) that dominates the posterior probability for a host association. The galaxy's measured properties are consistent with other FRB hosts with secure associations. △ Less

Submitted 29 May, 2022; originally announced May 2022.

Comments: 14 pages, 9 figures, 3 tables, accepted for publication in MNRAS

MeerTRAP: Twelve Galactic fast transients detected in a real-time, commensal MeerKAT survey

Authors: M. C. Bezuidenhout, E. Barr, M. Caleb, L. N. Driessen, F. Jankowski, M. Kramer, M. Malenta, V. Morello, K. Rajwade, S. Sanidas, B. W. Stappers, M. Surnis

Abstract: MeerTRAP is a real-time untargeted search project using the MeerKAT telescope to find single pulses from fast radio transients and pulsars. It is performed commensally with the MeerKAT large survey projects (LSPs), using data from up to 64 of MeerKAT's 13.96~m dishes to form hundreds of coherent beams on sky, each of which is processed in real time to search for millisecond-duration pulses. We pre… ▽ More MeerTRAP is a real-time untargeted search project using the MeerKAT telescope to find single pulses from fast radio transients and pulsars. It is performed commensally with the MeerKAT large survey projects (LSPs), using data from up to 64 of MeerKAT's 13.96~m dishes to form hundreds of coherent beams on sky, each of which is processed in real time to search for millisecond-duration pulses. We present the first twelve Galactic sources discovered by MeerTRAP, with DMs in the range of 33--381~pc~cm$^{-3}$. One source may be Galactic or extragalactic depending on the Galactic electron density model assumed. Follow-up observations performed with the MeerKAT, Lovell, and Parkes radio telescopes have detected repeat pulses from seven of the twelve sources. Pulse periods have been determined for four sources. Another four sources could be localised to the arcsecond-level using a novel implementation of the tied-array beam localisation method. △ Less

Submitted 4 March, 2022; v1 submitted 1 March, 2022; originally announced March 2022.

Comments: 16 pages, 14 figures

Modelling annual scintillation arc variations in PSR J1643-1224 using the Large European Array for Pulsars

Authors: G. Mall, R. A. Main, J. Antoniadis, C. G. Bassa, M. Burgay, S. Chen, I. Cognard, R. Concu, A. Corongiu, M. Gaikwad, H. Hu, G. H. Janssen, R. Karuppusamy, K. J. Lee, K. Liu, J. W. McKee, A. Melis, M. B. Mickaliger, D. Perrodin, M. Pilia, A. Possenti, D. J. Reardon, S. A. Sanidas, T. Sprenger, B. W. Stappers , et al. (3 additional authors not shown)

Abstract: In this work we study variations in the parabolic scintillation arcs of the binary millisecond pulsar PSR J1643-1224 over five years using the Large European Array for Pulsars (LEAP). The 2D power spectrum of scintillation, called the secondary spectrum, often shows a parabolic distribution of power, where the arc curvature encodes the relative velocities and distances of the pulsar, ionised inter… ▽ More In this work we study variations in the parabolic scintillation arcs of the binary millisecond pulsar PSR J1643-1224 over five years using the Large European Array for Pulsars (LEAP). The 2D power spectrum of scintillation, called the secondary spectrum, often shows a parabolic distribution of power, where the arc curvature encodes the relative velocities and distances of the pulsar, ionised interstellar medium (IISM), and Earth. We observe a clear parabolic scintillation arc which varies in curvature throughout the year. The distribution of power in the secondary spectra are inconsistent with a single scattering screen which is fully 1D, or entirely isotropic. We fit the observed arc curvature variations with two models; an isotropic scattering screen, and a model with two independent 1D screens. We measure the distance to the scattering screen to be in the range 114-223 pc, depending on the model, consistent with the known distance of the foreground large-diameter HII region Sh 2-27 (112+/-17 pc), suggesting that it is the dominant source of scattering. We obtain only weak constraints on the pulsar's orbital inclination and angle of periastron, since the scintillation pattern is not very sensitive to the pulsar's motion, since the screen is much closer to the Earth than the pulsar. More measurements of this kind - where scattering screens can be associated with foreground objects - will help to inform the origins and distribution of scattering screens within our galaxy. △ Less

Submitted 11 January, 2022; originally announced January 2022.

Comments: 11 pages, 6 figures. Accepted by MNRAS

The International Pulsar Timing Array second data release: Search for an isotropic Gravitational Wave Background

Authors: J. Antoniadis, Z. Arzoumanian, S. Babak, M. Bailes, A. -S. Bak Nielsen, P. T. Baker, C. G. Bassa, B. Becsy, A. Berthereau, M. Bonetti, A. Brazier, P. R. Brook, M. Burgay, S. Burke-Spolaor, R. N. Caballero, J. A. Casey-Clyde, A. Chalumeau, D. J. Champion, M. Charisi, S. Chatterjee, S. Chen, I. Cognard, J. M. Cordes, N. J. Cornish, F. Crawford , et al. (101 additional authors not shown)

Abstract: We searched for an isotropic stochastic gravitational wave background in the second data release of the International Pulsar Timing Array, a global collaboration synthesizing decadal-length pulsar-timing campaigns in North America, Europe, and Australia. In our reference search for a power law strain spectrum of the form $h_c = A(f/1\,\mathrm{yr}^{-1})^α$, we found strong evidence for a spectrally… ▽ More We searched for an isotropic stochastic gravitational wave background in the second data release of the International Pulsar Timing Array, a global collaboration synthesizing decadal-length pulsar-timing campaigns in North America, Europe, and Australia. In our reference search for a power law strain spectrum of the form $h_c = A(f/1\,\mathrm{yr}^{-1})^α$, we found strong evidence for a spectrally-similar low-frequency stochastic process of amplitude $A = 3.8^{+6.3}_{-2.5}\times10^{-15}$ and spectral index $α= -0.5 \pm 0.5$, where the uncertainties represent 95\% credible regions, using information from the auto- and cross-correlation terms between the pulsars in the array. For a spectral index of $α= -2/3$, as expected from a population of inspiralling supermassive black hole binaries, the recovered amplitude is $A = 2.8^{+1.2}_{-0.8}\times10^{-15}$. Nonetheless, no significant evidence of the Hellings-Downs correlations that would indicate a gravitational-wave origin was found. We also analyzed the constituent data from the individual pulsar timing arrays in a consistent way, and clearly demonstrate that the combined international data set is more sensitive. Furthermore, we demonstrate that this combined data set produces comparable constraints to recent single-array data sets which have more data than the constituent parts of the combination. Future international data releases will deliver increased sensitivity to gravitational wave radiation, and significantly increase the detection probability. △ Less

Submitted 11 January, 2022; originally announced January 2022.

Comments: 17 pages, 12 figures, accepted in MNRAS

A MeerKAT, e-MERLIN, H.E.S.S. and Swift search for persistent and transient emission associated with three localised FRBs

Authors: James O. Chibueze, M. Caleb, L. Spitler, H. Ashkar, F. Schussler, B. W. Stappers, C. Venter, I. Heywood, A. M. S. Richards, D. R. A. Williams, M. Kramer, R. Beswick, M. C. Bezuidenhout, R. P. Breton, L. N. Driessen, F. Jankowski, E. F. Keane, M. Malenta, M. Mickaliger, V. Morello, H. Qiu, K. Rajwade, S. Sanidas, M. Surnis, T. W. Scragg , et al. (134 additional authors not shown)

Abstract: We report on a search for persistent radio emission from the one-off Fast Radio Burst (FRB) 20190714A, as well as from two repeating FRBs, 20190711A and 20171019A, using the MeerKAT radio telescope. For FRB 20171019A we also conducted simultaneous observations with the High Energy Stereoscopic System (H.E.S.S.) in very high energy gamma rays and searched for signals in the ultraviolet, optical, an… ▽ More We report on a search for persistent radio emission from the one-off Fast Radio Burst (FRB) 20190714A, as well as from two repeating FRBs, 20190711A and 20171019A, using the MeerKAT radio telescope. For FRB 20171019A we also conducted simultaneous observations with the High Energy Stereoscopic System (H.E.S.S.) in very high energy gamma rays and searched for signals in the ultraviolet, optical, and X-ray bands. For this FRB, we obtain a UV flux upper limit of 1.39x10^-16 erg/cm^-2/s/Amstrong, X-ray limit of ~ 6.6x10^-14 erg/cm^-2/s and a limit on the very-high-energy gamma-ray flux (Phi) (E > 120 GeV) < 1.7 x 10^-12 erg/cm^-2/s. We obtain a radio upper limit of ~15 microJy/beam for persistent emission at the locations of both FRBs 20190711A and 20171019A, but detect diffuse radio emission with a peak brightness of ~53 microJy/beam associated with FRB 20190714A at z = 0.2365. This represents the first detection of the radio continuum emission potentially associated with the host (galaxy) of FRB 20190714A, and is only the third known FRB to have such an association. Given the possible association of a faint persistent source, FRB 20190714A may potentially be a repeating FRB whose age lies between that of FRB 20121102A and FRB 20180916A. A parallel search for repeat bursts from these FRBs revealed no new detections down to a fluence of 0.08 Jy ms for a 1 ms duration burst. △ Less

Submitted 31 December, 2021; originally announced January 2022.

Comments: 15 pages, 9 figures

Removal and replacement of interference in tied-array radio pulsar observations using the spectral kurtosis estimator

Authors: M. Purver, C. G. Bassa, I. Cognard, G. H. Janssen, R. Karuppusamy, M. Kramer, K. J. Lee, K. Liu, J. W. McKee, D. Perrodin, S. Sanidas, R. Smits, B. W. Stappers

Abstract: We describe how to implement the spectral kurtosis method of interference removal (zapping) on a digitized signal of averaged power values. Spectral kurtosis is a hypothesis test, analogous to the t-test, with a null hypothesis that the amplitudes from which power is formed belong to a `good' distribution -- typically Gaussian with zero mean -- where power values are zapped if the hypothesis is re… ▽ More We describe how to implement the spectral kurtosis method of interference removal (zapping) on a digitized signal of averaged power values. Spectral kurtosis is a hypothesis test, analogous to the t-test, with a null hypothesis that the amplitudes from which power is formed belong to a `good' distribution -- typically Gaussian with zero mean -- where power values are zapped if the hypothesis is rejected at a specified confidence level. We derive signal-to-noise ratios (SNRs) as a function of amount of zapping for folded radio pulsar observations consisting of a sum of signals from multiple telescopes in independent radio-frequency interference (RFI) environments, comparing four methods to compensate for lost data with coherent (tied-array) and incoherent summation. For coherently summed amplitudes, scaling amplitudes from non-zapped telescopes achieves a higher SNR than replacing zapped amplitudes with artificial noise. For incoherently summed power values, the highest SNR is given by scaling power from non-zapped telescopes to maintain a constant mean. We use spectral kurtosis to clean a tied-array radio pulsar observation by the Large European Array for Pulsars (LEAP): the signal from one telescope is zapped with time and frequency resolutions of 6.25 $μ$s and 0.16 MHz, removing interference along with 0.27 per cent of `good' data, giving an uncertainty of 0.25 $μ$s in pulse time of arrival (TOA) for PSR J1022+1001. We use a single-telescope observation to demonstrate recovery of the pulse profile shape, with 0.6 per cent of data zapped and a reduction from 1.22 to 0.70 $μ$s in TOA uncertainty. △ Less

Submitted 24 February, 2023; v1 submitted 28 November, 2021; originally announced November 2021.

Comments: 16 pages, 9 figures

Journal ref: MNRAS, 510, 1597 (2022)

Noise analysis in the European Pulsar Timing Array data release 2 and its implications on the gravitational-wave background search

Authors: A. Chalumeau, S. Babak, A. Petiteau, S. Chen, A. Samajdar, R. N. Caballero, G. Theureau, L. Guillemot, G. Desvignes, A. Parthasarathy, K. Liu, G. Shaifullah, H. Hu, E. van der Wateren, J. Antoniadis, A. -S. Bak Nielsen, C. G. Bassa, A. Berthereau, M. Burgay, D. J. Champion, I. Cognard, M. Falxa, R. D. Ferdman, P. C. C. Freire, J. R. Gair , et al. (27 additional authors not shown)

Abstract: The European Pulsar Timing Array (EPTA) collaboration has recently released an extended data set for six pulsars (DR2) and reported evidence for a common red noise signal. Here we present a noise analysis for each of the six pulsars. We consider several types of noise: (i) radio frequency independent, "achromatic", and time-correlated red noise; (ii) variations of dispersion measure and scattering… ▽ More The European Pulsar Timing Array (EPTA) collaboration has recently released an extended data set for six pulsars (DR2) and reported evidence for a common red noise signal. Here we present a noise analysis for each of the six pulsars. We consider several types of noise: (i) radio frequency independent, "achromatic", and time-correlated red noise; (ii) variations of dispersion measure and scattering; (iii) system and band noise; and (iv) deterministic signals (other than gravitational waves) that could be present in the PTA data. We perform Bayesian model selection to find the optimal combination of noise components for each pulsar. Using these custom models we revisit the presence of the common uncorrelated red noise signal previously reported in the EPTA DR2 and show that the data still supports it with a high statistical significance. Next, we confirm that there is no preference for or against the Hellings-Downs spatial correlations expected for the stochastic gravitational-wave background. The main conclusion of the EPTA DR2 paper remains unchanged despite a very significant change in the noise model of each pulsar. However, modelling the noise is essential for the robust detection of gravitational waves and its impact could be significant when analysing the next EPTA data release, which will include a larger number of pulsars and more precise measurements. △ Less

Submitted 9 November, 2021; originally announced November 2021.

Comments: 21 pages, 11 figures, 7 tables, 1 appendix figure and 1 appendix table, accepted for publication to MNRAS

Common-red-signal analysis with 24-yr high-precision timing of the European Pulsar Timing Array: Inferences in the stochastic gravitational-wave background search

Authors: S. Chen, R. N. Caballero, Y. J. Guo, A. Chalumeau, K. Liu, G. Shaifullah, K. J. Lee, S. Babak, G. Desvignes, A. Parthasarathy, H. Hu, E. van der Wateren, J. Antoniadis, A. -S. Bak Nielsen, C. G. Bassa, A. Berthereau, M. Burgay, D. J. Champion, I. Cognard, M. Falxa, R. D. Ferdman, P. C. C. Freire, J. R. Gair, E. Graikou, L. Guillemot , et al. (27 additional authors not shown)

Abstract: We present results from the search for a stochastic gravitational-wave background (GWB) as predicted by the theory of General Relativity using six radio millisecond pulsars from the Data Release 2 (DR2) of the European Pulsar Timing Array (EPTA) covering a timespan up to 24 years. A GWB manifests itself as a long-term low-frequency stochastic signal common to all pulsars, a common red signal (CRS)… ▽ More We present results from the search for a stochastic gravitational-wave background (GWB) as predicted by the theory of General Relativity using six radio millisecond pulsars from the Data Release 2 (DR2) of the European Pulsar Timing Array (EPTA) covering a timespan up to 24 years. A GWB manifests itself as a long-term low-frequency stochastic signal common to all pulsars, a common red signal (CRS), with the characteristic Hellings-Downs (HD) spatial correlation. Our analysis is performed with two independent pipelines, \eprise{} and \tn{}+\ftwo{}, which produce consistent results. A search for a CRS with simultaneous estimation of its spatial correlations yields spectral properties compatible with theoretical GWB predictions, but does not result in the required measurement of the HD correlation, as required for GWB detection. Further Bayesian model comparison between different types of CRSs, including a GWB, finds the most favoured model to be the common uncorrelated red noise described by a power-law with $A = 5.13_{-2.73}^{+4.20} \times 10^{-15}$ and $γ= 3.78_{-0.59}^{+0.69}$ (95\% credible regions). Fixing the spectral index to $γ=13/3$ as expected from the GWB by circular, inspiralling supermassive black-hole binaries results in an amplitude of $A =2.95_{-0.72}^{+0.89} \times 10^{-15}$. We implement three different models, BAYESEPHEM, LINIMOSS and EPHEMGP, to address possible Solar-system ephemeris (SSE) systematics and conclude that our results may only marginally depend on these effects. This work builds on the methods and models from the studies on the EPTA DR1. We show that under the same analysis framework the results remain consistent after the data set extension. △ Less

Submitted 25 October, 2021; originally announced October 2021.

Comments: 25 pages, 10 figures, 7 tables, 2 appendix tables and 1 appendix figure

Quasi-Simultaneous Radio/X-ray Observations of the Candidate Transitional Millisecond Pulsar 3FGL J1544.6-1125 During its Low-Luminosity Accretion-Disc State

Authors: Amruta Jaodand, Adam T. Deller, Nina Gusinskaia, Jason W. T. Hessels, James C. A. Miller-Jones, Anne M. Archibald, Slavko Bogdanov, Cees Bassa, Rudy Wijnands, Alessandro Patruno, Sotiris Sanidas

Abstract: 3FGL J1544.6-1125 is a candidate transitional millisecond pulsar (tMSP). Similar to the well-established tMSPs - PSR J1023+0038, IGR J18245-2452, and XSS J12270-4859 -- 3FGL J1544.6-1125 shows $γ$-ray emission and discrete X-ray "low" and "high" modes during its low-luminosity accretion-disk state. Coordinated radio/X-ray observations of PSR J1023+0038 in its current low-luminosity accretion-disk… ▽ More 3FGL J1544.6-1125 is a candidate transitional millisecond pulsar (tMSP). Similar to the well-established tMSPs - PSR J1023+0038, IGR J18245-2452, and XSS J12270-4859 -- 3FGL J1544.6-1125 shows $γ$-ray emission and discrete X-ray "low" and "high" modes during its low-luminosity accretion-disk state. Coordinated radio/X-ray observations of PSR J1023+0038 in its current low-luminosity accretion-disk state showed rapidly variable radio continuum emission-possibly originating from a compact, self-absorbed jet, the "propellering" of accretion material, and/or pulsar moding. 3FGL J1544.6-1125 is currently the only other (candidate) tMSP system in this state, and can be studied to see whether tMSPs are typically radio-loud compared to other neutron star binaries. In this work, we present a quasi-simultaneous Very Large Array and Swift radio/X-ray campaign on 3FGL J1544.6-1125. We detect 10 GHz radio emission varying in flux density from $47.7 \pm 6.0$ $μ$Jy down to $\sim$15 $μ$Jy (3$σ$ upper limit) at four epochs spanning three weeks. At the brightest epoch, the radio luminosity is $L_{5 GHz}$ $= (2.17 \pm 0.17) \times 10^{27}$ erg s$^{-1}$ for a quasi-simultaneous X-ray luminosity $L_{2-10 keV}$ $= (4.32 \pm 0.23) \times 10^{33}$ erg s$^{-1}$ (for an assumed distance of 3.8 kpc). These luminosities are close to those of PSR J1023+0038, and the results strengthen the case that 3FGL J1544.6-1125 is a tMSP showing similar phenomenology to PSR J1023+0038. △ Less

Submitted 20 October, 2021; originally announced October 2021.

Comments: Accepted for publication in the Astrophysical Journal

MeerTRAP in the era of multi-messenger astrophysics

Authors: Kaustubh Rajwade, Benjamin Stappers, Christopher Williams, Ewan Barr, Mechiel Christiaan Bezuidenhout, Manisha Caleb, Laura Driessen, Fabian Jankowski, Mateusz Malenta, Vincent Morello, Sotirios Sanidas, Mayuresh Surnis

Abstract: Real-time detections of transients and rapid multi-wavelength follow-up are at the core of modern multi-messenger astrophysics. MeerTRAP is one such instrument that has been deployed on the MeerKAT radio telescope in South Africa to search for fast radio transients in real-time. This, coupled with the ability to rapidly localize the transient in combination with optical co-pointing by the MeerLICH… ▽ More Real-time detections of transients and rapid multi-wavelength follow-up are at the core of modern multi-messenger astrophysics. MeerTRAP is one such instrument that has been deployed on the MeerKAT radio telescope in South Africa to search for fast radio transients in real-time. This, coupled with the ability to rapidly localize the transient in combination with optical co-pointing by the MeerLICHT telescope gives the instrument the edge in finding and identifying the nature of the transient on short timescales. The commensal nature of the project means that MeerTRAP will keep looking for transients even if the telescope is not being used specifically for that purpose. Here, we present a brief overview of the MeerTRAP project. We describe the overall design, specifications and the software stack required to implement such an undertaking. We conclude with some science highlights that have been enabled by this venture over the last 10 months of operation. △ Less

Submitted 15 March, 2021; originally announced March 2021.

Comments: 8 pages, 5 figures, published in the proceedings of SPIE Conference on Astronomy and Telescope Instrumentation, No. 114470J

Real-time triggering capabilities for Fast Radio Bursts at the MeerKAT telescope

Authors: F. Jankowski, M. Berezina, B. W. Stappers, E. D. Barr, M. C. Bezuidenhout, M. Caleb, L. N. Driessen, M. Malenta, V. Morello, K. M. Rajwade, S. Sanidas, M. P. Surnis

Abstract: Fast Radio Bursts (FRBs) are bright enigmatic radio pulses of roughly millisecond duration that come from extragalactic distances. As part of the MeerTRAP project, we use the MeerKAT telescope array in South Africa to search for and localise those bursts to high precision in real-time. We aim to pinpoint FRBs to their host galaxies and, thereby, to understand how they are created. However, the tra… ▽ More Fast Radio Bursts (FRBs) are bright enigmatic radio pulses of roughly millisecond duration that come from extragalactic distances. As part of the MeerTRAP project, we use the MeerKAT telescope array in South Africa to search for and localise those bursts to high precision in real-time. We aim to pinpoint FRBs to their host galaxies and, thereby, to understand how they are created. However, the transient nature of FRBs presents various challenges, e.g. in system design, raw compute power and real-time communication, where the real-time requirements are reasonably strict (a few tens of seconds). Rapid data processing is essential for us to be able to retain high-resolution data of the bursts, to localise them, and to minimise the delay for follow-up observations. We give a short overview of the data analysis pipeline, describe the challenges faced, and elaborate on our initial design and implementation of a real-time triggering infrastructure for FRBs at the MeerKAT telescope. △ Less

Submitted 9 December, 2020; originally announced December 2020.

Comments: 4 pages, 1 figure, proceedings of the ADASS XXX (2020) conference, to appear in ASP Conference Series

Measuring Interstellar Delays of PSR J0613-0200 over 7 years, using the Large European Array for Pulsars

Authors: R. A. Main, S. A. Sanidas, J. Antoniadis, C. Bassa, S. Chen, I. Cognard, M. Gaikwad, H. Hu, G. H. Janssen, R. Karuppusamy, M. Kramer, K. J. Lee, K. Liu, G. Mall, J. W. McKee, M. B. Mickaliger, D. Perrodin, B. W. Stappers, C. Tiburzi, O. Wucknitz, L. Wang, W. W. Zhu

Abstract: Using data from the Large European Array for Pulsars (LEAP), and the Effelsberg telescope, we study the scintillation parameters of the millisecond pulsar J0613-0200 over a 7 year timespan. The "secondary spectrum" -- the 2D power spectrum of scintillation -- presents the scattered power as a function of time delay, and contains the relative velocities of the pulsar, observer, and scattering mater… ▽ More Using data from the Large European Array for Pulsars (LEAP), and the Effelsberg telescope, we study the scintillation parameters of the millisecond pulsar J0613-0200 over a 7 year timespan. The "secondary spectrum" -- the 2D power spectrum of scintillation -- presents the scattered power as a function of time delay, and contains the relative velocities of the pulsar, observer, and scattering material. We detect a persistent parabolic scintillation arc, suggesting scattering is dominated by a thin, anisotropic region. The scattering is poorly described by a simple exponential tail, with excess power at high delays; we measure significant, detectable scattered power at times out to $\sim 5 μs$, and measure the bulk scattering delay to be between 50 to 200\,ns with particularly strong scattering throughout 2013. These delays are too small to detect a change of the pulse profile shape, yet they would change the times-of-arrival as measured through pulsar timing. The arc curvature varies annually, and is well fit by a one-dimensional scattering screen $\sim 40\%$ of the way towards the pulsar, with a changing orientation during the increased scattering in 2013. Effects of uncorrected scattering will introduce time delays correlated over time in individual pulsars, and may need to be considered in gravitational wave analyses. Pulsar timing programs would benefit from simultaneously recording in a way that scintillation can be resolved, in order to monitor the variable time delays caused by multipath propagation. △ Less

Submitted 3 November, 2020; v1 submitted 22 September, 2020; originally announced September 2020.

Comments: 12 pages, 10 figures. Published in MNRAS

Simultaneous multi-telescope observations of FRB 121102

Authors: M. Caleb, B. W. Stappers, T. D. Abbott, E. D. Barr, M. C. Bezuidenhout, S. J. Buchner, M. Burgay, W. Chen, I. Cognard, L. N. Driessen, R. Fender, G. H. Hilmarsson, J. Hoang, D. M. Horn, F. Jankowski, M. Kramer, D. R. Lorimer, M. Malenta, V. Morello, M. Pilia, E. Platts, A. Possenti, K. M. Rajwade, A. Ridolfi, L. Rhodes , et al. (7 additional authors not shown)

Abstract: We present 11 detections of FRB 121102 in ~3 hours of observations during its 'active' period on the 10th of September 2019. The detections were made using the newly deployed MeerTRAP system and single pulse detection pipeline at the MeerKAT radio telescope in South Africa. Fortuitously, the Nancay radio telescope observations on this day overlapped with the last hour of MeerKAT observations and r… ▽ More We present 11 detections of FRB 121102 in ~3 hours of observations during its 'active' period on the 10th of September 2019. The detections were made using the newly deployed MeerTRAP system and single pulse detection pipeline at the MeerKAT radio telescope in South Africa. Fortuitously, the Nancay radio telescope observations on this day overlapped with the last hour of MeerKAT observations and resulted in 4 simultaneous detections. The observations with MeerKAT's wide band receiver, which extends down to relatively low frequencies (900-1670 MHz usable L-band range), have allowed us to get a detailed look at the complex frequency structure, intensity variations and frequency-dependent sub-pulse drifting. The drift rates we measure for the full-band and sub-banded data are consistent with those published between 600-6500 MHz with a slope of -0.147 +/- 0.014 ms^-1. Two of the detected bursts exhibit fainter 'precursors' separated from the brighter main pulse by ~28 ms and ~34 ms. A follow-up multi-telescope campaign on the 6th and 8th October 2019 to better understand these frequency drifts and structures over a wide and continuous band was undertaken. No detections resulted, indicating that the source was 'inactive' over a broad frequency range during this time. △ Less

Submitted 15 June, 2020; originally announced June 2020.

Comments: Accepted for publication in MNRAS

The LOFAR Tied-Array All-Sky Survey: Timing of 21 pulsars including the first binary pulsar discovered with LOFAR

Authors: C. M. Tan, C. G. Bassa, S. Cooper, J. W. T. Hessels, V. I. Kondratiev, D. Michilli, S. Sanidas, B. W. Stappers, J. van Leeuwen, J. Y. Donner, J. -M. Grießmeier, M. Kramer, C. Tiburzi, P. Weltevrede, B. Ciardi, M. Hoeft, G. Mann, A. Miskolczi, D. J. Schwarz, C. Vocks, O. Wucknitz

Abstract: We report on the multi-frequency timing observations of 21 pulsars discovered in the LOFAR Tied-Array All-Sky Survey (LOTAAS). The timing data were taken at central frequencies of 149 MHz (LOFAR) as well as 334 and 1532 MHz (Lovell Telecope). The sample of pulsars includes 20 isolated pulsars and the first binary pulsar discovered by the survey, PSR J1658$+$3630. We modelled the timing properties… ▽ More We report on the multi-frequency timing observations of 21 pulsars discovered in the LOFAR Tied-Array All-Sky Survey (LOTAAS). The timing data were taken at central frequencies of 149 MHz (LOFAR) as well as 334 and 1532 MHz (Lovell Telecope). The sample of pulsars includes 20 isolated pulsars and the first binary pulsar discovered by the survey, PSR J1658$+$3630. We modelled the timing properties of the pulsars, which showed that they have, on average, larger characteristic ages. We present the pulse profiles of the pulsars across the three observing bands, where PSR J1643$+$1338 showed profile evolution that appears not to be well-described by the radius-to-frequency-mapping model. Furthermore, we modelled the spectra of the pulsars across the same observing bands, using a simple power law, and found an average spectral index of $-1.9 \pm 0.5$. Amongst the pulsars studied here, PSR J1657$+$3304 showed large flux density variations of a factor of 10 over 300 days, as well as mode changing and nulling on timescales of a few minutes. We modelled the rotational and orbital properties of PSR J1658$+$3630, which has a spin period of 33 ms in a binary orbit of 3.0 days with a companion of minimum mass of 0.87$M_{\odot}$, likely a Carbon-Oxygen or Oxygen-Neon-Magnesium type white dwarf. PSR J1658$+$3630 has a dispersion measure of 3.0 pc cm$^{-3}$, making it possibly one of the closest binary pulsars known. △ Less

Submitted 12 January, 2020; originally announced January 2020.

Comments: 21 pages, 12 figures, Accepted for publication in the Monthly Notices of the Royal Astronomical Society

A LOFAR radio search for single and periodic pulses from M31

Authors: Joeri van Leeuwen, Klim Mikhailov, Evan Keane, Thijs Coenen, Liam Connor, Vlad Kondratiev, Daniele Michilli, Sotiris Sanidas

Abstract: Bright, short radio bursts are emitted by sources at a large range of distances: from the nearby Crab pulsar to remote Fast Radio Bursts (FRBs). FRBs are likely to originate from distant neutron stars, but our knowledge of the radio pulsar population has been limited to the Galaxy and the Magellanic Clouds. In an attempt to increase our understanding of extragalactic pulsar populations, and its gi… ▽ More Bright, short radio bursts are emitted by sources at a large range of distances: from the nearby Crab pulsar to remote Fast Radio Bursts (FRBs). FRBs are likely to originate from distant neutron stars, but our knowledge of the radio pulsar population has been limited to the Galaxy and the Magellanic Clouds. In an attempt to increase our understanding of extragalactic pulsar populations, and its giant-pulse emission, we employed the low-frequency radio telescope LOFAR to search the Andromeda Galaxy (M31) for radio bursts emitted by young, Crab-like pulsars. For direct comparison we also present a LOFAR study on the low-frequency giant pulses from the Crab pulsar; their fluence distribution follows a power law with slope 3.04(3). A number of candidate signals were detected from M31 but none proved persistent. FRBs are sometimes thought of as Crab-like pulsars with exceedingly bright giant pulses -- given our sensitivity, we can rule out that M31 hosts pulsars more than an order of magnitude brighter than the Crab pulsar, assuming their pulse scattering follows that of the known FRBs. △ Less

Submitted 25 November, 2019; originally announced November 2019.

Comments: Accepted for publication in A&A. 6 pages with 4 nice figures

Journal ref: A&A 634, A3 (2020)

A pulsar-based timescale from the International Pulsar Timing Array

Authors: G. Hobbs, L. Guo, R. N. Caballero, W. Coles, K. J. Lee, R. N. Manchester, D. J. Reardon, D. Matsakis, M. L. Tong, Z. Arzoumanian, M. Bailes, C. G. Bassa, N. D. R. Bhat, A. Brazier, S. Burke-Spolaor, D. J. Champion, S. Chatterjee, I. Cognard, S. Dai, G. Desvignes, T. Dolch, R. D. Ferdman, E. Graikou, L. Guillemot, G. H. Janssen , et al. (34 additional authors not shown)

Abstract: We have constructed a new timescale, TT(IPTA16), based on observations of radio pulsars presented in the first data release from the International Pulsar Timing Array (IPTA). We used two analysis techniques with independent estimates of the noise models for the pulsar observations and different algorithms for obtaining the pulsar timescale. The two analyses agree within the estimated uncertainties… ▽ More We have constructed a new timescale, TT(IPTA16), based on observations of radio pulsars presented in the first data release from the International Pulsar Timing Array (IPTA). We used two analysis techniques with independent estimates of the noise models for the pulsar observations and different algorithms for obtaining the pulsar timescale. The two analyses agree within the estimated uncertainties and both agree with TT(BIPM17), a post-corrected timescale produced by the Bureau International des Poids et Mesures (BIPM). We show that both methods could detect significant errors in TT(BIPM17) if they were present. We estimate the stability of the atomic clocks from which TT(BIPM17) is derived using observations of four rubidium fountain clocks at the US Naval Observatory. Comparing the power spectrum of TT(IPTA16) with that of these fountain clocks suggests that pulsar-based timescales are unlikely to contribute to the stability of the best timescales over the next decade, but they will remain a valuable independent check on atomic timescales. We also find that the stability of the pulsar-based timescale is likely to be limited by our knowledge of solar-system dynamics, and that errors in TT(BIPM17) will not be a limiting factor for the primary goal of the IPTA, which is to search for the signatures of nano-Hertz gravitational waves. △ Less

Submitted 29 October, 2019; originally announced October 2019.

Comments: accepted by MNRAS

The LOFAR Tied-Array All-Sky Survey (LOTAAS): Characterization of 20 pulsar discoveries and their single-pulse behavior

Authors: D. Michilli, C. Bassa, S. Cooper, J. W. T. Hessels, V. I. Kondratiev, S. Sanidas, B. W. Stappers, C. M. Tan, J. van Leeuwen, I. Cognard, J. M. Griessmeier, A. G. Lyne, J. P. W. Verbiest, P. Weltevrede

Abstract: We are using the LOw-Frequency ARray (LOFAR) to perform the LOFAR Tied-Array All-Sky (LOTAAS) survey for pulsars and fast transients. Here we present the astrometric and rotational parameters of 20 pulsars discovered as part of LOTAAS. These pulsars have regularly been observed with LOFAR at 149 MHz and the Lovell telescope at 1532 MHz, supplemented by some observations with the Lovell telescope a… ▽ More We are using the LOw-Frequency ARray (LOFAR) to perform the LOFAR Tied-Array All-Sky (LOTAAS) survey for pulsars and fast transients. Here we present the astrometric and rotational parameters of 20 pulsars discovered as part of LOTAAS. These pulsars have regularly been observed with LOFAR at 149 MHz and the Lovell telescope at 1532 MHz, supplemented by some observations with the Lovell telescope at 334 MHz and the Nancay Radio Telescope at 1484 MHz. Timing models are calculated for the 20 pulsars, some of which are among the slowest-spinning pulsars known. PSR J1236-0159 rotates with a period P ~ 3.6 s, while 5 additional pulsars show P > 2 s. Also, the spin-down rates Pdot are, on average, low, with PSR J0815+4611 showing Pdot ~ 4E-18. Some of the pulse profiles, generically single-peaked, present complex shapes evolving with frequency. Multi-frequency flux measurements show that these pulsars have generically relatively steep spectra but exceptions are present, with values ranging between ~ -4 and -1. Among the pulsar sample, a large fraction shows large single-pulse variability, with 4 pulsars being undetectable more than 15% of the time and one tentatively classified as a Rotating Radio Transient. Two single-peaked pulsars show drifting sub-pulses. △ Less

Submitted 21 October, 2019; originally announced October 2019.

The International Pulsar Timing Array: Second data release

Authors: B. B. P. Perera, M. E. DeCesar, P. B. Demorest, M. Kerr, L. Lentati, D. J. Nice, S. Oslowski, S. M. Ransom, M. J. Keith, Z. Arzoumanian, M. Bailes, P. T. Baker, C. G. Bassa, N. D. R. Bhat, A. Brazier, M. Burgay, S. Burke-Spolaor, R. N. Caballero, D. J. Champion, S. Chatterjee, S. Chen, I. Cognard, J. M. Cordes, K. Crowter, S. Dai , et al. (50 additional authors not shown)

Abstract: In this paper, we describe the International Pulsar Timing Array second data release, which includes recent pulsar timing data obtained by three regional consortia: the European Pulsar Timing Array, the North American Nanohertz Observatory for Gravitational Waves, and the Parkes Pulsar Timing Array. We analyse and where possible combine high-precision timing data for 65 millisecond pulsars which a… ▽ More In this paper, we describe the International Pulsar Timing Array second data release, which includes recent pulsar timing data obtained by three regional consortia: the European Pulsar Timing Array, the North American Nanohertz Observatory for Gravitational Waves, and the Parkes Pulsar Timing Array. We analyse and where possible combine high-precision timing data for 65 millisecond pulsars which are regularly observed by these groups. A basic noise analysis, including the processes which are both correlated and uncorrelated in time, provides noise models and timing ephemerides for the pulsars. We find that the timing precisions of pulsars are generally improved compared to the previous data release, mainly due to the addition of new data in the combination. The main purpose of this work is to create the most up-to-date IPTA data release. These data are publicly available for searches for low-frequency gravitational waves and other pulsar science. △ Less

Submitted 10 September, 2019; originally announced September 2019.

Comments: Submitted to MNRAS and in review, 23 pages, 5 figures

Probing the gravitational wave background from cosmic strings with LISA

Authors: Pierre Auclair, Jose J. Blanco-Pillado, Daniel G. Figueroa, Alexander C. Jenkins, Marek Lewicki, Mairi Sakellariadou, Sotiris Sanidas, Lara Sousa, Daniele A. Steer, Jeremy M. Wachter, Sachiko Kuroyanagi

Abstract: Cosmic string networks offer one of the best prospects for detection of cosmological gravitational waves (GWs). The combined incoherent GW emission of a large number of string loops leads to a stochastic GW background (SGWB), which encodes the properties of the string network. In this paper we analyze the ability of the Laser Interferometer Space Antenna (LISA) to measure this background, consider… ▽ More Cosmic string networks offer one of the best prospects for detection of cosmological gravitational waves (GWs). The combined incoherent GW emission of a large number of string loops leads to a stochastic GW background (SGWB), which encodes the properties of the string network. In this paper we analyze the ability of the Laser Interferometer Space Antenna (LISA) to measure this background, considering leading models of the string networks. We find that LISA will be able to probe cosmic strings with tensions $Gμ\gtrsim \mathcal{O}(10^{-17})$, improving by about $6$ orders of magnitude current pulsar timing arrays (PTA) constraints, and potentially $3$ orders of magnitude with respect to expected constraints from next generation PTA observatories. We include in our analysis possible modifications of the SGWB spectrum due to different hypotheses regarding cosmic history and the underlying physics of the string network. These include possible modifications in the SGWB spectrum due to changes in the number of relativistic degrees of freedom in the early Universe, the presence of a non-standard equation of state before the onset of radiation domination, or changes to the network dynamics due to a string inter-commutation probability less than unity. In the event of a detection, LISA's frequency band is well-positioned to probe such cosmic events. Our results constitute a thorough exploration of the cosmic string science that will be accessible to LISA. △ Less

Submitted 2 September, 2019; originally announced September 2019.

Comments: 40 pages, 12 figures, 1 appendix

The LOFAR Tied-Array All-Sky Survey (LOTAAS): Survey overview and initial pulsar discoveries

Authors: S. Sanidas, S. Cooper, C. G. Bassa, J. W. T. Hessels, V. I. Kondratiev, D. Michilli, B. W. Stappers, C. M. Tan, J. van Leeuwen, L. Cerrigone, R. A. Fallows, M. Iacobelli, E. Orru, R. F. Pizzo, A. Shulevski, M. C. Toribio, S. ter Veen, P. Zucca, L. Bondonneau, J. -M. Griessmeier, A. Karastergiou, M. Kramer, C. Sobey

Abstract: We present an overview of the LOFAR Tied-Array All-Sky Survey (LOTAAS) for radio pulsars and fast transients. The survey uses the high-band antennas of the LOFAR Superterp, the dense inner part of the LOFAR core, to survey the northern sky (dec > 0 deg) at a central observing frequency of 135 MHz. A total of 219 tied-array beams (coherent summation of station signals, covering 12 square degrees),… ▽ More We present an overview of the LOFAR Tied-Array All-Sky Survey (LOTAAS) for radio pulsars and fast transients. The survey uses the high-band antennas of the LOFAR Superterp, the dense inner part of the LOFAR core, to survey the northern sky (dec > 0 deg) at a central observing frequency of 135 MHz. A total of 219 tied-array beams (coherent summation of station signals, covering 12 square degrees), as well as three incoherent beams (covering 67 square degrees) are formed in each survey pointing. For each ofthe 222 beams, total intensity is recorded at 491.52 us time resolution. Each observation integrates for 1 hr and covers 2592 channels from 119 to 151 MHz. This instrumental setup allows LOTAAS to reach a detection threshold of 1 to 5 mJy for periodic emission. Thus far, the LOTAAS survey has resulted in the discovery of 73 radio pulsars. Among these are two mildly recycled binary millisecond pulsars (P = 13 and 33 ms), as well as the slowest-spinning radio pulsar currently known (P = 23.5 s). The survey has thus far detected 311 known pulsars, with spin periods ranging from 4 ms to 5.0 s and dispersion measures from 3.0 to 217 pc/cc. Known pulsars are detected at flux densities consistent with literature values. We find that the LOTAAS pulsar discoveries have, on average, longer spin periods than the known pulsar population. This may reflect different selection biases between LOTAAS and previous surveys, though it is also possible that slower-spinning pulsars preferentially have steeper radio spectra. LOTAAS is the deepest all-sky pulsar survey using a digital aperture array; we discuss some of the lessons learned that can inform the approach for similar surveys using future radio telescopes such as the Square Kilometre Array. △ Less

Submitted 13 May, 2019; originally announced May 2019.

Comments: 22 pages, 9 figures, 3 tables. Accepted for publication in A&A

A detailed study of giant pulses from PSR B1937+21 using the Large European Array for Pulsars

Authors: J. W. McKee, B. W. Stappers, C. G. Bassa, S. Chen, I. Cognard, M. Gaikwad, G. H. Janssen, R. Karuppusamy, M. Kramer, K. J. Lee, K. Liu, D. Perrodin, S. A. Sanidas, R. Smits, L. Wang, W. W. Zhu

Abstract: We have studied 4265 giant pulses (GPs) from the millisecond pulsar B1937+21; the largest-ever sample gathered for this pulsar, in observations made with the Large European Array for Pulsars. The pulse energy distribution of GPs associated with the interpulse are well-described by a power law, with index $α= -3.99 \pm 0.04$, while those associated with the main pulse are best-described by a broken… ▽ More We have studied 4265 giant pulses (GPs) from the millisecond pulsar B1937+21; the largest-ever sample gathered for this pulsar, in observations made with the Large European Array for Pulsars. The pulse energy distribution of GPs associated with the interpulse are well-described by a power law, with index $α= -3.99 \pm 0.04$, while those associated with the main pulse are best-described by a broken power law, with the break occurring at $\sim7$ Jy $μ$s, with power law indices $α_{\text{low}} = -3.48 \pm 0.04$ and $α_{\text{high}} = -2.10 \pm 0.09$. The modulation indices of the GP emission are measured, which are found to vary by $\sim0.5$ at pulse phases close to the centre of the GP phase distributions. We find the frequency-resolved structure of GPs to vary significantly, and in a manner that cannot be attributed to the interstellar medium influence on the observed pulses. We examine the distribution of polarisation fractions of the GPs and find no correlation between GP emission phase and fractional polarisation. We use the GPs to time PSR B1937+21 and although the achievable time of arrival precision of the GPs is approximately a factor of two greater than that of the average pulse profile, there is a negligible difference in the precision of the overall timing solution when using the GPs. △ Less

Submitted 7 November, 2018; originally announced November 2018.

Comments: 20 pages, 21 figures, 3 tables, accepted for publication in MNRAS

Studying the solar system with the International Pulsar Timing Array

Authors: R. N. Caballero, Y. J. Guo, K. J. Lee, P. Lazarus, D. J. Champion, G. Desvignes, M. Kramer, K. Plant, Z. Arzoumanian, M. Bailes, C. G. Bassa, N. D. R. Bhat, A. Brazier, M. Burgay, S. Burke-Spolaor, S. J. Chamberlin, S. Chatterjee, I. Cognard, J. M. Cordes, S. Dai, P. Demorest, T. Dolch, R. D. Ferdman, E. Fonseca, J. R. Gair , et al. (55 additional authors not shown)

Abstract: Pulsar-timing analyses are sensitive to errors in the solar-system ephemerides (SSEs) that timing models utilise to estimate the location of the solar-system barycentre, the quasi-inertial reference frame to which all recorded pulse times-of-arrival are referred. Any error in the SSE will affect all pulsars, therefore pulsar timing arrays (PTAs) are a suitable tool to search for such errors and im… ▽ More Pulsar-timing analyses are sensitive to errors in the solar-system ephemerides (SSEs) that timing models utilise to estimate the location of the solar-system barycentre, the quasi-inertial reference frame to which all recorded pulse times-of-arrival are referred. Any error in the SSE will affect all pulsars, therefore pulsar timing arrays (PTAs) are a suitable tool to search for such errors and impose independent constraints on relevant physical parameters. We employ the first data release of the International Pulsar Timing Array to constrain the masses of the planet-moons systems and to search for possible unmodelled objects (UMOs) in the solar system. We employ ten SSEs from two independent research groups, derive and compare mass constraints of planetary systems, and derive the first PTA mass constraints on asteroid-belt objects. Constraints on planetary-system masses have been improved by factors of up to 20 from the previous relevant study using the same assumptions, with the mass of the Jovian system measured at 9.5479189(3)$\times10^{-4}$ $M_{\odot}$. The mass of the dwarf planet Ceres is measured at 4.7(4)$\times10^{-10}$ $M_{\odot}$. We also present the first sensitivity curves using real data that place generic limits on the masses of UMOs, which can also be used as upper limits on the mass of putative exotic objects. For example, upper limits on dark-matter clumps are comparable to published limits using independent methods. While the constraints on planetary masses derived with all employed SSEs are consistent, we note and discuss differences in the associated timing residuals and UMO sensitivity curves. △ Less

Submitted 27 September, 2018; originally announced September 2018.

Comments: Accepted for publication by the Monthly Notices of the Royal Astronomical Society

LOFAR discovery of a 23.5-second radio pulsar

Authors: C. M. Tan, C. G. Bassa, S. Cooper, T. J. Dijkema, P. Esposito, J. W. T. Hessels, V. I. Kondratiev, M. Kramer, D. Michilli, S. Sanidas, T. W. Shimwell, B. W. Stappers, J. van Leeuwen, I. Cognard, J. -M. Grießmeier, A. Karastergiou, E. F. Keane, C. Sobey, P. Weltevrede

Abstract: We present the discovery of PSR J0250+5854, a radio pulsar with a spin period of 23.5 s. This is the slowest-spinning radio pulsar known. PSR J0250+5854 was discovered by the LOFAR Tied-Array All-Sky Survey (LOTAAS), an all-Northern-sky survey for pulsars and fast transients at a central observing frequency of 135 MHz. We subsequently detected pulsations from the pulsar in the interferometric imag… ▽ More We present the discovery of PSR J0250+5854, a radio pulsar with a spin period of 23.5 s. This is the slowest-spinning radio pulsar known. PSR J0250+5854 was discovered by the LOFAR Tied-Array All-Sky Survey (LOTAAS), an all-Northern-sky survey for pulsars and fast transients at a central observing frequency of 135 MHz. We subsequently detected pulsations from the pulsar in the interferometric images of the LOFAR Two-metre Sky Survey, allowing for sub-arcsecond localization. This, along with a pre-discovery detection 2 years prior, allowed us to measure the spin-period derivative to be $\dot{P}=2.7 \times 10^{-14}$ s s$^{-1}$. The observed spin period derivative of PSR J0250+5854 indicates a surface magnetic field strength, characteristic age and spin-down luminosity of $2.6 \times 10^{13}$G, $13.7$ Myr and $8.2 \times 10^{28}$ erg s$^{-1}$ respectively, for a dipolar magnetic field configuration. This also places the pulsar beyond the conventional pulsar death line, where radio emission is expected to cease. The spin period of PSR J0250+5854 is similar to those of the high-energy-emitting magnetars and X-ray dim isolated neutron stars (XDINSs). However, the pulsar was not detected by the Swift/XRT in the energy band of 0.3-10 keV, placing a bolometric luminosity limit of $1.5 \times 10^{32}$ erg s$^{-1}$ for an assumed $N_{\rm H}=1.35\times10^{21}$ cm$^{-2}$ and a temperature of 85 eV (typical of XDINSs). We discuss the implications of the discovery for models of the pulsar death line as well as the prospect of finding more similarly long-period pulsars, including the advantages provided by LOTAAS for this. △ Less

Submitted 4 September, 2018; originally announced September 2018.

Comments: 16 pages, 9 figures, Accepted for publication in the Astrophysical Journal

Journal ref: ApJ 866 (2018) 54

Single-pulse classifier for the LOFAR Tied-Array All-sky Survey

Authors: D. Michilli, J. W. T. Hessels, R. J. Lyon, C. M. Tan, C. Bassa, S. Cooper, V. I. Kondratiev, S. Sanidas, B. W. Stappers, J. van Leeuwen

Abstract: Searches for millisecond-duration, dispersed single pulses have become a standard tool used during radio pulsar surveys in the last decade. They have enabled the discovery of two new classes of sources: rotating radio transients and fast radio bursts. However, we are now in a regime where the sensitivity to single pulses in radio surveys is often limited more by the strong background of radio freq… ▽ More Searches for millisecond-duration, dispersed single pulses have become a standard tool used during radio pulsar surveys in the last decade. They have enabled the discovery of two new classes of sources: rotating radio transients and fast radio bursts. However, we are now in a regime where the sensitivity to single pulses in radio surveys is often limited more by the strong background of radio frequency interference (RFI, which can greatly increase the false-positive rate) than by the sensitivity of the telescope itself. To mitigate this problem, we introduce the Single-pulse Searcher (SpS). This is a new machine-learning classifier designed to identify astrophysical signals in a strong RFI environment, and optimized to process the large data volumes produced by the new generation of aperture array telescopes. It has been specifically developed for the LOFAR Tied-Array All-Sky Survey (LOTAAS), an ongoing survey for pulsars and fast radio transients in the northern hemisphere. During its development, SpS discovered 7 new pulsars and blindly identified ~80 known sources. The modular design of the software offers the possibility to easily adapt it to other studies with different instruments and characteristics. Indeed, SpS has already been used in other projects, e.g. to identify pulses from the fast radio burst source FRB 121102. The software development is complete and SpS is now being used to re-process all LOTAAS data collected to date. △ Less

Submitted 16 August, 2018; originally announced August 2018.

Comments: Accepted in MNRAS

Improving timing sensitivity in the microhertz frequency regime: limits from PSR J1713$+$0747 on gravitational waves produced by super-massive black-hole binaries

Authors: B. B. P. Perera, B. W. Stappers, S. Babak, M. J. Keith, J. Antoniadis, C. G. Bassa, R. N. Caballero, D. J. Champion, I. Cognard, G. Desvignes, E. Graikou, L. Guillemot, G. H. Janssen, R. Karuppusamy, M. Kramer, P. Lazarus, L. Lentati, K. Liu, A. G. Lyne, J. W. McKee, S. Oslowski, D. Perrodin, S. A. Sanidas, A. Sesana, G. Shaifullah , et al. (3 additional authors not shown)

Abstract: We search for continuous gravitational waves (CGWs) produced by individual super-massive black-hole binaries (SMBHBs) in circular orbits using high-cadence timing observations of PSR J1713$+$0747. We observe this millisecond pulsar using the telescopes in the European Pulsar Timing Array (EPTA) with an average cadence of approximately 1.6 days over the period between April 2011 and July 2015, incl… ▽ More We search for continuous gravitational waves (CGWs) produced by individual super-massive black-hole binaries (SMBHBs) in circular orbits using high-cadence timing observations of PSR J1713$+$0747. We observe this millisecond pulsar using the telescopes in the European Pulsar Timing Array (EPTA) with an average cadence of approximately 1.6 days over the period between April 2011 and July 2015, including an approximately daily average between February 2013 and April 2014. The high-cadence observations are used to improve the pulsar timing sensitivity across the GW frequency range of $0.008-5$ $μ$Hz. We use two algorithms in the analysis, including a spectral fitting method and a Bayesian approach. For an independent comparison, we also use a previously published Bayesian algorithm. We find that the Bayesian approaches provide optimal results and the timing observations of the pulsar place a 95 per cent upper limit on the sky-averaged strain amplitude of CGWs to be $\lesssim3.5\times10^{-13}$ at a reference frequency of 1 $μ$Hz. We also find a 95 per cent upper limit on the sky-averaged strain amplitude of low-frequency CGWs to be $\lesssim1.4\times10^{-14}$ at a reference frequency of 20~nHz. △ Less

Submitted 27 April, 2018; originally announced April 2018.

Comments: 12 pages, 4 figures. Accepted for publication in MNRAS

Tests of Gravitational Symmetries with Pulsar Binary J1713+0747

Authors: W. W. Zhu, G. Desvignes, N. Wex, R. N. Caballero, D. J. Champion, P. B. Demorest, J. A. Ellis, G. H. Janssen, M. Kramer, A. Krieger, L. Lentati, D. J. Nice, S. M. Ransom, I. H. Stairs, B. W. Stappers, J. P. W. Verbiest, Z. Arzoumanian, C. G. Bassa, M. Burgay, I. Cognard, K. Crowter, T. Dolch, R. D. Ferdman, E. Fonseca, M. E. Gonzalez , et al. (28 additional authors not shown)

Abstract: Symmetries play an important role in modern theories of gravity. The strong equivalence principle (SEP) constitutes a collection of gravitational symmetries which are all implemented by general relativity. Alternative theories, however, are generally expected to violate some aspects of SEP. We test three aspects of SEP using observed change rates in the orbital period and eccentricity of binary pu… ▽ More Symmetries play an important role in modern theories of gravity. The strong equivalence principle (SEP) constitutes a collection of gravitational symmetries which are all implemented by general relativity. Alternative theories, however, are generally expected to violate some aspects of SEP. We test three aspects of SEP using observed change rates in the orbital period and eccentricity of binary pulsar J1713+0747: 1. the gravitational constant's constancy as part of locational invariance of gravitation; 2. the post-Newtonian parameter $\hatα_3$ in gravitational Lorentz invariance; 3. the universality of free fall (UFF) for strongly self-gravitating bodies. Based on the pulsar timing result of the combined dataset from the North American Nanohertz Gravitational Observatory (NANOGrav) and the European Pulsar Timing Array (EPTA), we find $\dot{G}/G = (-0.1 \pm 0.9) \times 10^{-12}\,{\rm yr}^{-1}$, which is weaker than Solar system limits, but applies for strongly self-gravitating objects. Furthermore, we obtain the constraints $|Δ|< 0.002$ for the UFF test and $-3\times10^{-20} < \hatα_3 < 4\times10^{-20}$ at 95% confidence. These are the first direct UFF and $\hatα_3$ tests based on pulsar binaries, and they overcome various limitations of previous tests. △ Less

Submitted 26 February, 2018; originally announced February 2018.

Comments: 11 pages, 6 figures, and 1 table; submitted to MNRAS

Low-frequency pulse profile variation in PSR B2217+47: evidence for echoes from the interstellar medium

Authors: D. Michilli, J. W. T. Hessels, J. Y. Donner, J. -M. Grießmeier, M. Serylak, B. Shaw, B. W. Stappers, J. P. W. Verbiest, A. T. Deller, L. N. Driessen, D. R. Stinebring, L. Bondonneau, M. Geyer, M. Hoeft, A. Karastergiou, M. Kramer, S. Osłowski, M. Pilia, S. Sanidas, P. Weltevrede

Abstract: We have observed a complex and continuous change in the integrated pulse profile of PSR B2217+47, manifested as additional components trailing the main peak. These transient components are detected over 6 years at $150$ MHz using the LOw Frequency ARray (LOFAR), but they are not seen in contemporaneous Lovell observations at $1.5$ GHz. We argue that propagation effects in the ionized interstellar… ▽ More We have observed a complex and continuous change in the integrated pulse profile of PSR B2217+47, manifested as additional components trailing the main peak. These transient components are detected over 6 years at $150$ MHz using the LOw Frequency ARray (LOFAR), but they are not seen in contemporaneous Lovell observations at $1.5$ GHz. We argue that propagation effects in the ionized interstellar medium (IISM) are the most likely cause. The putative structures in the IISM causing the profile variation are roughly half-way between the pulsar and the Earth and have transverse radii $R \sim 30$ AU. We consider different models for the structures. Under the assumption of spherical symmetry, their implied average electron density is $\overline{n}_e \sim 100$ cm$^{-3}$. Since PSR B2217+47 is more than an order of magnitude brighter than the average pulsar population visible to LOFAR, similar profile variations would not have been identified in most pulsars, suggesting that subtle profile variations in low-frequency profiles might be more common than we have observed to date. Systematic studies of these variations at low frequencies can provide a new tool to investigate the proprieties of the IISM and the limits to the precision of pulsar timing. △ Less

Submitted 9 February, 2018; originally announced February 2018.

Comments: to be published in MNRAS

Targeted millisecond pulsar surveys of Fermi gamma-ray sources with LOFAR

Authors: C. G. Bassa, Z. Pleunis, J. W. T. Hessels, E. C. Ferrara, V. I. Kondratiev, S. Sanidas, A. G. Lyne, B. W. Stappers, S. M. Ransom, the Fermi Pulsar Search Consortium

Abstract: We have used LOFAR to perform targeted millisecond pulsar surveys of Fermi gamma-ray sources. Operating at a center frequency of 135 MHz, the surveys use a novel semi-coherent dedispersion approach where coherently dedispersed trials at coarsely separated dispersion measures are incoherently dedispersed at finer steps. Three millisecond pulsars have been discovered as part of these surveys. We des… ▽ More We have used LOFAR to perform targeted millisecond pulsar surveys of Fermi gamma-ray sources. Operating at a center frequency of 135 MHz, the surveys use a novel semi-coherent dedispersion approach where coherently dedispersed trials at coarsely separated dispersion measures are incoherently dedispersed at finer steps. Three millisecond pulsars have been discovered as part of these surveys. We describe the LOFAR surveys and the properties of the newly discovered pulsars. △ Less

Submitted 14 December, 2017; originally announced December 2017.

Comments: 4 pages, 1 figure and 1 bet. Accepted for publication in the proceedings of the IAU Symposium 337 - Pulsar astrophysics: The Next Fifty Years

Pulsar science with data from the Large European Array for Pulsars

Authors: J. W. McKee, C. G. Bassa, S. Chen, M. Gaikwad, G. H. Janssen, R. Karuppusamy, M. Kramer, K. J. Lee, K. Liu, D. Perrodin, S. A. Sanidas, R. Smits, B. W. Stappers, L. Wang, W. W. Zhu

Abstract: The Large European Array for Pulsars (LEAP) is a European Pulsar Timing Array project that combines the Lovell, Effelsberg, Nançay, Sardinia, and Westerbork radio telescopes into a single tied-array, and makes monthly observations of a set of millisecond pulsars (MSPs). The overview of our experiment is presented in Bassa et al. (2016). Baseband data are recorded at a central frequency of 1396 MHz… ▽ More The Large European Array for Pulsars (LEAP) is a European Pulsar Timing Array project that combines the Lovell, Effelsberg, Nançay, Sardinia, and Westerbork radio telescopes into a single tied-array, and makes monthly observations of a set of millisecond pulsars (MSPs). The overview of our experiment is presented in Bassa et al. (2016). Baseband data are recorded at a central frequency of 1396 MHz and a bandwidth of 128 MHz at each telescope, and are correlated offline on a cluster at Jodrell Bank Observatory using a purpose-built correlator, detailed in Smits et al. (2017). LEAP offers a substantial increase in sensitivity over that of the individual telescopes, and can operate in timing and imaging modes (notably in observations of the galactic centre radio magnetar; Wucknitz 2015). To date, 4 years of observations have been reduced. Here, we report on the scientific projects which have made use of LEAP data. △ Less

Submitted 6 December, 2017; originally announced December 2017.

Comments: 2 pages, 1 figure, accepted for publication in the proceedings of IAU Symposium 337 - Pulsar Astrophysics: The Next Fifty Years

LOFAR discovery of the fastest-spinning millisecond pulsar in the Galactic field

Authors: C. G. Bassa, Z. Pleunis, J. W. T. Hessels, E. C. Ferrara, R. P. Breton, N. V. Gusinskaia, V. I. Kondratiev, S. Sanidas, L. Nieder, C. J. Clark, T. Li, A. S. van Amesfoort, T. H. Burnett, F. Camilo, P. F. Michelson, S. M. Ransom, P. S. Ray, K. Wood

Abstract: We report the discovery of PSR J0952$-$0607, a 707-Hz binary millisecond pulsar which is now the fastest-spinning neutron star known in the Galactic field (i.e., outside of a globular cluster). PSR J0952$-$0607 was found using LOFAR at a central observing frequency of 135 MHz, well below the 300 MHz to 3 GHz frequencies typically used in pulsar searches. The discovery is part of an ongoing LOFAR s… ▽ More We report the discovery of PSR J0952$-$0607, a 707-Hz binary millisecond pulsar which is now the fastest-spinning neutron star known in the Galactic field (i.e., outside of a globular cluster). PSR J0952$-$0607 was found using LOFAR at a central observing frequency of 135 MHz, well below the 300 MHz to 3 GHz frequencies typically used in pulsar searches. The discovery is part of an ongoing LOFAR survey targeting unassociated Fermi Large Area Telescope $γ$-ray sources. PSR J0952$-$0607 is in a 6.42-hr orbit around a very low-mass companion ($M_\mathrm{c}\gtrsim0.02$ M$_\odot$) and we identify a strongly variable optical source, modulated at the orbital period of the pulsar, as the binary companion. The light curve of the companion varies by 1.6 mag from $r^\prime=22.2$ at maximum to $r^\prime>23.8$, indicating that it is irradiated by the pulsar wind. Swift observations place a 3-$σ$ upper limit on the $0.3-10$ keV X-ray luminosity of $L_X < 1.1 \times 10^{31}$ erg s$^{-1}$ (using the 0.97 kpc distance inferred from the dispersion measure). Though no eclipses of the radio pulsar are observed, the properties of the system classify it as a black widow binary. The radio pulsed spectrum of PSR J0952$-$0607, as determined through flux density measurements at 150 and 350 MHz, is extremely steep with $α\sim-3$ (where $S \propto ν^α$). We discuss the growing evidence that the fastest-spinning radio pulsars have exceptionally steep radio spectra, as well as the prospects for finding more sources like PSR J0952$-$0607. △ Less

Submitted 5 September, 2017; originally announced September 2017.

Comments: 9 pages, 3 figures, 1 table, published in ApJ letters

Journal ref: Bassa, C. G., et al. 2017, ApJL, 846, L20

A millisecond pulsar discovery in a survey of unidentified Fermi $γ$-ray sources with LOFAR

Authors: Z. Pleunis, C. G. Bassa, J. W. T. Hessels, V. I. Kondratiev, F. Camilo, I. Cognard, J. -M. Griessmeier, B. W. Stappers, A. S. van Amesfoort, S. Sanidas

Abstract: Using LOFAR, we have performed a very-low-frequency (115-155 MHz) radio survey for millisecond pulsars (MSPs). The survey targeted 52 unidentified Fermi $γ$-ray sources. Employing a combination of coherent and incoherent dedispersion, we have mitigated the dispersive effects of the interstellar medium while maintaining sensitivity to fast-spinning pulsars. Toward 3FGL J1553.1+5437 we have found PS… ▽ More Using LOFAR, we have performed a very-low-frequency (115-155 MHz) radio survey for millisecond pulsars (MSPs). The survey targeted 52 unidentified Fermi $γ$-ray sources. Employing a combination of coherent and incoherent dedispersion, we have mitigated the dispersive effects of the interstellar medium while maintaining sensitivity to fast-spinning pulsars. Toward 3FGL J1553.1+5437 we have found PSR J1552+5437, the first MSP to be discovered (through its pulsations) at a radio frequency $<$ 200 MHz. PSR J1552+5437 is an isolated MSP with a 2.43 ms spin period and a dispersion measure of 22.9 pc cm$^{-3}$. The pulsar has a very steep radio spectral index ($α< -$2.8 $\pm$ 0.4). We obtain a phase-connected timing solution combining the 0.74 years of radio observations with $γ$-ray photon arrival times covering 7.5 years of Fermi observations. We find that the radio and $γ$-ray pulse profiles of PSR J1552+5437 appear to be nearly aligned. The very steep spectrum of PSR J1552+5437, along with other recent discoveries, hints at a population of radio MSPs that have been missed in surveys using higher observing frequencies. Detecting such steep spectrum sources is important for mapping the population of MSPs down to the shortest spin periods, understanding their emission in comparison to slow pulsars, and quantifying the prospects for future surveys with low-frequency radio telescopes like SKA-Low and its precursors. △ Less

Submitted 5 September, 2017; originally announced September 2017.

Comments: 10 pages, 4 figures, published in ApJL

Journal ref: 2017, ApJL, 846, L19

The beamformer and correlator for the Large European Array for Pulsars

Authors: R. Smits, C. G. Bassa, G. H. Janssen, R. Karuppusamy, M. Kramer, K. J. Lee, K. Liu, J. McKee, D. Perrodin, M. Purver, S. Sanidas, B. W. Stappers, W. W. Zhu

Abstract: The Large European Array for Pulsars combines Europe's largest radio telescopes to form a tied-array telescope that provides high signal-to-noise observations of millisecond pulsars (MSPs) with the objective to increase the sensitivity of detecting low-frequency gravitational waves. As part of this endeavor we have developed a software correlator and beamformer which enables the formation of a tie… ▽ More The Large European Array for Pulsars combines Europe's largest radio telescopes to form a tied-array telescope that provides high signal-to-noise observations of millisecond pulsars (MSPs) with the objective to increase the sensitivity of detecting low-frequency gravitational waves. As part of this endeavor we have developed a software correlator and beamformer which enables the formation of a tied-array beam from the raw voltages from each of telescopes. We explain the concepts and techniques involved in the process of adding the raw voltages coherently. We further present the software processing pipeline that is specifically designed to deal with data from widely spaced, inhomogeneous radio telescopes and describe the steps involved in preparing, correlating and creating the tied-array beam. This includes polarization calibration, bandpass correction, frequency dependent phase correction, interference mitigation and pulsar gating. A link is provided where the software can be obtained. △ Less

Submitted 19 March, 2017; originally announced March 2017.

Comments: 10 pages, 6 figures, accepted for publication in Astronomy and Computing