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 dataset. ...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, and place a Bayesian \(95\%\) upper limit on the level of broadband anisotropy such that \((C_{l>0} / C_{l=0}) < 20\%\). 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 dataset, and show that this dataset 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.
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-year period and low redshift ...($\sim0.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 on 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-year 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 eccentricity $e_0<0.5$
and the symmetric mass ratio $\eta>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\pm3.7$ ns for the Earth term-only and $81.74\pm0.86$
ns for the Earth+Pulsar term searches for $e_0<0.5$ and $\eta>0.1$. Similar 95%
upper limits on the chirp mass are $(1.98 \pm 0.05) \times 10^9\,M_{\odot}$ and
$(1.81 \pm 0.01) \times 10^9\,M_{\odot}$. These upper limits, while less
stringent than those calculated from a circular binary search in the NANOGrav
12.5-year data set, are consistent with the SMBHB model of 3C 66B developed
from electromagnetic observations.
The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical ...questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe.
Pulsar timing array collaborations, such as the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), are seeking to detect nanohertz gravitational waves emitted by supermassive ...black hole binaries formed in the aftermath of galaxy mergers. We have searched for continuous waves from individual circular supermassive black hole binaries using the NANOGrav's recent 12.5-year data set. We created new methods to accurately model the uncertainties on pulsar distances in our analysis, and we implemented new techniques to account for a common red noise process in pulsar timing array data sets while searching for deterministic gravitational wave signals, including continuous waves. As we found no evidence for continuous waves in our data, we placed 95\% upper limits on the strain amplitude of continuous waves emitted by these sources. At our most sensitive frequency of 7.65 nanohertz, we placed a sky-averaged limit of \(h_0 < \) \((6.82 \pm 0.35) \times 10^{-15}\), and \(h_0 <\) \((2.66 \pm 0.15) \times 10^{-15}\) in our most sensitive sky location. Finally, we placed a multi-messenger limit of \(\mathcal{M} <\) \((1.41 \pm 0.02) \times 10^9 M_\odot\) on the chirp mass of the supermassive black hole binary candidate 3C~66B.
We search for a first-order phase transition gravitational wave signal in 45 pulsars from the NANOGrav 12.5 year dataset. We find that the data can be modeled in terms of a strong first order phase ...transition taking place at temperatures below the electroweak scale. However, we do not observe any strong preference for a phase-transition interpretation of the signal over the standard astrophysical interpretation in terms of supermassive black holes mergers; but we expect to gain additional discriminating power with future datasets, improving the signal to noise ratio and extending the sensitivity window to lower frequencies. An interesting open question is how well gravitational wave observatories could separate such signals.
We search NANOGrav's 12.5-year data set for evidence of a gravitational wave background (GWB) with all the spatial correlations allowed by general metric theories of gravity. We find no substantial ...evidence in favor of the existence of such correlations in our data. We find that scalar-transverse (ST) correlations yield signal-to-noise ratios and Bayes factors that are higher than quadrupolar (tensor transverse, TT) correlations. Specifically, we find ST correlations with a signal-to-noise ratio of 2.8 that are preferred over TT correlations (Hellings and Downs correlations) with Bayesian odds of about 20:1. However, the significance of ST correlations is reduced dramatically when we include modeling of the Solar System ephemeris systematics and/or remove pulsar J0030\(+\)0451 entirely from consideration. Even taking the nominal signal-to-noise ratios at face value, analyses of simulated data sets show that such values are not extremely unlikely to be observed in cases where only the usual TT modes are present in the GWB. In the absence of a detection of any polarization mode of gravity, we place upper limits on their amplitudes for a spectral index of \(\gamma = 5\) and a reference frequency of \(f_\text{yr} = 1 \text{yr}^{-1}\). Among the upper limits for eight general families of metric theories of gravity, we find the values of \(A^{95\%}_{TT} = (9.7 \pm 0.4)\times 10^{-16}\) and \(A^{95\%}_{ST} = (1.4 \pm 0.03)\times 10^{-15}\) for the family of metric spacetime theories that contain both TT and ST modes.
Pulsar timing array (PTA) collaborations in North America, Australia, and Europe, have been exploiting the exquisite timing precision of millisecond pulsars over decades of observations to search for ...correlated timing deviations induced by gravitational waves (GWs). PTAs are sensitive to the frequency band ranging just below 1 nanohertz to a few tens of microhertz. The discovery space of this band is potentially rich with populations of inspiraling supermassive black-holes binaries, decaying cosmic string networks, relic post-inflation GWs, and even non-GW imprints of axionic dark matter. This article aims to provide an understanding of the exciting open science questions in cosmology, galaxy evolution, and fundamental physics that will be addressed by the detection and study of GWs through PTAs. The focus of the article is on providing an understanding of the mechanisms by which PTAs can address specific questions in these fields, and to outline some of the subtleties and difficulties in each case. The material included is weighted most heavily towards the questions which we expect will be answered in the near-term with PTAs; however, we have made efforts to include most currently anticipated applications of nanohertz GWs.
Supermassive black hole binaries (SMBHBs) should form frequently in galactic nuclei as a result of galaxy mergers. At sub-parsec separations, binaries become strong sources of low-frequency ...gravitational waves (GWs), targeted by Pulsar Timing Arrays (PTAs). We used recent upper limits on continuous GWs from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11yr dataset to place constraints on putative SMBHBs in nearby massive galaxies. We compiled a comprehensive catalog of ~44,000 galaxies in the local universe (up to redshift ~0.05) and populated them with hypothetical binaries, assuming that the total mass of the binary is equal to the SMBH mass derived from global scaling relations. Assuming circular equal-mass binaries emitting at NANOGrav's most sensitive frequency of 8nHz, we found that 216 galaxies are within NANOGrav's sensitivity volume. We ranked the potential SMBHBs based on GW detectability by calculating the total signal-to-noise ratio (S/N) such binaries would induce within the NANOGrav array. We placed constraints on the chirp mass and mass ratio of the 216 hypothetical binaries. For 19 galaxies, only very unequal-mass binaries are allowed, with the mass of the secondary less than 10 percent that of the primary, roughly comparable to constraints on a SMBHB in the Milky Way. Additionally, we were able to exclude binaries delivered by major mergers (mass ratio of at least 1/4) for several of these galaxies. We also derived the first limit on the density of binaries delivered by major mergers purely based on GW data.