More Like this

1. The NANOGrav 12.5 yr Data Set: A Computationally Efficient Eccentric Binary Search Pipeline and Constraints on an Eccentric Supermassive Binary Candidate in 3C 66B

   https://doi.org/10.3847/1538-4357/ad1f61  

   Agazie, Gabriella ; Arzoumanian, Zaven ; Baker, Paul T. ; Bécsy, Bence ; Blecha, Laura ; Blumer, Harsha ; Brazier, Adam ; Brook, Paul R. ; Burke-Spolaor, Sarah ; Casey-Clyde, J. Andrew ; et al ( March 2024 , The Astrophysical Journal)

   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 yr period and low redshift (∼0.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 to 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 yr 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 eccentricitye₀< 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 ± 3.7 ns for the Earth term-only and 81.74 ± 0.86 ns for the Earth+pulsar term searches fore₀< 0.5 andη> 0.1. Similar 95% upper limits on the chirp mass are (1.98 ± 0.05) × 10⁹and (1.81 ± 0.01) × 10⁹M_☉. These upper limits, while less stringent than those calculated from a circular binary search in the NANOGrav 12.5 yr data set, are consistent with the SMBHB model of 3C 66B developed from electromagnetic observations.

    

   more » « less
2. The NANOGrav 15 yr Data Set: Observations and Timing of 68 Millisecond Pulsars

   https://doi.org/10.3847/2041-8213/acda9a  

   Agazie, Gabriella ; Alam, Md Faisal ; Anumarlapudi, Akash ; Archibald, Anne M. ; Arzoumanian, Zaven ; Baker, Paul T. ; Blecha, Laura ; Bonidie, Victoria ; Brazier, Adam ; Brook, Paul R. ; et al ( June 2023 , The Astrophysical Journal Letters)

   Abstract We present observations and timing analyses of 68 millisecond pulsars (MSPs) comprising the 15 yr data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). NANOGrav is a pulsar timing array (PTA) experiment that is sensitive to low-frequency gravitational waves (GWs). This is NANOGrav’s fifth public data release, including both “narrowband” and “wideband” time-of-arrival (TOA) measurements and corresponding pulsar timing models. We have added 21 MSPs and extended our timing baselines by 3 yr, now spanning nearly 16 yr for some of our sources. The data were collected using the Arecibo Observatory, the Green Bank Telescope, and the Very Large Array between frequencies of 327 MHz and 3 GHz, with most sources observed approximately monthly. A number of notable methodological and procedural changes were made compared to our previous data sets. These improve the overall quality of the TOA data set and are part of the transition to new pulsar timing and PTA analysis software packages. For the first time, our data products are accompanied by a full suite of software to reproduce data reduction, analysis, and results. Our timing models include a variety of newly detected astrometric and binary pulsar parameters, including several significant improvements to pulsar mass constraints. We find that the time series of 23 pulsars contain detectable levels of red noise, 10 of which are new measurements. In this data set, we find evidence for a stochastic GW background. 

   more » « less
3. Pulsar observations at low frequencies: applications to pulsar timing and solar wind models

   https://doi.org/10.1093/mnras/stac316  

   Kumar, P. ; White, S. M. ; Stovall, K. ; Dowell, J. ; Taylor, G. B. ( February 2022 , Monthly Notices of the Royal Astronomical Society)

   Efforts are underway to use high-precision timing of pulsars in order to detect low-frequency gravitational waves. A limit to this technique is the timing noise generated by dispersion in the plasma along the line of sight to the pulsar, including the solar wind. The effects due to the solar wind vary with time, influenced by the change in solar activity on different time-scales, ranging up to ∼11 yr for a solar cycle. The solar wind contribution depends strongly on the angle between the pulsar line of sight and the solar disc, and is a dominant effect at small separations. Although solar wind models to mitigate these effects do exist, they do not account for all the effects of the solar wind and its temporal changes. Since low-frequency pulsar observations are most sensitive to these dispersive delays, they are most suited to test the efficacy of these models and identify alternative approaches. Here, we investigate the efficacy of some solar wind models commonly used in pulsar timing using long-term, high-cadence data on six pulsars taken with the Long Wavelength Array, and compare them with an operational solar wind model. Our results show that stationary models of the solar wind correction are insufficient to achieve the timing noise desired by pulsar timing experiments, and we need to use non-stationary models, which are informed by other solar wind observations, to obtain accurate timing residuals.

    

   more » « less
4. The NANOGrav 15 yr Data Set: Detector Characterization and Noise Budget

   https://doi.org/10.3847/2041-8213/acda88  

   Agazie, Gabriella ; Anumarlapudi, Akash ; Archibald, Anne M. ; Arzoumanian, Zaven ; Baker, Paul T. ; Bécsy, Bence ; Blecha, Laura ; Brazier, Adam ; Brook, Paul R. ; Burke-Spolaor, Sarah ; et al ( June 2023 , The Astrophysical Journal Letters)

   Pulsar timing arrays (PTAs) are galactic-scale gravitational wave (GW) detectors. Each individual arm, composed of a millisecond pulsar, a radio telescope, and a kiloparsecs-long path, differs in its properties but, in aggregate, can be used to extract low-frequency GW signals. We present a noise and sensitivity analysis to accompany the NANOGrav 15 yr data release and associated papers, along with an in-depth introduction to PTA noise models. As a first step in our analysis, we characterize each individual pulsar data set with three types of white-noise parameters and two red-noise parameters. These parameters, along with the timing model and, particularly, a piecewise-constant model for the time-variable dispersion measure, determine the sensitivity curve over the low-frequency GW band we are searching. We tabulate information for all of the pulsars in this data release and present some representative sensitivity curves. We then combine the individual pulsar sensitivities using a signal-to-noise ratio statistic to calculate the global sensitivity of the PTA to a stochastic background of GWs, obtaining a minimum noise characteristic strain of 7 × 10⁻¹⁵at 5 nHz. A power-law-integrated analysis shows rough agreement with the amplitudes recovered in NANOGrav’s 15 yr GW background analysis. While our phenomenological noise model does not model all known physical effects explicitly, it provides an accurate characterization of the noise in the data while preserving sensitivity to multiple classes of GW signals.

    

   more » « less
5. A New Model for Ionospheric Total Electron Content: The Impact of Solar Flux Proxies and Indices

   https://doi.org/10.1029/2020JA028466  

   Goncharenko, Larisa P. ; Tamburri, Cole A. ; Tobiska, W. Kent ; Schonfeld, Samuel J. ; Chamberlin, Phillip C. ; Woods, Thomas N. ; Didkovsky, Leonid ; Coster, Anthea J. ; Zhang, Shun‐Rong ( February 2021 , Journal of Geophysical Research: Space Physics)

   We present a new high‐resolution empirical model for the ionospheric total electron content (TEC). TEC data are obtained from the global navigation satellite system (GNSS) receivers with a 1° × 1° spatial resolution and 5‐min temporal resolution. The linear regression model is developed at 45°N, 0°E for the years 2000–2019 with 30‐min temporal resolution, unprecedented for typical empirical ionospheric models. The model describes dependency of TEC on solar flux, season, geomagnetic activity, and local time. Parameters describing solar and geomagnetic activity are evaluated. In particular, several options for solar flux input to the model are compared, including the 10.7 cm solar radio flux (F_10.7), the Mg II core‐to‐wing ratio, and formulations of the solar extreme ultraviolet flux (EUV). Ultimately, the extreme ultraviolet flux presented by the Flare Irradiance Spectral Model, integrated from 0.05 to 105.05 nm, best represents the solar flux input to the model. TEC time delays to this solar parameter on the order of several days as well as seasonal modulation of the solar flux terms are included. TheAp₃index and its history are used to reflect the influence of geomagnetic activity. The root mean squared error of the model (relative to the mean TEC observed in the 30‐min window) is 1.9539 TECu. A validation of this model for the first 3 months of 2020 shows excellent agreement with data. The new model shows significant improvement over the International Reference Ionosphere 2016 (IRI‐2016) when the two are compared during 2008 and 2012.

    

   more » « less