An official website of the United States government
Abstract Pulsar timing arrays (PTAs) are sensitive to low-frequency gravitational waves (GWs), which induce correlated changes in millisecond pulsars’ timing residuals. PTA collaborations around the world have recently announced evidence of a nanohertz gravitational wave background (GWB), which may be produced by a population of supermassive black hole binaries (SMBHBs). The GWB is often modeled as following a power-law power spectral density (PSD); however, a GWB produced by a cosmological population of SMBHBs is expected to have a more complex power spectrum due to the discrete nature of the sources. In this paper, we investigate using at-process PSD to model the GWB, which allows us to fit for both the underlying power-law amplitude and spectral index as well as deviations from that power law, which may be produced by individual nearby binaries. We create simulated data sets based on the properties of the NANOGrav 15 yr data set, and we demonstrate that thet-process PSD can accurately recover the PSD when deviations from a power law are present. With longer timed data sets and more pulsars, we expect the sensitivity of our PTAs to improve, which will allow us to precisely measure the PSD of the GWB and study the sources producing it.
more »
« less
Disentangling Multiple Stochastic Gravitational Wave Background Sources in PTA Data Sets
Abstract With strong evidence of a common-spectrum stochastic process in the most recent data sets from the NANOGrav Collaboration, the European Pulsar Timing Array (PTA), Parkes PTA, and the International PTA, it is crucial to assess the effects of the several astrophysical and cosmological sources that could contribute to the stochastic gravitational wave background (GWB). Using the same data set creation and injection techniques as in Pol et al., we assess the separability of multiple GWBs by creating single and multiple GWB source data sets. We search for these injected sources using Bayesian PTA analysis techniques to assess recovery and separability of multiple astrophysical and cosmological backgrounds. For a GWB due to supermassive black hole binaries and an underlying weaker background due to primordial gravitational waves with a GW energy-density ratio of ΩPGW/ΩSMBHB= 0.5, the Bayes’ factor for a second process exceeds unity at 17 yr, and increases with additional data. At 20 yr of data, we are able to constrain the spectral index and amplitude of the weaker GWB at this density ratio to a fractional uncertainty of 64% and 110%, respectively, using current PTA methods and techniques. Using these methods and findings, we outline a basic protocol to search for multiple backgrounds in future PTA data sets.
more »
« less
- PAR ID:
- 10376139
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 938
- Issue:
- 2
- ISSN:
- 0004-637X
- Format(s):
- Medium: X Size: Article No. 115
- Size(s):
- Article No. 115
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Pulsar timing arrays have found evidence for a low-frequency gravitational-wave background (GWB). Assuming that the GWB is produced by supermassive black hole binaries (SMBHBs), the next gravitational-wave (GW) signals astronomers anticipate are continuous waves (CWs) from single SMBHBs and their associated GWB anisotropy. The prospects for detecting CWs and anisotropy are highly dependent on the astrophysics of SMBHB populations. Thus, information from single sources can break degeneracies in astrophysical models and place much more stringent constraints than the GWB alone. We simulate and evolve SMBHB populations, model their GWs, and calculate their anisotropy and detectability. We investigate how varying components of our semianalytic model, including the galaxy stellar mass function, the SMBH–host galaxy relation (MBH–Mbulge), and the binary evolution prescription, impact the expected detections. The CW occurrence rate is greatest for few total binaries, high SMBHB masses, large scatter inMBH–Mbulge, and long hardening times. The occurrence rate depends most on the binary evolution parameters, implying that CWs offer a novel avenue to probe binary evolution. The most detectable CW sources are in the lowest frequency bin for a 16.03 yr PTA, have masses from ∼109to 1010M⊙, and are ∼1 Gpc away. The level of anisotropy increases with frequency, with the angular power spectrum over multipole modesℓvarying in low-frequencyCℓ>0/C0from ∼5 × 10−3to ∼2 × 10−1, depending on the model; typical values are near current upper limits. Observing this anisotropy would support SMBHB models for the GWB over cosmological models, which tend to be isotropic.more » « less
-
Abstract As pulsar timing arrays (PTAs) transition into the detection era of the stochastic gravitational wave background (GWB), it is important for PTA collaborations to review and possibly revise their observing campaigns. The detection of a “single source” would be a boon for gravitational astrophysics, as such a source would emit gravitational waves for millions of years in the PTA frequency band. Here we present generic methods for studying the effects of various observational strategies, taking advantage of detector sensitivity curves, i.e., noise-averaged, frequency-domain detection statistics. The statistical basis for these methods is presented along with myriad examples of how to tune a detector towards single, deterministic signals or a stochastic background. We demonstrate that trading observations of the worst pulsars for high cadence campaigns on the best pulsars increases sensitivity to single sources at high frequencies while hedging losses in GWB and single source sensitivity at low frequencies. We also find that sky-targeted observing campaigns yield minimal sensitivity improvements compared with other PTA tuning options. Lastly, we show the importance of the uncorrelated half of the GWB, i.e. the pulsar-term, as an increasingly prominent sources of noise and show the impact of this emerging noise source on various PTA configurations.more » « less
-
Abstract Recently, many different pulsar timing array (PTA) collaborations have reported strong evidence for a common stochastic process in their data sets. The reported amplitudes are in tension with previously computed upper limits. In this paper, we investigate how using a subset of a set of pulsars biases Bayesian upper limit recovery. We generate 500 simulated PTA data sets, based on the NANOGrav 11 yr data set with an injected stochastic gravitational-wave background (GWB). We then compute the upper limits by sampling the individual pulsar likelihoods, and combine them through a factorized version of the PTA likelihood to obtain upper limits on the GWB amplitude, using different numbers of pulsars. We find that it is possible to recover an upper limit (95% credible interval) below the injected value, and that it is significantly more likely for this to occur when using a subset of pulsars to compute the upper limit. When picking pulsars to induce the maximum possible bias, we find that the 95% Bayesian upper limit recovered is below the injected value in 10.6% of the realizations (53 of 500). Further, we find that if we choose a subset of pulsars in order to obtain a lower upper limit than when using the full set of pulsars, the distribution of the upper limits obtained from these 500 realizations is shifted to lower-amplitude values.more » « less
-
Abstract Analyses of pulsar timing data have provided evidence for a stochastic gravitational wave background in the nanohertz frequency band. The most plausible source of this background is the superposition of signals from millions of supermassive black hole binaries. The standard statistical techniques used to search for this background and assess its significance make several simplifying assumptions, namely (i) Gaussianity, (ii) isotropy, and most often, (iii) a power-law spectrum. However, a stochastic background from a finite collection of binaries does not exactly satisfy any of these assumptions. To understand the effect of these assumptions, we test standard analysis techniques on a large collection of realistic simulated data sets. The data-set length, observing schedule, and noise levels were chosen to emulate the NANOGrav 15 yr data set. Simulated signals from millions of binaries drawn from models based on the Illustris cosmological hydrodynamical simulation were added to the data. We find that the standard statistical methods perform remarkably well on these simulated data sets, even though their fundamental assumptions are not strictly met. They are able to achieve a confident detection of the background. However, even for a fixed set of astrophysical parameters, different realizations of the universe result in a large variance in the significance and recovered parameters of the background. We also find that the presence of loud individual binaries can bias the spectral recovery of the background if we do not account for them.more » « less
