Micromagnetic calculation of domain structure in thin magnetic film based on improved LaBonte method

A numerical method is developed to calculate equilibrium domain structures micromagnetically. The method is based on the LaBonte´s iteration which is improved in such a way that the demagnetizing field is calculated simultaneously in advance of each iteration cycle and that a torque term which represents the gyroscopic nature of the magnetization is included in the update rule. The former modification allows the use of the fast Fourier transform (FFT) to calculate the demagnetizing field in the LaBonte method. The torque term is added to prevent the LaBonte´s iteration from converging to a false solution different from that derived using the Landau-Lifshitz-Gilbert equation. The principle of calculation and the results of two-dimensional calculations which demonstrate the effectiveness of the developed method are presented