Towards More Accurate Molecular Dynamics Calculation of Thermal Conductivity. Case Study: GaN Bulk Crystals

Abstract

Significant differences exist among literature for thermal conductivity of various systems computed using molecular dynamics simulation. In some cases, unphysical results, for example, negative thermal conductivity, have been found. Using GaN as an example case and the direct non-equilibrium method, extensive molecular dynamics simulations and Monte Carlo analysis of the results have been carried out to quantify the uncertainty level of the molecular dynamics methods and to identify the conditions that can yield sufficiently accurate calculations of thermal conductivity. We found that the errors of the calculations are mainly due to the statistical thermal fluctuations. Extrapolating results to the limit of an infinite-size system tend to magnify the errors and occasionally lead to unphysical results. The error in bulk estimates can be reduced by performing longer time averages using properly selected systems over a range of sample lengths. If the errors in the conductivity estimates associated with each of the sample lengths are kept below a certain threshold, the likelihood of obtaining unphysical bulk values becomes insignificant. Using a Monte-Carlo approach developed here, we have determined the probability distributions for the bulk thermal conductivities obtained using the direct method. We also have observed a nonlinear effect that can become a source of significant errors. For the extremely accurate results presented here, we predict a [0001] GaN thermal conductivity of 185 $W/K\cdot \mathrm{m}$ at 300 K, 102 $W/K\cdot \mathrm{m}$ at 500 K, and 74 $W/K\cdot \mathrm{m}$ at 800 K. Using the insights obtained in the work, we have achieved a corresponding error level (standard deviation) for the bulk (infinite sample length) GaN thermal conductivity of less than 10 $W/K\cdot \mathrm{m}$, 5 $W/K\cdot \mathrm{m}$, and 15 $W/K\cdot \mathrm{m}$ at 300 K, 500 K, and 800 K respectively.

References in corpus (3)

• Molecular dynamics calculation of the thermal conductivity of vitreous silica
• Simulation of thermal conductivity and heat transport in solids
• Heat capacity of $\alpha$-GaN: Isotope Effects

Cited by in corpus (20)

• Generalization of Fourier's law into viscous heat equations
• Fast and slow thermal processes in harmonic scalar lattices
• Boltzmann Transport in Nanostructures as a Friction Effect
• Impact of Screw and Edge Dislocation on the Thermal Conductivity of Nanowires and Bulk GaN
• Size Effects in Thermal Conduction by Phonons
• Effects of Cutoff Functions of Tersoff Potentials on Molecular Dynamics Simulations of Thermal Transport
• Ballistic thermophoresis of adsorbates on free-standing graphene
• Molecular Dynamics Prediction of Thermal Conductivity of GaN Films and Wires at Realistic Length Scales
• Temperature in a Peierls-Boltzmann Treatment of Nonlocal Phonon Heat Transport
• Transport Coefficients from Large Deviation Functions
• An analytical law for size effects on thermal conductivity of nanostructures
• Decorated dislocations against phonon propagation for thermal management
• Effects of nano-void density, size, and spatial population on thermal conductivity: a case study of GaN crystal
• Analysis of nonlocal phonon thermal conductivity simulations showing the ballistic to diffusive crossover
• Strong Reduction of Thermal Conductivity of WSe2 with Introduction of Atomic Defects
• Thermal Boundary Characteristics of Homo-/Heterogeneous Interfaces
• Non-local linear-response functions for thermal transport computed with equilibrium molecular-dynamics simulation
• Size and Stoichiometric Dependence of Thermal Conductivities of InxGa1-xN: A Molecular Dynamics Study
• The influence of defects on the thermal conductivity of compressed LiF
• Vibronics of multi-material nanopillared membranes and impact on the thermal conductivity