“…Among analytical and numerical methods available to solve the BTE [1,2,4,23,26,28], Monte Carlo (MC) simulations are proven to be the most flexible and accurate, yet they can be slow and expensive in terms of computational resources depending on the levels of physics included in the simulation. Many researchers have used MC simulations for phonon transport at nanoscales because of its flexibility in accounting complicated geometries and the correct phonon dispersion relation and different polarization branches [14,[16][17][18][19][20]25,[29][30][31]. While the simulation has been successfully used for predicting thermal conductivities of nanostructures such as nanowires and nanofilms [14,[18][19][20][31][32][33][34], there is plenty of room for improvement in the algorithm, especially for treating the phonon-phonon scattering mechanisms.…”