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.