Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNumerical integration of the stochastic Landau-Lifshitz-Gilbert equation in generic time-discretisation schemes
717
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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.
We investigate in details the inertial dynamics of a uniform magnetization in the ferromagnetic resonance (FMR) context. Analytical predictions and numerical simulations of the complete equations within the Inertial Landau-Lifshitz-Gilbert (ILLG) model are presented. In addition to the usual precession resonance, the inertial model gives a second resonance peak associated to the nutation dynamics provided that the damping is not too large. The analytical resolution of the equations of motion yields both the precession and nutation angular frequencies. They are function of the inertial dynamics characteristic time $tau$, the dimensionless damping $alpha$ and the static magnetic field $H$. A scaling function with respect to $alphataugamma H$ is found for the nutation angular frequency, also valid for the precession angular frequency when $alphataugamma Hgg 1$. Beyond the direct measurement of the nutation resonance peak, we show that the inertial dynamics of the magnetization has measurable effects on both the width and the angular frequency of the precession resonance peak when varying the applied static field. These predictions could be used to experimentally identify the inertial dynamics of the magnetization proposed in the ILLG model.
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 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.
We present a comparison between finite differences schemes and a pseudospectral method applied to the numerical integration of stochastic partial differential equations that model surface growth. We have studied, in 1+1 dimensions, the Kardar, Parisi and Zhang model (KPZ) and the Lai, Das Sarma and Villain model (LDV). The pseudospectral method appears to be the most stable for a given time step for both models. This means that the time up to which we can follow the temporal evolution of a given system is larger for the pseudospectral method. Moreover, for the KPZ model, a pseudospectral scheme gives results closer to the predictions of the continuum model than those obtained through finite difference methods. On the other hand, some numerical instabilities appearing with finite difference methods for the LDV model are absent when a pseudospectral integration is performed. These numerical instabilities give rise to an approximate multiscaling observed in the numerical simulations. With the pseudospectral approach no multiscaling is seen in agreement with the continuum model.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
