How accurate are DFT forces? Unexpectedly large uncertainties in molecular datasets

Training of general-purpose machine learning interatomic potentials (MLIPs) relies on large datasets with properties usually computed with density functional theory (DFT). A prerequisite for accurate MLIPs is that the DFT data are well converged to minimize numerical errors. A possible symptom of errors in DFT force components is nonzero net force. Here, we consider net forces in datasets including SPICE, Transition1x, ANI-1x, ANI-1xbb, AIMNet2, QCML, and OMol25. Several of these datasets suffer from significant nonzero DFT net forces. We also quantify individual force component errors by comparison to recomputed forces using more reliable DFT settings at the same level of theory, and we find significant discrepancies in force components averaging from 1.7 meV/Å in the SPICE dataset to 33.2 meV/Å in the ANI-1x dataset. These findings underscore the importance of well converged DFT data as increasingly accurate MLIP architectures become available.