More Like this

1. The NANOGrav Search for Signals from New Physics: MCMC chains

   https://doi.org/10.5281/zenodo.8083620  

   NANOGrav, The Collaboration ( January 2023 , Zenodo)

   MCMC chains for the GWB analyses performed in the paper "The NANOGrav 15 yr Data Set: Search for Signals from New Physics". 

   The data is provided in pickle format. Each file contains a NumPy array with the MCMC chain (with burn-in already removed), and a dictionary with the model parameters' names as keys and their priors as values. You can load them as

   with open ('path/to/file.pkl', 'rb') as pick: temp = pickle.load(pick) params = temp[0] chain = temp[1]

   The naming convention for the files is the following:

   • igw: inflationary Gravitational Waves (GWs)
   • sigw: scalar-induced GWs
     • sigw_box: assumes a box-like feature in the primordial power spectrum.
     • sigw_delta: assumes a delta-like feature in the primordial power spectrum.
     • sigw_gauss: assumes a Gaussian peak feature in the primordial power spectrum.
   • pt: cosmological phase transitions
     • pt_bubble: assumes that the dominant contribution to the GW productions comes from bubble collisions.
     • pt_sound: assumes that the dominant contribution to the GW productions comes from sound waves.
   • stable: stable cosmic strings
     • stable-c: stable strings emitting GWs only in the form of GW bursts from cusps on closed loops.
     • stable-k: stable strings emitting GWs only in the form of GW bursts from kinks on closed loops.
     • stable-m: stable strings emitting monochromatic GW at the fundamental frequency.
     • stable-n: stable strings described by numerical simulations including GWs from cusps and kinks.
   • meta: metastable cosmic strings
     • meta-l: metastable strings with GW emission from loops only.
     • meta-ls metastable strings with GW emission from loops and segments.
   • super: cosmic superstrings.
   • dw: domain walls
     • dw-sm: domain walls decaying into Standard Model particles.
     • dw-dr: domain walls decaying into dark radiation.

   For each model, we provide four files. One for the run where the new-physics signal is assumed to be the only GWB source. One for the run where the new-physics signal is superimposed to the signal from Supermassive Black Hole Binaries (SMBHB), for these files "_bhb" will be appended to the model name. Then, for both these scenarios, in the "compare" folder we provide the files for the hypermodel runs that were used to derive the Bayes' factors.

   In addition to chains for the stochastic models, we also provide data for the two deterministic models considered in the paper (ULDM and DM substructures). For the ULDM model, the naming convention of the files is the following (all the ULDM signals are superimposed to the SMBHB signal, see the discussion in the paper for more details)

   • uldm_e: ULDM Earth signal.
   • uldm_p: ULDM pulsar signal
     • uldm_p_cor: correlated limit
     • uldm_p_unc: uncorrelated limit
   • uldm_c: ULDM combined Earth + pulsar signal direct coupling 
     • uldm_c_cor: correlated limit
     • uldm_c_unc: uncorrelated limit
   • uldm_vecB: vector ULDM coupled to the baryon number
     • uldm_vecB_cor: correlated limit
     • uldm_vecB_unc: uncorrelated limit 
   • uldm_vecBL: vector ULDM coupled to B-L
     • uldm_vecBL_cor: correlated limit
     • uldm_vecBL_unc: uncorrelated limit
   • uldm_c_grav: ULDM combined Earth + pulsar signal for gravitational-only coupling
     • uldm_c_grav_cor: correlated limit
       • uldm_c_cor_grav_low: low mass region  
       • uldm_c_cor_grav_mon: monopole region
       • uldm_c_cor_grav_low: high mass region
     • uldm_c_unc: uncorrelated limit
       • uldm_c_unc_grav_low: low mass region  
       • uldm_c_unc_grav_mon: monopole region
       • uldm_c_unc_grav_low: high mass region

   For the substructure (static) model, we provide the chain for the marginalized distribution (as for the ULDM signal, the substructure signal is always superimposed to the SMBHB signal)

    

   more » « less
2. 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
3. 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
4. The NANOGrav 12.5 yr Data Set: Search for Gravitational Wave Memory

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

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

   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. We find a Bayes factor of 2.8 in favor of a model that includes a memory signal and common spatially uncorrelated red noise (CURN) compared to a model including only a CURN. However, further investigation shows that a disproportionate amount of support for the memory signal comes from three dubious pulsars. Using a more flexible red-noise model in these pulsars reduces the Bayes factor to 1.3. Having found no compelling evidence, 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 × 10⁻¹⁴. 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 that 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.

    

   more » « less
5. Disentangling the primordial nature of stochastic gravitational wave backgrounds with CMB spectral distortions

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

   Cyr, Bryce ; Kite, Thomas ; Chluba, Jens ; Hill, J. Colin ; Jeong, Donghui ; Acharya, Sandeep Kumar ; Bolliet, Boris ; Patil, Subodh P. ( December 2023 , Monthly Notices of the Royal Astronomical Society)

   The recent detection of a stochastic gravitational wave background (SGWB) at nanohertz frequencies by pulsar timing arrays (PTAs) has sparked a flurry of interest. Beyond the standard interpretation that the progenitor is a network of supermassive black hole binaries, many exotic models have also been proposed, some of which can potentially offer a better fit to the data. We explore how the various connections between gravitational waves (GWs) and cosmic microwave background (CMB) spectral distortions (SDs) can be leveraged to help determine whether an SGWB was generated primordially or astrophysically. To this end, we present updated k-space window functions that can be used for distortion parameter estimation on enhancements to the primordial scalar power spectrum. These same enhancements can also source GWs directly at second order in perturbation theory, so-called scalar-induced GWs (SIGWs), and indirectly through the formation of primordial black holes (PBHs). We perform a mapping of scalar power spectrum constraints into limits on the GW parameter space of SIGWs for δ-function features. We highlight that broader features in the scalar spectrum can explain the PTA results while simultaneously producing an SD within reach of future experiments. We additionally update PBH constraints from μ- and y-type SDs. Refined treatments of the distortion window functions widen existing SD constraints, and we find that a future CMB spectrometer could play a pivotal role in unravelling the origin of GWs imprinted at or below CMB anisotropy scales.

    

   more » « less