Thermal conductivity of crystalline quartz from classical simulations

We calculate the thermal conductivity of crystalline α- and β-quartz in the high-temperature range (5000.3em0exK to 11000.3em0exK) using nonequilibrium molecular dynamics simulations and an empirical interatomic potential. We find that finite-size effects associated with the nonequilibrium dynamics are not negligible, which implies that reliable results for the bulk thermal conductivity of quartz must be obtained by extrapolation to infinite sizes. The calculated thermal conductivity is nearly temperature independent over a wide range of temperature, in agreement with experiment.