Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userStable and fast semi-implicit integration of the stochastic Landau-Lifshitz equation
111
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We propose new semi-implicit numerical methods for the integration of the stochastic Landau-Lifshitz equation with built-in angular momentum conservation. The performance of the proposed integrators is tested on the 1D Heisenberg chain. For this system, our schemes show better stability properties and allow us to use considerably larger time steps than standard explicit methods. At the same time, these semi-implicit schemes are also of comparable accuracy to and computationally much cheaper than the standard midpoint implicit method. The results are of key importance for atomistic spin dynamics simulations and the study of spin dynamics beyond the macro spin approximation.
rate research
Read More
We introduce a numerical method to integrate the stochastic Landau-Lifshitz-Gilbert equation in spherical coordinates for generic discretization schemes. This method conserves the magnetization modulus and ensures the approach to equilibrium under the expected conditions. We test the algorithm on a benchmark problem: the dynamics of a uniformly magnetized ellipsoid. We investigate the influence of various parameters, and in particular, we analyze the efficiency of the numerical integration, in terms of the number of steps needed to reach a chosen long time with a given accuracy.
In this paper, we consider homogenization of the Landau-Lifshitz equation with a highly oscillatory material coefficient with period $varepsilon$ modeling a ferromagnetic composite. We derive equations for the homogenized solution to the problem and the corresponding correctors and obtain estimates for the difference between the exact and homogenized solution as well as corrected approximations to the solution. Convergence rates in $varepsilon$ over times $O(varepsilon^sigma)$ with $0leq sigmaleq 2$ are given in the Sobolev norm $H^q$, where $q$ is limited by the regularity of the solution to the detailed Landau-Lifshitz equation and the homogenized equation. The rates depend on $q$, $sigma$ and the the number of correctors.
Recent theoretical and experimental advances show that the inertia of magnetization emerges at sub-picoseconds and contributes to the ultrafast magnetization dynamics which cannot be captured intrinsically by the LLG equation. Therefore, as a generalization, the inertial Landau-Lifshitz-Gilbert (iLLG) equation is proposed to model the ultrafast magnetization dynamics. Mathematically, the LLG equation is a nonlinear system of parabolic type with (possible) degeneracy. However, the iLLG equation is a nonlinear system of mixed hyperbolic-parabolic type with degeneracy, and exhibits more complicated structures. It behaves like a hyperbolic system at the sub-picosecond scale while behaves like a parabolic system at larger timescales. Such hybrid behaviors impose additional difficulties on designing numerical methods for the iLLG equation. In this work, we propose a second-order semi-implicit scheme to solve the iLLG equation. The second temporal derivative of magnetization is approximated by the standard centered difference scheme and the first derivative is approximated by the midpoint scheme involving three time steps. The nonlinear terms are treated semi-implicitly using one-sided interpolation with the second-order accuracy. At each step, the unconditionally unique solvability of the unsymmetric linear system of equations in the proposed method is proved with a detailed discussion on the condition number. Numerically, the second-order accuracy in both time and space is verified. Using the proposed method, the inertial effect of ferromagnetics is observed in micromagnetics simulations at small timescales, in consistency with the hyperbolic property of the model at sub-picoseconds. For long time simulations, the results of the iLLG model are in nice agreements with those of the LLG model, in consistency with the parabolic feature of the iLLG model at larger timescales.
Magnetization dynamics in magnetic materials is modeled by the Landau-Lifshitz-Gilbert (LLG) equation. In the LLG equation, the length of magnetization is conserved and the system energy is dissipative. Implicit and semi-implicit schemes have been used in micromagnetics simulations due to their unconditional numerical stability. In more details, implicit schemes preserve the properties of the LLG equation, but solve a nonlinear system of equations per time step. In contrast, semi-implicit schemes only solve a linear system of equations, while additional operations are needed to preserve the length of magnetization. It still remains unclear which one shall be used if both implicit and semi-implicit schemes are available. In this work, using the implicit Crank-Nicolson (ICN) scheme as a benchmark, we propose to make this implicit scheme semi-implicit. It can be proved that both schemes are second-order accurate in space and time. For the unique solvability of nonlinear systems of equations in the ICN scheme, we require that the temporal step size scales quadratically with the spatial mesh size. It is numerically verified that the convergence of the nonlinear solver becomes slower for larger temporal step size and multiple magnetization profiles are obtained for different initial guesses. The linear systems of equations in the semi-implicit CN (SICN) scheme are unconditionally uniquely solvable, and the condition that the temporal step size scales linearly with the spatial mesh size is needed in the convergence of the SICN scheme. In terms of numerical efficiency, the SICN scheme achieves the same accuracy as the ICN scheme with less computational time. Based on these results, we conclude that a semi-implicit scheme is superior to its implicit analog both theoretically and numerically, and we recommend the semi-implicit scheme in micromagnetics simulations if both methods are available.
In this article, we study the persistence of properties of a given classical deter-ministic dierential equation under a stochastic perturbation of two distinct forms: external and internal. The rst case corresponds to add a noise term to a given equation using the framework of It^o or Stratonovich stochastic dierential equations. The second case corresponds to consider a parameters dependent dierential equations and to add a stochastic dynamics on the parameters using the framework of random ordinary dierential equations. Our main concerns for the preservation of properties is stability/instability of equilibrium points and symplectic/Poisson Hamiltonian structures. We formulate persistence theorem in these two cases and prove that the cases of external and internal stochastic perturbations are drastically dierent. We then apply our results to develop a stochastic version of the Landau-Lifshitz equation. We discuss in particular previous results obtain by Etore and al. in [P. Etore, S.Labbe, J. Lelong, Long time behaviour of a stochastic nanoparticle, J. Differential Equations 257 (2014), 2115-2135] and we nally propose a new family of stochastic Landau-Lifshitz equations.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
