Clustering of primordial black holes with non-Gaussian initial
  fluctuations
  We formulate the two-point correlation function of primordial black holes
(PBHs) at their formation time, based on the functional integration approach
which has often been used in the context of halo clustering. We find that PBH
clustering on super-Hubble scales could never be induced in the case where the
initial primordial fluctuations are Gaussian, while it can be enhanced by the
so-called local-type trispectrum (four-point correlation function) of the
primordial curvature perturbations.
