Figure and Table demonstrate our calculated thermal conductivity in this workflow is in remarkable agreement with the experimental values (R 2 = 0.98 for LTC data). The benchmarking results for both CTE and LTC as downstream outputs validate the accuracy of the anharmonic force constants produced in our workflow using the selected parameters.