Search for: All records

Abstract While supermassive black hole (SMBH) binaries are not the only viable source for the low-frequency gravitational wave background (GWB) signal evidenced by the most recent pulsar timing array (PTA) data sets, they are expected to be the most likely. Thus, connecting the measured PTA GWB spectrum and the underlying physics governing the demographics and dynamics of SMBH binaries is extremely important. Previously, Gaussian processes (GPs) and dense neural networks have been used to make such a connection by being built as conditional emulators; their input is some selected evolution or environmental SMBH binary parameters and their output is the emulated mean and standard deviation of the GWB strain ensemble distribution over many Universes. In this paper, we use a normalizing flow (NF) emulator that is trained on the entirety of the GWB strain ensemble distribution, rather than only mean and standard deviation. As a result, we can predict strain distributions that mirror underlying simulations very closely while also capturing frequency covariances in the strain distributions as well as statistical complexities such as tails, non-Gaussianities, and multimodalities that are otherwise not learnable by existing techniques. In particular, we feature various comparisons between the NF-based emulator and the GP approach used extensively in past efforts. Our analyses conclude that the NF-based emulator not only outperforms GPs in the ease and computational cost of training but also outperforms in the fidelity of the emulated GWB strain ensemble distributions. 

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 (M_BH–M_bulge), and the binary evolution prescription, impact the expected detections. The CW occurrence rate is greatest for few total binaries, high SMBHB masses, large scatter inM_BH–M_bulge, 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 ∼10⁹to 10¹⁰M_⊙, 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/C₀from ∼5 × 10⁻³to ∼2 × 10⁻¹, 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. 

Abstract Pulsar timing array observations have found evidence for an isotropic gravitational-wave background with the Hellings–Downs angular correlations between pulsar pairs. This interpretation hinges on the measured shape of the angular correlations, which is predominantly quadrupolar under general relativity. Here we explore a more flexible parameterization: we expand the angular correlations into a sum of Legendre polynomials and use a Bayesian analysis to constrain their coefficients with the 15 yr pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). When including Legendre polynomials with multipolesℓ≥ 2, we only find a significant signal in the quadrupole with an amplitude consistent with general relativity and nonzero at the ∼95% confidence level and a Bayes factor of 200. When we include multipolesℓ≤ 1, the Bayes factor evidence for quadrupole correlations decreases by more than an order of magnitude due to evidence for a monopolar signal at approximately 4 nHz, which has also been noted in previous analyses of the NANOGrav 15 yr data. Further work needs to be done in order to better characterize the properties of this monopolar signal and its effect on the evidence for quadrupolar angular correlations. 

Abstract Evidence has emerged for a stochastic signal correlated among 67 pulsars within the 15 yr pulsar-timing data set compiled by the NANOGrav collaboration. Similar signals have been found in data from the European, Indian, Parkes, and Chinese pulsar timing arrays. This signal has been interpreted as indicative of the presence of a nanohertz stochastic gravitational-wave background (GWB). To explore the internal consistency of this result, we investigate how the recovered signal strength changes as we remove the pulsars one by one from the data set. We calculate the signal strength using the (noise-marginalized) optimal statistic, a frequentist metric designed to measure the correlated excess power in the residuals of the arrival times of the radio pulses. We identify several features emerging from this analysis that were initially unexpected. The significance of these features, however, can only be assessed by comparing the real data to synthetic data sets. After conducting identical analyses on simulated data sets, we do not find anything inconsistent with the presence of a stochastic GWB in the NANOGrav 15 yr data. The methodologies developed here can offer additional tools for application to future, more sensitive data sets. While this analysis provides an internal consistency check of the NANOGrav results, it does not eliminate the necessity for additional investigations that could identify potential systematics or uncover unmodeled physical phenomena in the data. 

Abstract The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has reported evidence for the presence of an isotropic nanohertz gravitational-wave background (GWB) in its 15 yr data set. However, if the GWB is produced by a population of inspiraling supermassive black hole binary (SMBHB) systems, then the background is predicted to be anisotropic, depending on the distribution of these systems in the local Universe and the statistical properties of the SMBHB population. In this work, we search for anisotropy in the GWB using multiple methods and bases to describe the distribution of the GWB power on the sky. We do not find significant evidence of anisotropy. By modeling the angular power distribution as a sum over spherical harmonics (where the coefficients are not bound to always generate positive power everywhere), we find that the Bayesian 95% upper limit on the level of dipole anisotropy is (C_l=1/C_l=0) < 27%. This is similar to the upper limit derived under the constraint of positive power everywhere, indicating that the dipole may be close to the data-informed regime. By contrast, the constraints on anisotropy at higher spherical-harmonic multipoles are strongly prior dominated. We also derive conservative estimates on the anisotropy expected from a random distribution of SMBHB systems using astrophysical simulations conditioned on the isotropic GWB inferred in the 15 yr data set and show that this data set has sufficient sensitivity to probe a large fraction of the predicted level of anisotropy. We end by highlighting the opportunities and challenges in searching for anisotropy in pulsar timing array data. 

Abstract The NANOGrav 15 yr data set shows evidence for the presence of a low-frequency gravitational-wave background (GWB). While many physical processes can source such low-frequency gravitational waves, here we analyze the signal as coming from a population of supermassive black hole (SMBH) binaries distributed throughout the Universe. We show that astrophysically motivated models of SMBH binary populations are able to reproduce both the amplitude and shape of the observed low-frequency gravitational-wave spectrum. While multiple model variations are able to reproduce the GWB spectrum at our current measurement precision, our results highlight the importance of accurately modeling binary evolution for producing realistic GWB spectra. Additionally, while reasonable parameters are able to reproduce the 15 yr observations, the implied GWB amplitude necessitates either a large number of parameters to be at the edges of expected values or a small number of parameters to be notably different from standard expectations. While we are not yet able to definitively establish the origin of the inferred GWB signal, the consistency of the signal with astrophysical expectations offers a tantalizing prospect for confirming that SMBH binaries are able to form, reach subparsec separations, and eventually coalesce. As the significance grows over time, higher-order features of the GWB spectrum will definitively determine the nature of the GWB and allow for novel constraints on SMBH populations.