From each binned category, 50 structures were chosen at random for a total of 150 pairs for each H-bonding and stacking set of interactions and a DFT EDA analysis was run on each structure. For the multiplets, because of the near cubic scaling of DFT calculations, only 50 total computations were able to be computed and chosen at random from all multiplets found.