For the hematite surface, we benchmarked four different force fields and again observed that, whilst all force fields could accurately predict the density, the predicted thermal conductivity varied significantly from experimental results. The most accurate force field we tested was the Morse potential. The thermal conductivity calculated using this model was in excellent agreement with experimental data obtained in all three crystallographic directions.