Abstract:We consider the numerical approximation of a continuum model of antiferromagnetic and ferrimagnetic materials. The state of the material is described in terms of two unit-length vector fields, which can be interpreted as the magnetizations averaging the spins of two sublattices. For the static setting, which requires the solution of a constrained energy minimization problem, we introduce a discretization based on first-order finite elements and prove its $\Gamma$-convergence. Then, we propose and analyze two iterative algorithms for the computation of low-energy stationary points. The algorithms are obtained from (semi-)implicit time discretizations of gradient flows of the energy. Finally, we extend the algorithms to the dynamic setting, which consists of a nonlinear system of two Landau-Lifshitz-Gilbert equations solved by the two fields, and we prove unconditional stability and convergence of the finite element approximations toward a weak solution of the problem. Numerical experiments assess the performance of the algorithms and demonstrate their applicability for the simulation of physical processes involving antiferromagnetic and ferrimagnetic materials.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop an effective numerical method to simulate static and dynamic magnetic processes in antiferromagnetic (AFM) and ferrimagnetic (FiM) materials. Specifically, the paper proposes a finite - element method for solving the continuous models of the magnetization states in these materials and proves the mathematical validity of these methods.
### Background of the Paper
Antiferromagnetic and ferrimagnetic materials show great potential in improving the functions of spintronic devices due to their unique magnetic properties. The magnetic moments of these materials tend to be arranged in an antiparallel manner, which is in contrast to the parallel arrangement in ferromagnetic materials. However, due to the complexity of AFM and FiM materials, their numerical simulations face many challenges.
### Main Contributions
1. **Mathematical Model**:
- The paper proposes a continuous model for describing the magnetization states of AFM and FiM materials. This model uses two vector fields \( \mathbf{m}_1 \) and \( \mathbf{m}_2 \) with unit length to represent the average magnetization of two sub - lattices.
- The energy functional \( E[\mathbf{m}_1,\mathbf{m}_2] \) contains multiple terms, and each term represents a specific physical effect, such as inhomogeneous lattice exchange energy, inhomogeneous inter - lattice exchange energy and homogeneous inter - lattice exchange energy.
2. **Static Problem**:
- The goal of the static problem is to minimize the energy functional \( E[\mathbf{m}_1,\mathbf{m}_2] \). The paper proposes a discretization method based on the first - order finite element and proves that the discrete energy functional converges to the continuous energy functional in the sense of Γ - convergence.
- Two iterative algorithms (Algorithm 4.4 and Algorithm 4.5) are proposed, which are respectively based on the time discretization of the gradient flow. Both of these algorithms can ensure the monotonic decrease of energy during the iteration process and can converge to the stable points of energy under certain conditions.
3. **Dynamic Problem**:
- The dynamic problem is described by two coupled Landau - Lifshitz - Gilbert (LLG) equations. The paper extends the algorithm with the best performance (Algorithm 5.1) to handle the dynamic problem.
- It is proved that this algorithm is well - posed, stable, and the approximate solutions generated unconditionally converge to the weak solutions of the coupled LLG equations.
### Numerical Experiments
The paper verifies the effectiveness of the proposed algorithms through numerical experiments. The experimental results show that the algorithms can find the minimum value of energy under different gradient flow metrics and have good performance.
### Conclusion
This paper provides a mathematically rigorous and effective framework for the numerical simulation of AFM and FiM materials. By introducing and analyzing new numerical methods, the paper fills the research gap in this field and provides a theoretical basis for future applications.