-
A Gaussian-processes approach to fitting for time-variable spherical solar wind in pulsar timing data
Authors:
Iuliana C. Niţu,
Michael J. Keith,
Caterina Tiburzi,
Marcus Brüggen,
David J. Champion,
Siyuan Chen,
Ismaël Cognard,
Gregory Desvignes,
Ralf-Jürgen Dettmar,
Jean-Mathias Grießmeier,
Lucas Guillemot,
Yanjun Guo,
Matthias Hoeft,
Huanchen Hu,
Jiwoong Jang,
Gemma H. Janssen,
Jedrzej Jawor,
Ramesh Karuppusamy,
Evan F. Keane,
Michael Kramer,
Jörn Künsemöller,
Kristen Lackeos,
Kuo Liu,
Robert A. Main,
James W. McKee
, et al. (4 additional authors not shown)
Abstract:
Propagation effects are one of the main sources of noise in high-precision pulsar timing. For pulsars below an ecliptic latitude of $5^\circ$, the ionised plasma in the solar wind can introduce dispersive delays of order 100 microseconds around solar conjunction at an observing frequency of 300 MHz. A common approach to mitigate this assumes a spherical solar wind with a time-constant amplitude. H…
▽ More
Propagation effects are one of the main sources of noise in high-precision pulsar timing. For pulsars below an ecliptic latitude of $5^\circ$, the ionised plasma in the solar wind can introduce dispersive delays of order 100 microseconds around solar conjunction at an observing frequency of 300 MHz. A common approach to mitigate this assumes a spherical solar wind with a time-constant amplitude. However, this has been shown to be insufficient to describe the solar wind. We present a linear, Gaussian-process piecewise Bayesian approach to fit a spherical solar wind of time-variable amplitude, which has been implemented in the pulsar software run_enterprise. Through simulations, we find that the current EPTA+InPTA data combination is not sensitive to such variations; however, solar wind variations will become important in the near future with the addition of new InPTA data and data collected with the low-frequency LOFAR telescope. We also compare our results for different high-precision timing datasets (EPTA+InPTA, PPTA, and LOFAR) of three millisecond pulsars (J0030$+$0451, J1022$+$1001, J2145$-$0450), and find that the solar-wind amplitudes are generally consistent for any individual pulsar, but they can vary from pulsar to pulsar. Finally, we compare our results with those of an independent method on the same LOFAR data of the three millisecond pulsars. We find that differences between the results of the two methods can be mainly attributed to the modelling of dispersion variations in the interstellar medium, rather than the solar wind modelling.
△ Less
Submitted 15 January, 2024;
originally announced January 2024.
-
The LOFAR Tied-Array All-Sky Survey: Timing of 35 radio pulsars and an overview of the properties of the LOFAR pulsar discoveries
Authors:
E. van der Wateren,
C. G. Bassa,
S. Cooper,
J. -M. Grießmeier,
B. W. Stappers,
J. W. T. Hessels,
V. I. Kondratiev,
D. Michilli,
C. M. Tan,
C. Tiburzi,
P. Weltevrede,
A. -S. Bak Nielsen,
T. D. Carozzi,
B. Ciardi,
I. Cognard,
R. -J. Dettmar,
A. Karastergiou,
M. Kramer,
J. Künsemöller,
S. Osłowski,
M. Serylak,
C. Vocks,
O. Wucknitz
Abstract:
The LOFAR Tied-Array All-Sky Survey (LOTAAS) is the most sensitive untargeted radio pulsar survey performed at low radio frequencies (119--151\,MHz) to date and has discovered 76 new radio pulsars, among which the 23.5-s pulsar J0250+5854, up until recently the slowest-spinning radio pulsar known. Here, we report on the timing solutions of 35 pulsars discovered by LOTAAS, which include a nulling p…
▽ More
The LOFAR Tied-Array All-Sky Survey (LOTAAS) is the most sensitive untargeted radio pulsar survey performed at low radio frequencies (119--151\,MHz) to date and has discovered 76 new radio pulsars, among which the 23.5-s pulsar J0250+5854, up until recently the slowest-spinning radio pulsar known. Here, we report on the timing solutions of 35 pulsars discovered by LOTAAS, which include a nulling pulsar and a mildly recycled pulsar, and thereby complete the full timing analysis of the LOTAAS pulsar discoveries. We give an overview of the findings from the full LOTAAS sample of 76 pulsars, discussing their pulse profiles, radio spectra and timing parameters. We found that the pulse profiles of some of the pulsars show profile variations in time or frequency and while some pulsars show signs of scattering, a large majority display no pulse broadening. The LOTAAS discoveries have on average steeper radio spectra and have longer spin periods ($1.4\times$) as well as lower spin-down rates ($3.1\times$) compared to the known pulsar population. We discuss the cause of these differences, and attribute them to a combination of selection effects of the LOTAAS survey as well as previous pulsar surveys, though can not rule out that older pulsars tend to have steeper radio spectra.
△ Less
Submitted 20 November, 2022;
originally announced November 2022.
-
Pulsar scintillation studies with LOFAR. I. The census
Authors:
Ziwei Wu,
Joris P. W. Verbiest,
Robert A. Main,
Jean-Mathias Grießmeier,
Yulan Liu,
Stefan Osłowski,
Krishnakumar Moochickal Ambalappat,
Ann-Sofie Bak Nielsen,
Jörn Künsemöller,
Julian Y. Donner,
Caterina Tiburzi,
Nataliya Porayko,
Maciej Serylak,
Lars Künkel,
Marcus Brüggen,
Christian Vocks
Abstract:
Context. Interstellar scintillation (ISS) of pulsar emission can be used both as a probe of the ionised interstellar medium (IISM) and cause corruptions in pulsar timing experiments. Of particular interest are so-called scintillation arcs which can be used to measure time-variable interstellar scattering delays directly, potentially allowing high-precision improvements to timing precision.
Aims.…
▽ More
Context. Interstellar scintillation (ISS) of pulsar emission can be used both as a probe of the ionised interstellar medium (IISM) and cause corruptions in pulsar timing experiments. Of particular interest are so-called scintillation arcs which can be used to measure time-variable interstellar scattering delays directly, potentially allowing high-precision improvements to timing precision.
Aims. The primary aim of this study is to carry out the first sizeable and self-consistent census of diffractive pulsar scintillation and scintillation-arc detectability at low frequencies, as a primer for larger-scale IISM studies and pulsar-timing related propagation studies with the LOw-Frequency ARray (LOFAR) High Band Antennae (HBA).
Results. In this initial set of 31 sources, 15 allow full determination of the scintillation properties; nine of these show detectable scintillation arcs at 120-180 MHz. Eight of the observed sources show unresolved scintillation; and the final eight don't display diffractive scintillation. Some correlation between scintillation detectability and pulsar brightness and dispersion measure is apparent, although no clear cut-off values can be determined. Our measurements across a large fractional bandwidth allow a meaningful test of the frequency scaling of scintillation parameters, uncorrupted by influences from refractive scintillation variations.
Conclusions. Our results indicate the powerful advantage and great potential of ISS studies at low frequencies and the complex dependence of scintillation detectability on parameters like pulsar brightness and interstellar dispersion. This work provides the first installment of a larger-scale census and longer-term monitoring of interstellar scintillation effects at low frequencies.
△ Less
Submitted 19 March, 2022;
originally announced March 2022.
-
A broadband radio study of PSR J0250+5854: the slowest-spinning radio pulsar known
Authors:
C. H. Agar,
P. Weltevrede,
L. Bondonneau,
J. -M. Grießmeier,
J. W. T. Hessels,
W. J. Huang,
A. Karastergiou,
M. J. Keith,
V. I. Kondratiev,
J. Künsemöller,
D. Li,
B. Peng,
C. Sobey,
B. W. Stappers,
C. M. Tan,
G. Theureau,
H. G. Wang,
C. M. Zhang,
B. Cecconi,
J. N. Girard,
A. Loh,
P. Zarka
Abstract:
We present radio observations of the most slowly rotating known radio pulsar PSR J0250+5854. With a 23.5 s period, it is close, or even beyond, the $P$-$\dot{P}$ diagram region thought to be occupied by active pulsars. The simultaneous observations with FAST, the Chilbolton and Effelsberg LOFAR international stations, and NenuFAR represent a five-fold increase in the spectral coverage of this obje…
▽ More
We present radio observations of the most slowly rotating known radio pulsar PSR J0250+5854. With a 23.5 s period, it is close, or even beyond, the $P$-$\dot{P}$ diagram region thought to be occupied by active pulsars. The simultaneous observations with FAST, the Chilbolton and Effelsberg LOFAR international stations, and NenuFAR represent a five-fold increase in the spectral coverage of this object, with the detections at 1250 MHz (FAST) and 57 MHz (NenuFAR) being the highest- and lowest-frequency published respectively to date. We measure a flux density of $4\pm2$ $μ$Jy at 1250 MHz and an exceptionally steep spectral index of $-3.5^{+0.2}_{-1.5}$, with a turnover below $\sim$95 MHz. In conjunction with observations of this pulsar with the GBT and the LOFAR Core, we show that the intrinsic profile width increases drastically towards higher frequencies, contrary to the predictions of conventional radius-to-frequency mapping. We examine polarimetric data from FAST and the LOFAR Core and conclude that its polar cap radio emission is produced at an absolute height of several hundreds of kilometres around 1.5 GHz, similar to other rotation-powered pulsars across the population. Its beam is significantly underfilled at lower frequencies, or it narrows because of the disappearance of conal outriders. Finally, the results for PSR J0250+5854 and other slowly spinning rotation-powered pulsars are contrasted with the radio-detected magnetars. We conclude that magnetars have intrinsically wider radio beams than the slow rotation-powered pulsars, and that consequently the latter's lower beaming fraction is what makes objects such as PSR J0250+5854 so scarce.
△ Less
Submitted 1 September, 2021;
originally announced September 2021.
-
The impact of Solar wind variability on pulsar timing
Authors:
C. Tiburzi,
G. M. Shaifullah,
C. G. Bassa,
P. Zucca,
J. P. W. Verbiest,
N. K. Porayko,
E. van der Wateren,
R. A. Fallows,
R. A. Main,
G. H. Janssen,
J. M. Anderson,
A-. S. Bak Nielsen,
J. Y. Donner,
E. F. Keane,
J. Künsemöller,
S. Osłowski,
J-. M. Grießmeier,
M. Serylak,
M. Brüggen,
B. Ciardi,
R. -J. Dettmar,
M. Hoeft,
M. Kramer,
G. Mann,
C. Vocks
Abstract:
High-precision pulsar timing requires accurate corrections for dispersive delays of radio waves, parametrized by the dispersion measure (DM), particularly if these delays are variable in time. In a previous paper we studied the Solar-wind (SW) models used in pulsar timing to mitigate the excess of DM annually induced by the SW, and found these to be insufficient for high-precision pulsar timing. H…
▽ More
High-precision pulsar timing requires accurate corrections for dispersive delays of radio waves, parametrized by the dispersion measure (DM), particularly if these delays are variable in time. In a previous paper we studied the Solar-wind (SW) models used in pulsar timing to mitigate the excess of DM annually induced by the SW, and found these to be insufficient for high-precision pulsar timing. Here we analyze additional pulsar datasets to further investigate which aspects of the SW models currently used in pulsar timing can be readily improved, and at what levels of timing precision SW mitigation is possible. Our goals are to verify: a) whether the data are better described by a spherical model of the SW with a time-variable amplitude rather than a time-invariant one as suggested in literature, b) whether a temporal trend of such a model's amplitudes can be detected. We use the pulsar-timing technique on low-frequency pulsar observations to estimate the DM and quantify how this value changes as the Earth moves around the Sun. Specifically, we monitor the DM in weekly to monthly observations of 14 pulsars taken with LOFAR across time spans of up to 6 years. We develop an informed algorithm to separate the interstellar variations in DM from those caused by the SW and demonstrate the functionality of this algorithm with extensive simulations. Assuming a spherically symmetric model for the SW density, we derive the amplitude of this model for each year of observations. We show that a spherical model with time-variable amplitude models the observations better than a spherical model with constant amplitude, but that both approaches leave significant SW induced delays uncorrected in a number of pulsars in the sample. The amplitude of the spherical model is found to be variable in time, as opposed to what has been previously suggested.
△ Less
Submitted 21 December, 2020;
originally announced December 2020.
-
Dispersion measure variability for 36 millisecond pulsars at 150MHz with LOFAR
Authors:
J. Y. Donner,
J. P. W. Verbiest,
C. Tiburzi,
S. Osłowski,
J. Künsemöller,
A. -S. Bak Nielsen,
J. -M. Grießmeier,
M. Serylak,
M. Kramer,
J. M. Anderson,
O. Wucknitz,
E. Keane,
V. Kondratiev,
C. Sobey,
J. W. McKee,
A. V. Bilous,
R. P. Breton,
M. Brüggen,
B. Ciardi,
M. Hoeft,
J. van Leeuwen,
C. Vocks
Abstract:
Radio pulses from pulsars are affected by plasma dispersion, which results in a frequency-dependent propagation delay. Variations in the magnitude of this effect lead to an additional source of red noise in pulsar timing experiments, including pulsar timing arrays that aim to detect nanohertz gravitational waves.
We aim to quantify the time-variable dispersion with much improved precision and ch…
▽ More
Radio pulses from pulsars are affected by plasma dispersion, which results in a frequency-dependent propagation delay. Variations in the magnitude of this effect lead to an additional source of red noise in pulsar timing experiments, including pulsar timing arrays that aim to detect nanohertz gravitational waves.
We aim to quantify the time-variable dispersion with much improved precision and characterise the spectrum of these variations.
We use the pulsar timing technique to obtain highly precise dispersion measure (DM) time series. Our dataset consists of observations of 36 millisecond pulsars, which were observed for up to 7.1 years with the LOFAR telescope at a centre frequency of ~150 MHz. Seventeen of these sources were observed with a weekly cadence, while the rest were observed at monthly cadence.
We achieve a median DM precision of the order of 10^-5 cm^-3 pc for a significant fraction of our sources. We detect significant variations of the DM in all pulsars with a median DM uncertainty of less than 2x10^-4 cm^-3 pc. The noise contribution to pulsar timing experiments at higher frequencies is calculated to be at a level of 0.1-10 us at 1.4 GHz over a timespan of a few years, which is in many cases larger than the typical timing precision of 1 us or better that PTAs aim for. We found no evidence for a dependence of DM on radio frequency for any of the sources in our sample.
The DM time series we obtained using LOFAR could in principle be used to correct higher-frequency data for the variations of the dispersive delay. However, there is currently the practical restriction that pulsars tend to provide either highly precise times of arrival (ToAs) at 1.4 GHz or a high DM precision at low frequencies, but not both, due to spectral properties. Combining the higher-frequency ToAs with those from LOFAR to measure the infinite-frequency ToA and DM would improve the result.
△ Less
Submitted 27 November, 2020;
originally announced November 2020.
-
On the usefulness of existing Solar-wind models for pulsar timing corrections
Authors:
C. Tiburzi,
J. P. W. Verbiest,
G. M. Shaifullah,
G. H. Janssen,
J. M. Anderson,
A. Horneffer,
J. Kuensemoeller,
S. Oslowski,
J. Y. Donner,
M. Kramer,
A. Kumari,
N. K. Porayko,
P. Zucca,
B. Ciardi,
R. -J. Dettmar,
J. -M. Griessmeier,
M. Hoeft,
M. Serylak
Abstract:
Dispersive delays due to the Solar wind introduce excess noise in high-precision pulsar timing experiments, and must be removed in order to achieve the accuracy needed to detect, e.g., low-frequency gravitational waves. In current pulsar timing experiments, this delay is usually removed by approximating the electron density distribution in the Solar wind either as spherically symmetric, or with a…
▽ More
Dispersive delays due to the Solar wind introduce excess noise in high-precision pulsar timing experiments, and must be removed in order to achieve the accuracy needed to detect, e.g., low-frequency gravitational waves. In current pulsar timing experiments, this delay is usually removed by approximating the electron density distribution in the Solar wind either as spherically symmetric, or with a two-phase model that describes the contributions from both high- and low-speed phases of the Solar wind. However, no dataset has previously been available to test the performance and limitations of these models over extended timescales and with sufficient sensitivity. Here we present the results of such a test with an optimal dataset of observations of pulsar J0034-0534, taken with the German stations of LOFAR. We conclude that the spherical approximation performs systematically better than the two-phase model at almost all angular distances, with a residual root-mean-square (rms) given by the two-phase model being up to 28% larger than the result obtained with the spherical approximation. Nevertheless, the spherical approximation remains insufficiently accurate in modelling the Solar-wind delay (especially within 20 degrees of angular distance from the Sun), as it leaves timing residuals with rms values that reach the equivalent of 0.3 microseconds at 1400 MHz. This is because a spherical model ignores the large daily variations in electron density observed in the Solar wind. In the short term, broadband observations or simultaneous observations at low frequencies are the most promising way forward to correct for Solar-wind induced delay variations.
△ Less
Submitted 8 May, 2019;
originally announced May 2019.
-
First detection of frequency-dependent, time-variable dispersion measures
Authors:
J. Y. Donner,
J. P. W. Verbiest,
C. Tiburzi,
S. Osłowski,
D. Michilli,
M. Serylak,
J. M. Anderson,
A. Horneffer,
M. Kramer,
J. -M. Grießmeier,
J. Künsemöller,
J. W. T. Hessels,
M. Hoeft,
A. Miskolczi
Abstract:
Context. High-precision pulsar-timing experiments are affected by temporal variations of the Dispersion Measure (DM), which are related to spatial variations in the interstellar electron content. Correcting for DM variations relies on the cold-plasma dispersion law which states that the dispersive delay varies with the squared inverse of the observing frequency. This may however give incorrect mea…
▽ More
Context. High-precision pulsar-timing experiments are affected by temporal variations of the Dispersion Measure (DM), which are related to spatial variations in the interstellar electron content. Correcting for DM variations relies on the cold-plasma dispersion law which states that the dispersive delay varies with the squared inverse of the observing frequency. This may however give incorrect measurements if the probed electron content (and therefore the DM) varies with observing frequency, as is predicted theoretically.
Aims. We study small-scale density variations in the ionised interstellar medium. These structures may lead to frequency-dependent DMs in pulsar signals and could inhibit the use of lower-frequency pulsar observations to correct time-variable interstellar dispersion in higher-frequency pulsar-timing data.
Methods. We used high-cadence, low-frequency observations with three stations from the German LOng-Wavelength (GLOW) consortium, which are part of the LOw Frequency ARray (LOFAR). Specifically, 3.5 years of weekly observations of PSR J2219+4754 are presented.
Results. We present the first detection of frequency-dependent DMs towards any interstellar object and a precise multi-year time-series of the time- and frequency-dependence of the measured DMs. The observed DM variability is significant and may be caused by extreme scattering events. Potential causes for frequency-dependent DMs are quantified and evaluated.
Conclusions. We conclude that frequency-dependence of DMs has been reliably detected and is caused by small-scale (up to 10s of AUs) but steep density variations in the interstellar electron content. We find that long-term trends in DM variability equally affect DMs measured at both ends of our frequency band and hence the negative impact on long-term high-precision timing projects is expected to be limited.
△ Less
Submitted 15 February, 2019; v1 submitted 11 February, 2019;
originally announced February 2019.
-
Testing the accuracy of the ionospheric Faraday rotation corrections through LOFAR observations of bright northern pulsars
Authors:
N. K. Porayko,
A. Noutsos,
C. Tiburzi,
J. P. W. Verbiest,
A. Horneffer,
J. Künsemöller,
S. Osłowski,
M. Kramer,
D. H. F. M. Schnitzeler,
J. M. Anderson,
M. Brüggen,
J. -M. Grießmeier,
M. Hoeft,
D. J. Schwarz,
M. Serylak,
O. Wucknitz
Abstract:
Faraday rotation of polarized emission from pulsars measured at radio frequencies provides a powerful tool to investigate the interstellar and interplanetary magnetic fields. However, besides being sensitive to the astrophysical media, pulsar observations in radio are affected by the highly time-variable ionosphere. In this article, the amount of ionospheric Faraday rotation has been computed by a…
▽ More
Faraday rotation of polarized emission from pulsars measured at radio frequencies provides a powerful tool to investigate the interstellar and interplanetary magnetic fields. However, besides being sensitive to the astrophysical media, pulsar observations in radio are affected by the highly time-variable ionosphere. In this article, the amount of ionospheric Faraday rotation has been computed by assuming a thin layer model. For this aim, ionospheric maps of the free electron density (based on Global Positioning System data) and semi-empirical geomagnetic models are needed. Through the data of five highly polarized pulsars observed with the individual German LOw-Frequency ARray stations, we investigate the performances of the ionospheric modelling. In addition, we estimate the parameters of the systematics and the correlated noise generated by the residual unmodelled ionospheric effects, and show the comparison of the different free-electron density maps. For the best ionospheric maps, we have found that the rotation measure corrections on one-year timescales after subtraction of diurnal periodicity are accurate to $\sim$ 0.06--0.07 rad m$^{-2}$.
△ Less
Submitted 4 December, 2018;
originally announced December 2018.
-
Multifrequency behaviour of the anomalous events of PSR J0922+0638
Authors:
G. Shaifullah,
C. Tiburzi,
S. Osłowski,
J. P. W. Verbiest,
A. Szary,
J. Künsemöller,
A. Horneffer,
J. Anderson,
M. Kramer,
D. J. Schwarz,
G. Mann,
M. Steinmetz,
C. Vocks
Abstract:
PSR J0922+0638 (B0919+06) shows unexplained anomalous variations in the on-pulse phase, where the pulse appears to episodically move to an earlier longitude for a few tens of rotations before reverting to the usual phase for approximately several hundred to more than a thousand rotations. These events, where the pulse moves in phase by up to 5$^{\circ}$, have been previously detected in observatio…
▽ More
PSR J0922+0638 (B0919+06) shows unexplained anomalous variations in the on-pulse phase, where the pulse appears to episodically move to an earlier longitude for a few tens of rotations before reverting to the usual phase for approximately several hundred to more than a thousand rotations. These events, where the pulse moves in phase by up to 5$^{\circ}$, have been previously detected in observations from $\sim$300 to 2000 MHz. We present simultaneous observations from the Effelsberg 100-m radio telescope at 1350 MHz and the Bornim (Potsdam) station of the LOw Frequency ARray at 150 MHz. Our observations present the first evidence for an absence of the anomalous phase-shifting behaviour at 150 MHz. Instead, the observed intensity at the usual pulse-phase typically decreases, often showing a pseudo-nulling feature corresponding to the times when phase shifts are observed at 1350 MHz. The presence of weak emission at the usual pulse-phase supports the theory that these shifts may result from processes similar to the 'profile-absorption' expected to operate for PSR J0814+7429 (B0809+74). A possible mechanism for this could be intrinsic variations of the emission within the pulsar's beam combined with absorption by expanding shells of electrons in the line of sight.
△ Less
Submitted 23 March, 2018;
originally announced March 2018.