Accurate calculation of interatomic forces is critical for estimating reliable thermal and thermodynamic properties of crystalline and amorphous materials through large-scale molecular dynamics simulations. First-principles methods based on density functional theory (DFT) offer the most reliable approaches for calculating interatomic interactions, but the increased computational cost has been a limiting factor that hinders the study of realistic systems with defects and anharmonicity at scale. An alternative approach, classical force fields (FFs) constructed with parameterized functions, is a computationally efficient solution, but its accuracy in evaluating interatomic forces is limited compared to DFT.