Numerical technique for integration of the Landau–Lifshitz equation

In the article, a finite difference scheme for the numerical integration of the Landau–Lifshitz equation is presented. The scheme is based on the application of the midpoint rule coupled with a suitable extrapolation formula. The important properties of the scheme are the conservation of magnetization magnitude, its linearity, its second order truncation error accuracy, and the absence of spatial coupling. The accuracy of the scheme has been extensively tested by comparing numerical solutions with exact analytical solutions and by applying the scheme to the analysis of magnetization dynamics in conducting thin films.