Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSolving the stochastic Landau-Lifshitz-Gilbert-Slonczewski equation for monodomain nanomagnets : A survey and analysis of numerical techniques
71
0
0.0
(
0
)
Added by
Nikhil Rangarajan
Publication date
2016
fields
Informatics Engineering
and research's language is
English
Authors
Sebastian Ament
Ask ChatGPT about the research
No Arabic abstract
The stochastic Landau-Lifshitz-Gilbert-Slonczewski (s-LLGS) equation is widely used to study the temporal evolution of the macrospin subject to spin torque and thermal noise. The numerical simulation of the s-LLGS equation requires an appropriate choice of stochastic calculus and numerical integration scheme. In this paper, we comprehensively evaluate the accuracy and complexity of various numerical techniques to solve the s-LLGS equation. We focus on implicit midpoint, Heun, and Euler-Heun methods that converge to the Stratonovich solution of the s-LLGS equation. By performing numerical tests for both strong (path-wise) and weak (statistical) convergence, we quantify the accuracy of various numerical schemes used to solve the s-LLGS equation. We demonstrate a new method intended to solve Stochastic Differential Equations (SDEs) with small noise (RK4-Heun), and test its capability to handle the s-LLGS equation. We also discuss the circuit implementation of nanomagnets for large-scale SPICE-based simulations. We evaluate the efficacy of SPICE in handling the stochastic dynamics of the multiplicative noise in the s-LLGS equation. Numerical schemes such as Euler and Gear, typically used by SPICE-based circuit simulators do not yield the expected outcome when solving the Stratonovich s-LLGS equation. While the trapezoidal method in SPICE does solve for the Stratonovich solution, its accuracy is limited by the minimum time step of integration in SPICE. We implement the s-LLGS equation in both its cartesian and spherical coordinates form in SPICE and compare the stability and accuracy of the two implementations. The results in this paper will serve as guidelines for researchers to understand the tradeoffs between accuracy and complexity of various numerical methods and the choice of appropriate calculus to solve the s-LLGS equation.
rate research
Read More
We consider the numerical approximation of the inertial Landau-Lifshitz-Gilbert (iLLG) equation, which describes the dynamics of the magnetization in ferromagnetic materials at subpicosecond time scales. We propose and analyze two fully discrete numerical schemes based on two different approaches: The first method is based on a reformulation of the problem as a linear constrained variational formulation for the time derivative of the magnetization. The second method exploits a reformulation of the problem as a first order system in time for the magnetization and the angular momentum. Both schemes are implicit, based on first-order finite elements, and the constructed numerical approximations satisfy the inherent unit-length constraint of iLLG at the vertices of the underlying mesh. For both schemes, we establish a discrete energy law and prove convergence of the approximations towards a weak solution of the problem. Numerical experiments validate the theoretical results and show the applicability of the methods for the simulation of ultrafast magnetic processes.
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.
We construct a path-integral representation of the generating functional for the dissipative dynamics of a classical magnetic moment as described by the stochastic generalization of the Landau-Lifshitz-Gilbert equation proposed by Brown, with the possible addition of spin-torque terms. In the process of constructing this functional in the Cartesian coordinate system, we critically revisit this stochastic equation. We present it in a form that accommodates for any discretization scheme thanks to the inclusion of a drift term. The generalized equation ensures the conservation of the magnetization modulus and the approach to the Gibbs-Boltzmann equilibrium in the absence of non-potential and time-dependent forces. The drift term vanishes only if the mid-point Stratonovich prescription is used. We next reset the problem in the more natural spherical coordinate system. We show that the noise transforms non-trivially to spherical coordinates acquiring a non-vanishing mean value in this coordinate system, a fact that has been often overlooked in the literature. We next construct the generating functional formalism in this system of coordinates for any discretization prescription. The functional formalism in Cartesian or spherical coordinates should serve as a starting point to study different aspects of the out-of-equilibrium dynamics of magnets. Extensions to colored noise, micro-magnetism and disordered problems are straightforward.
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.
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
