skip to main content


Search for: All records

Creators/Authors contains: "Vigeland, Sarah J."

Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

  1. Abstract

    Pulsar timing arrays (PTAs) are Galactic-scale gravitational wave (GW) detectors consisting of precisely timed pulsars distributed across the sky. Within the decade, PTAs are expected to detect nanohertz GWs emitted by close-separation supermassive black hole binaries (SMBHBs), thereby opening up the low-frequency end of the GW spectrum for science. Individual SMBHBs which power active galactic nuclei are also promising multi-messenger sources; they may be identified via theoretically predicted electromagnetic (EM) signatures and be followed up by PTAs for GW observations. In this work, we study the detection and parameter estimation prospects of a PTA which targets EM-selected SMBHBs. Adopting a simulated Galactic millisecond pulsar population, we envisage three different pulsar timing campaigns which observe three mock sources at different sky locations. We find that an all-sky PTA which times the best pulsars is an optimal and feasible approach to observe EM-selected SMBHBs and measure their source parameters to high precision (i.e., comparable to or better than conventional EM measurements). We discuss the implications of our findings in the context of future PTA experiments with the planned Deep Synoptic Array-2000 and the multi-messenger studies of SMBHBs such as the well-known binary candidate OJ 287.

     
    more » « less
  2. 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
  3. Abstract Pulsar timing array (PTA) experiments are becoming increasingly sensitive to gravitational waves (GWs) in the nanohertz frequency range, where the main astrophysical sources are supermassive black hole binaries (SMBHBs), which are expected to form following galaxy mergers. Some of these individual SMBHBs may power active galactic nuclei, and thus their binary parameters could be obtained electromagnetically, which makes it possible to apply electromagnetic (EM) information to aid the search for a GW signal in PTA data. In this work, we investigate the effects of such an EM-informed search on binary detection and parameter estimation by performing mock data analyses on simulated PTA data sets. We find that by applying EM priors, the Bayes factor of some injected signals with originally marginal or sub-threshold detectability (i.e., Bayes factor ∼1) can increase by a factor of a few to an order of magnitude, and thus an EM-informed targeted search is able to find hints of a signal when an uninformed search fails to find any. Additionally, by combining EM and GW data, one can achieve an overall improvement in parameter estimation, regardless of the source’s sky location or GW frequency. We discuss the implications for the multi-messenger studies of SMBHBs with PTAs. 
    more » « less
  4. 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 ΩPGWSMBHB= 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
  5. null (Ed.)
  6. 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 eccentricitye0< 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 fore0< 0.5 andη> 0.1. Similar 95% upper limits on the chirp mass are (1.98 ± 0.05) × 109and (1.81 ± 0.01) × 109M. 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.

     
    more » « less
  7. 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 (Cl=1/Cl=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.

     
    more » « less
    Free, publicly-accessible full text available October 1, 2024
  8. Abstract

    We present the results of a Bayesian search for gravitational wave (GW) memory in the NANOGrav 12.5 yr data set. We find no convincing evidence for any gravitational wave memory signals in this data set. We find a Bayes factor of 2.8 in favor of a model that includes a memory signal and common spatially uncorrelated red noise (CURN) compared to a model including only a CURN. However, further investigation shows that a disproportionate amount of support for the memory signal comes from three dubious pulsars. Using a more flexible red-noise model in these pulsars reduces the Bayes factor to 1.3. Having found no compelling evidence, we go on to place upper limits on the strain amplitude of GW memory events as a function of sky location and event epoch. These upper limits are computed using a signal model that assumes the existence of a common, spatially uncorrelated red noise in addition to a GW memory signal. The median strain upper limit as a function of sky position is approximately 3.3 × 10−14. We also find that there are some differences in the upper limits as a function of sky position centered around PSR J0613−0200. This suggests that this pulsar has some excess noise that can be confounded with GW memory. Finally, the upper limits as a function of burst epoch continue to improve at later epochs. This improvement is attributable to the continued growth of the pulsar timing array.

     
    more » « less
  9. Abstract Evidence for a low-frequency stochastic gravitational-wave background has recently been reported based on analyses of pulsar timing array data. The most likely source of such a background is a population of supermassive black hole binaries, the loudest of which may be individually detected in these data sets. Here we present the search for individual supermassive black hole binaries in the NANOGrav 15 yr data set. We introduce several new techniques, which enhance the efficiency and modeling accuracy of the analysis. The search uncovered weak evidence for two candidate signals, one with a gravitational-wave frequency of ∼4 nHz, and another at ∼170 nHz. The significance of the low-frequency candidate was greatly diminished when Hellings–Downs correlations were included in the background model. The high-frequency candidate was discounted due to the lack of a plausible host galaxy, the unlikely astrophysical prior odds of finding such a source, and since most of its support comes from a single pulsar with a commensurate binary period. Finding no compelling evidence for signals from individual binary systems, we place upper limits on the strain amplitude of gravitational waves emitted by such systems. At our most sensitive frequency of 6 nHz, we place a sky-averaged 95% upper limit of 8 × 10 −15 on the strain amplitude. We also calculate an exclusion volume and a corresponding effective radius, within which we can rule out the presence of black hole binaries emitting at a given frequency. 
    more » « less
    Free, publicly-accessible full text available July 1, 2024