The values of R min and R cut used for this work are shown in table of the supporting information. The a i parameters were fitted using the Levenberg-Marquardt algorithm to reproduce the DFT energies for each crystal system. The degree m of the polynomial was chosen to be 8 and the c i coefficients were optimized through a least-squares fit.