Gravitational waves from binary supermassive black holes missing in pulsar observations

Abstract

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 will modulate the arrival times of pulses from radio pulsars. Using observations of millisecond pulsars obtained with the Parkes radio telescope, we constrain the characteristic amplitude of this background, A c,yr , to be < 1.0×10-15 with 95% confidence. This limit excludes predicted ranges for A c,yr from current models with 91-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 result in an increased sensitivity to gravitational waves.Studies of the dynamics of stars and gas in nearby galaxies provide strong evidence for the ubiquity of supermassive (> 10 6 solar mass) black holes (SMBHs) (1). Observations of luminous quasars indicate that SMBHs are hosted by galaxies throughout the history of the universe (2) and affect global properties of the host galaxies (3). The prevailing dark energycold dark matter cosmological paradigm predicts that large galaxies are assembled through the hierarchical merging of smaller galaxies. The remnants of mergers can host gravitationally bound binary SMBHs with orbits decaying through the emission of gravitational waves (GWs) (4).Gravitational waves from binary SMBHs, with periods between ~ 0.1 and 30 yr (5), can be detected or constrained by monitoring, for years to decades, a set of rapidly rotating millisecond pulsars (MSPs) distributed throughout our galaxy. Radio emission beams from MSPs are observed as pulses that can be time-tagged with as small as 20 ns precision (6). When traveling across the pulsar-Earth line of sight, GWs induce variations in the arrival times of the pulses (7).The superposition of GWs from the binary SMBH population is a stochastic background (GWB), which is typically characterized by the strain-amplitude spectrum h c (f)=A c,yr [f/(1 yr -1 )] -2/3 , where f is the GW frequency, A c,yr is the characteristic amplitude of the GWB measured at f = 1 yr -1 , predicted to be A c,yr > 10 -15 (5,(8)(9)(10)(11)(12), and -2/3 is the predicted spectral index (5,(8)(9)(10)(11)(12). The GWB will add low-frequency perturbations to pulse arrival times. While the detection of the GWB would confirm the presence of a cosmological population of binary SMBHs, limits on its amplitude constrain models of galaxy and SMBH evolution (8).As part of the Parkes Pulsar Timing Array project to detect GWs (6), we have been monitoring 24 pulsars with the 64-m Parkes radio telescope. We have produced a new data set, using observations taken at a central wavelength of 10 cm and previously reported methods (6,8), that spans 11 yr, which is 3 yr longer than previous data sets analyzed at this wavelength. In addition to having greater sensitivity to the GWB because of the longer duration, the data set was improved by identifying and correc...

Selection bias in dynamically measured supermassive black hole samples: consequences for pulsar timing arrays

“…The expected characteristic amplitude ranges obtained in the two cases are compared in figure 2, together with current upper limits placed by the three pulsar timing arrays: EPTA (Lentati et al 2015), NANOGrav (Arzoumanian et al 2015) and PPTA (Shannon et al 2015). Including the S16 scaling relations in the model lowers the median value of the expected signal at f = 1yr −1 by a factor of three, from 1.3 × 10 −15 to 4 × 10 −16 .…”

“…PTA limits on a stochastic GWB are usually quoted in terms of A (Lentati et al 2015;Arzoumanian et al 2015;Shannon et al 2015).…”

“…In recent years, improvements in SMBH dynamical mass measurements, together with the discovery of few "overmassive" black holes in brightest cluster galaxies (McConnell & Ma 2013a) resulted in a upward revision of the SMBH-host galaxy relations (Kormendy & Ho 2013, hereinafter KH13), implying a median expected GWB signal with strain amplitude A ∼ 10 −15 at f = 1/yr (although values of A < 10 −15 are not excluded). Therefore, recent PTA non detections of the GWB at the A ∼ 10 −15 level (Shannon et al 2015) has been interpreted as being somewhat in tension with currently favoured SMBH assembly scenarios, pointing to a possible important role of SMBH binary eccentricity or environmental coupling (Arzoumanian et al 2015).…”

Production of high stellar-mass primordial black holes in trapped inflation

“…The peak is associated with the production of about 14M ⊙ PBHs. The figure also shows the current upper limits on gravitational wave background made by the pulsar timing array experiments EPTA [43], NANOGrav [44], and PPTA [45], and the projected sensitivity of SKA radio telescope [46]. The gravitational waves associated with the production of 14M ⊙ PBHs are about an order of magnitude below the current pulsar timing array sensitivity.…”

Prospects for high-precision pulsar timing with the new Effelsberg PSRIX backend

“…Desvignes et al, submitted; Manchester et al 2013;Arzoumanian et al 2015;Shannon et al 2015). It is important to note that observing campaigns will only benefit from high-frequency detections of pulsars that are sufficiently bright to be able to take advantage of the reduced ISM effects.…”