Extensive calculations of potential energy surfaces for parallel-displaced configurations of pyrrole-pyrrole systems have been carried out by the use of a dispersion-corrected density functional. System geometries associated with the energy minima have been found. The minimum interaction energy has been calculated as -5.38 kcal/mol. However, bonding boundaries appeared to be relatively broad, and stacking interactions can be binding even for ring centroid distances larger than 6 Å . Though the contribution of the correlation energy to intermolecular interaction in pyrrole dimers appeared to be relatively small (around 1.6 smaller than it is in a benzenebenzene system), this system's minimum interaction energy is lower than those calculated for benzene-benzene, benzene-pyridine and even pyridine-pyridine configurations. The calculation of the charges and energy decomposition analysis revealed that the specific charge distribution in a pyrrole molecule and its relatively high polarization are the significant source of the intermolecular interaction in pyrrole dimer systems.