Equilibration of energies in a two-dimensional harmonic graphene lattice

We study dynamical phenomena in a harmonic graphene (honeycomb) lattice, consisting of equal particles connected by linear and angular springs. Equations of in-plane motion for the lattice are derived. Initial conditions typical for molecular dynamic modelling are considered. Particles have random initial velocities and zero displacements. In this case, the lattice is far from thermal equilibrium. In particular, initial kinetic and potential energies are not equal. Moreover, initial kinetic energies (and temperatures), corresponding to degrees of freedom of the unit cell, are generally different. The motion of particles leads to equilibration of kinetic and potential energies and redistribution of kinetic energy among degrees of freedom. During equilibration, the kinetic energy performs decaying high-frequency oscillations. We show that these oscillations are accurately described by an integral depending on dispersion relation and polarization matrix of the lattice. At large times, kinetic and potential energies tend to equal values. Kinetic energy is partially redistributed among degrees of freedom of the unit cell. Equilibrium distribution of the kinetic energies is accurately predicted by the non-equipartition theorem. Presented results may serve for better understanding of the approach to thermal equilibrium in graphene. This article is part of the theme issue ‘Modelling of dynamic phenomena and localization in structured media (part 2)’.