  Primordial black holes (PBH) can efficiently form black hole binaries in the
early universe. We update the resulting constraints on PBH abundance using data
from the third observational run (O3) of LIGO-Virgo-KAGRA. To capture a wide
range of PBH scenarios, we consider a variety of mass functions, including
critical collapse in the QCD epoch in the presence of non-Gaussianities.
Applying hierarchical Bayesian analysis to a population binaries consisting of
primordial and astrophysical black holes, we find that, in every scenario, the
PBHs can make up at most $f_{\rm PBH} \lesssim 10^{-3}$ of dark matter in the
mass range $1-200~M_\odot$. The shape and strength of the constraints are
insensitive to the type of non-Gaussianities, the modifications to the mass
function during the QCD epoch, or the modelling of the astrophysical PBH
population.
