Perturbation theory for LCDM and Modified Gravity theories using "fk"-Kernels.
fkpt is a code that computes the 1-loop redshift space power spectrum for tracers.
- Alejandro Aviles (ICF-UNAM, Mexico), avilescervantes@gmail.com
- Mario A. Rodriguez-Meza (ININ, Mexico), marioalberto.rodriguezmeza@gmail.com
Special thanks to Sebastien Fromenteau and Hernán E. Noriega for helping with the Python interface.
Hu-Sawicky f(R) is the only modified gravity model implemented so far. It is straightforward to do it for other models, hope I can find the time to do it soon. Contact me if you have questions about how to do it.
If you run the LCDM model, it will use the exact time dependence for the kernels.
First use git clone:
git clone https://github.com/alejandroaviles/fkpt.git
Compile:
/fkpt$ make
Run:
/fkpt$ ./fkpt
This will compute the LCDM redshift space multipoles of the power spectrum
For help:
/fkpt$ ./fkpt -help
In help you can see how to change parameters, in the form [option]=[value], for example:
/fkpt$ ./fkpt Om=0.3 h=0.7 model=HS fR0=1.0e-6 suffix=_F6z05 zout=0.5 fnamePS=pklin.dat
computes Hu-Sawicky f_R0 = -10^-6 with background cosmology h=0.7, Omega_m = 0.3, at z=0.5, for the input real space linear power spectrum pklin.dat. The output files will have a suffix _F6z05.
The input linear power spectrum should be the LCDM one and given in a two column (k,pkl) file in Mpc/h units at the desire output redshift.
Dependences
The code employs the standard libraries:
Once compiled with everything ready, check the Jupyter Notebook which contains a helpful example.
AP is implemented for Python only
Please cite https://arxiv.org/abs/2312.10510 if you use this code for you research.