Showing 1–50 of 145 results for author: Crawford, F

Detection of 49 Weak Dispersed Radio Pulses in a Parkes Observation of the X-ray Pulsar PSR J0537$-$6910

Authors: Fronefield Crawford

Abstract: I conducted a new search for dispersed radio pulses from the X-ray pulsar PSR J0537$-$6910 in the Large Magellanic Cloud in a long (11.6 hr) archival 1.4 GHz Parkes search observation. I searched dispersion measures (DMs) between 0 and 10000 pc cm$^{-3}$ and detected 49 pulses with a signal-to-noise ratio (S/N) greater than 7 at a wide range of DMs using the HEIMDALL and FETCH pulse detection and… ▽ More I conducted a new search for dispersed radio pulses from the X-ray pulsar PSR J0537$-$6910 in the Large Magellanic Cloud in a long (11.6 hr) archival 1.4 GHz Parkes search observation. I searched dispersion measures (DMs) between 0 and 10000 pc cm$^{-3}$ and detected 49 pulses with a signal-to-noise ratio (S/N) greater than 7 at a wide range of DMs using the HEIMDALL and FETCH pulse detection and classification packages. All of the pulses were weak, with none having a S/N above 8.5. There was a significant excess of pulses observed in the DM range of the known pulsar population in the LMC, suggesting that these pulses may originate from LMC pulsars. Three repeat pulses, each having widths $\lesssim 1$ ms, were detected in a single DM trial of 103.412 pc cm$^{-3}$, which is in the LMC DM range. This is unlikely to occur by chance in a single DM trial in this search at the (marginally significant) 4.3$σ$ level. It remains unclear whether any of the detected pulses in the sample are from PSR J0537$-$6910 itself. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: Accepted by the Astrophysical Journal

The NANOGrav 15 yr Data Set: Looking for Signs of Discreteness in the Gravitational-wave Background

Authors: Gabriella Agazie, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Lucas Brown, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Heling Deng, Timothy Dolch, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Gabriel E. Freedman, Nate Garver-Daniels , et al. (58 additional authors not shown)

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 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 semi-analyt… ▽ More 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 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 semi-analytic SMBHB population model, fit to NANOGrav's 15 yr GWB amplitude, and with 1,000 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 \; \mathrm{nHz}$, is below our GWB realizations with $p$-value significance $p = 0.05$ to $0.06$ ($\approx 1.8 σ- 1.9 σ$). The second, at $16 \; \mathrm{nHz}$, is above our GWB realizations with $p = 0.04$ to $0.15$ ($\approx 1.4 σ- 2.1 σ$). We explore the properties of a loud SMBHB which could cause such an excursion. Our simulations also show that the expected number of SMBHBs decreases by three orders of magnitude, from $\sim 10^6$ to $\sim 10^3$, between $2\; \mathrm{nHz}$ and $20 \; \mathrm{nHz}$. This causes a break in the strain spectrum as the stochasticity of the background breaks down at $26^{+28}_{-19} \; \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. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 10 pages, 8 figures, 1 appendix, submitted to ApJ

A search for millisecond radio bursts with Messier 82

Authors: S. Paine, T. Hawkins, D. R. Lorimer, J. Stanley, J. Kania, F. Crawford, N. Fairfield

Abstract: Fast radio bursts (FRBs) are short-duration radio pulses of cosmological origin. Among the most common sources predicted to explain this phenomenon are bright pulses from a class of extremely highly magnetized neutron stars known as magnetars. Motivated by the discovery of an FRB-like pulse from the Galactic magnetar SGR 1935+2154, we searched for similar events in Messier 82 (M82). With a star fo… ▽ More Fast radio bursts (FRBs) are short-duration radio pulses of cosmological origin. Among the most common sources predicted to explain this phenomenon are bright pulses from a class of extremely highly magnetized neutron stars known as magnetars. Motivated by the discovery of an FRB-like pulse from the Galactic magnetar SGR 1935+2154, we searched for similar events in Messier 82 (M82). With a star formation rate 40 times that of the Milky Way, one might expect that the implied rate of events similar to that seen from SGR 1935+2154 from M82 should be 40 times higher than that of the Milky Way. We observed M82 at 1.4 GHz with the 20-m telescope at the Green Bank Observatory for 34.8 days. While we found many candidate events, none had a signal-to-noise ratio greater than 8. We also show that there are insufficient numbers of repeating low-significance events at similar dispersion measures to constitute a statistically significant detection. From these results, we place an upper bound for the rate of radio pulses from M82 to be 30 per year above a fluence limit of 8.5 Jy ms. While this is less than 9 times the rate of radio bursts from magnetars in the Milky Way inferred from the previous radio detections of SGR 1935+2154, it is possible that propagation effects from interstellar scattering are currently limiting our ability to detect sources in M82. Further searches of M82 and other nearby galaxies are encouraged to probe this putative FRB population. △ Less

Submitted 30 January, 2024; originally announced January 2024.

No Dispersed Single Radio Pulses Detected in Archival Parkes Pulsar Observations Targeting Supernova Remnants and Anomalous X-ray Pulsars

Authors: Fronefield Crawford

Abstract: Four supernova remnants and four anomalous X-ray pulsars were previously observed with the Parkes telescope in a campaign to detect pulsed radio emission from associated neutron stars. No signals were detected in the original searches of these data. I have reprocessed the data with the more recently developed HEIMDALL and FETCH software packages, which are optimized for single-pulse detection and… ▽ More Four supernova remnants and four anomalous X-ray pulsars were previously observed with the Parkes telescope in a campaign to detect pulsed radio emission from associated neutron stars. No signals were detected in the original searches of these data. I have reprocessed the data with the more recently developed HEIMDALL and FETCH software packages, which are optimized for single-pulse detection and classification. In this new analysis, no astrophysical pulses were detected having a signal-to-noise ratio above 7 from any of the targets at dispersion measures ranging from 0 to $10^{4}$ pc cm$^{-3}$. I include calculated fluence limits on single radio pulses from these targets. △ Less

Submitted 2 November, 2023; originally announced November 2023.

The NANOGrav 15-year data set: Search for Transverse Polarization Modes in the Gravitational-Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Jeremy Baier, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Dallas DeGan, Paul B. Demorest , et al. (74 additional authors not shown)

Abstract: Recently we found compelling evidence for a gravitational wave background with Hellings and Downs (HD) correlations in our 15-year data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes which produce different interpulsar correl… ▽ More Recently we found compelling evidence for a gravitational wave background with Hellings and Downs (HD) correlations in our 15-year data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes which produce different interpulsar correlations. In this work we search the NANOGrav 15-year data set for evidence of a gravitational wave background with quadrupolar Hellings and Downs (HD) and Scalar Transverse (ST) correlations. We find that HD correlations are the best fit to the data, and no significant evidence in favor of ST correlations. While Bayes factors show strong evidence for a correlated signal, the data does not strongly prefer either correlation signature, with Bayes factors $\sim 2$ when comparing HD to ST correlations, and $\sim 1$ for HD plus ST correlations to HD correlations alone. However, when modeled alongside HD correlations, the amplitude and spectral index posteriors for ST correlations are uninformative, with the HD process accounting for the vast majority of the total signal. Using the optimal statistic, a frequentist technique that focuses on the pulsar-pair cross-correlations, we find median signal-to-noise-ratios of 5.0 for HD and 4.6 for ST correlations when fit for separately, and median signal-to-noise-ratios of 3.5 for HD and 3.0 for ST correlations when fit for simultaneously. While the signal-to-noise-ratios for each of the correlations are comparable, the estimated amplitude and spectral index for HD are a significantly better fit to the total signal, in agreement with our Bayesian analysis. △ Less

Submitted 18 October, 2023; originally announced October 2023.

Comments: 11 pages, 5 figures

The NANOGrav 12.5-year data set: A computationally efficient eccentric binary search pipeline and constraints on an eccentric supermassive binary candidate in 3C 66B

Authors: Gabriella Agazie, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Harsha Blumer, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Belinda D. Cheeseboro, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Lankeswar Dey, Timothy Dolch, Justin A. Ellis, Robert D. Ferdman, Elizabeth C. Ferrara , et al. (63 additional authors not shown)

Abstract: The radio galaxy 3C 66B has been hypothesized to host a supermassive black hole binary (SMBHB) at its center based on electromagnetic observations. Its apparent 1.05-year period and low redshift ($\sim0.02$) make it an interesting testbed to search for low-frequency gravitational waves (GWs) using Pulsar Timing Array (PTA) experiments. This source has been subjected to multiple searches for contin… ▽ More The radio galaxy 3C 66B has been hypothesized to host a supermassive black hole binary (SMBHB) at its center based on electromagnetic observations. Its apparent 1.05-year period and low redshift ($\sim0.02$) make it an interesting testbed to search for low-frequency gravitational waves (GWs) using Pulsar Timing Array (PTA) experiments. This source has been subjected to multiple searches for continuous GWs from a circular SMBHB, resulting in progressively more stringent constraints on its GW amplitude and chirp mass. In this paper, we develop a pipeline for performing Bayesian targeted searches for eccentric SMBHBs in PTA data sets, and test its efficacy by applying it on simulated data sets with varying injected signal strengths. We also search for a realistic eccentric SMBHB source in 3C 66B using the NANOGrav 12.5-year data set employing PTA signal models containing Earth term-only as well as Earth+Pulsar term contributions using this pipeline. Due to limitations in our PTA signal model, we get meaningful results only when the initial eccentricity $e_0<0.5$ and the symmetric mass ratio $η>0.1$. We find no evidence for an eccentric SMBHB signal in our data, and therefore place 95% upper limits on the PTA signal amplitude of $88.1\pm3.7$ ns for the Earth term-only and $81.74\pm0.86$ ns for the Earth+Pulsar term searches for $e_0<0.5$ and $η>0.1$. Similar 95% upper limits on the chirp mass are $(1.98 \pm 0.05) \times 10^9\,M_{\odot}$ and $(1.81 \pm 0.01) \times 10^9\,M_{\odot}$. These upper limits, while less stringent than those calculated from a circular binary search in the NANOGrav 12.5-year data set, are consistent with the SMBHB model of 3C 66B developed from electromagnetic observations. △ Less

Submitted 15 January, 2024; v1 submitted 29 September, 2023; originally announced September 2023.

Comments: 27 Pages, 10 Figures, 1 Table, Accepted for publication in ApJ

How to Detect an Astrophysical Nanohertz Gravitational-Wave Background

Authors: Bence Bécsy, Neil J. Cornish, Patrick M. Meyers, Luke Zoltan Kelley, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Tyler Cohen, James M. Cordes, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch , et al. (71 additional authors not shown)

Abstract: Analysis of pulsar timing data have provided evidence for a stochastic gravitational wave background in the nHz frequency band. The most plausible source of such a background is the superposition of signals from millions of supermassive black hole binaries. The standard statistical techniques used to search for such a background and assess its significance make several simplifying assumptions, nam… ▽ More Analysis of pulsar timing data have provided evidence for a stochastic gravitational wave background in the nHz frequency band. The most plausible source of such a background is the superposition of signals from millions of supermassive black hole binaries. The standard statistical techniques used to search for such a background and assess its significance make several simplifying assumptions, namely: i) Gaussianity; ii) isotropy; and most often iii) a power-law spectrum. However, a stochastic background from a finite collection of binaries does not exactly satisfy any of these assumptions. To understand the effect of these assumptions, we test standard analysis techniques on a large collection of realistic simulated datasets. The dataset length, observing schedule, and noise levels were chosen to emulate the NANOGrav 15-year dataset. Simulated signals from millions of binaries drawn from models based on the Illustris cosmological hydrodynamical simulation were added to the data. We find that the standard statistical methods perform remarkably well on these simulated datasets, despite their fundamental assumptions not being strictly met. They are able to achieve a confident detection of the background. However, even for a fixed set of astrophysical parameters, different realizations of the universe result in a large variance in the significance and recovered parameters of the background. We also find that the presence of loud individual binaries can bias the spectral recovery of the background if we do not account for them. △ Less

Submitted 1 December, 2023; v1 submitted 8 September, 2023; originally announced September 2023.

Comments: 14 pages, 8 figures, version matching published paper

Journal ref: ApJ 959 9 (2023)

Comparing recent PTA results on the nanohertz stochastic gravitational wave background

Authors: The International Pulsar Timing Array Collaboration, G. Agazie, J. Antoniadis, A. Anumarlapudi, A. M. Archibald, P. Arumugam, S. Arumugam, Z. Arzoumanian, J. Askew, S. Babak, M. Bagchi, M. Bailes, A. -S. Bak Nielsen, P. T. Baker, C. G. Bassa, A. Bathula, B. Bécsy, A. Berthereau, N. D. R. Bhat, L. Blecha, M. Bonetti, E. Bortolas, A. Brazier, P. R. Brook, M. Burgay , et al. (220 additional authors not shown)

Abstract: The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTA… ▽ More The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTAs that constitute the International Pulsar Timing Array (IPTA). We show that despite making different modeling choices, there is no significant difference in the GWB parameters that are measured by the different PTAs, agreeing within $1σ$. The pulsar noise parameters are also consistent between different PTAs for the majority of the pulsars included in these analyses. We bridge the differences in modeling choices by adopting a standardized noise model for all pulsars and PTAs, finding that under this model there is a reduction in the tension in the pulsar noise parameters. As part of this reanalysis, we "extended" each PTA's data set by adding extra pulsars that were not timed by that PTA. Under these extensions, we find better constraints on the GWB amplitude and a higher signal-to-noise ratio for the Hellings and Downs correlations. These extensions serve as a prelude to the benefits offered by a full combination of data across all pulsars in the IPTA, i.e., the IPTA's Data Release 3, which will involve not just adding in additional pulsars, but also including data from all three PTAs where any given pulsar is timed by more than as single PTA. △ Less

Submitted 1 September, 2023; originally announced September 2023.

Comments: 21 pages, 9 figures, submitted to ApJ

The NANOGrav 12.5-year Data Set: Search for Gravitational Wave Memory

Authors: Gabriella Agazie, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Harsha Blumer, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Dallas DeGan, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Justin A. Ellis , et al. (65 additional authors not shown)

Abstract: We present the results of a Bayesian search for gravitational wave (GW) memory in the NANOGrav 12.5-yr data set. We find no convincing evidence for any gravitational wave memory signals in this data set (Bayes factor = 2.8). As such, we go on to place upper limits on the strain amplitude of GW memory events as a function of sky location and event epoch. These upper limits are computed using a sign… ▽ More We present the results of a Bayesian search for gravitational wave (GW) memory in the NANOGrav 12.5-yr data set. We find no convincing evidence for any gravitational wave memory signals in this data set (Bayes factor = 2.8). As such, we go on to place upper limits on the strain amplitude of GW memory events as a function of sky location and event epoch. These upper limits are computed using a signal model that assumes the existence of a common, spatially uncorrelated red noise in addition to a GW memory signal. The median strain upper limit as a function of sky position is approximately $3.3 \times 10^{-14}$. We also find that there are some differences in the upper limits as a function of sky position centered around PSR J0613$-$0200. This suggests that this pulsar has some excess noise which can be confounded with GW memory. Finally, the upper limits as a function of burst epoch continue to improve at later epochs. This improvement is attributable to the continued growth of the pulsar timing array. △ Less

Submitted 25 July, 2023; originally announced July 2023.

Comments: 29 pages, 5 figures

The NANOGrav 15-year Gravitational-Wave Background Analysis Pipeline

Authors: Aaron D. Johnson, Patrick M. Meyers, Paul T. Baker, Neil J. Cornish, Jeffrey S. Hazboun, Tyson B. Littenberg, Joseph D. Romano, Stephen R. Taylor, Michele Vallisneri, Sarah J. Vigeland, Ken D. Olum, Xavier Siemens, Justin A. Ellis, Rutger van Haasteren, Sophie Hourihane, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Bence Bécsy, J. Andrew Casey-Clyde , et al. (71 additional authors not shown)

Abstract: This paper presents rigorous tests of pulsar timing array methods and software, examining their consistency across a wide range of injected parameters and signal strength. We discuss updates to the 15-year isotropic gravitational-wave background analyses and their corresponding code representations. Descriptions of the internal structure of the flagship algorithms \texttt{Enterprise} and \texttt{P… ▽ More This paper presents rigorous tests of pulsar timing array methods and software, examining their consistency across a wide range of injected parameters and signal strength. We discuss updates to the 15-year isotropic gravitational-wave background analyses and their corresponding code representations. Descriptions of the internal structure of the flagship algorithms \texttt{Enterprise} and \texttt{PTMCMCSampler} are given to facilitate understanding of the PTA likelihood structure, how models are built, and what methods are currently used in sampling the high-dimensional PTA parameter space. We introduce a novel version of the PTA likelihood that uses a two-step marginalization procedure that performs much faster when the white noise parameters remain fixed. We perform stringent tests of consistency and correctness of the Bayesian and frequentist analysis software. For the Bayesian analysis, we test prior recovery, injection recovery, and Bayes factors. For the frequentist analysis, we test that the cross-correlation-based optimal statistic, when modified to account for a non-negligible gravitational-wave background, accurately recovers the amplitude of the background. We also summarize recent advances and tests performed on the optimal statistic in the literature from both GWB detection and parameter estimation perspectives. The tests presented here validate current and future analyses of PTA data. △ Less

Submitted 7 July, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 30 pages, 10 figures, 1 table; Companion paper to "The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background"; For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year Data Set: Bayesian Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Robin Case, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan DeCesar, Paul B. Demorest, Matthew C. Digman, Timothy Dolch, Brendan Drachler , et al. (74 additional authors not shown)

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 datasets. Here we present the search for individual supermassive black hole binaries in the NANOGrav… ▽ More 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 datasets. Here we present the search for individual supermassive black hole binaries in the NANOGrav 15-year dataset. 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 $\sim$4 nHz, and another at $\sim$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. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 23 pages, 13 figures, 2 tables. Accepted for publication in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year Data Set: Search for Anisotropy in the Gravitational-Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Elizabeth C. Ferrara, William Fiore , et al. (68 additional authors not shown)

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 dataset. 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… ▽ More 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 dataset. 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, and place a Bayesian $95\%$ upper limit on the level of broadband anisotropy such that $(C_{l>0} / C_{l=0}) < 20\%$. 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 dataset, and show that this dataset 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. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 19 pages, 11 figures; submitted to Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year Data Set: Constraints on Supermassive Black Hole Binaries from the Gravitational Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Paul T. Baker, Bence Bécsy, Laura Blecha, Alexander Bonilla, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Belinda D. Cheeseboro, Siyuan Chen, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Curt J. Cutler , et al. (89 additional authors not shown)

Abstract: The NANOGrav 15-year 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 popul… ▽ More The NANOGrav 15-year 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-year 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 sub-parsec 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. △ Less

Submitted 18 July, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: Accepted by Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org. Edited to fix two equation typos (Eq.13 & 21), and minor text typos

The NANOGrav 15-year Data Set: Search for Signals from New Physics

Authors: Adeela Afzal, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Jose Juan Blanco-Pillado, Laura Blecha, Kimberly K. Boddy, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Belinda D. Cheeseboro, Siyuan Chen, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie , et al. (98 additional authors not shown)

Abstract: The 15-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) shows positive evidence for the presence of a low-frequency gravitational-wave (GW) background. In this paper, we investigate potential cosmological interpretations of this signal, specifically cosmic inflation, scalar-induced GWs, first-order phase transitions, cosmic string… ▽ More The 15-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) shows positive evidence for the presence of a low-frequency gravitational-wave (GW) background. In this paper, we investigate potential cosmological interpretations of this signal, specifically cosmic inflation, scalar-induced GWs, first-order phase transitions, cosmic strings, and domain walls. We find that, with the exception of stable cosmic strings of field theory origin, all these models can reproduce the observed signal. When compared to the standard interpretation in terms of inspiraling supermassive black hole binaries (SMBHBs), many cosmological models seem to provide a better fit resulting in Bayes factors in the range from 10 to 100. However, these results strongly depend on modeling assumptions about the cosmic SMBHB population and, at this stage, should not be regarded as evidence for new physics. Furthermore, we identify excluded parameter regions where the predicted GW signal from cosmological sources significantly exceeds the NANOGrav signal. These parameter constraints are independent of the origin of the NANOGrav signal and illustrate how pulsar timing data provide a new way to constrain the parameter space of these models. Finally, we search for deterministic signals produced by models of ultralight dark matter (ULDM) and dark matter substructures in the Milky Way. We find no evidence for either of these signals and thus report updated constraints on these models. In the case of ULDM, these constraints outperform torsion balance and atomic clock constraints for ULDM coupled to electrons, muons, or gluons. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 74 pages, 31 figures, 4 tables; published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-Year Data Set: Detector Characterization and Noise Budget

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. Decesar, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca , et al. (66 additional authors not shown)

Abstract: Pulsar timing arrays (PTAs) are galactic-scale gravitational wave 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 gravitational wave (GW) signals. We present a noise and sensitivity analysis to accompany the NANOGrav 15-year data release and associated… ▽ More Pulsar timing arrays (PTAs) are galactic-scale gravitational wave 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 gravitational wave (GW) signals. We present a noise and sensitivity analysis to accompany the NANOGrav 15-year 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\times 10^{-15}$ at 5 nHz. A power law-integrated analysis shows rough agreement with the amplitudes recovered in NANOGrav's 15-year 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. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 67 pages, 73 figures, 3 tables; published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year Data Set: Observations and Timing of 68 Millisecond Pulsars

Authors: Gabriella Agazie, Md Faisal Alam, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Laura Blecha, Victoria Bonidie, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Bence Bécsy, Christopher Chapman, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Brendan Drachler , et al. (75 additional authors not shown)

Abstract: We present observations and timing analyses of 68 millisecond pulsars (MSPs) comprising the 15-year 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. This is NANOGrav's fifth public data release, including both "narrowband" and "wideband" time-of-arrival… ▽ More We present observations and timing analyses of 68 millisecond pulsars (MSPs) comprising the 15-year 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. 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 three years, now spanning nearly 16 years 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 gravitational-wave background. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 90 pages, 74 figures, 6 tables; published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Becsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Belinda D. Cheeseboro, Siyuan Chen, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Curt J. Cutler, Megan E. DeCesar , et al. (89 additional authors not shown)

Abstract: We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15-year 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-spectr… ▽ More We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15-year 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 these latter Bayes factors using a method that removes inter-pulsar correlations from our data set, finding $p = 10^{-3}$ (approx. $3σ$) for the observed Bayes factors in the null no-correlation scenario. A frequentist test statistic built directly as a weighted sum of inter-pulsar correlations yields $p = 5 \times 10^{-5} - 1.9 \times 10^{-4}$ (approx. $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.7}_{-0.6} \times 10^{-15}$ (median + 90% credible interval) at a reference frequency of 1/(1 yr). 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. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 30 pages, 18 figures. Published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

Measurements of the Crab Pulsar's Giant Radio Pulse Amplitude Power-Law Index Using Low-Frequency Arecibo and Green Bank Telescope Observations

Authors: F. Crawford, T. J. W. Lazio, A. McEwen, J. S. Deneva, J. M. Cordes, L. Spitler, R. F. Trainor

Abstract: We report two low-frequency measurements of the power-law index for the amplitudes of giant radio pulses from the Crab pulsar. The two observations were taken with the Arecibo and Green Bank radio telescopes at center frequencies of 327 MHz and 350 MHz, respectively. We find best-fit values for the differential power-law index $β$ (where $dN/dS \propto S^β$ and $S$ is pulse amplitude) of… ▽ More We report two low-frequency measurements of the power-law index for the amplitudes of giant radio pulses from the Crab pulsar. The two observations were taken with the Arecibo and Green Bank radio telescopes at center frequencies of 327 MHz and 350 MHz, respectively. We find best-fit values for the differential power-law index $β$ (where $dN/dS \propto S^β$ and $S$ is pulse amplitude) of $-2.63 \pm 0.05$ and $-3.6 \pm 0.5$ from the Arecibo and Green Bank data sets, respectively. Both values are broadly consistent with other values previously measured for the Crab pulsar at low radio frequencies. These reported values may be useful in future giant pulse studies of the Crab pulsar. △ Less

Submitted 5 April, 2023; originally announced April 2023.

Comments: 7 pages with 1 table and 1 figure. Accepted by The Astrophysical Journal

Searching for continuous Gravitational Waves in the second data release of the International Pulsar Timing Array

Authors: M. Falxa, S. Babak, P. T. Baker, B. Bécsy, A. Chalumeau, S. Chen, Z. Chen, N. J. Cornish, L. Guillemot, J. S. Hazboun, C. M. F. Mingarelli, A. Parthasarathy, A. Petiteau, N. S. Pol, A. Sesana, S. B. Spolaor, S. R. Taylor, G. Theureau, M. Vallisneri, S. J. Vigeland, C. A. Witt, X. Zhu, J. Antoniadis, Z. Arzoumanian, M. Bailes , et al. (102 additional authors not shown)

Abstract: The International Pulsar Timing Array 2nd data release is the combination of datasets from worldwide collaborations. In this study, we search for continuous waves: gravitational wave signals produced by individual supermassive black hole binaries in the local universe. We consider binaries on circular orbits and neglect the evolution of orbital frequency over the observational span. We find no evi… ▽ More The International Pulsar Timing Array 2nd data release is the combination of datasets from worldwide collaborations. In this study, we search for continuous waves: gravitational wave signals produced by individual supermassive black hole binaries in the local universe. We consider binaries on circular orbits and neglect the evolution of orbital frequency over the observational span. We find no evidence for such signals and set sky averaged 95% upper limits on their amplitude h 95 . The most sensitive frequency is 10nHz with h 95 = 9.1 10-15 . We achieved the best upper limit to date at low and high frequencies of the PTA band thanks to improved effective cadence of observations. In our analysis, we have taken into account the recently discovered common red noise process, which has an impact at low frequencies. We also find that the peculiar noise features present in some pulsars data must be taken into account to reduce the false alarm. We show that using custom noise models is essential in searching for continuous gravitational wave signals and setting the upper limit. △ Less

Submitted 19 March, 2023; originally announced March 2023.

The NANOGrav 12.5-year Data Set: Bayesian Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries

Authors: Zaven Arzoumanian, Paul T. Baker, Laura Blecha, Harsha Blumer, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Bence Bécsy, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Siyuan Chen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Justin A. Ellis, E. C. Ferrara, William Fiore, Emmanuel Fonseca, Gabriel E. Freedman , et al. (53 additional authors not shown)

Abstract: Pulsar timing array collaborations, such as the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), are seeking to detect nanohertz gravitational waves emitted by supermassive black hole binaries formed in the aftermath of galaxy mergers. We have searched for continuous waves from individual circular supermassive black hole binaries using the NANOGrav's recent 12.5-year data s… ▽ More Pulsar timing array collaborations, such as the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), are seeking to detect nanohertz gravitational waves emitted by supermassive black hole binaries formed in the aftermath of galaxy mergers. We have searched for continuous waves from individual circular supermassive black hole binaries using the NANOGrav's recent 12.5-year data set. We created new methods to accurately model the uncertainties on pulsar distances in our analysis, and we implemented new techniques to account for a common red noise process in pulsar timing array data sets while searching for deterministic gravitational wave signals, including continuous waves. As we found no evidence for continuous waves in our data, we placed 95\% upper limits on the strain amplitude of continuous waves emitted by these sources. At our most sensitive frequency of 7.65 nanohertz, we placed a sky-averaged limit of $h_0 < $ $(6.82 \pm 0.35) \times 10^{-15}$, and $h_0 <$ $(2.66 \pm 0.15) \times 10^{-15}$ in our most sensitive sky location. Finally, we placed a multi-messenger limit of $\mathcal{M} <$ $(1.41 \pm 0.02) \times 10^9 M_\odot$ on the chirp mass of the supermassive black hole binary candidate 3C~66B. △ Less

Submitted 6 June, 2023; v1 submitted 9 January, 2023; originally announced January 2023.

Comments: 22 pages, 12 figures. Accepted by ApJL

An unusual pulse shape change event in PSR J1713+0747 observed with the Green Bank Telescope and CHIME

Authors: Ross J. Jennings, James M. Cordes, Shami Chatterjee, Maura A. McLaughlin, Paul B. Demorest, Zaven Arzoumanian, Paul T. Baker, Harsha Blumer, Paul R. Brook, Tyler Cohen, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Timothy Dolch, Elizabeth C. Ferrara, Emmanuel Fonseca, Deborah C. Good, Jeffrey S. Hazboun, Megan L. Jones, David L. Kaplan, Michael T. Lam, T. Joseph W. Lazio, Duncan R. Lorimer, Jing Luo, Ryan S. Lynch , et al. (19 additional authors not shown)

Abstract: The millisecond pulsar J1713+0747 underwent a sudden and significant pulse shape change between April 16 and 17, 2021 (MJDs 59320 and 59321). Subsequently, the pulse shape gradually recovered over the course of several months. We report the results of continued multi-frequency radio observations of the pulsar made using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) and the 100-meter G… ▽ More The millisecond pulsar J1713+0747 underwent a sudden and significant pulse shape change between April 16 and 17, 2021 (MJDs 59320 and 59321). Subsequently, the pulse shape gradually recovered over the course of several months. We report the results of continued multi-frequency radio observations of the pulsar made using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) and the 100-meter Green Bank Telescope (GBT) in a three-year period encompassing the shape change event, between February 2020 and February 2023. As of February 2023, the pulse shape had returned to a state similar to that seen before the event, but with measurable changes remaining. The amplitude of the shape change and the accompanying TOA residuals display a strong non-monotonic dependence on radio frequency, demonstrating that the event is neither a glitch (the effects of which should be independent of radio frequency, $ν$) nor a change in dispersion measure (DM) alone (which would produce a delay proportional to $ν^{-2}$). However, it does bear some resemblance to the two previous "chromatic timing events" observed in J1713+0747 (Demorest et al. 2013; Lam et al. 2016), as well as to a similar event observed in PSR J1643-1224 in 2015 (Shannon et al. 2016). △ Less

Submitted 31 January, 2024; v1 submitted 21 October, 2022; originally announced October 2022.

Comments: 18 pages, 8 figures. Submitted to ApJ. Data available at https://doi.org/10.5281/zenodo.7236459

Four New Fast Radio Bursts Discovered in the Parkes 70-cm Pulsar Survey Archive

Authors: F. Crawford, S. Hisano, M. Golden, T. Kikunaga, A. Laity, D. Zoeller

Abstract: We present four new fast radio bursts discovered in a search of the Parkes 70-cm pulsar survey data archive for dispersed single pulses and bursts. We searched dispersion measures (DMs) ranging between 0 to 5000 pc cm$^{-3}$ with the HEIMDALL and FETCH detection and classification algorithms. All four of the FRBs discovered have significantly larger widths ($> 50$ ms) than almost all of the FRBs d… ▽ More We present four new fast radio bursts discovered in a search of the Parkes 70-cm pulsar survey data archive for dispersed single pulses and bursts. We searched dispersion measures (DMs) ranging between 0 to 5000 pc cm$^{-3}$ with the HEIMDALL and FETCH detection and classification algorithms. All four of the FRBs discovered have significantly larger widths ($> 50$ ms) than almost all of the FRBs detected and cataloged to date. The large pulse widths are not dominated by interstellar scattering or dispersive smearing within channels. One of the FRBs has a DM of 3338 pc cm$^{3}$, the largest measured for any FRB to date. These are also the first FRBs detected by any radio telescope so far, predating the Lorimer Burst by almost a decade. Our results suggest that pulsar survey archives remain important sources of previously undetected FRBs and that searches for FRBs on time scales extending beyond $\sim 100$ ms may reveal the presence of a larger population of wide-pulse FRBs. △ Less

Submitted 25 July, 2022; originally announced July 2022.

Comments: Accepted for publication in MNRAS

Discovery of PSR J0523-7125 as a Circularly Polarized Variable Radio Source in the Large Magellanic Cloud

Authors: Yuanming Wang, Tara Murphy, David L. Kaplan, Teresa Klinner-Teo, Alessandro Ridolfi, Matthew Bailes, Fronefield Crawford, Shi Dai, Dougal Dobie, B. M. Gaensler, Vanessa Graber, Ian Heywood, Emil Lenc, Duncan R. Lorimer, Maura A. McLaughlin, Andrew O'Brien, Sergio Pintaldi, Joshua Pritchard, Nanda Rea, Joshua P. Ridley, Michele Ronchi, Ryan M. Shannon, Gregory R. Sivakoff, Adam Stewart, Ziteng Wang , et al. (1 additional authors not shown)

Abstract: We report the discovery of a highly circularly polarized, variable, steep-spectrum pulsar in the Australian Square Kilometre Array Pathfinder (ASKAP) Variables and Slow Transients (VAST) survey. The pulsar is located about $1^\circ$ from the center of the Large Magellanic Cloud, and has a significant fractional circular polarization of $\sim$20%. We discovered pulsations with a period of 322.5 ms,… ▽ More We report the discovery of a highly circularly polarized, variable, steep-spectrum pulsar in the Australian Square Kilometre Array Pathfinder (ASKAP) Variables and Slow Transients (VAST) survey. The pulsar is located about $1^\circ$ from the center of the Large Magellanic Cloud, and has a significant fractional circular polarization of $\sim$20%. We discovered pulsations with a period of 322.5 ms, dispersion measure (DM) of 157.5 pc cm$^{-3}$, and rotation measure (RM) of $+456$ rad m$^{-2}$ using observations from the MeerKAT and the Parkes telescopes. This DM firmly places the source, PSR J0523$-$7125, in the Large Magellanic Cloud (LMC). This RM is extreme compared to other pulsars in the LMC (more than twice that of the largest previously reported one). The average flux density of $\sim$1 mJy at 1400 MHz and $\sim$25 mJy at 400 MHz places it among the most luminous radio pulsars known. It likely evaded previous discovery because of its very steep radio spectrum (spectral index $α\approx -3$, where $S_ν\propto ν^α$) and broad pulse profile (duty cycle $\gtrsim35$%). We discuss implications for searches for unusual radio sources in continuum images, as well as extragalactic pulsars in the Magellanic Clouds and beyond. Our result highlighted the possibility of identifying pulsars, especially extreme pulsars, from radio continuum images. Future large-scale radio surveys will give us an unprecedented opportunity to discover more pulsars and potentially the most distant pulsars beyond the Magellanic Clouds. △ Less

Submitted 3 May, 2022; v1 submitted 1 May, 2022; originally announced May 2022.

Comments: 17 pages, 7 figures; accepted by ApJ

A Parkes "Murriyang" Search for Pulsars and Transients in the Large Magellanic Cloud

Authors: Shinnosuke Hisano, Fronefield Crawford, Victoria Bonidie, Md F. Alam, Keitaro Takahashi, Duncan R. Lorimer, Josh P. Ridley, Maura M. McLaughlin, Benetge B. P. Perera

Abstract: We report on the latest results of a Parkes multibeam survey for pulsars and dispersed radio bursts in the Large Magellanic Cloud (LMC). We conducted both periodicity and single-pulse searches at a much larger range of trial dispersion measures (DMs) than previously searched. We detected 229 single pulses with signal-to-noise ratio ($\rm{S/N}) > 7$ that were classified by the deep learning network… ▽ More We report on the latest results of a Parkes multibeam survey for pulsars and dispersed radio bursts in the Large Magellanic Cloud (LMC). We conducted both periodicity and single-pulse searches at a much larger range of trial dispersion measures (DMs) than previously searched. We detected 229 single pulses with signal-to-noise ratio ($\rm{S/N}) > 7$ that were classified by the deep learning network FETCH as being real (with $> 90$% likelihood), of which 9 were from the known giant-pulse-emitting pulsar PSR B0540$-$69. Two possibly repeating sources were detected with DMs suggesting that they lie within the LMC, but these require confirmation. Only 3 of the 220 unknown pulses had S/N greater than 8, and the DM distribution for these detected pulses follows an exponential falloff with increasing DM and does not show any excess of signals at DM values expected for the LMC. These features suggest that the detected pulses are not likely to be real, although they are visually compelling. We also report the discovery of a new pulsar (PSR J0556$-$67) in our periodicity search. This pulsar has a spin period of 791 ms, a DM of 71 cm$^{-3}$ pc, an estimated 1400 MHz flux density of $\sim 0.12$ mJy, and shows no evidence of binary motion. Future observations may be able to confirm whether any of the weak but promising candidates in our single pulse and periodicity searches, including our two possible repeaters, are real or not. △ Less

Submitted 23 February, 2022; v1 submitted 22 February, 2022; originally announced February 2022.

Comments: 17 pages, 11 figures

The International Pulsar Timing Array second data release: Search for an isotropic Gravitational Wave Background

Authors: J. Antoniadis, Z. Arzoumanian, S. Babak, M. Bailes, A. -S. Bak Nielsen, P. T. Baker, C. G. Bassa, B. Becsy, A. Berthereau, M. Bonetti, A. Brazier, P. R. Brook, M. Burgay, S. Burke-Spolaor, R. N. Caballero, J. A. Casey-Clyde, A. Chalumeau, D. J. Champion, M. Charisi, S. Chatterjee, S. Chen, I. Cognard, J. M. Cordes, N. J. Cornish, F. Crawford , et al. (101 additional authors not shown)

Abstract: We searched for an isotropic stochastic gravitational wave background in the second data release of the International Pulsar Timing Array, a global collaboration synthesizing decadal-length pulsar-timing campaigns in North America, Europe, and Australia. In our reference search for a power law strain spectrum of the form $h_c = A(f/1\,\mathrm{yr}^{-1})^α$, we found strong evidence for a spectrally… ▽ More We searched for an isotropic stochastic gravitational wave background in the second data release of the International Pulsar Timing Array, a global collaboration synthesizing decadal-length pulsar-timing campaigns in North America, Europe, and Australia. In our reference search for a power law strain spectrum of the form $h_c = A(f/1\,\mathrm{yr}^{-1})^α$, we found strong evidence for a spectrally-similar low-frequency stochastic process of amplitude $A = 3.8^{+6.3}_{-2.5}\times10^{-15}$ and spectral index $α= -0.5 \pm 0.5$, where the uncertainties represent 95\% credible regions, using information from the auto- and cross-correlation terms between the pulsars in the array. For a spectral index of $α= -2/3$, as expected from a population of inspiralling supermassive black hole binaries, the recovered amplitude is $A = 2.8^{+1.2}_{-0.8}\times10^{-15}$. Nonetheless, no significant evidence of the Hellings-Downs correlations that would indicate a gravitational-wave origin was found. We also analyzed the constituent data from the individual pulsar timing arrays in a consistent way, and clearly demonstrate that the combined international data set is more sensitive. Furthermore, we demonstrate that this combined data set produces comparable constraints to recent single-array data sets which have more data than the constituent parts of the combination. Future international data releases will deliver increased sensitivity to gravitational wave radiation, and significantly increase the detection probability. △ Less

Submitted 11 January, 2022; originally announced January 2022.

Comments: 17 pages, 12 figures, accepted in MNRAS

The NANOGrav 12.5-year data set: Search for Non-Einsteinian Polarization Modes in theGravitational-Wave Background

Authors: Zaven Arzoumanian, Paul T. Baker, Harsha Blumer, Bence Becsy, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Maria Charisi, Shami Chatterjee, Siyuan Chen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Dallas M. DeGan, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Justin A. Ellis, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Nathan Garver-Daniels, Peter A. Gentile , et al. (46 additional authors not shown)

Abstract: We search NANOGrav's 12.5-year data set for evidence of a gravitational wave background (GWB) with all the spatial correlations allowed by general metric theories of gravity. We find no substantial evidence in favor of the existence of such correlations in our data. We find that scalar-transverse (ST) correlations yield signal-to-noise ratios and Bayes factors that are higher than quadrupolar (ten… ▽ More We search NANOGrav's 12.5-year data set for evidence of a gravitational wave background (GWB) with all the spatial correlations allowed by general metric theories of gravity. We find no substantial evidence in favor of the existence of such correlations in our data. We find that scalar-transverse (ST) correlations yield signal-to-noise ratios and Bayes factors that are higher than quadrupolar (tensor transverse, TT) correlations. Specifically, we find ST correlations with a signal-to-noise ratio of 2.8 that are preferred over TT correlations (Hellings and Downs correlations) with Bayesian odds of about 20:1. However, the significance of ST correlations is reduced dramatically when we include modeling of the Solar System ephemeris systematics and/or remove pulsar J0030$+$0451 entirely from consideration. Even taking the nominal signal-to-noise ratios at face value, analyses of simulated data sets show that such values are not extremely unlikely to be observed in cases where only the usual TT modes are present in the GWB. In the absence of a detection of any polarization mode of gravity, we place upper limits on their amplitudes for a spectral index of $γ= 5$ and a reference frequency of $f_\text{yr} = 1 \text{yr}^{-1}$. Among the upper limits for eight general families of metric theories of gravity, we find the values of $A^{95\%}_{TT} = (9.7 \pm 0.4)\times 10^{-16}$ and $A^{95\%}_{ST} = (1.4 \pm 0.03)\times 10^{-15}$ for the family of metric spacetime theories that contain both TT and ST modes. △ Less

Submitted 29 September, 2021; originally announced September 2021.

Comments: 24 pages, 18 figures, 3 appendices. Please send any comments/questions to Nima Laal (laaln@oregonstate.edu)

Journal ref: The Astrophysical Journal Letters, vol. 923, no. 2, p. L22, Dec. 2021

Study of 72 pulsars discovered in the PALFA survey: Timing analysis, glitch activity, emission variability, and a pulsar in an eccentric binary

Authors: E. Parent, H. Sewalls, P. C. C. Freire, T. Matheny, A. G. Lyne, B. B. P. Perera, F. Cardoso, M. A. McLaughlin, B. Allen, A. Brazier, F. Camilo, S. Chatterjee, J. M. Cordes, F. Crawford, J. S. Deneva, F. A. Dong, R. D. Ferdman, E. Fonseca, J. W. T. Hessels, V. M. Kaspi, B. Knispel, J. van Leeuwen, R. S. Lynch, B. M. Meyers, J. W. McKee , et al. (9 additional authors not shown)

Abstract: We present new discoveries and results from long-term timing of 72 pulsars discovered in the Arecibo PALFA survey, including precise determination of astrometric and spin parameters, and flux density and scatter broadening measurements at 1.4 GHz. Notable discoveries include two young pulsars (characteristic ages $\sim$30 kyr) with no apparent supernova remnant associations, three mode changing, 1… ▽ More We present new discoveries and results from long-term timing of 72 pulsars discovered in the Arecibo PALFA survey, including precise determination of astrometric and spin parameters, and flux density and scatter broadening measurements at 1.4 GHz. Notable discoveries include two young pulsars (characteristic ages $\sim$30 kyr) with no apparent supernova remnant associations, three mode changing, 12 nulling and two intermittent pulsars. We detected eight glitches in five pulsars. Among them is PSR J1939+2609, an apparently old pulsar (characteristic age $\sim$1 Gy), and PSR J1954+2529, which likely belongs to a newly-emerging class of binary pulsars. The latter is the only pulsar among the 72 that is clearly not isolated: a non-recycled neutron star with a 931-ms spin period in an eccentric ($e\,=\,0.114$) wide ($P_b\,=\,82.7\,$d) orbit with a companion of undetermined nature having a minimum mass of $\sim0.6\,M_{\odot}$. Since operations at Arecibo ceased in 2020 August, we give a final tally of PALFA sky coverage, and compare its 207 pulsar discoveries to the known population. On average, they are 50% more distant than other Galactic plane radio pulsars; PALFA millisecond pulsars (MSP) have twice the dispersion measure per unit spin period than the known population of MSP in the Plane. The four intermittent pulsars discovered by PALFA more than double the population of such objects, which should help to improve our understanding of pulsar magnetosphere physics. The statistics for these, RRATS, and nulling pulsars suggest that there are many more of these objects in the Galaxy than was previously thought. △ Less

Submitted 17 November, 2021; v1 submitted 4 August, 2021; originally announced August 2021.

Comments: 39 pages, 18 figures, 9 tables. Accepted for publication in ApJ

An Arecibo Search for Fast Radio Transients from M87

Authors: Akshay Suresh, Shami Chatterjee, James M. Cordes, Fronefield Crawford

Abstract: The possible origin of millisecond bursts from the giant elliptical galaxy M87 has been scrutinized since the earliest searches for extragalactic fast radio transients undertaken in the late 1970s. Motivated by rapid technological advancements in recent years, we conducted $\rm \simeq 10~hours$ of L-band ($\rm 1.15-1.75~GHz$) observations of the core of M87 with the Arecibo radio telescope in 2019… ▽ More The possible origin of millisecond bursts from the giant elliptical galaxy M87 has been scrutinized since the earliest searches for extragalactic fast radio transients undertaken in the late 1970s. Motivated by rapid technological advancements in recent years, we conducted $\rm \simeq 10~hours$ of L-band ($\rm 1.15-1.75~GHz$) observations of the core of M87 with the Arecibo radio telescope in 2019. Adopting a matched filtering approach, we searched our data for single pulses using trial dispersion measures up to $\rm 5500~pc~cm^{-3}$ and burst durations between $\rm 0.3-123~ms$. We find no evidence of astrophysical bursts in our data above a 7$σ$ detection threshold. Our observations thus constrain the burst rate from M87 to $\rm \lesssim 0.1~bursts~hr^{-1}$ above $\rm 1.4~Jy~ms$, the most stringent upper limit obtained to date. Our non-detection of radio bursts is consistent with expectations of giant pulse emission from a Crab-like young neutron star population in M87. However, the dense, strongly magnetized interstellar medium surrounding the central $\sim 10^9 \ M_{\odot}$ supermassive black hole of M87 may potentially harbor magnetars that can emit detectable radio bursts during their flaring states. △ Less

Submitted 11 October, 2021; v1 submitted 20 July, 2021; originally announced July 2021.

Comments: 15 pages, 4 figures, published in ApJ

Journal ref: 2021 ApJ 920 16

An Arecibo 327 MHz Search for Radio Pulsars and Bursts in the Dwarf Irregular Galaxies Leo A and T

Authors: Fronefield Crawford, Kohei Hayashi

Abstract: We have conducted an Arecibo 327 MHz search of two dwarf irregular galaxies in the Local Group, Leo A and T, for radio pulsars and single pulses from fast radio bursts and other giant pulse emitters. We detected no astrophysical signals in this search, and we estimate flux density limits on both periodic and burst emission. Our derived luminosity limits indicate that only the most luminous radio p… ▽ More We have conducted an Arecibo 327 MHz search of two dwarf irregular galaxies in the Local Group, Leo A and T, for radio pulsars and single pulses from fast radio bursts and other giant pulse emitters. We detected no astrophysical signals in this search, and we estimate flux density limits on both periodic and burst emission. Our derived luminosity limits indicate that only the most luminous radio pulsars known in our Galaxy and in the Magellanic Clouds (MCs) would have been detectable in our search if they were at the distances of Leo A and T. Given the much smaller stellar mass content and star formation rates of Leo A and T compared to the Milky Way and the MCs, there are likely to be few (if any) extremely luminous pulsars in these galaxies. It is therefore not surprising that we detected no pulsars in our search. △ Less

Submitted 16 June, 2021; originally announced June 2021.

Searching For Gravitational Waves From Cosmological Phase Transitions With The NANOGrav 12.5-year dataset

Authors: Zaven Arzoumanian, Paul T. Baker, Harsha Blumer, Bence Bécsy, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Maria Charisi, Shami Chatterjee, Siyuan Chen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Justin A. Ellis, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Nathan Garver-Daniels, Peter A. Gentile, Deborah C. Good, Jeffrey S. Hazboun , et al. (40 additional authors not shown)

Abstract: We search for a first-order phase transition gravitational wave signal in 45 pulsars from the NANOGrav 12.5 year dataset. We find that the data can be modeled in terms of a strong first order phase transition taking place at temperatures below the electroweak scale. However, we do not observe any strong preference for a phase-transition interpretation of the signal over the standard astrophysical… ▽ More We search for a first-order phase transition gravitational wave signal in 45 pulsars from the NANOGrav 12.5 year dataset. We find that the data can be modeled in terms of a strong first order phase transition taking place at temperatures below the electroweak scale. However, we do not observe any strong preference for a phase-transition interpretation of the signal over the standard astrophysical interpretation in terms of supermassive black holes mergers; but we expect to gain additional discriminating power with future datasets, improving the signal to noise ratio and extending the sensitivity window to lower frequencies. An interesting open question is how well gravitational wave observatories could separate such signals. △ Less

Submitted 11 January, 2022; v1 submitted 28 April, 2021; originally announced April 2021.

Comments: 13 pages, 4 figures. v2: updated to match published version

Journal ref: Phys.Rev.Lett. 127 (2021) 25, 251302

Reprocessing of a Green Bank 43-meter Telescope Survey of Unidentified Bright Radio Sources for Pulsars and Radio Bursts

Authors: F. Crawford, J. Margeson, B. Nguyen, T. Saigal, O. Young, D. Agarwal, K. Aggarwal

Abstract: We have reprocessed a set of observations of 75 bright, unidentified, steep-spectrum polarized radio sources taken with the Green Bank 43-m telescope to find previously undetected sub-millisecond pulsars and radio bursts. The (null) results of the first search of these data were reported by Schmidt et al. Our reprocessing searched for single pulses out to a dispersion measure (DM) of 1000 pc cm… ▽ More We have reprocessed a set of observations of 75 bright, unidentified, steep-spectrum polarized radio sources taken with the Green Bank 43-m telescope to find previously undetected sub-millisecond pulsars and radio bursts. The (null) results of the first search of these data were reported by Schmidt et al. Our reprocessing searched for single pulses out to a dispersion measure (DM) of 1000 pc cm$^{-3}$ which were classified using the Deep Learning based classifier FETCH. We also searched for periodicities at a wider range of DMs and accelerations. Our search was sensitive to highly accelerated, rapidly rotating pulsars (including sub-millisecond pulsars) in compact binary systems as well as to highly-dispersed impulsive signals, such as fast radio bursts. No pulsars or astrophysical burst signals were found in the reprocessing. △ Less

Submitted 1 February, 2021; originally announced February 2021.

Journal ref: Research Notes of the American Astronomical Society, 5, 21 (2021)

The NANOGrav 11yr Data Set: Limits on Supermassive Black Hole Binaries in Galaxies within 500Mpc

Authors: Zaven Arzoumanian, Paul T. Baker, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Bence Becsy, Maria Charisi, Shami Chatterjee, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Rodney D. Elliott, Justin A. Ellis, Elizabeth C. Ferrara, Emmanuel Fonseca, Nathan Garver-Daniels, Peter A. Gentile, Deborah C. Good, Jeffrey S. Hazboun, Kristina Islo, Ross J. Jennings , et al. (32 additional authors not shown)

Abstract: Supermassive black hole binaries (SMBHBs) should form frequently in galactic nuclei as a result of galaxy mergers. At sub-parsec separations, binaries become strong sources of low-frequency gravitational waves (GWs), targeted by Pulsar Timing Arrays (PTAs). We used recent upper limits on continuous GWs from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11yr dataset to… ▽ More Supermassive black hole binaries (SMBHBs) should form frequently in galactic nuclei as a result of galaxy mergers. At sub-parsec separations, binaries become strong sources of low-frequency gravitational waves (GWs), targeted by Pulsar Timing Arrays (PTAs). We used recent upper limits on continuous GWs from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11yr dataset to place constraints on putative SMBHBs in nearby massive galaxies. We compiled a comprehensive catalog of ~44,000 galaxies in the local universe (up to redshift ~0.05) and populated them with hypothetical binaries, assuming that the total mass of the binary is equal to the SMBH mass derived from global scaling relations. Assuming circular equal-mass binaries emitting at NANOGrav's most sensitive frequency of 8nHz, we found that 216 galaxies are within NANOGrav's sensitivity volume. We ranked the potential SMBHBs based on GW detectability by calculating the total signal-to-noise ratio (S/N) such binaries would induce within the NANOGrav array. We placed constraints on the chirp mass and mass ratio of the 216 hypothetical binaries. For 19 galaxies, only very unequal-mass binaries are allowed, with the mass of the secondary less than 10 percent that of the primary, roughly comparable to constraints on a SMBHB in the Milky Way. Additionally, we were able to exclude binaries delivered by major mergers (mass ratio of at least 1/4) for several of these galaxies. We also derived the first limit on the density of binaries delivered by major mergers purely based on GW data. △ Less

Submitted 7 January, 2021; originally announced January 2021.

Comments: Submitted to ApJ. Send comments to Maria Charisi (maria.charisi@nanograv.org)

Astrophysics Milestones For Pulsar Timing Array Gravitational Wave Detection

Authors: Nihan S. Pol, Stephen R. Taylor, Luke Zoltan Kelley, Sarah J. Vigeland, Joseph Simon, Siyuan Chen, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Shami Chatterjee, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Nathan Garver-Daniels, Deborah C. Good , et al. (27 additional authors not shown)

Abstract: The NANOGrav Collaboration reported strong Bayesian evidence for a common-spectrum stochastic process in its 12.5-yr pulsar timing array dataset, with median characteristic strain amplitude at periods of a year of $A_{\rm yr} = 1.92^{+0.75}_{-0.55} \times 10^{-15}$. However, evidence for the quadrupolar Hellings \& Downs interpulsar correlations, which are characteristic of gravitational wave sign… ▽ More The NANOGrav Collaboration reported strong Bayesian evidence for a common-spectrum stochastic process in its 12.5-yr pulsar timing array dataset, with median characteristic strain amplitude at periods of a year of $A_{\rm yr} = 1.92^{+0.75}_{-0.55} \times 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 dataset, 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 yrs of data, i.e., an additional 2--5 yrs from the 12.5 yr dataset; (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~yrs 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 datasets. Most notably, by combining the data of individual arrays into the International Pulsar Timing Array, all of these milestones can be reached significantly earlier. △ Less

Submitted 24 March, 2021; v1 submitted 22 October, 2020; originally announced October 2020.

Comments: 15 pages, 7 figures

The NANOGrav 12.5-year Data Set: Search For An Isotropic Stochastic Gravitational-Wave Background

Authors: Zaven Arzoumanian, Paul T. Baker, Harsha Blumer, Bence Becsy, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Shami Chatterjee, Siyuan Chen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Justin A. Ellis, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Nathan Garver-Daniels, Peter A. Gentile, Deborah C. Good, Jeffrey S. Hazboun, A. Miguel Holgado , et al. (36 additional authors not shown)

Abstract: We search for an isotropic stochastic gravitational-wave background (GWB) in the $12.5$-year 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. The Bayesian posterior of the amplitude for an $f^{-2/3}$ power-… ▽ More We search for an isotropic stochastic gravitational-wave background (GWB) in the $12.5$-year 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. The Bayesian posterior of the amplitude for an $f^{-2/3}$ power-law spectrum, expressed as the characteristic GW strain, has median $1.92 \times 10^{-15}$ and $5\%$--$95\%$ quantiles of $1.37$--$2.67 \times 10^{-15}$ at a reference frequency of $f_\mathrm{yr} = 1 ~\mathrm{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. △ Less

Submitted 8 January, 2021; v1 submitted 9 September, 2020; originally announced September 2020.

Comments: 25 pages, 14 figures, 5 tables, 3 appendices. Published in The Astrophysical Journal Letters. Please send any comments/questions to Joseph Simon (joe.simon@nanograv.org). Jupyter notebook tutorials and some MCMC chain files are available at https://github.com/nanograv/12p5yr_stochastic_analysis

Journal ref: The Astrophysical Journal Letters, Volume 905, Number 2 (2020)

Asymmetric mass ratios for bright double neutron-star mergers

Authors: R. D. Ferdman, P. C. C. Freire, B. B. P. Perera, N. Pol, F. Camilo, S. Chatterjee, J. M. Cordes, F. Crawford, J. W. T. Hessels, V. M. Kaspi, M. A. McLaughlin, E. Parent, I. H. Stairs, J. van Leeuwen

Abstract: The discovery of a radioactively powered kilonova associated with the binary neutron star merger GW170817 was the first - and still only - confirmed electromagnetic counterpart to a gravitational-wave event. However, observations of late-time electromagnetic emission are in tension with the expectations from standard neutron-star merger models. Although the large measured ejecta mass is potentiall… ▽ More The discovery of a radioactively powered kilonova associated with the binary neutron star merger GW170817 was the first - and still only - confirmed electromagnetic counterpart to a gravitational-wave event. However, observations of late-time electromagnetic emission are in tension with the expectations from standard neutron-star merger models. Although the large measured ejecta mass is potentially explained by a progenitor system that is asymmetric in terms of the stellar component masses, i.e. with a mass ratio $q$ of 0.7-0.8, the known Galactic population of merging double neutron star (DNS) systems (i.e. those that will coalesce within billions of years or less) has, until now, only consisted of nearly equal-mass ($q > 0.9$) binaries. PSR J1913+1102 is a DNS system in a 5-hour, low-eccentricity ($e = 0.09$) orbit, implying an orbital separation of 1.8 solar radii, with the two neutron stars predicted to coalesce in 470 million years due to gravitational-wave emission. Here we report that the masses of the two neutron stars, as measured by a dedicated pulsar timing campaign, are $1.62 \pm 0.03$ and $1.27 \pm 0.03$ solar masses for the pulsar and companion neutron star, respectively; with a measured mass ratio $q = 0.78 \pm 0.03$, it is the most asymmetric DNS among known merging systems. Based on this detection, our population synthesis analysis implies that such asymmetric binaries represent between 2 and 30% (90% confidence) of the total population of merging DNS binaries. The coalescence of a member of this population offers a possible explanation for the anomalous properties of GW170817, including the observed kilonova emission from that event. △ Less

Submitted 8 July, 2020; originally announced July 2020.

Comments: 17 pages, 3 figures, published in Nature on 9 July 2020 (https://www.nature.com/articles/s41586-020-2439-x); this is the authors' final submitted version before review. Includes methods and supplementary information

Multi-Messenger Gravitational Wave Searches with Pulsar Timing Arrays: Application to 3C66B Using the NANOGrav 11-year Data Set

Authors: Zaven Arzoumanian, Paul T. Baker, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Bence Becsy, Maria Charisi, Shami Chatterjee, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Rodney D. Elliott, Justin A. Ellis, Robert D. Ferdman, Elizabeth C. Ferrara, Emmanuel Fonseca, Nathan Garver-Daniels, Peter A. Gentile, Deborah C. Good, Jeffrey S. Hazboun , et al. (34 additional authors not shown)

Abstract: When galaxies merge, the supermassive black holes in their centers may form binaries and, during the process of merger, emit low-frequency gravitational radiation in the process. In this paper we consider the galaxy 3C66B, which was used as the target of the first multi-messenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of th… ▽ More When galaxies merge, the supermassive black holes in their centers may form binaries and, during the process of merger, emit low-frequency gravitational radiation in the process. In this paper we consider the galaxy 3C66B, which was used as the target of the first multi-messenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C66B to less than $(1.65\pm0.02) \times 10^9~{M_\odot}$ using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits, and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data to `blind' pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences. △ Less

Submitted 12 August, 2020; v1 submitted 14 May, 2020; originally announced May 2020.

Comments: 14 pages, 6 figures. Accepted by ApJ

The NANOGrav 12.5-year Data Set: Wideband Timing of 47 Millisecond Pulsars

Authors: Md F. Alam, Zaven Arzoumanian, Paul T. Baker, Harsha Blumer, Keith E. Bohler, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Keeisi Caballero, Richard S. Camuccio, Rachel L. Chamberlain, Shami Chatterjee, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Justin A. Ellis, Robert D. Ferdman, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Yhamil Garcia , et al. (45 additional authors not shown)

Abstract: We present a new analysis of the profile data from the 47 millisecond pulsars comprising the 12.5-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), which is presented in a parallel paper (Alam et al. 2021a; NG12.5). Our reprocessing is performed using "wideband" timing methods, which use frequency-dependent template profiles, simultaneous time-of-arrival… ▽ More We present a new analysis of the profile data from the 47 millisecond pulsars comprising the 12.5-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), which is presented in a parallel paper (Alam et al. 2021a; NG12.5). Our reprocessing is performed using "wideband" timing methods, which use frequency-dependent template profiles, simultaneous time-of-arrival (TOA) and dispersion measure (DM) measurements from broadband observations, and novel analysis techniques. In particular, the wideband DM measurements are used to constrain the DM portion of the timing model. We compare the ensemble timing results to NG12.5 by examining the timing residuals, timing models, and noise model components. There is a remarkable level of agreement across all metrics considered. Our best-timed pulsars produce encouragingly similar results to those from NG12.5. In certain cases, such as high-DM pulsars with profile broadening, or sources that are weak and scintillating, wideband timing techniques prove to be beneficial, leading to more precise timing model parameters by 10-15%. The high-precision, multi-band measurements of several pulsars indicate frequency-dependent DMs. Compared to the narrowband analysis in NG12.5, the TOA volume is reduced by a factor of 33, which may ultimately facilitate computational speed-ups for complex pulsar timing array analyses. This first wideband pulsar timing data set is a stepping stone, and its consistent results with NG12.5 assure us that such data sets are appropriate for gravitational wave analyses. △ Less

Submitted 18 December, 2020; v1 submitted 13 May, 2020; originally announced May 2020.

Comments: 62 pages, 55 figures, 5 tables, 3 appendices. Data available at http://nanograv.org/data/ and via DOI 10.5281/zenodo.4312887

Journal ref: The Astrophysical Journal Supplement Series, 252, 5 (2021)

The NANOGrav 12.5 yr Data Set: Observations and Narrowband Timing of 47 Millisecond Pulsars

Authors: Md F. Alam, Zaven Arzoumanian, Paul T. Baker, Harsha Blumer, Keith E. Bohler, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Keeisi Caballero, Richard S. Camuccio, Rachel L. Chamberlain, Shami Chatterjee, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Justin A. Ellis, Robert D. Ferdman, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Yhamil Garcia , et al. (45 additional authors not shown)

Abstract: We present time-of-arrival (TOA) measurements and timing models of 47 millisecond pulsars (MSPs) observed from 2004 to 2017 at the Arecibo Observatory and the Green Bank Telescope by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). The observing cadence was three to four weeks for most pulsars over most of this time span, with weekly observations of six sources. These d… ▽ More We present time-of-arrival (TOA) measurements and timing models of 47 millisecond pulsars (MSPs) observed from 2004 to 2017 at the Arecibo Observatory and the Green Bank Telescope by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). The observing cadence was three to four weeks for most pulsars over most of this time span, with weekly observations of six sources. These data were collected for use in low-frequency gravitational wave searches and for other astrophysical purposes. We detail our observational methods and present a set of TOA measurements, based on "narrowband" analysis, in which many TOAs are calculated within narrow radio-frequency bands for data collected simultaneously across a wide bandwidth. A separate set of "wideband" TOAs will be presented in a companion paper. We detail a number of methodological changes compared to our previous work which yield a cleaner and more uniformly processed data set. Our timing models include several new astrometric and binary pulsar measurements, including previously unpublished values for the parallaxes of PSRs J1832-0836 and J2322+2057, the secular derivatives of the projected semi-major orbital axes of PSRs J0613-0200 and J2229+2643, and the first detection of the Shapiro delay in PSR J2145-0750. We report detectable levels of red noise in the time series for 14 pulsars. As a check on timing model reliability, we investigate the stability of astrometric parameters across data sets of different lengths. We report flux density measurements for all pulsars observed. Searches for stochastic and continuous gravitational waves using these data will be subjects of forthcoming publications. △ Less

Submitted 23 December, 2020; v1 submitted 13 May, 2020; originally announced May 2020.

Comments: 54 pages, 52 figures, 7 tables, 1 appendix. Data are available at http://nanograv.org/data/ and via the DOI 10.5281/zenodo.4312297

Journal ref: The Astrophysical Journal Supplement Series, 252, 4 (2021)

Analysis of quasar magnitudes

Authors: David F. Crawford

Abstract: Since their discovery the analysis of quasar magnitudes has generally required some form of evolution. Assuming that quasars do not have evolution this paper shows that they have a well-defined intrinsic magnitude distribution that is independent of cosmological models. However the average apparent magnitudes are essentially constant which means that the only cosmological information they contain… ▽ More Since their discovery the analysis of quasar magnitudes has generally required some form of evolution. Assuming that quasars do not have evolution this paper shows that they have a well-defined intrinsic magnitude distribution that is independent of cosmological models. However the average apparent magnitudes are essentially constant which means that the only cosmological information they contain is that the absolute flux density has a power law distribution. Thus quasar magnitudes, by themselves, are essentially useless for cosmological investigations. △ Less

Submitted 29 January, 2020; originally announced January 2020.

Comments: 4 pages, 2 figures

Modeling the uncertainties of solar-system ephemerides for robust gravitational-wave searches with pulsar timing arrays

Authors: M. Vallisneri, S. R. Taylor, J. Simon, W. M. Folkner, R. S. Park, C. Cutler, J. A. Ellis, T. J. W. Lazio, S. J. Vigeland, K. Aggarwal, Z. Arzoumanian, P. T. Baker, A. Brazier, P. R. Brook, S. Burke-Spolaor, S. Chatterjee, J. M. Cordes, N. J. Cornish, F. Crawford, H. T. Cromartie, K. Crowter, M. DeCesar, P. B. Demorest, T. Dolch, R. D. Ferdman , et al. (39 additional authors not shown)

Abstract: 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 Earth's 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 perform… ▽ More 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 Earth's 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 (Jupiter's orbital elements) most relevant to gravitational-wave searches with pulsar timing arrays. This model, BayesEphem, was used to derive ephemeris-robust results in NANOGrav's 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. △ Less

Submitted 6 January, 2020; v1 submitted 2 January, 2020; originally announced January 2020.

Comments: Fixed typo in author list. Code that supports all calculations and figures is available at github.com/nanograv/11yr_stochastic_analysis/tree/master/bayesephem

The NANOGrav 11-Year Data Set: Limits on Gravitational Wave Memory

Authors: K. Aggarwal, Z. Arzoumanian, P. T. Baker, A. Brazier, P. R. Brook, S. Burke-Spolaor, S. Chatterjee, J. M. Cordes, N. J. Cornish, F. Crawford, H. T. Cromartie, K. Crowter, M. Decesar, P. B. Demorest, T. Dolch, J. A. Ellis, R. D. Ferdman, E. C. Ferrara, E. Fonseca, N. Garver-Daniels, P. Gentile, D. Good, J. S. Hazboun, A. M. Holgado, E. A. Huerta , et al. (36 additional authors not shown)

Abstract: The mergers of supermassive black hole binaries (SMBHBs) promise to be incredible sources of gravitational waves (GWs). While the oscillatory part of the merger gravitational waveform will be outside the frequency sensitivity range of pulsar timing arrays (PTAs), the non-oscillatory GW memory effect is detectable. Further, any burst of gravitational waves will produce GW memory, making memory a us… ▽ More The mergers of supermassive black hole binaries (SMBHBs) promise to be incredible sources of gravitational waves (GWs). While the oscillatory part of the merger gravitational waveform will be outside the frequency sensitivity range of pulsar timing arrays (PTAs), the non-oscillatory GW memory effect is detectable. Further, any burst of gravitational waves will produce GW memory, making memory a useful probe of unmodeled exotic sources and new physics. We searched the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GW memory. This dataset is sensitive to very low frequency GWs of $\sim3$ to $400$ nHz (periods of $\sim11$ yr $-$ $1$ mon). Finding no evidence for GWs, we placed limits on the strain amplitude of GW memory events during the observation period. We then used the strain upper limits to place limits on the rate of GW memory causing events. At a strain of $2.5\times10^{-14}$, corresponding to the median upper limit as a function of source sky position, we set a limit on the rate of GW memory events at $<0.4$ yr$^{-1}$. That strain corresponds to a SMBHB merger with reduced mass of $ηM \sim 2\times10^{10}M_\odot$ and inclination of $ι=π/3$ at a distance of 1 Gpc. As a test of our analysis, we analyzed the NANOGrav 9-year data set as well. This analysis found an anomolous signal, which does not appear in the 11-year data set. This signal is not a GW, and its origin remains unknown. △ Less

Submitted 6 December, 2019; v1 submitted 19 November, 2019; originally announced November 2019.

Comments: 10 pages, 6 figures, submitted to ApJ

The NANOGrav 11-Year Data Set: Evolution of Gravitational Wave Background Statistics

Authors: J. S. Hazboun, J. Simon, S. R. Taylor, M. T. Lam, S. J. Vigeland, K. Islo, J. S. Key, Z. Arzoumanian, P. T. Baker, A. Brazier, P. R. Brook, S. Burke-Spolaor, S. Chatterjee, J. M. Cordes, N. J. Cornish, F. Crawford, K. Crowter, H. T. Cromartie, M. DeCesar, P. B. Demorest, T. Dolch, J. A. Ellis, R. D. Ferdman, E. Ferrara, E. Fonseca , et al. (38 additional authors not shown)

Abstract: An ensemble of inspiraling supermassive black hole binaries should produce a stochastic background of very low frequency gravitational waves. This stochastic background is predicted to be a power law, with a spectral index of -2/3, and it should be detectable by a network of precisely timed millisecond pulsars, widely distributed on the sky. This paper reports a new "time slicing" analysis of the… ▽ More An ensemble of inspiraling supermassive black hole binaries should produce a stochastic background of very low frequency gravitational waves. This stochastic background is predicted to be a power law, with a spectral index of -2/3, and it should be detectable by a network of precisely timed millisecond pulsars, widely distributed on the sky. This paper reports a new "time slicing" analysis of the 11-year data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) using 34 millisecond pulsars. Methods to flag potential "false positive" signatures are developed, including techniques to identify responsible pulsars. Mitigation strategies are then presented. We demonstrate how an incorrect noise model can lead to spurious signals, and show how independently modeling noise across 30 Fourier components, spanning NANOGrav's frequency range, effectively diagnoses and absorbs the excess power in gravitational-wave searches. This results in a nominal, and expected, progression of our gravitational-wave statistics. Additionally we show that the first interstellar medium event in PSR J1713+0747 pollutes the common red noise process with low-spectral index noise, and use a tailored noise model to remove these effects. △ Less

Submitted 20 September, 2019; v1 submitted 18 September, 2019; originally announced September 2019.

Comments: 14 pages, 13 figures, fixed typo in abstract of earlier version

Eight Millisecond Pulsars Discovered in the Arecibo PALFA Survey

Authors: E. Parent, V. M. Kaspi, S. M. Ransom, P. C. C. Freire, A. Brazier, F. Camilo, S. Chatterjee, J. M. Cordes, F. Crawford, J. S. Deneva, R. D. Ferdman, J. W. T. Hessels, J. van Leeuwen, A. G. Lyne, E. C. Madsen, M. A. McLaughlin, C. Patel, P. Scholz, I. H. Stairs, B. W. Stappers, W. W. Zhu

Abstract: We report on eight millisecond pulsars (MSPs) in binary systems discovered with the Arecibo PALFA survey. Phase-coherent timing solutions derived from 2.5 to 5 years of observations carried out at Arecibo and Jodrell Bank observatories are provided. PSR J1921+1929 is a 2.65-ms pulsar in a 39.6-day orbit for which we detect $γ$-ray pulsations in archival Fermi data. PSR J1928+1245 is a very low-mas… ▽ More We report on eight millisecond pulsars (MSPs) in binary systems discovered with the Arecibo PALFA survey. Phase-coherent timing solutions derived from 2.5 to 5 years of observations carried out at Arecibo and Jodrell Bank observatories are provided. PSR J1921+1929 is a 2.65-ms pulsar in a 39.6-day orbit for which we detect $γ$-ray pulsations in archival Fermi data. PSR J1928+1245 is a very low-mass-function system with an orbital period of 3.3 hours that belongs to the non-eclipsing black widow population. We also present PSR J1932+1756, the longest-orbital-period (41.5 days) intermediate-mass binary pulsar known to date. In light of the numerous discoveries of binary MSPs over the past years, we characterize the Galactic distribution of known MSP binaries in terms of binary class. Our results support and strengthen previous claims that the scatter in the Galactic scale height distribution correlates inversely with the binary mass function. We provide evidence of observational biases against detecting the most recycled pulsars near the Galactic plane, which overestimates the scale height of lighter systems. A possible bimodality in the mass function of MSPs with massive white dwarfs is also reported. △ Less

Submitted 26 August, 2019; originally announced August 2019.

Comments: Submitted to ApJ, 21 pages, 10 figures

NANOGrav Education and Outreach: Growing a Diverse and Inclusive Collaboration for Low-Frequency Gravitational Wave Astronomy

Authors: The NANOGrav Collaboration, P. T. Baker, H. Blumer, A. Brazier, S. Chatterjee, B. Christy, F. Crawford, M. E. DeCesar, T. Dolch, N. E. Garver-Daniels, J. S. Hazboun, K. Holley-Bockelmann, D. L. Kaplan, J . S. Key, T. C. Klein, M. T. Lam, N. Lewandowska, D. R. Lorimer, R. S. Lynch, M. A. McLaughlin, N. McMann, J. Page, N. T. Palliyaguru, J. D. Romano, X. Siemens , et al. (3 additional authors not shown)

Abstract: The new field of gravitational wave astrophysics requires a growing pool of students and researchers with unique, interdisciplinary skill sets. It also offers an opportunity to build a diverse, inclusive astronomy community from the ground up. We describe the efforts used by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) NSF Physics Frontiers Center to foster such grow… ▽ More The new field of gravitational wave astrophysics requires a growing pool of students and researchers with unique, interdisciplinary skill sets. It also offers an opportunity to build a diverse, inclusive astronomy community from the ground up. We describe the efforts used by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) NSF Physics Frontiers Center to foster such growth by involving students at all levels in low-frequency gravitational wave astrophysics with pulsar timing arrays (PTAs) and establishing collaboration policies that ensure broad participation by diverse groups. We describe and illustrate the impact of these techniques on our collaboration as a case study for other distributed collaborations. △ Less

Submitted 17 July, 2019; originally announced July 2019.

Comments: Submitted as a State of the Profession White Paper for the Astro2020 Decadal Survey (10 pages + cover and references, 6 figures)

Mass measurements for two binary pulsars discovered in the PALFA survey

Authors: W. W. Zhu, P. C. C. Freire, B. Knispel, B. Allen, B. W. Stappers, A. G. Lyne, S. Chatterjee, J. M. Cordes, F. Crawford, J. S. Deneva, R. D. Ferdman, J. W. T. Hessels, V. M. Kaspi, P. Lazarus, R. Lynch, S. M. Ransom, K. Stovall, J. Y. Donner

Abstract: In this paper, we present the results of timing observations of PSRs J1949+3106 and J1950+2414, two binary millisecond pulsars discovered in data from the Arecibo ALFA pulsar survey (PALFA). The timing parameters include precise measurements of the proper motions of both pulsars, which show that PSR J1949+3106 has a transversal motion very similar to that of an object in the local standard of rest… ▽ More In this paper, we present the results of timing observations of PSRs J1949+3106 and J1950+2414, two binary millisecond pulsars discovered in data from the Arecibo ALFA pulsar survey (PALFA). The timing parameters include precise measurements of the proper motions of both pulsars, which show that PSR J1949+3106 has a transversal motion very similar to that of an object in the local standard of rest. The timing also includes measurements of the Shapiro delay and the rate of advance of periastron for both systems. Assuming general relativity, these allow estimates of the masses of the components of the two systems; for PSR J1949+3106, the pulsar mass is $M_p \, = \, 1.34^{+0.17}_{-0.15} \, M_{\odot}$ and the companion mass $M_c \, = \, 0.81^{+0.06}_{-0.05}\, M_{\odot}$; for PSR J1950+2414 $M_p \, = \, 1.496 \, \pm \, 0.023\, M_{\odot}$ and $M_c \, = \, 0.280^{+0.005}_{-0.004}\, M_{\odot}$ (all values 68.3 % confidence limits). We use these masses and proper motions to investigate the evolutionary history of both systems: PSR J1949+3106 is likely the product of a low-kick supernova; PSR J1950+2414 is a member of a new class of eccentric millisecond pulsar binaries with an unknown formation mechanism. We discuss the proposed hypotheses for the formations of these systems in light of our new mass measurements. △ Less

Submitted 11 July, 2019; originally announced July 2019.

Comments: 12 pages, 5 figures, accepted for publication in the Astrophysical Journal

Unidentified FRBs in archival data

Authors: E. F. Keane, D. R. Lorimer, F. Crawford

Abstract: Recently Zhang et al. (2019) reported the discovery of FRB 010312 in a dataset previously searched for FRBs. Here we explain how and why this FRB was initially missed, and highlight several caveats relevant to FRB search completeness. It is possible that up to $\sim40$% of discoverable FRBs remain undiscovered in some existing public domain archival data sets. The situation could be even more pron… ▽ More Recently Zhang et al. (2019) reported the discovery of FRB 010312 in a dataset previously searched for FRBs. Here we explain how and why this FRB was initially missed, and highlight several caveats relevant to FRB search completeness. It is possible that up to $\sim40$% of discoverable FRBs remain undiscovered in some existing public domain archival data sets. The situation could be even more pronounced given that FRB search pipelines generally look for temporally-symmetric broadband flat-spectrum pulses; most FRBs do not look like that. △ Less

Submitted 1 March, 2019; originally announced March 2019.

Comments: 2 pages. A follow-up note on arXiv:1902.06009. Published in RNAAS, submitted 2019/02/21, accepted 2019/02/23, published 2019/02/27

Journal ref: RNAAS, 2019, volume 3, number 2, article id 41

The NANOGrav 11-Year Data Set: Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries

Authors: K. Aggarwal, Z. Arzoumanian, P. T. Baker, A. Brazier, M. R. Brinson, P. R. Brook, S. Burke-Spolaor, S. Chatterjee, J. M. Cordes, N. J. Cornish, F. Crawford, K. Crowter, H. T. Cromartie, M. DeCesar, P. B. Demorest, T. Dolch, J. A. Ellis, R. D. Ferdman, E. Ferrara, E. Fonseca, N. Garver-Daniels, P. Gentile, J. S. Hazboun, A. M. Holgado, E. A. Huerta , et al. (38 additional authors not shown)

Abstract: Observations indicate that nearly all galaxies contain supermassive black holes (SMBHs) 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 (PTAs). We have searched the recently-released North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 1… ▽ More Observations indicate that nearly all galaxies contain supermassive black holes (SMBHs) 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 (PTAs). We have searched the recently-released North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year 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 as a function of GW frequency and sky location. We placed a sky-averaged upper limit on the GW strain of $h_0 < 7.3(3) \times 10^{-15}$ at $f_\mathrm{gw}= 8$ nHz. 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 in measurements from individual pulsars. We used our upper limits on the GW strain to place lower limits on the distances to individual SMBHBs. At the most-sensitive sky location, we ruled out SMBHBs emitting GWs with $f_\mathrm{gw}= 8$ nHz within 120 Mpc for $\mathcal{M} = 10^9 \, M_\odot$, and within 5.5 Gpc for $\mathcal{M} = 10^{10} \, M_\odot$. We also determined that there are no SMBHBs with $\mathcal{M} > 1.6 \times 10^9 \, M_\odot$ emitting GWs in the Virgo Cluster. Finally, we estimated the number of potentially detectable sources given our current strain upper limits based on galaxies in Two Micron All-Sky Survey (2MASS) and merger rates from the Illustris cosmological simulation project. Only 34 out of 75,000 realizations of the local Universe contained a detectable source, from which we concluded it was unsurprising that we did not detect any individual sources given our current sensitivity to GWs. △ Less

Submitted 21 May, 2019; v1 submitted 30 December, 2018; originally announced December 2018.

Comments: 10 pages, 11 figures. Accepted by Astrophysical Journal. Please send any comments/questions to S. J. Vigeland (vigeland@uwm.edu)

Journal ref: Astrophys. J. 880, 2 (2019)

FRB 121102 Bursts Show Complex Time-Frequency Structure

Authors: J. W. T. Hessels, L. G. Spitler, A. D. Seymour, J. M. Cordes, D. Michilli, R. S. Lynch, K. Gourdji, A. M. Archibald, C. G. Bassa, G. C. Bower, S. Chatterjee, L. Connor, F. Crawford, J. S. Deneva, V. Gajjar, V. M. Kaspi, A. Keimpema, C. J. Law, B. Marcote, M. A. McLaughlin, Z. Paragi, E. Petroff, S. M. Ransom, P. Scholz, B. W. Stappers , et al. (1 additional authors not shown)

Abstract: FRB 121102 is the only known repeating fast radio burst source. Here we analyze a wide-frequency-range (1-8 GHz) sample of high-signal-to-noise, coherently dedispersed bursts detected using the Arecibo and Green Bank telescopes. These bursts reveal complex time-frequency structures that include sub-bursts with finite bandwidths. The frequency-dependent burst structure complicates the determination… ▽ More FRB 121102 is the only known repeating fast radio burst source. Here we analyze a wide-frequency-range (1-8 GHz) sample of high-signal-to-noise, coherently dedispersed bursts detected using the Arecibo and Green Bank telescopes. These bursts reveal complex time-frequency structures that include sub-bursts with finite bandwidths. The frequency-dependent burst structure complicates the determination of a dispersion measure (DM); we argue that it is appropriate to use a DM metric that maximizes frequency-averaged pulse structure, as opposed to peak signal-to-noise, and find DM = 560.57 +/- 0.07 pc/cc at MJD 57644. After correcting for dispersive delay, we find that the sub-bursts have characteristic frequencies that typically drift lower at later times in the total burst envelope. In the 1.1-1.7 GHz band, the ~ 0.5-1-ms sub-bursts have typical bandwidths ranging from 100-400 MHz, and a characteristic drift rate of ~ 200 MHz/ms towards lower frequencies. At higher radio frequencies, the sub-burst bandwidths and drift rate are larger, on average. While these features could be intrinsic to the burst emission mechanism, they could also be imparted by propagation effects in the medium local to the source. Comparison of the burst DMs with previous values in the literature suggests an increase of Delta(DM) ~ 1-3 pc/cc in 4 years, though this could be a stochastic variation as opposed to a secular trend. This implies changes in the local medium or an additional source of frequency-dependent delay. Overall, the results are consistent with previously proposed scenarios in which FRB 121102 is embedded in a dense nebula. △ Less

Submitted 26 November, 2018; originally announced November 2018.

Comments: Submitted to ApJ; comments welcome

PALFA Single-Pulse Pipeline: New Pulsars, Rotating Radio Transients and a Candidate Fast Radio Burst

Authors: C. Patel, D. Agarwal, M. Bhardwaj, M. M. Boyce, A. Brazier, S. Chatterjee, P. Chawla, V. M. Kaspi, D. R. Lorimer, M. M. McLaughlin, E. Parent, Z. Pleunis, S. M. Ransom, P. Scholz, R. S. Wharton, W. W. Zhu, M. Alam, K. Caballero Valdez, F. Camilo, J. M. Cordes, F. Crawford, J. S. Deneva, R. D. Ferdman, P. C. C. Freire, J. W. T. Hessels , et al. (4 additional authors not shown)

Abstract: We present a newly implemented single-pulse pipeline for the PALFA survey to efficiently identify single radio pulses from pulsars, Rotating Radio Transients (RRATs) and Fast Radio Bursts (FRBs). We have conducted a sensitivity analysis of this new pipeline in which multiple single pulses with a wide range of parameters were injected into PALFA data sets and run through the pipeline. Based on the… ▽ More We present a newly implemented single-pulse pipeline for the PALFA survey to efficiently identify single radio pulses from pulsars, Rotating Radio Transients (RRATs) and Fast Radio Bursts (FRBs). We have conducted a sensitivity analysis of this new pipeline in which multiple single pulses with a wide range of parameters were injected into PALFA data sets and run through the pipeline. Based on the recovered pulses, we find that for pulse widths $\rm < 5\ ms$ the sensitivity of the PALFA pipeline is at most a factor of $\rm \sim 2$ less sensitive to single pulses than our theoretical predictions. For pulse widths $\rm > 10\ ms$, as the $\rm DM$ decreases, the degradation in sensitivity gets worse and can increase up to a factor of $\rm \sim 4.5$. Using this pipeline, we have thus far discovered 7 pulsars and 2 RRATs and identified 3 candidate RRATs and 1 candidate FRB. The confirmed pulsars and RRATs have DMs ranging from 133 to 386 pc cm$^{-3}$ and flux densities ranging from 20 to 160 mJy. The pulsar periods range from 0.4 to 2.1 s. We report on candidate FRB 141113, which we argue is likely astrophysical and extragalactic, having $\rm DM \simeq 400\ pc~cm^{-3}$, which represents an excess over the Galactic maximum along this line of sight of $\rm \sim$ 100 - 200 pc cm$^{-3}$. We consider implications for the FRB population and show via simulations that if FRB 141113 is real and extragalactic, the slope $α$ of the distribution of integral source counts as a function of flux density ($N (>S) \propto S^{-α}$) is $1.4 \pm 0.5$ (95% confidence range). However this conclusion is dependent on several assumptions that require verification. △ Less

Submitted 10 August, 2018; originally announced August 2018.

Comments: 21 pages, 12 figures

The implementation of a Fast-Folding pipeline for long-period pulsar searching in the PALFA survey

Authors: E. Parent, V. M. Kaspi, S. M. Ransom, M. Krasteva, C. Patel, P. Scholz, A. Brazier, M. A. McLaughlin, M. Boyce, W. W. Zhu, Z. Pleunis, B. Allen, S. Bogdanov, K. Caballero, F. Camilo, R. Camuccio, S. Chatterjee, J. M. Cordes, F. Crawford, J. S. Deneva, R. Ferdman, P. C. C. Freire, J. W. T. Hessels, F. A. Jenet, B. Knispel , et al. (9 additional authors not shown)

Abstract: The Pulsar Arecibo L-Band Feed Array (PALFA) survey, the most sensitive blind search for radio pulsars yet conducted, is ongoing at the Arecibo Observatory in Puerto Rico. The vast majority of the 180 pulsars discovered by PALFA have spin periods shorter than 2 seconds. Pulsar surveys may miss long-period radio pulsars due to the summing of a finite number of harmonic components in conventional Fo… ▽ More The Pulsar Arecibo L-Band Feed Array (PALFA) survey, the most sensitive blind search for radio pulsars yet conducted, is ongoing at the Arecibo Observatory in Puerto Rico. The vast majority of the 180 pulsars discovered by PALFA have spin periods shorter than 2 seconds. Pulsar surveys may miss long-period radio pulsars due to the summing of a finite number of harmonic components in conventional Fourier analyses (typically $\sim$16), or due to the strong effect of red noise at low modulation frequencies. We address this reduction in sensitivity by using a time-domain search technique: the Fast-Folding Algorithm (FFA). We designed a program that implements a FFA-based search in the PALFA processing pipeline, and tested the efficiency of the algorithm by performing tests under both ideal, white noise conditions, as well as with real PALFA observational data. In the two scenarios, we show that the time-domain algorithm has the ability to outperform the FFT-based periodicity search implemented in the survey. We perform simulations to compare the previously reported PALFA sensitivity with that obtained using our new FFA implementation. These simulations show that for a pulsar having a pulse duty cycle of roughly 3%, the performance of our FFA pipeline exceeds that of our FFT pipeline for pulses with DM $\lesssim$ 40 pc cm$^{-3}$ and for periods as short as $\sim$500 ms, and that the survey sensitivity is improved by at least a factor of two for periods $\gtrsim$ 6 sec. Discoveries from the implementation of the algorithm in PALFA are also presented in this paper. △ Less

Submitted 21 May, 2018; originally announced May 2018.

Comments: Accepted for publication in The Astrophysical Journal; 19 pages, 9 figures, 2 tables