A recent article  presents a semianalytical method to solve the Landau Lifshitz (LL) equation. Spin motion is computed analytically as precession about the effective field H, where H is assumed fixed over the time step. However, the exchange field dominates at short range and varies at the time scale of neighbor spin precessions, undermining the fixed field assumption. We present an axis corrected version of this algorithm. We add a scalar multiple of m to H (preserving torque and hence the LL solution) to produce a more stable precession axis parallel to the cross product of the torques m H at two closely spaced time steps. We build a predictorcorrector solver on this foundation. The second order convergence of the solver enables calculation of adjustable time steps to meets a desired error magnitude.
Citation: Journal of Applied Physics
Pub Type: Journals