The potential energy and atomic forces of each molecular cluster in the dataset were calculated with Gaussian 16 (59) at the MN15/6-31G** level . The MN15 functional was chosen because it is specifically designed to have broad accuracy for both multireference and single-reference systems. The stability of wavefunction was checked during self-consistent field calculation, and if instability was found, further optimization was performed until a stable solution was achieved.