Gravitational-wave Background (GWB): A Breakthrough in Direct Detection

Pulsar timing arrays (PTAs) have revealed galaxy-size gravitational waves (GWs) in the form of a stochastic gravitational wave background (SGWB), correlating the radio pulses emitted by millisecond pulsars. This discovery naturally leads to the question of the origin and the nature of the SGWB; the latter is synonymous to testing how quadrupolar the inter-pulsar spatial correlation is. In this paper, we investigate the nature of the SGWB by considering correlations beyond the Hellings–Downs (HD) curve of Einstein’s general relativity. We scrutinize the HD and non-Einsteinian GW correlations with the North American Nanohertz Observatory for Gravitational Waves and the Chinese PTA data, and find that both data sets allow a graviton mass of mg ≲ 1.04 × 10−22 eV/c2 and subluminal traveling waves. We discuss gravitational physics scenarios beyond general relativity that could host non-Einsteinian GW correlations in the SGWB and highlight the importance of the cosmic variance inherited from stochastic variations across realizations in interpreting PTA observations.

We search for an isotropic stochastic gravitational-wave background (GWB) in the newly released 11 year data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). While we find no evidence for a GWB, we place constraints on a population of inspiraling supermassive black hole (SMBH) binaries, a network of decaying cosmic strings, and a primordial GWB. For the first time, we find that the GWB constraints are sensitive to the solar system ephemeris (SSE) model used and that SSE errors can mimic a GWB signal. We developed an approach that bridges systematic SSE differences, producing the first pulsar-timing array (PTA) constraints that are robust against SSE errors. We thus place a 95% upper limit on the GW-strain amplitude of A GWB < 1.45 × 10−15 at a frequency of f = 1 yr−1 for a fiducial f −2/3 power-law spectrum and with interpulsar correlations modeled. This is a factor of ∼2 improvement over the NANOGrav nine-year limit calculated using the same procedure. Previous PTA upper limits on the GWB (as well as their astrophysical and cosmological interpretations) will need revision in light of SSE systematic errors. We use our constraints to characterize the combined influence on the GWB of the stellar mass density in galactic cores, the eccentricity of SMBH binaries, and SMBH–galactic-bulge scaling relationships. We constrain the cosmic-string tension using recent simulations, yielding an SSE-marginalized 95% upper limit of Gμ < 5.3 × 10−11—a factor of ∼2 better than the published NANOGrav nine-year constraints. Our SSE-marginalized 95% upper limit on the energy density of a primordial GWB (for a radiation-dominated post-inflation universe) is ΩGWB(f) h 2 < 3.4 × 10−10.

We present observations and timing analyses of 68 millisecond pulsars (MSPs) comprising the 15 yr data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). NANOGrav is a pulsar timing array (PTA) experiment that is sensitive to low-frequency gravitational waves (GWs). This is NANOGrav’s fifth public data release, including both “narrowband” and “wideband” time-of-arrival (TOA) measurements and corresponding pulsar timing models. We have added 21 MSPs and extended our timing baselines by 3 yr, now spanning nearly 16 yr for some of our sources. The data were collected using the Arecibo Observatory, the Green Bank Telescope, and the Very Large Array between frequencies of 327 MHz and 3 GHz, with most sources observed approximately monthly. A number of notable methodological and procedural changes were made compared to our previous data sets. These improve the overall quality of the TOA data set and are part of the transition to new pulsar timing and PTA analysis software packages. For the first time, our data products are accompanied by a full suite of software to reproduce data reduction, analysis, and results. Our timing models include a variety of newly detected astrometric and binary pulsar parameters, including several significant improvements to pulsar mass constraints. We find that the time series of 23 pulsars contain detectable levels of red noise, 10 of which are new measurements. In this data set, we find evidence for a stochastic GW background.

Primordial black holes (PBHs) can constitute the predominant fraction of dark matter (DM) if PBHs reside in the currently unconstrained “sublunar” mass range. PBHs originating from scalar perturbations generated during inflation can naturally appear with a broad spectrum in a class of models. The resulting stochastic gravitational wave (GW) background generated from such PBH production can account for the recently reported North American Nanohertz Observatory for Gravitational Waves (NANOGrav) pulsar timing array data signal, and will be testable in future GW observations by interferometer-type experiments such as Laser Interferometer Space Antenna (LISA). We show that the broad mass function of such PBH DM is already being probed by Subaru Hyper Suprime-Cam (HSC) microlensing data and is consistent with a detected candidate event. Upcoming observations of HSC will be able to provide an independent definitive test of the stochastic GW signals originating from such PBH DM production scenarios.

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) conducts long-term timing campaigns on large numbers of millisecond pulsars (MSPs). Wideband (WB) timing of the NANOGrav datasets was first done in the 12.5yr dataset, and proved to be extremely helpful for measuring the time-varying dispersion measure, which can significantly limit timing precision. However, we are limited by the fact that NANOGrav only observes most pulsars once per month. In contrast, the Canadian Hydrogen Intensity Mapping Experiment (CHIME) Telescope currently observes most NANOGrav MSPs at nearly daily cadence at 600 MHz with a 400 MHz bandwidth. Using WB timing techniques, we combined one year of data taken at CHIME with data over the same time period from the NANOGrav 15yr Dataset for three MSPs.

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.

Gravitational waves from a variety of sources are predicted to superpose to create a stochastic background. This background is expected to contain unique information from throughout the history of the Universe that is unavailable through standard electromagnetic observations, making its study of fundamental importance to understanding the evolution of the Universe. We carry out a search for the stochastic background with the latest data from the LIGO and Virgo detectors. Consistent with predictions from most stochastic gravitational-wave background models, the data display no evidence of a stochastic gravitational-wave signal. Assuming a gravitational-wave spectrum of Ω_{GW}(f)=Ω_{α}(f/f_{ref})^{α}, we place 95% confidence level upper limits on the energy density of the background in each of four frequency bands spanning 41.5-1726Hz. In the frequency band of 41.5-169.25Hz for a spectral index of α=0, we constrain the energy density of the stochastic background to be Ω_{GW}(f)<5.6×10^{-6}. For the 600-1000Hz band, Ω_{GW}(f)<0.14(f/900 Hz)^{3}, a factor of 2.5 lower than the best previously reported upper limits. We find Ω_{GW}(f)<1.8×10^{-4} using a spectral index of zero for 170-600Hz and Ω_{GW}(f)<1.0(f/1300 Hz)^{3} for 1000-1726Hz, bands in which no previous direct limits have been placed. The limits in these four bands are the lowest direct measurements to date on the stochastic background. We discuss the implications of these results in light of the recent claim by the BICEP2 experiment of the possible evidence for inflationary gravitational waves.

We performed a simultaneous observational run with prototypes of Torsion-bar Antenna (TOBA) and searched for a stochastic gravitational waves (GW) background. TOBA is a new type of GW detector which measures differential rotation of two test-mass bars caused by tidal force from GWs. It fundamentally has a good sensitivity at lower frequencies, such as 0.1–1.0 Hz. The prototype has a 20-cm test mass bar which is levitated by the pinning effect of a superconductor. The data was taken from 1:00 am to 10:00 am on March 11th 2011, at Tokyo and Kyoto in Japan. As a result, we did not detect a stochastic GW background with false alarm rate of 5 %, and set an upper limit on a stochastic GW background. Our 95 % confidence upper limit is Ωgwh20 < 1.2 × 1019 at 0.06 – 0.9 Hz, where Ωgw is the GW energy density per logarithmic frequency interval in units of the critical density and ho is the Hubble constant per 100 km/sec/Mpc. We had established the simultaneous observation and the analysis pipeline with two TOBAs, and set an upper limit at a wider frequency band.

The stochastic gravitational wave (GW) background provides a fascinating window to the physics of the very early universe. Beyond the nearly scale-invariant primordial GW spectrum produced during inflation, a spectrum with a much richer structure is typically generated during the preheating phase after inflation (or after some other phase transition at lower energies). This raises the question of what one can learn from a future observation of the stochastic gravitational wave background spectrum about the underlying physics during preheating. Recently, it has been shown that during preheating non-perturbative quasi-stable objects like oscillons can act as strong sources for GW, leading to characteristic features such as distinct peaks in the spectrum. In this paper, we study the GW production from oscillons using semi-analytical techniques. In particular, we discuss how the GW spectrum is affected by the parameters that characterise a given oscillon system, e.g. by the background cosmology, the asymmetry of the oscillons and the evolution of the number density of the oscillons. We compare our semi-analytic results with numerical lattice simulations for a hilltop inflation model and a KKLT scenario, which differ strongly in some of these characteristics, and find very good agreement.

We study the co-evolution of supermassive black holes (SMBHs) with galaxies by means of semi-analytic model (SAM) of galaxy formation based on sub-halo merger trees built from Millennium and Millennium-II simulation. We utilize the simulation results from Guo 2013 and Henriques 2015 to study two aspects of the co-evolution, \emph{i.e.} the stochastic gravitational wave (GW) background generated by SMBH merger and the SMBH/galaxy clustering. The characteristic strain amplitude of GW background predicted by Guo 2013 and Henriques 2015 models are $A_{yr^{-1}}=5.00\times10^{-16}$ and $A_{yr^{-1}}=9.42\times10^{-17}$, respectively. We find the GW amplitude is very sensitive to the galaxy merger rate. The difference in the galaxy merger rate between Guo 2013 and Henriques 2015, results in a factor $5$ deviation in the GW strain amplitude. For clusterings, we calculate the spatially isotropic two point auto- and cross-correlation functions (2PCFs) for both SMBHs and galaxies by using the mock catalogs generated from Guo 2013 model. We find that all 2PCFs have positive dependence on both SMBH and galaxy mass. And there exist a significant time evolution in 2PCFs, namely, the clustering effect is enhanced at lower redshifts. Interestingly, this result is not reported in the active galactic nuclei samples in SDSS. Our analysis also shows that, roughly, SMBHs and galaxies, with galaxy mass $10^2\sim10^3$ larger than SMBH mass, have similar pattern of clustering, which is a reflection of the co-evolution of SMBH and galaxy. Finally, we calculate the first ten multiples of the angular power spectrum of the energy density of GW background. We find the amplitude of angular power spectrum of the first ten multiples is about $10\%$ to $60\%$ of the monopole component in the whole frequency range.

The observations from pulsar timing arrays (PTAs), led by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), have provided opportunities to constrain primordial gravitational waves at low frequencies. In this paper, we analyze the best-fit parameter values for different Gauss-Bonnet Inflationary Gravitational Wave (GB-IGW) models with the PTArcade program, and we compare the results with the observations of NANOGrav, European Pulsar Timing Array (EPTA), Parkes Pulsar Timing Array (PPTA), and International Pulsar Timing Array (IPTA). We find the potential parameter n derived from the GB-IGW model is not necessarily positive. Instead, PTA data suggest the possibility of new parameter ranges. Meanwhile, the other parameters' results, such as the spectral indices of the scalar and tensor perturbation also have differences with traditional Cosmic Microwave Background analyses from Planck 2018, showing the reason for the different phenomena of Gravitational Wave power spectra expected. Additionally, the GB-IGW models we calculated are closer to the central values of multiple PTA datasets compared to the standard inflation model, making the GB-IGW model more likely to be detected by future space-based gravitational wave observatories.

The North American Nanohertz Observatory for Gravitational Waves has recently reported strong evidence for a stochastic common-spectrum process affecting the pulsar timing residuals in its 12.5-year data set. We demonstrate that this process admits an interpretation in terms of a stochastic gravitational-wave background emitted by a cosmic-string network in the early Universe. We study stable Nambu-Goto strings in dependence of their tension Gμ and loop size α and show that the entire viable parameter space will be probed by an array of future experiments.

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 = 8 nHz, we placed a sky-averaged upper limit of h 0 &lt; 7.3(3) × 10−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 with = 8 nHz within 120 Mpc for , and within 5.5 Gpc for at our most sensitive sky location. We also determined that there are no SMBHBs with emitting GWs with = 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.

This paper reports a comprehensive study on the gravitational wave (GW) background from compact binary coalescences. We consider in our calculations newly available observation-based neutron star and black hole mass distributions and complete analytical waveforms that include post-Newtonian amplitude corrections. Our results show that: (i) post-Newtonian effects cause a small reduction in the GW background signal; (ii) below 100 Hz the background depends primarily on the local coalescence rate $r_0$ and the average chirp mass and is independent of the chirp mass distribution; (iii) the effects of cosmic star formation rates and delay times between the formation and merger of binaries are linear below 100 Hz and can be represented by a single parameter within a factor of ~ 2; (iv) a simple power law model of the energy density parameter $\Omega_{GW}(f) ~ f^{2/3}$ up to 50-100 Hz is sufficient to be used as a search template for ground-based interferometers. In terms of the detection prospects of the background signal, we show that: (i) detection (a signal-to-noise ratio of 3) within one year of observation by the Advanced LIGO detectors (H1-L1) requires a coalescence rate of $r_0 = 3 (0.2) Mpc^{-3} Myr^{-1}$ for binary neutron stars (binary black holes); (ii) this limit on $r_0$ could be reduced 3-fold for two co-located detectors, whereas the currently proposed worldwide network of advanced instruments gives only ~ 30% improvement in detectability; (iii) the improved sensitivity of the planned Einstein Telescope allows not only confident detection of the background but also the high frequency components of the spectrum to be measured. Finally we show that sub-threshold binary neutron star merger events produce a strong foreground, which could be an issue for future terrestrial stochastic searches of primordial GWs.

We investigate the formation and growth of massive black hole (BH) seeds in dusty star-forming galaxies, relying and extending the framework proposed by [1]. Specifically, the latter envisages the migration of stellar compact remnants (neutron stars and stellar-mass black holes) via gaseous dynamical friction towards the galaxy nuclear region, and their subsequent merging to grow a massive central BH seed. In this paper we add two relevant ingredients: (i) we include primordial BHs, that could constitute a fraction f pBH of the dark matter, as an additional component participating in the seed growth; (ii) we predict the stochastic gravitational wave background originated during the seed growth, both from stellar compact remnant and from primordial BH mergers. We find that the latter events contribute most to the initial growth of the central seed during a timescale of 106–107 yr, before stellar compact remnant mergers and gas accretion take over. In addition, if the fraction of primordial BHs f pBH is large enough, gravitational waves emitted by their mergers in the nuclear galactic regions could be detected by future interferometers like Einsten Telescope, DECIGO and LISA. As for the associated stochastic gravitational wave background, we predict that it extends over the wide frequency band 10-6 ≲ f[Hz] ≲ 10, which is very different from the typical range originated by mergers of isolated binary compact objects. On the one hand, the detection of such a background could be a smoking gun to test the proposed seed growth mechanism; on the other hand, it constitutes a relevant contaminant from astrophysical sources to be characterized and subtracted, in the challenging search for a primordial background of cosmological origin.