In practice, one of the most accurate routes to theoretically describing lattice vibration at the level of densityfunctional theory (DFT) is ab initio molecular dynamics (AIMD). Correction for quantum nuclear motion may further be done using path-integral molecular dynamics . Due to the immense computational cost of running AIMD, however, it is unwieldy to be applied in high-throughput for a large-scale coverage of the material space. Molecular dynamics using machine-learned force fields have emerged in the past few years and have demonstrated near-AIMD accuracy on several occasions . However, force-field training itself still requires a considerable amount of ab initio calculations (at least hundreds of trajectory structures). A much more cost-effective way to describe lattice dynamics is in terms of interatomic force constants (IFC). They are defined within the Taylor expansion of the total energy with respect to atomic displacements (u), whose derivative yields interatomic forces (F a i ) of atom a (at its equilibrium site) in direction i as: