pbhstat: A Python package for calculating the primordial black hole abundance
arXiv:2509.14326
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.
13 pages, 2 figures plus references. Accepted in Universe as part of the Special Issue: Primordial Black Holes from Inflation