Supplementary Material

Figures depicting the evolution of loss function components during training, the cross-plane thermal conductivity of bulk h-BN calculated by the NEP and DFT, and the SED spectra calculated by MD simulations based on the NEP for h-BN at 300 K, and definitions of hyperparameters for the NEP training.