High-frequency thermal processes in harmonic crystals

We consider two high frequency thermal processes in uniformly heated harmonic crystals relaxing towards equilibrium: (i) equilibration of kinetic and potential energies and (ii) redistribution of energy among spatial directions. Equation describing these processes with deterministic initial conditions is derived. Solution of the equation shows that characteristic time of these processes is of the order of ten periods of atomic vibrations. After that time the system practically reaches the stationary state. It is shown analytically that in harmonic crystals temperature tensor is not isotropic even in the stationary state. As an example, harmonic triangular lattice is considered. Simple formula relating the stationary value of the temperature tensor and initial conditions is derived. The function describing equilibration of kinetic and potential energies is obtained. It is shown that the difference between the energies (Lagrangian) oscillates around zero. Amplitude of these oscillations decays inversely proportional to time. Analytical results are in a good agreement with numerical simulations.