The reference QM calculations were performed using ωb97xD/6-31G*, as implemented in Gaussian 16. To avoid unphysical structures that could result in not converged ISA-Pol calculations, we incorporated an structural filter within the active learning workflow. This filter eliminates structures with bond lengths less than 0.8 or more than 1.4 times their equilibrium values.