These notebooks reproduce the PTA forecasts from:

Xin, Mingarelli and Hazboun, 2020. "Multimessenger pulsar timing array constraints on supermassive black hole binaries traced by periodic light curves"

They use the Python package Hasasia developed to produce sensitivity curves and skymaps for pulsar timing arrays. The notebook that makes the various forecasted PTA products can be accessed in Google Colab, via the badge above.

The simulated PTAs are made by pulling parameters from various empirical distributions based on current and forecasted PTA propertires, hence they are realization dependent. The exact versions used in the paper are included in the repo, except for the IPTA 2025 forecasted PTA which can be downloaded via this link. Of particular interest are the SKA1 and SKA2 sky maps which are based on a PTA using the first generation of SKA millisecond pulsars and the long baseline datasets of current PTA experiments forecasted into the SKA era.