Since intermolecular interactions in AC(1)-QM/MM are solely described at a hybrid QM/MM level, devia-tions might be expected for short contacts and potentially accumulate, given their large numbers. To improve upon the single molecule QM model, interactions with neighbouring molecules need to be corrected with QM based descriptions. The latter can be calculated using a supramolecular two-body interaction scheme. A two-body region R within the cluster is selected. All molecules that are within a cutoff distance from the "central molecule" are included. The QM/MM interactions with these molecules need to be removed and QM based interactions need to be added pairwise.