Semiquantum molecular dynamics simulation of thermal properties and heat transport in low-dimensional nanostructures

We present a detailed description of semiquantum molecular dynamics simulation of stochastic dynamics of a system of interacting particles. Within this approach, the dynamics of the system is described with the use of classical Newtonian equations of motion in which the effects of phonon quantum statistics are introduced through random Langevin-like forces with a specific power spectral density (the color noise). The color noise describes the interaction of the molecular system with the thermostat. We apply this technique to the simulation of thermal properties and heat transport in different low-dimensional nanostructures. We describe the determination of temperature in quantum lattice systems, to which the equipartition limit is not applied. We show that one can determine the temperature of such a system from the measured power spectrum and temperature-and relaxation-rate-independent density of vibrational (phonon) states. We simulate the specific heat and heat transport in carbon nanotubes, as well as the heat transport in molecular nanoribbons with perfect (atomically smooth) and rough (porous) edges, and in nanoribbons with strongly anharmonic periodic interatomic potentials. We show that the effects of quantum statistics of phonons are essential for the carbon nanotube in the whole temperature range T < 500 K, in which the values of the specific heat and thermal conductivity of the nanotube are considerably less than that obtained within the description based on classical statistics of phonons. This conclusion is also applicable to other carbon-based materials and systems with high Debye temperature like graphene, graphene nanoribbons, fullerene, diamond, diamond nanowires, etc. We show that the existence of rough edges and quantum statistics of phonons change drastically the low-temperature thermal conductivity of the nanoribbon in comparison with that of the nanoribbon with perfect edges and classical phonon dynamics and statistics. The semiquantum molecular dynamics approach allows us to model the transition in the rough-edge nanoribbons from the thermal-insulator-like behavior at high temperature, when the thermal conductivity decreases with the conductor length, to the ballistic-conductor-like behavior at low temperature, when the thermal conductivity increases with the conductor length. We also show how the combination of strong nonlinearity of periodic interatomic potentials with the quantum statistics of phonons changes completely the low-temperature thermal conductivity of the system.