ABSTRACT We discovered fine-scale structure within the scattering disk of PSR B0329+54 in observations with the RadioAstron ground-space radio interferometer. Here we describe this phenomenon, ...characterize it with averages and correlation functions, and interpret it as the result of decorrelation of the impulse-response function of interstellar scattering between the widely separated antennas. This instrument included the 10 m Space Radio Telescope, the 110 m Green Bank Telescope, the 14 × 25 m Westerbork Synthesis Radio Telescope, and the 64 m Kalyazin Radio Telescope. The observations were performed at 324 MHz on baselines of up to 235,000 km in 2012 November and 2014 January. In the delay domain, on long baselines the interferometric visibility consists of many discrete spikes within a limited range of delays. On short baselines it consists of a sharp spike surrounded by lower spikes. The average envelope of correlations of the visibility function shows two exponential scales, with characteristic delays of τ 1 = 4.1 0.3 s and τ 2 = 23 3 s , indicating the presence of two scales of scattering in the interstellar medium. These two scales are present in the pulse-broadening function. The longer scale contains 0.38 times the scattered power of the shorter one. We suggest that the longer tail arises from highly scattered paths, possibly from anisotropic scattering or from substructure at large angles.
We present temporal scattering measurements of single pulses and average profiles of PSR J1745-2900, a magnetar recently discovered only 3 arcsec away from Sagittarius A* (Sgr A*), from 1.2 to 18.95 ...GHz using the Effelsberg 100 m Radio Telescope, the Nancay Decimetric Radio Telescope, and the Jodrell Bank Lovell Telescope. Single pulse analysis shows that the integrated pulse profile above 2 GHz is dominated by pulse jitter, while below 2 GHz the pulse profile shape is dominated by scattering. This is the first object in the Galactic center (GC) with both pulse broadening and angular broadening measurements. We measure a pulse broadening time scale at 1 GHz of tau sub(1GHz) = 1.3 + or - 0.2 and pulse broadening spectral index of alpha = -3.8 + or - 0.2, which is several orders of magnitude lower than predicted by the NE2001 model (Cordes & Lazio 2002). If this scattering time scale is representative of the GC as a whole, then previous surveys should have detected many pulsars. The lack of detections implies either our understanding of scattering in the GC is incomplete or there are fewer pulsars in the GC than previously predicted. Given that magnetars are a rare class of radio pulsar, there are likely many canonical and millisecond pulsars in the GC, and not surprisingly, scattering in the GC is spatially complex.
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 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 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: the signal from one telescope is zapped with time and frequency resolutions of $6.25\, \mathrm{ms}$ and $0.16\, \mathrm{MHz}$, removing interference, along with 0.27 per cent of ‘good’ data, giving an uncertainty of $0.25\, \mathrm{\mu \mathrm{ 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\, \mathrm{\mu \mathrm{ s}}$ in TOA uncertainty.
ABSTRACT
We present the first large sample of scintillation arcs in millisecond pulsars (MSPs), 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 time-scales over a 10 yr time span. 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 multiscreen 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. Unmodelled 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 very long baseline interferometry (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.
We report on 22 yr of radio timing observations of the millisecond pulsar J1024−0719 by the telescopes participating in the European Pulsar Timing Array (EPTA). These observations reveal a ...significant second derivative of the pulsar spin frequency and confirm the discrepancy between the parallax and Shklovskii distances that has been reported earlier. We also present optical astrometry, photometry and spectroscopy of 2MASS J10243869−0719190. We find that it is a low-metallicity main-sequence star (K7V spectral type, M/H = −1.0, T
eff = 4050 ± 50 K) and that its position, proper motion and distance are consistent with those of PSR J1024−0719. We conclude that PSR J1024−0719 and 2MASS J10243869−0719190 form a common proper motion pair and are gravitationally bound. The gravitational interaction between the main-sequence star and the pulsar accounts for the spin frequency derivatives, which in turn resolves the distance discrepancy. Our observations suggest that the pulsar and main-sequence star are in an extremely wide (P
b > 200 yr) orbit. Combining the radial velocity of the companion and proper motion of the pulsar, we find that the binary system has a high spatial velocity of 384 ± 45 km s−1 with respect to the local standard of rest and has a Galactic orbit consistent with halo objects. Since the observed main-sequence companion star cannot have recycled the pulsar to millisecond spin periods, an exotic formation scenario is required. We demonstrate that this extremely wide-orbit binary could have evolved from a triple system that underwent an asymmetric supernova explosion, though find that significant fine-tuning during the explosion is required. Finally, we discuss the implications of the long period orbit on the timing stability of PSR J1024−0719 in light of its inclusion in pulsar timing arrays.
In 2015, Torne et al. reported detections of the magnetar SGR J1745-2900 up to 225 GHz (1.33 mm), which was the highest radio frequency detection of pulsar emission at that time. In this work, we ...present the results of new observations of the same magnetar with detections up to 291 GHz (1.03 mm), together with evidence of linear polarization in its millimetre emission. SGR J1745-2900 continues to show variability and is, on average, a factor ~4 brighter in the millimetre band than in our observations of 2014 July. The new measured spectrum is slightly inverted, with ... However, the spectrum does not seem to be well described by a single power law, which might be due to the intrinsic variability of the source, or perhaps a turn-up somewhere between 8.35 and 87 GHz. These results may help us to improve our still incomplete model of pulsar emission and, in addition, they further support the search for and study of pulsars located at the Galactic Centre using millimetre wavelengths. (ProQuest: ... denotes formulae/symbols omitted.)
Abstract
We search for continuous gravitational waves (CGWs) produced by individual supermassive black hole binaries 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 with an average cadence of approximately 1.6 d over the period between 2011 April and 2015 July, including an approximately daily average between 2013 February and 2014 April. The high-cadence observations are used to improve the pulsar timing sensitivity across the gravitational wave frequency range of 0.008−5$\mu$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 ≲3.5 × 10−13 at a reference frequency of 1 $\mu$Hz. We also find a 95 per cent upper limit on the sky-averaged strain amplitude of low-frequency CGWs to be ≲1.4 × 10−14 at a reference frequency of 20 nHz.
Single pulses preserve information about the pulsar radio emission and propagation in the pulsar magnetosphere, and understanding the behaviour of their variability is essential for estimating the ...fundamental limit on the achievable pulsar timing precision. Here we report the findings of our analysis of single pulses from PSR J1713+0747 with data collected by the Large European Array for Pulsars (LEAP). We present statistical studies of the pulse properties that include distributions of their energy, phase and width. Two modes of systematic sub-pulse drifting have been detected, with a periodicity of seven and three pulse periods. The two modes appear at different ranges of pulse longitude but overlap under the main peak of the integrated profile. No evidence for pulse micro-structure is seen with a time resolution down to 140 ns. In addition, we show that the fractional polarization of single pulses increases with their pulse peak flux density. By mapping the probability density of linear polarization position angle with pulse longitude, we reveal the existence of two orthogonal polarization modes. Finally, we find that the resulting phase jitter of integrated profiles caused by single pulse variability can be described by a Gaussian probability distribution only when at least 100 pulses are used for integration. Pulses of different flux densities and widths contribute approximately equally to the phase jitter, and no improvement on timing precision is achieved by using a sub-set of pulses with a specific range of flux density or width.
Context.
Extracting precise pulse times of arrival (TOAs) and their uncertainties is the first and most fundamental step in high-precision pulsar timing. In the classical method, TOAs are derived ...from total intensity pulse profiles of pulsars via cross-correlation with an idealised 1D template of that profile. While a number of results have been presented in the literature that rely on the ever increasing sensitivity of these pulsar timing experiments, there is no consensus on the most reliable methods for creating TOAs, and, more importantly, on the associated TOA uncertainties for each scheme.
Aims.
We present a comprehensive comparison of TOA determination practices. We focus on creating timing templates, TOA determination methods, and the most useful TOA bandwidth. The aim is to present a possible approach towards TOA optimisation, the (partial) identification of an optimal TOA-creation scheme, and the demonstration of optimisation differences between pulsars and data sets.
Methods.
We compared the values of data-derived template profiles with analytic profiles and evaluated the three most commonly used template-matching methods. Finally, we studied the relation between timing precision and TOA bandwidth to identify any potential breaks in this relation. As a practical demonstration, we applied our selected methods to European Pulsar Timing Array data on three test pulsars, PSRs J0218+4232, J1713+0747, and J2145−0750.
Results.
Our demonstration shows that data-derived and smoothed templates are typically preferred to some more commonly applied alternatives. The template-matching method called Fourier domain with Markov chain Monte Carlo is generally superior to or competitive with other methods. While the optimal TOA bandwidth is strongly dependent on pulsar brightness, telescope sensitivity, and scintillation properties, some significant frequency averaging seems required for the data we investigated.