Our workflow for developing GAP models and calculating thermal properties is depicted in Figure . The procedure for calculating thermal properties of 2D materials is similar to our earlier studies as shown on the right side of Figure . The key difference here is that we use the thermal property calculations to tune the hyperparameters of MLIP training and to validate the accuracy of the generated potentials.