Midpoint numerical technique for stochastic Landau-Lifshitz-Gilbert dynamics

The implicit midpoint time-integration technique is applied to the stochastic Landau-Lifshitz-Gilbert (LLG) equation. The numerical scheme converges to the Stratonovich solution in the limit of vanishing time step. It preserves the magnetization magnitude and the main energy balance properties of the LLG equation independently of the time step. The numerical technique is then applied to the study of superparamagnetic state in a small spheroidal particle, and the numerical results are compared with the theory.