Large European Array For Pulsars Research Articles

The second data release from the European Pulsar Timing Array

Pulsar timing arrays offer a probe of the low-frequency gravitational wave spectrum (1–100 nHz), 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 yr. 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 yr. 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.

Variable scintillation arcs of millisecond pulsars observed with the Large European Array for Pulsars

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.

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

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 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 polarized. In PSR J1022+1001, the linear polarization 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 explain 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.

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

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 two-dimensional (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, ionized interstellar medium, and Earth. We observe a clear parabolic scintillation arc, which varies in curvature throughout the year. The distribution of power in the secondary spectra is 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 H ii 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 longitude of ascending node, since the scintillation pattern is not very sensitive to the pulsar’s motion and 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.

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

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.

Measuring interstellar delays of PSR J0613−0200 over 7 yr, using the Large European Array for Pulsars

ABSTRACT Using data from the Large European Array for Pulsars, and the Effelsberg telescope, we study the scintillation parameters of the millisecond pulsar PSR J0613−0200 over a 7 yr 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 \, \mu {\rm 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 fitted by a one-dimensional scattering screen ${\sim}40{{\ \rm per\ cent}}$ 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 programmes 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.

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

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 $\alpha = -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 $\mu$s, with power law indices $\alpha_{\text{low}} = -3.48 \pm 0.04$ and $\alpha_{\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.

The Sardinia Radio Telescope

[Abridged] The Sardinia Radio Telescope (SRT) is the new 64-m dish operated by INAF (Italy). Its active surface will allow us to observe at frequencies of up to 116 GHz. At the moment, three receivers, one per focal position, have been installed and tested. The SRT was officially opened in October 2013, upon completion of its technical commissioning phase. In this paper, we provide an overview of the main science drivers for the SRT, describe the main outcomes from the scientific commissioning of the telescope, and discuss a set of observations demonstrating the SRT's scientific capabilities. One of the main objectives of scientific commissioning was the identification of deficiencies in the instrumentation and/or in the telescope sub-systems for further optimization. As a result, the overall telescope performance has been significantly improved. As part of the scientific commissioning activities, different observing modes were tested and validated, and first astronomical observations were carried out to demonstrate the science capabilities of the SRT. In addition, we developed astronomer-oriented software tools, to support future observers on-site. The astronomical validation activities were prioritized based on technical readiness and scientific impact. The highest priority was to make the SRT available for joint observations as part of European networks. As a result, the SRT started to participate (in shared-risk mode) in EVN (European VLBI Network) and LEAP (Large European Array for Pulsars) observing sessions in early 2014. The validation of single-dish operations for the suite of SRT first light receivers and backends continued in the following years, and was concluded with the first call for shared-risk/early-science observations issued at the end of 2015.

Pulsar science with data from the Large European Array for Pulsars

AbstractThe 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.

Pulsar science at the Sardinia Radio Telescope

AbstractThe Sardinia Radio Telescope (SRT) is a modern, fully-steerable 64-m dish located in San Basilio, Sardinia (Italy). It is characterized by an active surface that allows it to cover a wide range of radio frequencies (300 MHz to 100 GHz). During SRT’s commissioning phase, we installed the hardware and software needed for pulsar observations. Since then, SRT has taken part in Large European Array for Pulsars and European Pulsar Timing Array observations for the purpose of gravitational wave detection. We have installed a new S-band receiver that will allow us to search for pulsars in the Galactic Center. We also plan to combine our efforts to search for Extraterrestrial Intelligence (SETI) with the search for pulsars and Fast Radio Bursts.

The beamformer and correlator for the Large European Array for Pulsars

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.

Variability, polarimetry, and timing properties of single pulses from PSR J1713+0747 using the Large European Array for Pulsars

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 7 and 3 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 polarisation of single pulses increases with their pulse peak flux density. By mapping the probability density of linear polarisation position angle with pulse longitude, we reveal the existence of two orthogonal polarisation 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.

LEAP: the Large European Array for Pulsars

The Large European Array for Pulsars (LEAP) is an experiment that harvests the collective power of Europe's largest radio telescopes in order to increase the sensitivity of high-precision pulsar timing. As part of the ongoing effort of the European Pulsar Timing Array (EPTA), LEAP aims to go beyond the sensitivity threshold needed to deliver the first direct detection of gravitational waves. The five telescopes presently included in LEAP are: the Effelsberg telescope, the Lovell telescope at Jodrell Bank, the Nan\c cay radio telescope, the Sardinia Radio Telescope and the Westerbork Synthesis Radio Telescope. Dual polarization, Nyquist-sampled time-series of the incoming radio waves are recorded and processed offline to form the coherent sum, resulting in a tied-array telescope with an effective aperture equivalent to a 195-m diameter circular dish. All observations are performed using a bandwidth of 128 MHz centered at a frequency of 1396 MHz. In this paper, we present the design of the LEAP experiment, the instrumentation, the storage and transfer of data, and the processing hardware and software. In particular, we present the software pipeline that was designed to process the Nyquist-sampled time-series, measure the phase and time delays between each individual telescope and a reference telescope and apply these delays to form the tied-array coherent addition. The pipeline includes polarization calibration and interference mitigation. We also present the first results from LEAP and demonstrate the resulting increase in sensitivity, which leads to an improvement in the pulse arrival times.

The European Pulsar Timing Array and the Large European Array for Pulsars

The European Pulsar Timing Array (EPTA) and the Large European Array for Pulsars (LEAP) play crucial roles in the global effort to detect gravitational waves (GWs) with a Pulsar Timing Array (PTA) experiment. While the EPTA uses five of the world’s largest cm-radio telescopes, LEAP harvests their combined power to synthesize a 194 m equivalent dish to provide high-precision PTA data for most of the sky. The EPTA has already produced a large variety of results, including astrophysical studies of individual pulsars, tests of theories of gravity, stringent limits on a GW background produced by super-massive binary black holes or the vibration of cosmic strings. It has also undertaken the development of new analysis methods and techniques, and studies of the astrophysics and population of expected GW background sources. This review gives an overview of the EPTA and LEAP set-ups and corresponding activities.

Constraints on cosmic string tension imposed by the limit on the stochastic gravitational wave background from the European Pulsar Timing Array

We investigate the constraints that can be placed on the cosmic string tension by using the current Pulsar Timing Array limits on the stochastic gravitational wave background (SGWB). We have developed a code to compute the spectrum of gravitational waves (GWs) based on the widely accepted one-scale model. In its simplest form the one-scale model allows one to vary: (i) the string tension, G\mu/c^2; (ii) the size of cosmic string loops relative to the horizon at birth, \alpha; (iii) the spectral index of the emission spectrum, q; (iv) the cut-off in the emission spectrum, n_*; and (v) the intercommutation probability, p. The amplitude and slope of the spectrum in the nHz frequency range is very sensitive to these unknown parameters. We have also investigated the impact of more complicated scenarios with multiple initial loop sizes, in particular the 2-\alpha models proposed in the literature and a log-normal distribution for \alpha. We have computed the constraint on G\mu/c^2 due to the limit on a SGWB imposed by data from the European Pulsar Timing Array. Taking into account all the possible uncertainties in the parameters we find a conservative upper limit of G\mu/c^2<5.3x 10^{-7} which typically occurs when the loop production scale is close to the gravitational backreaction scale, \alpha\approx\Gamma G\mu/c^2. Stronger limits are possible for specific values of the parameters which typically correspond to the extremal cases \alpha\ll \Gamma G\mu/c^2 and \alpha\gg \Gamma G\mu/c^2. This limit is less stringent than the previously published limits which are based on cusp emission, an approach which does not necessarily model all the possible uncertainties. We discuss the prospects for lowering this limit by two orders of magnitude, or even a detection of the SGWB, in the very near future in the context of the Large European Array for Pulsars and the Square Kilometre Array.

The European Pulsar Timing Array: current efforts and a LEAP toward the future

The European Pulsar Timing Array (EPTA) is a multi-institutional, multi-telescope collaboration, with the goal of using high-precision pulsar timing to directly detect gravitational waves. In this paper we discuss the EPTA member telescopes, current achieved timing precision and near-future goals. We report a preliminary upper limit to the amplitude of a gravitational wave background. We also discuss the Large European Array for Pulsars, in which the five major European telescopes involved in pulsar timing will be combined to provide a coherent array that will give similar sensitivity to the Arecibo radio telescope, and larger sky coverage.