In recent years, instrumentation enabling pulsar observations with unprecedentedly high fractional bandwidth has been under development which can be used to substantially improve the precision of ...pulsar timing experiments. The traditional template-matching method used to calculate pulse times of arrival (ToAs) may not function effectively on these broad-band data due to a variety of effects such as diffractive scintillation in the interstellar medium, profile variation as a function of frequency, dispersion measure (DM) evolution, and so forth. In this paper, we describe the channelized discrete Fourier transform method that can greatly mitigate the influence of the aforementioned effects when measuring ToAs from broad-band timing data. The method is tested on simulated data, and its potential in improving timing precision is shown. We further apply the method to PSR J1909−3744 data collected at the Nançay Radio Telescope with the Nançay Ultimate Pulsar Processing Instrument. We demonstrate removal of systematics due to the scintillation effect as well as improvement on ToA measurement uncertainties. Our method also determines temporal variations in DM, which are consistent with multichannel timing approaches used earlier.
We report on the high-precision timing of 42 radio millisecond pulsars (MSPs) observed by the European Pulsar Timing Array (EPTA). This EPTA Data Release 1.0 extends up to mid-2014 and baselines ...range from 7–18 yr. It forms the basis for the stochastic gravitational-wave background, anisotropic background, and continuous-wave limits recently presented by the EPTA elsewhere. The Bayesian timing analysis performed with temponest
yields the detection of several new parameters: seven parallaxes, nine proper motions and, in the case of six binary pulsars, an apparent change of the semimajor axis. We find the NE2001 Galactic electron density model to be a better match to our parallax distances (after correction from the Lutz–Kelker bias) than the M2 and M3 models by Schnitzeler. However, we measure an average uncertainty of 80 per cent (fractional) for NE2001, three times larger than what is typically assumed in the literature. We revisit the transverse velocity distribution for a set of 19 isolated and 57 binary MSPs and find no statistical difference between these two populations. We detect Shapiro delay in the timing residuals of PSRs J1600−3053 and J1918−0642, implying pulsar and companion masses
$m_{\rm p}=1.22_{-0.35}^{+0.5}\ {\rm M}_{{\odot }}$
,
$m_{\rm c} = 0.21_{-0.04}^{+0.06}\ {\rm M}_{{\odot } }$
and
$m_{\rm p}=1.25_{-0.4}^{+0.6}\ {\rm M}_{{\odot }}$
,
$m_{\rm c} = 0.23_{-0.05}^{+0.07}\ {\rm M}_{{\odot } }$
, respectively. Finally, we use the measurement of the orbital period derivative to set a stringent constraint on the distance to PSRs J1012+5307 and J1909−3744, and set limits on the longitude of ascending node through the search of the annual-orbital parallax for PSRs J1600−3053 and J1909−3744.
PSR J1903+0327, a millisecond pulsar in an eccentric (e= 0.44) 95-d orbit with an ∼1 M⊙ companion poses a challenge to our understanding of stellar evolution in binary and multiple-star systems. Here ...we describe optical and radio observations which rule out most of the scenarios proposed to explain formation of this system. Radio timing measurements of three post-Keplerian effects yield the most precise measurement of the mass of a millisecond pulsar to date: 1.667 ± 0.021 solar masses (99.7 per cent confidence limit). This rules out some equations of state for superdense matter; furthermore, it is consistent with the spin-up of the pulsar by mass accretion, as suggested by its short spin period and low magnetic field. Optical spectroscopy of a proposed main-sequence counterpart shows that its orbital motion mirrors the pulsar's 95-d orbit; being therefore its binary companion. This finding rules out a previously suggested scenario which proposes that the system is presently a hierarchical triple. Conventional binary evolution scenarios predict that, after recycling a neutron star into a millisecond pulsar, the binary companion should become a white dwarf and its orbit should be nearly circular. This suggests that if PSR J1903+0327 was recycled, its present companion was not responsible for it. The optical detection also provides a measurement of the systemic radial velocity of the binary; this and the proper motion measured from pulsar timing allow the determination of the systemic 3D velocity in the Galaxy. We find that the system is always within 270 pc of the plane of the Galaxy, but always more than 3 kpc away from the Galactic Centre. Thus an exchange interaction in a dense stellar environment (like a globular cluster or the Galactic Centre) is not likely to be the origin of this system. We suggest that after the supernova that formed it, the neutron star was in a tight orbit with a main-sequence star and the present companion was a tertiary farther out. The neutron star then accreted matter from its evolving inner companion, forming a millisecond pulsar. The inner companion then disappeared, either due to a chaotic three-body interaction with the outer star (caused by the expansion of the inner orbit that necessarily results from mass transfer), or in the case of a very compact inner system, due to ablation/accretion by the newly formed millisecond pulsar. We discuss in detail the possible evolution of such a system before the supernova.
We present new limits on an isotropic stochastic gravitational-wave background (GWB) using a six pulsar data set spanning 18 yr of observations from the 2015 European Pulsar Timing Array data ...release. Performing a Bayesian analysis, we fit simultaneously for the intrinsic noise parameters for each pulsar, along with common correlated signals including clock, and Solar system ephemeris errors, obtaining a robust 95 per cent upper limit on the dimensionless strain amplitude A of the background of A < 3.0 × 10−15 at a reference frequency of 1 yr−1 and a spectral index of 13/3, corresponding to a background from inspiralling supermassive black hole binaries, constraining the GW energy density to Ωgw(f)h
2 < 1.1 × 10−9 at 2.8 nHz. We also present limits on the correlated power spectrum at a series of discrete frequencies, and show that our sensitivity to a fiducial isotropic GWB is highest at a frequency of ∼5 × 10−9 Hz. Finally, we discuss the implications of our analysis for the astrophysics of supermassive black hole binaries, and present 95 per cent upper limits on the string tension, Gμ/c
2, characterizing a background produced by a cosmic string network for a set of possible scenarios, and for a stochastic relic GWB. For a Nambu–Goto field theory cosmic string network, we set a limit Gμ/c
2 < 1.3 × 10−7, identical to that set by the Planck Collaboration, when combining Planck and high-ℓ cosmic microwave background data from other experiments. For a stochastic relic background, we set a limit of
$\Omega ^\mathrm{relic}_\mathrm{gw}(f)h^2<1.2 \times 10^{-9}$
, a factor of 9 improvement over the most stringent limits previously set by a pulsar timing array.
The PSRIX backend is the primary pulsar timing instrument of the Effelsberg 100 m radio telescope since early 2011. This new ROACH-based system enables bandwidths up to 500 MHz to be recorded, ...significantly more than what was possible with its predecessor, the Effelsberg–Berkeley Pulsar Processor (EBPP). We review the first four years of PSRIX timing data for 33 pulsars collected as part of the monthly European Pulsar Timing Array (EPTA) observations. We describe the automated data analysis pipeline, coastguard, that we developed to reduce these observations. We also introduce toaster, the EPTA timing data base, used to store timing results, processing information and observation metadata. Using these new tools, we measure the phase-averaged flux densities at 1.4 GHz of all 33 pulsars. For seven of these pulsars, our flux density measurements are the first values ever reported. For the other 26 pulsars, we compare our flux density measurements with previously published values. By comparing PSRIX data with EBPP data, we find an improvement of ∼2–5 times in signal-to-noise ratio, which translates to an increase of ∼2–5 times in pulse time-of-arrival (TOA) precision. We show that such an improvement in TOA precision will improve the sensitivity to the stochastic gravitational wave background. Finally, we showcase the flexibility of the new PSRIX backend by observing several millisecond-period pulsars (MSPs) at 5 and 9 GHz. Motivated by our detections, we discuss the potential for complementing existing pulsar timing array data sets with MSP monitoring campaigns at these higher frequencies.
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 yr. 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, ENTERPRISE, and TEMPONEST+FORTYTWO, 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 $\gamma = 3.78_{-0.59}^{+0.69}$ (95 per cent 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.
We have searched for continuous gravitational wave (CGW) signals produced by individually resolvable, circular supermassive black hole binaries (SMBHBs) in the latest European Pulsar Timing Array ...(EPTA) data set, which consists of ultraprecise timing data on 41-ms pulsars. We develop frequentist and Bayesian detection algorithms to search both for monochromatic and frequency-evolving systems. None of the adopted algorithms show evidence for the presence of such a CGW signal, indicating that the data are best described by pulsar and radiometer noise only. Depending on the adopted detection algorithm, the 95 per cent upper limit on the sky-averaged strain amplitude lies in the range ... This limit varies by a factor of five, depending on the assumed source position and the most constraining limit is achieved towards the positions of the most sensitive pulsars in the timing array. The most robust upper limit - obtained via a full Bayesian analysis searching simultaneously over the signal and pulsar noise on the subset of ours six best pulsars -- is ... These limits, the most stringent to date at f < 10...nHz, exclude the presence of sub-centiparsec binaries with chirp mass ... out to a distance of about 25 Mpc, and with ... out to a distance of about 1Gpc (...). We show that state-of-the-art SMBHB population models predict <1 per cent probability of detecting a CGW with the current EPTA data set, consistent with the reported non-detection. We stress, however, that PTA limits on individual CGW have improved by almost an order of magnitude in the last five years. The continuing advances in pulsar timing data acquisition and analysis techniques will allow for strong astrophysical constraints on the population of nearby SMBHBs in the coming years. (ProQuest: ... denotes formulae/symbols omitted.)
A glitch in the millisecond pulsar J0613−0200 McKee, J. W; Janssen, G. H; Stappers, B. W ...
Monthly notices of the Royal Astronomical Society,
09/2016, Volume:
461, Issue:
3
Journal Article
Peer reviewed
Open access
We present evidence for a small glitch in the spin evolution of the millisecond pulsar J0613−0200, using the EPTA Data Release 1.0, combined with Jodrell Bank analogue filterbank times of arrival ...(TOAs) recorded with the Lovell telescope and Effelsberg Pulsar Observing System TOAs. A spin frequency step of 0.82(3) nHz and frequency derivative step of −1.6(39) × 10−19 Hz s−1 are measured at the epoch of MJD 50888(30). After PSR B1821−24A, this is only the second glitch ever observed in a millisecond pulsar, with a fractional size in frequency of Δν/ν = 2.5(1) × 10−12, which is several times smaller than the previous smallest glitch. PSR J0613−0200 is used in gravitational wave searches with pulsar timing arrays, and is to date only the second such pulsar to have experienced a glitch in a combined 886 pulsar-years of observations. We find that accurately modelling the glitch does not impact the timing precision for pulsar timing array applications. We estimate that for the current set of millisecond pulsars included in the International Pulsar Timing Array, there is a probability of ∼50 per cent that another glitch will be observed in a timing array pulsar within 10 years.
The analysis of high-precision timing observations of an array of ∼20 millisecond pulsars (a so-called ‘timing array’) may ultimately result in the detection of a stochastic gravitational-wave ...background. The feasibility of such a detection and the required duration of this type of experiment are determined by the achievable rms of the timing residuals and the timing stability of the pulsars involved. We present results of the first long-term, high-precision timing campaign on a large sample of millisecond pulsars used in gravitational-wave detection projects. We show that the timing residuals of most pulsars in our sample do not contain significant low-frequency noise that could limit the use of these pulsars for decade-long gravitational-wave detection efforts. For our most precisely timed pulsars, intrinsic instabilities of the pulsars or the observing system are shown to contribute to timing irregularities on a 5-year time-scale below the 100 ns level. Based on those results, realistic sensitivity curves for planned and ongoing timing array efforts are determined. We conclude that prospects for detection of a gravitational-wave background through pulsar timing array efforts within 5 years to a decade are good.
The pulsar spectral index distribution Bates, S. D; Lorimer, D. R; Verbiest, J. P. W
Monthly notices of the Royal Astronomical Society,
05/2013, Volume:
431, Issue:
2
Journal Article
Peer reviewed
Open access
The flux-density spectra of radio pulsars are known to be steep and, to first order, described by a power-law relationship of the form S
ν ∝ να, where S
ν is the flux density at some frequency ν and ...α is the spectral index. Although measurements of α have been made over the years for several hundred pulsars, a study of the intrinsic distribution of pulsar spectra has not been carried out. From the result of pulsar surveys carried out at three different radio frequencies, we use population synthesis techniques and a likelihood analysis to deduce what underlying spectral index distribution is required to replicate the results of these surveys. We find that in general the results of the surveys can be modelled by a Gaussian distribution of spectral indices with a mean of −1.4 and unit standard deviation. We also consider the impact of the so-called gigahertz-peaked spectrum pulsars proposed by Kijak et al. The fraction of peaked-spectrum sources in the population with any significant turnover at low frequencies appears to be at most 10 per cent. We demonstrate that high-frequency (>2 GHz) surveys preferentially select flatter spectrum pulsars and the converse is true for lower frequency (<1 GHz) surveys. This implies that any correlations between α and other pulsar parameters (for example age or magnetic field) need to carefully account for selection biases in pulsar surveys. We also expect that many known pulsars which have been detected at high frequencies will have shallow, or positive, spectral indices. The majority of pulsars do not have recorded flux density measurements over a wide frequency range, making it impossible to constrain their spectral shapes. We also suggest that such measurements would allow an improved description of any populations of pulsars with 'non-standard' spectra. Further refinements to this picture will soon be possible from the results of surveys with the Green Bank Telescope and LOFAR.