Search for: All records

Abstract Massive black hole binaries (MBHBs) produce gravitational waves (GWs) that are detectable with pulsar timing arrays. We determine the properties of the host galaxies of simulated MBHBs at the time they are producing detectable GW signals. The population of MBHB systems we evaluate is from theIllustriscosmological simulations taken in tandem with post processing semi-analytic models of environmental factors in the evolution of binaries. Upon evolving to the GW frequency regime accessible by pulsar timing arrays, we calculate the detection probability of each system using a variety of different values for pulsar noise characteristics in a plausible near-future International Pulsar Timing Array dataset. We find that detectable systems have host galaxies that are clearly distinct from the overall binary population and from most galaxies in general. With conservative noise factors, we find that host stellar metallicity, for example, peaks at $\sim 2Z_{\odot }$as opposed to the total population of galaxies which peaks at $\sim 0.6Z_{\odot }$. Additionally, the most detectable systems are much brighter in magnitude and more red in color than the overall population, indicating their likely identity as large ellipticals with diminished star formation. These results can be used to develop effective search strategies for identifying host galaxies and electromagnetic counterparts following GW detection by pulsar timing arrays. 

Abstract Supermassive black hole binaries (SMBHBs) are thought to form in galaxy mergers, possessing the potential to produce electromagnetic (EM) radiation as well as gravitational waves (GWs) detectable with pulsar timing arrays (PTAs). Once GWs from individually resolved SMBHBs are detected, the identification of the host galaxy will be a major challenge due to the ambiguity in possible EM signatures and the poor localization capability of PTAs. To aid EM observations in choosing follow-up sources, we use NANOGrav’s galaxy catalog to quantify the number of plausible hosts in both realistic and idealistic scenarios. We outline a host identification pipeline that injects a single-source GW signal into a simulated PTA data set, recovers the signal using production-level techniques, quantifies the localization region and number of galaxies contained therein, and finally imposes cuts on the galaxies using parameter estimates from the GW search. In an ideal case, the 90% credible areas span 29–241 deg², containing about 14–341 galaxies. After cuts, the number of galaxies remaining ranges from 22 at worst to one true host at best. In a realistic case, these areas range from 287 to 530 deg²and enclose about 285–1238 galaxies. After cuts, the number of galaxies is 397 at worst and 27 at best. While the signal-to-noise ratio is the primary determinant of the localization area of a given source, we find that the area is also influenced by the proximity to nearby pulsars on the sky and the binary chirp mass. 

Abstract A population of compact object binaries emitting gravitational waves that are not individually resolvable will form a stochastic gravitational-wave signal. While the expected spectrum over population realizations is well known from Phinney, its higher-order moments have not been fully studied before or computed in the case of arbitrary binary evolution. We calculate analytic scaling relationships as a function of gravitational-wave frequency for the statistical variance, skewness, and kurtosis of a stochastic gravitational-wave signal over population realizations due to finite source effects. If the time derivative of the binary orbital frequency can be expressed as a power law in frequency, we find that these moment quantities also take the form of power-law relationships. We also develop a numerical population synthesis framework against which we compare our analytic results, finding excellent agreement. These new scaling relationships provide physical context to understanding spectral fluctuations in a gravitational-wave background signal and may provide additional information that can aid in explaining the origin of the nanohertz-frequency signal observed by pulsar timing array campaigns. 

Abstract The cosmic merger history of supermassive black hole binaries (SMBHBs) is expected to produce a low-frequency gravitational wave background (GWB). Here we investigate how signs of the discrete nature of this GWB can manifest in pulsar timing arrays (PTAs) through excursions from, and breaks in, the expected $f_{\mathrm{GW}}^{-2/3}$power law of the GWB strain spectrum. To do this, we create a semianalytic SMBHB population model, fit to North American Nanohertz Observatory for Gravitational Waves (NANOGrav’s) 15 yr GWB amplitude, and with 1000 realizations, we study the populations’ characteristic strain and residual spectra. Comparing our models to the NANOGrav 15 yr spectrum, we find two interesting excursions from the power law. The first, at 2 nHz, is below our GWB realizations with ap-value significancep= 0.05–0.06 (≈1.8σ–1.9σ). The second, at 16 nHz, is above our GWB realizations withp= 0.04–0.15 (≈1.4σ–2.1σ). We explore the properties of a loud SMBHB that could cause such an excursion. Our simulations also show that the expected number of SMBHBs decreases by 3 orders of magnitude, from ∼10⁶to ∼10³, between 2 and 20 nHz. This causes a break in the strain spectrum as the stochasticity of the background breaks down at ${26}_{-19}^{+28}\mathrm{nHz}$, consistent with predictions pre-dating GWB measurements. The diminished GWB signal from SMBHBs at frequencies above the 26 nHz break opens a window for PTAs to detect continuous GWs from individual SMBHBs or GWs from the early Universe. 

Abstract The radio galaxy 3C 66B has been hypothesized to host a supermassive black hole binary (SMBHB) at its center based on electromagnetic observations. Its apparent 1.05 yr period and low redshift (∼0.02) make it an interesting testbed to search for low-frequency gravitational waves (GWs) using pulsar timing array (PTA) experiments. This source has been subjected to multiple searches for continuous GWs from a circular SMBHB, resulting in progressively more stringent constraints on its GW amplitude and chirp mass. In this paper, we develop a pipeline for performing Bayesian targeted searches for eccentric SMBHBs in PTA data sets, and test its efficacy by applying it to simulated data sets with varying injected signal strengths. We also search for a realistic eccentric SMBHB source in 3C 66B using the NANOGrav 12.5 yr data set employing PTA signal models containing Earth term-only as well as Earth+pulsar term contributions using this pipeline. Due to limitations in our PTA signal model, we get meaningful results only when the initial eccentricitye₀< 0.5 and the symmetric mass ratioη> 0.1. We find no evidence for an eccentric SMBHB signal in our data, and therefore place 95% upper limits on the PTA signal amplitude of 88.1 ± 3.7 ns for the Earth term-only and 81.74 ± 0.86 ns for the Earth+pulsar term searches fore₀< 0.5 andη> 0.1. Similar 95% upper limits on the chirp mass are (1.98 ± 0.05) × 10⁹and (1.81 ± 0.01) × 10⁹M_☉. These upper limits, while less stringent than those calculated from a circular binary search in the NANOGrav 12.5 yr data set, are consistent with the SMBHB model of 3C 66B developed from electromagnetic observations.