Search for: All records

Abstract Pulsar timing array (PTA) searches for gravitational waves (GWs) aim to detect a characteristic correlation pattern in the timing residuals of galactic millisecond pulsars. This pattern is described by the PTA overlap reduction function (ORF) ${\mathrm{\Gamma }}_{ab}\left({\xi }_{ab}\right)$, which is known as the Hellings–Downs (HD) curve in general relativity (GR). In theories of modified gravity, the HD curve often receives corrections. Assuming, e.g. a subluminal GW phase velocity, one finds a drastically enhanced ORF in the limit of small angular separations between pulsaraand pulsarbin the sky, ${\xi }_{ab}\to 0$. In particular, working in harmonic space and performing an approximate resummation of all multipole contributions, the auto correlation coefficientΓ_aaseems to diverge. In this paper, we confirm that this divergence is unphysical and provide an exact and analytical expression forΓ_aain dependence of the pulsar distanceL_aand the GW phase velocity $v_{\mathrm{ph}}$. In the GR limit and assuming a large pulsar distance, our expression reduces to ${\mathrm{\Gamma }}_{aa}=1$. In the case of subluminal phase velocity, we show that the regularization of the naive divergent result is a finite-distance effect, meaning thatΓ_aascales linearly withfL_a, wherefis the GW frequency. For superluminal phase velocity (subluminal group velocity), which is relevant in the case of massive gravity, we correct an earlier analytical result forΓ_ab. Our results pave the way for fitting modified-gravity theories with nonstandard phase velocity to PTA data, which requires a proper understanding of the auto correlation coefficientΓ_aa. 

Metric perturbations induced by ultralight dark matter (ULDM) fields have long been identified as a potential target for pulsar timing array (PTA) observations. Previous works have focused on the coherent oscillation of metric perturbations at the characteristic frequency set by the ULDM mass. In this work, we show that ULDM fields source low-frequency stochastic metric fluctuations and that these low-frequency fluctuations can produce distinctive detectable signals in PTA data. Using the NANOGrav 12.5-yr dataset and synthetic datasets mimicking present and future PTA capabilities, we show that the current and future PTA observations provide the strongest probe of ULDM density within the Solar System for masses in the range of ${10}^{-18}\mathrm{eV}-{10}^{-16}\mathrm{eV}$. Published by the American Physical Society2024 

Abstract Pulsar timing arrays have found evidence for a low-frequency gravitational-wave background (GWB). Assuming that the GWB is produced by supermassive black hole binaries (SMBHBs), the next gravitational-wave (GW) signals astronomers anticipate are continuous waves (CWs) from single SMBHBs and their associated GWB anisotropy. The prospects for detecting CWs and anisotropy are highly dependent on the astrophysics of SMBHB populations. Thus, information from single sources can break degeneracies in astrophysical models and place much more stringent constraints than the GWB alone. We simulate and evolve SMBHB populations, model their GWs, and calculate their anisotropy and detectability. We investigate how varying components of our semianalytic model, including the galaxy stellar mass function, the SMBH–host galaxy relation (M_BH–M_bulge), and the binary evolution prescription, impact the expected detections. The CW occurrence rate is greatest for few total binaries, high SMBHB masses, large scatter inM_BH–M_bulge, and long hardening times. The occurrence rate depends most on the binary evolution parameters, implying that CWs offer a novel avenue to probe binary evolution. The most detectable CW sources are in the lowest frequency bin for a 16.03 yr PTA, have masses from ∼10⁹to 10¹⁰M_⊙, and are ∼1 Gpc away. The level of anisotropy increases with frequency, with the angular power spectrum over multipole modesℓvarying in low-frequencyC_ℓ>0/C₀from ∼5 × 10⁻³to ∼2 × 10⁻¹, depending on the model; typical values are near current upper limits. Observing this anisotropy would support SMBHB models for the GWB over cosmological models, which tend to be isotropic. 

Abstract 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. 

Abstract Evidence has emerged for a stochastic signal correlated among 67 pulsars within the 15 yr pulsar-timing data set compiled by the NANOGrav collaboration. Similar signals have been found in data from the European, Indian, Parkes, and Chinese pulsar timing arrays. This signal has been interpreted as indicative of the presence of a nanohertz stochastic gravitational-wave background (GWB). To explore the internal consistency of this result, we investigate how the recovered signal strength changes as we remove the pulsars one by one from the data set. We calculate the signal strength using the (noise-marginalized) optimal statistic, a frequentist metric designed to measure the correlated excess power in the residuals of the arrival times of the radio pulses. We identify several features emerging from this analysis that were initially unexpected. The significance of these features, however, can only be assessed by comparing the real data to synthetic data sets. After conducting identical analyses on simulated data sets, we do not find anything inconsistent with the presence of a stochastic GWB in the NANOGrav 15 yr data. The methodologies developed here can offer additional tools for application to future, more sensitive data sets. While this analysis provides an internal consistency check of the NANOGrav results, it does not eliminate the necessity for additional investigations that could identify potential systematics or uncover unmodeled physical phenomena in the data. 

Abstract The cosmic merger history of supermassive black hole binaries (SMBHBs) is expected to produce a low-frequency gravitational wave background (GWB). Here we investigate how signs of the discrete nature of this GWB can manifest in pulsar timing arrays (PTAs) through excursions from, and breaks in, the expected $f_{\mathrm{GW}}^{-2/3}$power law of the GWB strain spectrum. To do this, we create a semianalytic SMBHB population model, fit to North American Nanohertz Observatory for Gravitational Waves (NANOGrav’s) 15 yr GWB amplitude, and with 1000 realizations, we study the populations’ characteristic strain and residual spectra. Comparing our models to the NANOGrav 15 yr spectrum, we find two interesting excursions from the power law. The first, at 2 nHz, is below our GWB realizations with ap-value significancep= 0.05–0.06 (≈1.8σ–1.9σ). The second, at 16 nHz, is above our GWB realizations withp= 0.04–0.15 (≈1.4σ–2.1σ). We explore the properties of a loud SMBHB that could cause such an excursion. Our simulations also show that the expected number of SMBHBs decreases by 3 orders of magnitude, from ∼10⁶to ∼10³, between 2 and 20 nHz. This causes a break in the strain spectrum as the stochasticity of the background breaks down at ${26}_{-19}^{+28}\mathrm{nHz}$, consistent with predictions pre-dating GWB measurements. The diminished GWB signal from SMBHBs at frequencies above the 26 nHz break opens a window for PTAs to detect continuous GWs from individual SMBHBs or GWs from the early Universe.