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.
Microstructural and geochronological analysis of shocked zircon has greatly advanced understanding the formation and evolution of impact structures. However, fundamental aspects of shock-produced ...planar microstructures in zircon remain poorly known, such as their deformation mechanisms, crystallographic orientations, and how planar microstructures visible at the grain scale by scanning electron microscopy correlate to microstructures visible at sub-micrometer scales by transmission electron microscopy and electron backscatter diffraction (EBSD). To unify observations of planar microstructures in zircon made at different scales into a consistent framework, we integrate the results of: (1) three-dimensional crystallographic modeling of planar microstructure orientations, with (2) 360° external prism backscattered electron imaging at the grain scale, and (3) polished section cathodoluminescence and EBSD analysis at the sub-micrometer scale for a suite of detrital shocked zircons eroded from the Vredefort Dome in South Africa. Our combined approach resulted in the documentation of seven planar microstructure orientations that can be correlated from grain to sub-micrometer scales of observation: (010), (100), (112), (112), (112), (112), and (011). All orientations of planar microstructures exhibit minor variations in style, however all are considered to be fractures; no amorphous ZrSiO4 lamellae were identified. We therefore favor the usage of "planar fracture" (PF) over "planar deformation feature" (PDF) for describing the observed planar microstructures in zircon based broadly on the nomenclature developed for shocked quartz. Some {112} PFs visible at the grain scale contain impact microtwins detectable by EBSD, and are the first report of polysynthetic twinning in zircon. The microtwins consist of parallel sets of thin lamellae of zircon oriented 65° about and occur in multiple crosscutting {112} orientations within single grains. Curviplanar fractures and injected melt are additional impact-related microstructures associated with PF formation. Crosscutting relations of shock microstructures reveal the following chronology: (1) Early development of c-axis parallel PFs in (010) and (100) orientations; (2) the development of up to four {112} PFs, including some with microtwins; (3) the development of curviplanar fractures and the injection of impact derived melt; (4) the development of (011) PFs associated with compressional deformation; and (5) grain-scale non-discrete crystal plastic deformation. Experimental constraints for the onset of PFs, together with the absence of reidite, suggest formation conditions from 20 to 40 GPa for all of the planar microstructures described here.
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 report results from continued timing observations of PSR J0740+6620, a high-mass, 2.8 ms radio pulsar in orbit with a likely ultracool white dwarf companion. Our data set consists of ...combined pulse arrival-time measurements made with the 100 m Green Bank Telescope and the Canadian Hydrogen Intensity Mapping Experiment telescope. We explore the significance of timing-based phenomena arising from general relativistic dynamics and variations in pulse dispersion. When using various statistical methods, we find that combining ∼1.5 yr of additional, high-cadence timing data with previous measurements confirms and improves on previous estimates of relativistic effects within the PSR J0740+6620 system, with the pulsar mass
m
p
=
2.08
−
0.07
+
0.07
M
⊙
(68.3% credibility) determined by the relativistic Shapiro time delay. For the first time, we measure secular variation in the orbital period and argue that this effect arises from apparent acceleration due to significant transverse motion. After incorporating contributions from Galactic differential rotation and off-plane acceleration in the Galactic potential, we obtain a model-dependent distance of
d
=
1.14
−
0.15
+
0.17
kpc (68.3% credibility). This improved distance confirms the ultracool nature of the white dwarf companion determined from recent optical observations. We discuss the prospects for future observations with next-generation facilities, which will likely improve the precision on
m
p
for J0740+6620 by an order of magnitude within the next few years.
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.
Since the discovery of superconductivity, there has been a drive to understand the mechanisms by which it occurs. The BCS (Bardeen-Cooper-Schrieffer) model successfully treats the electrons in ...conventional superconductors as pairs coupled by phonons (vibrational modes of oscillation) moving through the material, but there is as yet no accepted model for high-transition-temperature, organic or 'heavy fermion' superconductivity. Experiments that reveal unusual properties of those superconductors could therefore point the way to a deeper understanding of the underlying physics. In particular, the response of a material to a magnetic field can be revealing, because this usually reduces or quenches superconductivity. Here we report measurements of the heat capacity and magnetization that show that, for particular orientations of an external magnetic field, superconductivity in the heavy-fermion material CeCoIn5 is enhanced through the magnetic moments (spins) of individual electrons. This enhancement occurs by fundamentally altering how the superconducting state forms, resulting in regions of superconductivity alternating with walls of spin-polarized unpaired electrons; this configuration lowers the free energy and allows superconductivity to remain stable. The large magnetic susceptibility of this material leads to an unusually strong coupling of the field to the electron spins, which dominates over the coupling to the electron orbits.
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.
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.