This case report concerns a woman with intellectual disability, schizoaffective disorder, and avoidance and refusal of having her blood drawn. She required but refused blood draws to properly monitor ...the therapeutic dose of a necessary psychotropic medication. During intervention at a community-based habilitation setting and under simulated conditions, direct-care and supervisory staff (a) gradually exposed her to steps constituting a blood draw, (b) reinforced completed steps (praise and tokens), and (c) offered a monetary incentive to participate in an actual blood draw. Evaluated in a changing criterion design, the woman demonstrated increased compliance and eventually completed blood draws immediately following intervention, months later, and at 12- to 24-month follow-ups. We discuss implications from the case and the advantages of in vivo intervention for overcoming medical fears among people who have intellectual disability.
PINT is a pure-Python framework for high-precision pulsar timing developed on top of widely used and well-tested Python libraries, supporting both interactive and programmatic data analysis ...workflows. We present a new frequentist framework within PINT to characterize the single-pulsar noise processes present in pulsar timing data sets. This framework enables parameter estimation for both uncorrelated and correlated noise processes, as well as model comparison between different timing and noise models in a computationally inexpensive way. We demonstrate the efficacy of the new framework by applying it to simulated data sets as well as a real data set of PSR B1855+09. We also describe the new features implemented in PINT since it was first described in the literature.
Abstract
We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15 yr pulsar timing data set collected by the North American Nanohertz Observatory ...for Gravitational Waves. The correlations follow the Hellings–Downs pattern expected for a stochastic gravitational-wave background. The presence of such a gravitational-wave background with a power-law spectrum is favored over a model with only independent pulsar noises with a Bayes factor in excess of 10
14
, and this same model is favored over an uncorrelated common power-law spectrum model with Bayes factors of 200–1000, depending on spectral modeling choices. We have built a statistical background distribution for the latter Bayes factors using a method that removes interpulsar correlations from our data set, finding
p
= 10
−3
(≈3
σ
) for the observed Bayes factors in the null no-correlation scenario. A frequentist test statistic built directly as a weighted sum of interpulsar correlations yields
p
= 5 × 10
−5
to 1.9 × 10
−4
(≈3.5
σ
–4
σ
). Assuming a fiducial
f
−2/3
characteristic strain spectrum, as appropriate for an ensemble of binary supermassive black hole inspirals, the strain amplitude is
2.4
−
0.6
+
0.7
×
10
−
15
(median + 90% credible interval) at a reference frequency of 1 yr
−1
. The inferred gravitational-wave background amplitude and spectrum are consistent with astrophysical expectations for a signal from a population of supermassive black hole binaries, although more exotic cosmological and astrophysical sources cannot be excluded. The observation of Hellings–Downs correlations points to the gravitational-wave origin of this signal.
Abstract
The NANOGrav 15 yr data set shows evidence for the presence of a low-frequency gravitational-wave background (GWB). While many physical processes can source such low-frequency gravitational ...waves, here we analyze the signal as coming from a population of supermassive black hole (SMBH) binaries distributed throughout the Universe. We show that astrophysically motivated models of SMBH binary populations are able to reproduce both the amplitude and shape of the observed low-frequency gravitational-wave spectrum. While multiple model variations are able to reproduce the GWB spectrum at our current measurement precision, our results highlight the importance of accurately modeling binary evolution for producing realistic GWB spectra. Additionally, while reasonable parameters are able to reproduce the 15 yr observations, the implied GWB amplitude necessitates either a large number of parameters to be at the edges of expected values or a small number of parameters to be notably different from standard expectations. While we are not yet able to definitively establish the origin of the inferred GWB signal, the consistency of the signal with astrophysical expectations offers a tantalizing prospect for confirming that SMBH binaries are able to form, reach subparsec separations, and eventually coalesce. As the significance grows over time, higher-order features of the GWB spectrum will definitively determine the nature of the GWB and allow for novel constraints on SMBH populations.
Abstract
We search for an isotropic stochastic gravitational-wave background (GWB) in the 12.5 yr pulsar-timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. ...Our analysis finds strong evidence of a stochastic process, modeled as a power law, with common amplitude and spectral slope across pulsars. Under our fiducial model, the Bayesian posterior of the amplitude for an
f
−2/3
power-law spectrum, expressed as the characteristic GW strain, has median 1.92 × 10
−15
and 5%–95% quantiles of 1.37–2.67 × 10
−15
at a reference frequency of
f
yr
=
1
yr
−
1
;
the Bayes factor in favor of the common-spectrum process versus independent red-noise processes in each pulsar exceeds 10,000. However, we find no statistically significant evidence that this process has quadrupolar spatial correlations, which we would consider necessary to claim a GWB detection consistent with general relativity. We find that the process has neither monopolar nor dipolar correlations, which may arise from, for example, reference clock or solar system ephemeris systematics, respectively. The amplitude posterior has significant support above previously reported upper limits; we explain this in terms of the Bayesian priors assumed for intrinsic pulsar red noise. We examine potential implications for the supermassive black hole binary population under the hypothesis that the signal is indeed astrophysical in nature.
Abstract
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.
Abstract
Pulsar timing arrays (PTAs) are galactic-scale gravitational wave (GW) detectors. Each individual arm, composed of a millisecond pulsar, a radio telescope, and a kiloparsecs-long path, ...differs in its properties but, in aggregate, can be used to extract low-frequency GW signals. We present a noise and sensitivity analysis to accompany the NANOGrav 15 yr data release and associated papers, along with an in-depth introduction to PTA noise models. As a first step in our analysis, we characterize each individual pulsar data set with three types of white-noise parameters and two red-noise parameters. These parameters, along with the timing model and, particularly, a piecewise-constant model for the time-variable dispersion measure, determine the sensitivity curve over the low-frequency GW band we are searching. We tabulate information for all of the pulsars in this data release and present some representative sensitivity curves. We then combine the individual pulsar sensitivities using a signal-to-noise ratio statistic to calculate the global sensitivity of the PTA to a stochastic background of GWs, obtaining a minimum noise characteristic strain of 7 × 10
−15
at 5 nHz. A power-law-integrated analysis shows rough agreement with the amplitudes recovered in NANOGrav’s 15 yr GW background analysis. While our phenomenological noise model does not model all known physical effects explicitly, it provides an accurate characterization of the noise in the data while preserving sensitivity to multiple classes of GW signals.
The NANOGrav Collaboration reported strong Bayesian evidence for a common-spectrum stochastic process in its12.5 yr pulsar timing array data set, with median characteristic strain amplitude at ...periods of a year of A(yr) = 1.92(+0.75,-0.55) x 10^(-15). However, evidence for the quadrupolar Hellings & Downs interpulsar correlations, which are characteristic of gravitational-wave signals, was not yet significant. We emulate and extend the NANOGrav data set, injecting a wide range of stochastic gravitational-wave background(GWB)signals that encompass a variety of amplitudes and spectral shapes, and quantify three key milestones.(I)Given the amplitude measured in the 12.5 yr analysis and assuming this signal is a GWB, we expect to accumulate robust evidence of an interpulsar-correlated GWB signal with 15–17 yr of data, i.e., an additional 2–5 yr from the 12.5 yr data set.(II)At the initial detection, we expect a fractional uncertainty of 40% on the power-law strain spectrum slope, which is sufficient to distinguish a GWB of supermassive black hole binary origin from some models predicting more exotic origins.(III)Similarly, the measured GWB amplitude will have an uncertainty of 44% upon initial detection, allowing us to arbitrate between some population models of supermassive black hole binaries. In addition, power-law models are distinguishable from those having low-frequency spectral turnovers once 20 yr of data are reached. Even though our study is based on the NANOGrav data, we also derive relations that allow for a generalization to other pulsar timing array data sets. Most notably, by combining the data of individual arrays into the International Pulsar Timing Array, all of these milestones can be reached significantly earlier.
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.
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 (
C
l
=1
/
C
l
=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.