The observed difference in force directions has potential repercussions for MLFF methods, where molecular data is routinely optimized at one level of theory, method or functional or dispersion correction, and computed at a different one to serve as input or reference for energies and forces. Unfortunately, while for E int we have highly accurate LNO-CCSD(T) reference data and even the ability to achieve a "platinum standard" confidence by comparing with QMC, obtaining forces at benchmark ab initio level is prohibitively expensive. Current research in developing gradients for LNO-CCSD(T) and QMC could facilitate a reference benchmark for the vdW forces in the future and provide a clearer picture for the accuracy of the DFT methods for MD of ligands binding to a pocket.