Accessing the thermal conductivities of Sb 2 Te 3 and Bi 2 Te 3 /Sb 2 Te 3 superlattices by molecular dynamics simulations with a deep neural network potential

Phonon thermal transport is a key feature for the operation of thermoelectric materials, but it is challenging to accurately calculate the thermal conductivity of materials with strong anharmonicity or large cells. In this work, a deep neural network potential (NNP) is developed using a dataset based on density functional theory (DFT) and applied to describe the lattice dynamics of Sb Te and Bi Te /Sb Te superlattices. The lattice thermal conductivities of Sb Te are first predicted using equilibrium molecular dynamics (EMD) simulations combined with an NNP and the results match well with experimental values. Then, through further exploration of weighted phase spaces and the Grüneisen parameter, we find that there is a stronger anharmonicity in the out-of-plane direction in Sb Te , which is the reason why the thermal conductivities are overestimated more in the out-of-plane direction than in the in-plane direction by solving the phonon Boltzmann transport equation (BTE) with only three-phonon scattering processes being considered. More importantly, the lattice thermal conductivities of Bi Te /Sb Te superlattices with different periods are accurately predicted using non-equilibrium molecular dynamics (NEMD) simulations together with an NNP, which serves as a good example to explore the thermal transport physics of superlattices using a deep neural network potential.