Gravitational waves are expected to be radiated by supermassive black hole binaries formed during galaxy mergers. A stochastic superposition of gravitational waves from all such binary systems would ...modulate the arrival times of pulses from radio pulsars. Using observations of millisecond pulsars obtained with the Parkes radio telescope, we constrained the characteristic amplitude of this background, Ac,yr, to be <1.0 × 10–15 with 95% confidence. This limit excludes predicted ranges for Ac,yr from current models with 91 to 99.7% probability. We conclude that binary evolution is either stalled or dramatically accelerated by galactic-center environments and that higher-cadence and shorter-wavelength observations would be more sensitive to gravitational waves.
We present timing models for 20 millisecond pulsars in the Parkes Pulsar Timing Array. The precision of the parameter measurements in these models has been improved over earlier results by using ...longer data sets and modelling the non-stationary noise. We describe a new noise modelling procedure and demonstrate its effectiveness using simulated data. Our methodology includes the addition of annual dispersion measure (DM) variations to the timing models of some pulsars. We present the first significant parallax measurements for PSRs J1024−0719, J1045−4509, J1600−3053, J1603−7202, and J1730−2304, as well as the first significant measurements of some post-Keplerian orbital parameters in six binary pulsars, caused by kinematic effects. Improved Shapiro delay measurements have resulted in much improved pulsar mass measurements, particularly for PSRs J0437−4715 and J1909−3744 with M
p = 1.44 ± 0.07 and 1.47 ± 0.03 M⊙, respectively. The improved orbital period-derivative measurement for PSR J0437−4715 results in a derived distance measurement at the 0.16 per cent level of precision, D = 156.79 ± 0.25 pc, one of the most fractionally precise distance measurements of any star to date.
We present results of an all-sky search in the Parkes Pulsar Timing Array (PPTA) Data Release 1 data set for continuous gravitational waves (GWs) in the frequency range from 5 × 10−9 to 2 × 10−7 Hz. ...Such signals could be produced by individual supermassive binary black hole systems in the early stage of coalescence. We phase up the pulsar timing array data set to form, for each position on the sky, two data streams that correspond to the two GW polarizations and then carry out an optimal search for GW signals on these data streams. Since no statistically significant GWs were detected, we place upper limits on the intrinsic GW strain amplitude h
0 for a range of GW frequencies. For example, at 10−8 Hz our analysis has excluded with 95 per cent confidence the presence of signals with h
0 ≥ 1.7 × 10−14. Our new limits are about a factor of 4 more stringent than those of Yardley et al. based on an earlier PPTA data set and a factor of 2 better than those reported in the recent Arzoumanian et al. paper. We also present PPTA directional sensitivity curves and find that for the most sensitive region on the sky, the current data set is sensitive to GWs from circular supermassive binary black holes with chirp masses of 109 M⊙ out to a luminosity distance of about 100 Mpc. Finally, we set an upper limit of 4 × 10−3 Mpc−3 Gyr−1 at 95 per cent confidence on the coalescence rate of nearby (z ≲ 0.1) supermassive binary black holes in circular orbits with chirp masses of 1010 M⊙.
High-sensitivity radio-frequency observations of millisecond pulsars usually show stochastic, broad-band, pulse-shape variations intrinsic to the pulsar emission process. These variations induce ...jitter noise in pulsar timing observations; understanding the properties of this noise is of particular importance for the effort to detect gravitational waves with pulsar timing arrays. We assess the short-term profile and timing stability of 22 millisecond pulsars that are part of the Parkes Pulsar Timing Array sample by examining intraobservation arrival time variability and single-pulse phenomenology. In 7 of the 22 pulsars, in the band centred at approximately 1400 MHz, we find that the brightest observations are limited by intrinsic jitter. We find consistent results, either detections or upper limits, for jitter noise in other frequency bands. PSR J1909−3744 shows the lowest levels of jitter noise, which we estimate to contribute ∼10 ns root mean square error to the arrival times for hour-duration observations. Larger levels of jitter noise are found in pulsars with wider pulses and distributions of pulse intensities. The jitter noise in PSR J0437−4715 decorrelates over a bandwidth of ∼2 GHz. We show that the uncertainties associated with timing pulsar models can be improved by including physically motivated jitter uncertainties. Pulse-shape variations will limit the timing precision at future, more sensitive, telescopes; it is imperative to account for this noise when designing instrumentation and timing campaigns for these facilities.
Anisotropic bursts of gravitational radiation produced by events such as supermassive black hole mergers leave permanent imprints on space. Such gravitational wave 'memory' (GWM) signals are, in ...principle, detectable through pulsar timing as sudden changes in the apparent pulse frequency of a pulsar. If an array of pulsars is monitored as a GWM signal passes over the Earth, the pulsars would simultaneously appear to change pulse frequency by an amount that varies with their sky position in a quadrupolar fashion. Here, we describe a search algorithm for such events and apply the algorithm to approximately six years of data from the Parkes Pulsar Timing Array. We find no GWM events and set an upper bound on the rate for events which could have been detected. We show, using simple models of black hole coalescence rates, that this non-detection is not unexpected.
We have constructed a new time-scale, 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 time-scale. The two analyses agree within the estimated uncertainties and both agree with TT(BIPM17), a post-corrected time-scale 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 time-scales are unlikely to contribute to the stability of the best time-scales over the next decade, but they will remain a valuable independent check on atomic time-scales. We also find that the stability of the pulsar-based time-scale 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.
Pulsar-timing analyses are sensitive to errors in the Solar-system ephemerides (SSEs) that timing models utilize 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 10 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) × 10^−4 M_⊙. The mass of the dwarf planet Ceres is measured at 4.7(4) × 10^−10 M_⊙. 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.
Despite the widely known fact that mantle flow in and around subduction zones produces the development of considerable seismic anisotropy, most P‐wave tomography efforts still rely on the assumption ...of isotropy. In this study, we explore the potential effects of erroneous assumption on tomographic images and explore an alternative approach. We conduct a series of synthetic tomography tests based on a geodynamic simulation of subduction and rollback. The simulation results provide a self‐consistent distribution of isotropic (thermal) anomalies and seismic anisotropy which we use to calculate synthetic delay times for a number of realistic and hypothetical event distributions. We find that anisotropy‐induced artifacts are abundant and significant for teleseismic, local and mixed event distributions. The occurrence of artifacts is not reduced, and indeed can be exacerbated, by increasing richness in ray‐path azimuths and incidence angles. The artifacts that we observe are, in all cases, important enough to significantly impact the interpretation of the images. We test an approach based on prescribing the anisotropy field as an a priori constraint and find that even coarse approximations to the true anisotropy field produce useful results. Using approximate anisotropy, fields can result in reduced RMS misfit to the travel time delays and reduced abundance and severity of imaging artifacts. We propose that the use of anisotropy fields derived from geodynamic modeling and constrained by seismic observables may constitute a viable alternative to isotropic tomography that does not require the inversion for anisotropy parameters in each node of the model.
Key Points:
Isotropic tomography images of subduction zones are likely contaminated by artifacts
This problem is not solved by increasing the azimuth and incidence angle diversity of the ray coverage
Prescribing an approximate anisotropy field as an a priori constraint may be a useful approach
We present high signal-to-noise ratio, multifrequency polarization pulse profiles for 24 millisecond pulsars that are being observed as part of the Parkes Pulsar Timing Array project. The pulsars are ...observed in three bands, centred close to 730, 1400 and 3100 MHz, using a dual-band 10 cm/50 cm receiver and the central beam of the 20-cm multibeam receiver. Observations spanning approximately six years have been carefully calibrated and summed to produce high S/N profiles. This allows us to study the individual profile components and in particular how they evolve with frequency. We also identify previously undetected profile features. For many pulsars we show that pulsed emission extends across almost the entire pulse profile. The pulse component widths and component separations follow a complex evolution with frequency; in some cases these parameters increase and in other cases they decrease with increasing frequency. The evolution with frequency of the polarization properties of the profile is also non-trivial. We provide evidence that the pre- and post-cursors generally have higher fractional linear polarization than the main pulse. We have obtained the spectral index and rotation measure for each pulsar by fitting across all three observing bands. For the majority of pulsars, the spectra follow a single power-law and the position angles follow a ... relation, as expected. However, clear deviations are seen for some pulsars. We also present phase-resolved measurements of the spectral index, fractional linear polarization and rotation measure. All these properties are shown to vary systematically over the pulse profile. (ProQuest: ... denotes formulae/symbols omitted.)