COSMOS: A numerical relativity code specialized for PBH formation

Abstract

Primordial black holes (PBHs) are black holes generated in the early universe without having gone through stellar evolution. In the standard formation process, PBHs are formed from super-horizon primordial fluctuations with non-linearly large initial amplitude. In order to simulate the non-linear gravitational dynamics of PBH formation, one has to rely on numerical relativity solvers to approximate the solution of the Einstein equations. COSMOS is a C++ package for solving the Einstein equations in 3+1 dimensions, providing simple tools for the simulation of PBH formation. In order to resolve the collapsing region, non-Cartesian scale-up coordinates and a fixed mesh-refinement procedure are implemented. In COSMOS, a massless scalar field and a perfect fluid with a linear equation of state are implemented as matter fields. To achieve a practically acceptable computational speed, OpenMP is used for the parallelization. COSMOS has no other dependencies, which makes for an easier installation.