Forecasting the sensitivity of pulsar timing arrays to gravitational wave backgrounds

Pulsar timing array (PTA) observations hinted towards the existence of a stochastic gravitational wave background (SGWB) in the nHz frequency band. Still, the nature of the SGWB signal cannot be confidently inferred from current data, and the leading explanation invokes mergers of supermassive black holes. If confirmed, such discovery would not only represent a turning point in our understanding of astrophysics, but it may severely limit the capability of searching for additional cosmological sources in the nHz frequency range. In this work, we build a simple framework to forecast the sensitivity of future pulsar timing array configurations and assess the parameter estimation of SGWB, which could consist of several contributions. We release the python code fastpta, implementing this framework and ready to use.