Another relevant property for understanding NCIs in QUID dimers is the atomic forces, which are widely used to parameterize MLFFs for (bio)molecular systems. Indeed, the molecular conformational sampling at a given temperature strongly depends on the accuracy of the chosen computational method in adequately describing the forces acting on the atoms in the molecule. Accordingly, we have analyzed the atomic force distributions using a selection of DFT methods: the best performing ones PBE0+MBD and PBE0+D4, and the well-performing pairwise PBE0+XDM method. Since these methods share the DFT functional and our primary interest lies in NCIs, the focus will be on the vdW components of the atomic forces.