Search

The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTAs that constitute the International Pulsar Timing Array (IPTA). We show that despite making different modeling choices, there is no significant difference in the GWB parameters that are measured by the different PTAs, agreeing within $1\sigma$. The pulsar noise parameters are also consistent between different PTAs for the majority of the pulsars included in these analyses. We bridge the differences in modeling choices by adopting a standardized noise model for all pulsars and PTAs, finding that under this model there is a reduction in the tension in the pulsar noise parameters. As part of this reanalysis, we "extended" each PTA's data set by adding extra pulsars that were not timed by that PTA. Under these extensions, we find better constraints on the GWB amplitude and a higher signal-to-noise ratio for the Hellings and Downs correlations. These extensions serve as a prelude to the benefits offered by a full combination of data across all pulsars in the IPTA, i.e., the IPTA's Data Release 3, which will involve not just adding in additional pulsars, but also including data from all three PTAs where any given pulsar is timed by more than as single PTA., Comment: 21 pages, 9 figures, submitted to ApJ

The International Pulsar Timing Array 2nd data release is the combination of datasets from worldwide collaborations. In this study, we search for continuous waves: gravitational wave signals produced by individual supermassive black hole binaries in the local universe. We consider binaries on circular orbits and neglect the evolution of orbital frequency over the observational span. We find no evidence for such signals and set sky averaged 95% upper limits on their amplitude h 95 . The most sensitive frequency is 10nHz with h 95 = 9.1 10-15 . We achieved the best upper limit to date at low and high frequencies of the PTA band thanks to improved effective cadence of observations. In our analysis, we have taken into account the recently discovered common red noise process, which has an impact at low frequencies. We also find that the peculiar noise features present in some pulsars data must be taken into account to reduce the false alarm. We show that using custom noise models is essential in searching for continuous gravitational wave signals and setting the upper limit.

We present extensive multi-frequency VLA and VLBA observations of the radio-bright supernova (SN) IIb SN 2004C that span $\sim(40-2793)$ days post-explosion. We interpret the temporal evolution of the radio spectral energy distribution (SED) in the context of synchrotron self-absorbed (SSA) emission from the explosion's forward shock as it expands in the circumstellar medium (CSM) previously sculpted by the mass-loss history of the stellar progenitor. VLBA observations and modeling of the VLA data point to a blastwave with average velocity $\sim0.06c$ that carries an energy of $\sim 10^{49}$ erg. Our modeling further reveals a flat CSM density profile $\rho_{\rm{CSM}} \propto R^{-0.03 \pm0.22}$ up to a break radius $R_{br} \approx (1.96 \pm 0.10) \times 10^{16}$ cm, with a steep density gradient following $\rho_{\rm{CSM}} \propto R^{-2.3 \pm 0.5}$ at larger radii. We infer that the flat part of the density profile corresponds to a CSM shell with mass $\sim0.021 M_{\odot}$, and that the progenitor's effective mass-loss rate varied with time over the range $(50-500) \times 10^{-5} M_{\odot} \rm{yr}^{-1}$ for an adopted wind velocity $v_w =1000$ km $s^{-1}$ and shock microphysical parameters ${\epsilon}_e = 0.1, {\epsilon}_B = 0.01$. These results add to the mounting observational evidence for departures from the traditional single-wind mass-loss scenarios in evolved, massive stars in the centuries leading up to core collapse. Potentially viable scenarios include mass loss powered by gravity waves and/or interaction with a binary companion., Comment: 32 pages, 10 figures, submitted to ApJ

We searched for an isotropic stochastic gravitational wave background in the second data release of the International Pulsar Timing Array, a global collaboration synthesizing decadal-length pulsar-timing campaigns in North America, Europe, and Australia. In our reference search for a power law strain spectrum of the form $h_c = A(f/1\,\mathrm{yr}^{-1})^{\alpha}$, we found strong evidence for a spectrally-similar low-frequency stochastic process of amplitude $A = 3.8^{+6.3}_{-2.5}\times10^{-15}$ and spectral index $\alpha = -0.5 \pm 0.5$, where the uncertainties represent 95\% credible regions, using information from the auto- and cross-correlation terms between the pulsars in the array. For a spectral index of $\alpha = -2/3$, as expected from a population of inspiralling supermassive black hole binaries, the recovered amplitude is $A = 2.8^{+1.2}_{-0.8}\times10^{-15}$. Nonetheless, no significant evidence of the Hellings-Downs correlations that would indicate a gravitational-wave origin was found. We also analyzed the constituent data from the individual pulsar timing arrays in a consistent way, and clearly demonstrate that the combined international data set is more sensitive. Furthermore, we demonstrate that this combined data set produces comparable constraints to recent single-array data sets which have more data than the constituent parts of the combination. Future international data releases will deliver increased sensitivity to gravitational wave radiation, and significantly increase the detection probability., Comment: 17 pages, 12 figures, accepted in MNRAS

The regularity of pulsar emissions becomes apparent once we reference the pulses' times of arrivals to the inertial rest frame of the solar system. It follows that errors in the determination of Earth's position with respect to the solar-system barycenter can appear as a time-correlated bias in pulsar-timing residual time series, affecting the searches for low-frequency gravitational waves performed with pulsar timing arrays. Indeed, recent array datasets yield different gravitational-wave background upper limits and detection statistics when analyzed with different solar-system ephemerides. Crucially, the ephemerides do not generally provide usable error representations. In this article we describe the motivation, construction, and application of a physical model of solar-system ephemeris uncertainties, which focuses on the degrees of freedom (Jupiter's orbital elements) most relevant to gravitational-wave searches with pulsar timing arrays. This model, BayesEphem, was used to derive ephemeris-robust results in NANOGrav's 11-yr stochastic-background search, and it provides a foundation for future searches by NANOGrav and other consortia. The analysis and simulations reported here suggest that ephemeris modeling reduces the gravitational-wave sensitivity of the 11-yr dataset; and that this degeneracy will vanish with improved ephemerides and with the longer pulsar timing datasets that will become available in the near future., Comment: Fixed typo in author list. Code that supports all calculations and figures is available at github.com/nanograv/11yr_stochastic_analysis/tree/master/bayesephem

The mergers of supermassive black hole binaries (SMBHBs) promise to be incredible sources of gravitational waves (GWs). While the oscillatory part of the merger gravitational waveform will be outside the frequency sensitivity range of pulsar timing arrays (PTAs), the non-oscillatory GW memory effect is detectable. Further, any burst of gravitational waves will produce GW memory, making memory a useful probe of unmodeled exotic sources and new physics. We searched the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GW memory. This dataset is sensitive to very low frequency GWs of $\sim3$ to $400$ nHz (periods of $\sim11$ yr $-$ $1$ mon). Finding no evidence for GWs, we placed limits on the strain amplitude of GW memory events during the observation period. We then used the strain upper limits to place limits on the rate of GW memory causing events. At a strain of $2.5\times10^{-14}$, corresponding to the median upper limit as a function of source sky position, we set a limit on the rate of GW memory events at $<0.4$ yr$^{-1}$. That strain corresponds to a SMBHB merger with reduced mass of $\eta M \sim 2\times10^{10}M_\odot$ and inclination of $\iota=\pi/3$ at a distance of 1 Gpc. As a test of our analysis, we analyzed the NANOGrav 9-year data set as well. This analysis found an anomolous signal, which does not appear in the 11-year data set. This signal is not a GW, and its origin remains unknown., Comment: 10 pages, 6 figures, submitted to ApJ

An ensemble of inspiraling supermassive black hole binaries should produce a stochastic background of very low frequency gravitational waves. This stochastic background is predicted to be a power law, with a spectral index of -2/3, and it should be detectable by a network of precisely timed millisecond pulsars, widely distributed on the sky. This paper reports a new "time slicing" analysis of the 11-year data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) using 34 millisecond pulsars. Methods to flag potential "false positive" signatures are developed, including techniques to identify responsible pulsars. Mitigation strategies are then presented. We demonstrate how an incorrect noise model can lead to spurious signals, and show how independently modeling noise across 30 Fourier components, spanning NANOGrav's frequency range, effectively diagnoses and absorbs the excess power in gravitational-wave searches. This results in a nominal, and expected, progression of our gravitational-wave statistics. Additionally we show that the first interstellar medium event in PSR J1713+0747 pollutes the common red noise process with low-spectral index noise, and use a tailored noise model to remove these effects., Comment: 14 pages, 13 figures, fixed typo in abstract of earlier version

Observations indicate that nearly all galaxies contain supermassive black holes (SMBHs) at their centers. When galaxies merge, their component black holes form SMBH binaries (SMBHBs), which emit low-frequency gravitational waves (GWs) that can be detected by pulsar timing arrays (PTAs). We have searched the recently-released North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GWs from individual SMBHBs in circular orbits. As we did not find strong evidence for GWs in our data, we placed 95\% upper limits on the strength of GWs from such sources as a function of GW frequency and sky location. We placed a sky-averaged upper limit on the GW strain of $h_0 < 7.3(3) \times 10^{-15}$ at $f_\mathrm{gw}= 8$ nHz. We also developed a technique to determine the significance of a particular signal in each pulsar using ``dropout' parameters as a way of identifying spurious signals in measurements from individual pulsars. We used our upper limits on the GW strain to place lower limits on the distances to individual SMBHBs. At the most-sensitive sky location, we ruled out SMBHBs emitting GWs with $f_\mathrm{gw}= 8$ nHz within 120 Mpc for $\mathcal{M} = 10^9 \, M_\odot$, and within 5.5 Gpc for $\mathcal{M} = 10^{10} \, M_\odot$. We also determined that there are no SMBHBs with $\mathcal{M} > 1.6 \times 10^9 \, M_\odot$ emitting GWs in the Virgo Cluster. Finally, we estimated the number of potentially detectable sources given our current strain upper limits based on galaxies in Two Micron All-Sky Survey (2MASS) and merger rates from the Illustris cosmological simulation project. Only 34 out of 75,000 realizations of the local Universe contained a detectable source, from which we concluded it was unsurprising that we did not detect any individual sources given our current sensitivity to GWs., Comment: 10 pages, 11 figures. Accepted by Astrophysical Journal. Please send any comments/questions to S. J. Vigeland (vigeland@uwm.edu)

The advent of time domain astronomy is revolutionizing our understanding of the Universe. Programs such as the Catalina Real-time Transient Survey (CRTS) or the Palomar Transient Factory (PTF) surveyed millions of objects for several years, allowing variability studies on large statistical samples. The inspection of $\approx$250k quasars in CRTS resulted in a catalogue of 111 potentially periodic sources, put forward as supermassive black hole binary (SMBHB) candidates. A similar investigation on PTF data yielded 33 candidates from a sample of $\approx$35k quasars. Working under the SMBHB hypothesis, we compute the implied SMBHB merger rate and we use it to construct the expected gravitational wave background (GWB) at nano-Hz frequencies, probed by pulsar timing arrays (PTAs). After correcting for incompleteness and assuming virial mass estimates, we find that the GWB implied by the CRTS sample exceeds the current most stringent PTA upper limits by almost an order of magnitude. After further correcting for the implicit bias in virial mass measurements, the implied GWB drops significantly but is still in tension with the most stringent PTA upper limits. Similar results hold for the PTF sample. Bayesian model selection shows that the null hypothesis (whereby the candidates are false positives) is preferred over the binary hypothesis at about $2.3\sigma$ and $3.6\sigma$ for the CRTS and PTF samples respectively. Although not decisive, our analysis highlights the potential of PTAs as astrophysical probes of individual SMBHB candidates and indicates that the CRTS and PTF samples are likely contaminated by several false positives., Comment: 14 pages, 11 figures, 3 tables. Resubmitted to the Astrophysical Journal after some major revision of the results including a proper estimate of the intrinsic mass of the binary candidates

Observations indicate that nearly all galaxies contain supermassive black holes at their centers. When galaxies merge, their component black holes form SMBH binaries (SMBHBs), which emit low-frequency gravitational waves (GWs) that can be detected by pulsar timing arrays. We have searched the North American Nanohertz Observatory for Gravitational Waves 11 yr data set for GWs from individual SMBHBs in circular orbits. As we did not find strong evidence for GWs in our data, we placed 95% upper limits on the strength of GWs from such sources. At f(gw) = 8 nHz, we placed a sky-averaged upper limit of h(0) < 7.3(3) × 10(exp −15). We also developed a technique to determine the significance of a particular signal in each pulsar using "dropout" parameters as a way of identifying spurious signals. From these upper limits, we ruled out SMBHBs emitting GWs f(gw) = 8 nHz within 120 Mpc for M = 10(exp 9) Solar Mass, and within 5.5 Gpc for M= 10(exp 10) Solar Mass at our most sensitive sky location. We also determined that there are no SMBHBs with M > 1.6 x 10(exp 9) Solar Mass emitting GWs with f(gw) = 2.8–317.8 nHz in the Virgo Cluster. Finally, we compared our strain upper limits to simulated populations of SMBHBs, based on galaxies in the Two Micron All-Sky Survey and merger rates from the Illustris cosmological simulation project, and found that only 34 out of 75,000 realizations of the local universe contained a detectable source.

Pulsar timing arrays (PTAs) are on the verge of detecting low-frequency gravitational waves (GWs)from supermassive black hole binaries (SMBHBs). With continued observations of a large sampleof millisecond pulsars, PTAs will reach this major milestone within the next decade. Already,SMBHB candidates are being identied by electromagnetic surveys in ever-increasing numbers;upcoming surveys will enhance our ability to detect and verify candidates, and will be instrumentalin identifying the host galaxies of GW sources. Multi-messenger (GW and electromagnetic) obser-vations of SMBHBs will revolutionize our understanding of the co-evolution of SMBHs with theirhost galaxies, the dynamical interactions between binaries and their galactic environments, and thefundamental physics of accretion. Multi-messenger observations can also make SMBHBs `standardsirens' for cosmological distance measurements out to z ~ 0.5 LIGO has already ushered in break-through insights in our knowledge of black holes. The multi-messenger detection of SMBHBs withPTAs will be a breakthrough in the years 2020-2030 and beyond, and prepare us for LISA to helpcomplete our views of black hole demographics and evolution at higher redshifts.