Li cubic spline potential fit to the mlearn Li training set using KLIFF. The training set contained ground-state crystal, strained, and slab configurations, as well as configurations taken from NVT AIMD simulations at 300 K and 0.5×, 0.9×, 1.5×, and 2.0× the melting point, with and without a single vacancy.