An Efficient Energy Stable Structure Preserving Method for The Landau-Lifshitz Equation

TL;DR

This paper introduces a new first-order numerical method for simulating the Landau-Lifshitz equation that maintains the magnetization norm, enhances stability, and avoids complex coupled system solutions.

Contribution

It proposes a structure-preserving, first-order numerical scheme combining Gauss-Seidel, double diffusion, and Crank-Nicolson iterations for improved micromagnetic simulations.

Findings

Ensures norm preservation at discrete and continuous levels.

Provides stability without severe time step restrictions.

Avoids solving complex coupled systems.

Abstract

One of the main difficulties in micromagnetics simulation is the norm preserving constraints $\|\mathbf{m}\|=1$ at the continuous or the discrete level. Another difficulty is the stability with the time step constraint. Using standard explicit integrators leads to a physical time step of sub-pico seconds, which is often two orders of magnitude smaller than the fastest physical time scales. Direct implicit integrators require solving complicated, coupled systems. Another major difficulty with the projection method in this field is the lack of rigorous theoretical guarantees regarding its stability of the projection step. In this paper, we introduce a first order method. Such a method is structure preserving based on a combination of a Gauss-Seidel iteration, a double diffusion iteration and a Crank-Nicolson iteration to preserve the norm constraints.