Monte Carlo (MC) simulations for the Doppler dark matter signals in pulsar timing arrays (PTA). The main module computes the SNR for each of these searches. There is also a supplementary module that computes the signal shape for a realization.
Use a (mini)conda to create an environment with all the necessary dependencies. If you do not have conda installed, install it from here
Create the conda environment with:
conda create --name <env> --file requirements.txt
conda activate <env>
The MC is parallelized using Open MPI. Check the documentation of your system for instruction to load this module. On HPC system it can usually be loaded with
module load openmpi/4.0.4
Once in the conda environment, from the main folder run
mpirun -np 100 python dm-pta-mc.py input.txt
The main module is parallelized over realizations. The argument after '-np' is the number of processors used for the simulation. For optimal performace it should be no more than the number of universes. The last argument passed to the code is the filename of the input parameter file. Some examples are contained in the examples foler:
- pbh_input.txt: exmple input file for PBH (monochromatic halo mass function, point mass)
- AX_input.txt: example input file for axion minicluster
Example data for the axion minicluster model can be found in /examples/data/.
Details of each parameter can be found in the example input files.
The output is saved in the /data folder.
The module that computes the signal shape is /src/signals.py
The MC was authored by Vincent Lee and Tanner Trickle. Contributors include Andrea Mitridate.
ArXiv link to the paper associated to this code will be added shortly.