Computer-Aided Series Expansion for Phonon Self-Energy

In this paper we discuss an algorithm for calculation of the temperature-dependent anharmonic correction to the phonon spectrum in atomic and molecular crystals. We show how the equation of motion method can be used to compute corrections of arbitrary perturbation order to the phonon self-energy. Complete analytical expressions up to order λ6 are obtained as an application of the method.