Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNoise-marginalized optimal statistic: A robust hybrid frequentist-Bayesian statistic for the stochastic gravitational-wave background in pulsar timing arrays
83
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Observations have revealed that nearly all galaxies contain supermassive black holes (SMBHs) at their centers. When galaxies merge, these SMBHs form SMBH binaries (SMBHBs) that emit low-frequency gravitational waves (GWs). The incoherent superposition of these sources produce a stochastic GW background (GWB) that can be observed by pulsar timing arrays (PTAs). The optimal statistic is a frequentist estimator of the amplitude of the GWB that specifically looks for the spatial correlations between pulsars induced by the GWB. In this paper, we introduce an improved method for computing the optimal statistic that marginalizes over the red noise in individual pulsars. We use simulations to demonstrate that this method more accurately determines the strength of the GWB, and we use the noise-marginalized optimal statistic to compare the significance of monopole, dipole, and Hellings-Downs (HD) spatial correlations and perform sky scrambles.
rate research
Read More
Supermassive black hole binaries, cosmic strings, relic gravitational waves from inflation, and first order phase transitions in the early universe are expected to contribute to a stochastic background of gravitational waves in the 10^(-9) Hz-10^(-7) Hz frequency band. Pulsar timing arrays (PTAs) exploit the high precision timing of radio pulsars to detect signals at such frequencies. Here we present a time-domain implementation of the optimal cross-correlation statistic for stochastic background searches in PTA data. Due to the irregular sampling typical of PTA data as well as the use of a timing model to predict the times-of-arrival of radio pulses, time-domain methods are better suited for gravitational wave data analysis of such data. We present a derivation of the optimal cross-correlation statistic starting from the likelihood function, a method to produce simulated stochastic background signals, and a rigorous derivation of the scaling laws for the signal-to-noise ratio of the cross-correlation statistic in the two relevant PTA regimes: the weak signal limit where instrumental noise dominates over the gravitational wave signal at all frequencies, and a second regime where the gravitational wave signal dominates at the lowest frequencies.
We have begun an exciting era for gravitational wave detection, as several world-leading experiments are breaching the threshold of anticipated signal strengths. Pulsar timing arrays (PTAs) are pan-Galactic gravitational wave detectors that are already cutting into the expected strength of gravitational waves from cosmic strings and binary supermassive black holes in the nHz-$mu$Hz gravitational wave band. These limits are leading to constraints on the evolutionary state of the Universe. Here, we provide a broad review of this field, from how pulsars are used as tools for detection, to astrophysical sources of uncertainty in the signals PTAs aim to see, to the primary current challenge areas for PTA work. This review aims to provide an up-to-date reference point for new parties interested in the field of gravitational wave detection via pulsar timing.
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 Earths 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 (Jupiters orbital elements) most relevant to gravitational-wave searches with pulsar timing arrays. This model, BayesEphem, was used to derive ephemeris-robust results in NANOGravs 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.
We derive scaling laws for the signal-to-noise ratio of the optimal cross-correlation statistic, and show that the large power-law increase of the signal-to-noise ratio as a function of the the observation time $T$ that is usually assumed holds only at early times. After enough time has elapsed, pulsar timing arrays enter a new regime where the signal to noise only scales as $sqrt{T}$. In addition, in this regime the quality of the pulsar timing data and the cadence become relatively un-important. This occurs because the lowest frequencies of the pulsar timing residuals become gravitational-wave dominated. Pulsar timing arrays enter this regime more quickly than one might naively suspect. For T=10 yr observations and typical stochastic background amplitudes, pulsars with residual RMSs of less than about $1,mu$s are already in that regime. The best strategy to increase the detectability of the background in this regime is to increase the number of pulsars in the array. We also perform realistic simulations of the NANOGrav pulsar timing array, which through an aggressive pulsar survey campaign adds new millisecond pulsars regularly to its array, and show that a detection is possible within a decade, and could occur as early as 2016.
We search for an isotropic stochastic gravitational-wave background (GWB) in the newly released $11$-year dataset from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). While we find no significant evidence for a GWB, we place constraints on a GWB from a population of supermassive black-hole binaries, cosmic strings, and a primordial GWB. For the first time, we find that the GWB upper limits and detection statistics 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 PTA constraints that are robust against SSE uncertainties. We thus place a $95%$ upper limit on the GW strain amplitude of $A_mathrm{GWB}<1.45times 10^{-15}$ at a frequency of $f=1$ yr$^{-1}$ for a fiducial $f^{-2/3}$ power-law spectrum, and with inter-pulsar correlations modeled. This is a factor of $sim 2$ improvement over the NANOGrav $9$-year limit, calculated using the same procedure. Previous PTA upper limits on the GWB will need revision in light of SSE systematic uncertainties. 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 cosmic-string tension using recent simulations, yielding an SSE-marginalized $95%$ upper limit on the cosmic string tension of $Gmu < 5.3times 10^{-11}$---a factor of $sim 2$ better than the published NANOGrav $9$-year constraints. Our SSE-marginalized $95%$ upper limit on the energy density of a primordial GWB (for a radiation-dominated post-inflation Universe) is $Omega_mathrm{GWB}(f)h^2<3.4times10^{-10}$.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
