pbhstat: A Python package for calculating the primordial black hole abundance
Cosmology and Nongalactic Astrophysics
2025-09-19 v1
Abstract
We present pbhstat, a publicly available Python package designed to compute the mass function and total abundance of primordial black holes (PBHs) from a given primordial power spectrum. The package offers a modular framework using multiple statistical approaches, including Press-Schechter theory, peaks theory, and formalisms based on the non-linear compaction function. Currently, the implementation is limited to scenarios with nearly Gaussian initial conditions.
Keywords
Cite
@article{arxiv.2509.14326,
title = {pbhstat: A Python package for calculating the primordial black hole abundance},
author = {Philippa S. Cole and Jacopo Fumagalli},
journal= {arXiv preprint arXiv:2509.14326},
year = {2025}
}
Comments
13 pages, 2 figures plus references. Accepted in Universe as part of the Special Issue: Primordial Black Holes from Inflation