Skip to main content

Workshop 2025: High order finite difference method for exchange field computation

Session Information

Location: Lecture room F3213 - 18
Day: Wednesday, 14 May
Time: 18:00 - 19:00
Chairperson: -

Presentation Details

Presentation Type: Poster presentation
Title: High order finite difference method for exchange field computation
Abstract: Many micromagnetic simulation codes employ a nearest neighbour approach with the central finite difference method to compute magnetization derivatives for the exchange field and DMI. While effective for obtaining results, this method requires a finer mesh to achieve a specific accuracy as compared to higher order finite difference approaches. In addition, for time integration this approach needs smaller time steps for accurate time domain ù evaluation.
The Petaspin [https://www.petaspin.com] [A. Giordano et al., J. Appl. Phys. 111, 07D112 (2012)] solver has been already used to describe several experiments and predict new phenomena in spintronics and micromagnetics. It’s aimed at simulating the spin dynamic of ferromagnetic materials by numerically solving the Landau-Lifshitz-Gilbert (LLG) equation. It incorporates various energy terms influencing the material's effective magnetic field, such as exchange interaction, material anisotropy, the applied external fields, the demagnetizing field and magnetoelastic interactions. The user can also choose between different time integrators, like the 5th-order Runge-Kutta method method, the Heun method and the Adams2-Bashforth3 method.
In our latest development , Petaspin now supports high order (4th) neighbour central finite difference method for derivative computation in order to consider high order exchange and DMI interactions allowing the study of frustration and stabilization of three-dimensional solitons such as hopfions. From a numerical point of view, this advancement allows us to use fewer discretization cells and larger time steps while. All computations are natively optimized with NVIDIA's CUDA platform to fully harness GPU parallel computing capabilities, offering superior performance compared to CPU-based solvers.
Here we will show results achieved in the study of the static properties of spin textures, like bubbles, skyrmions, hopfions and vortices, and on new dynamical properties exploiting high order effects.

Presenter

Mr Edoardo Piccolo
Politecnico di Bari, Bari, Italy | Italy

Authors

1. Piccolo Edoardo | Department of Electrical and Information Engineering, Politecnico di Bari, 70125 Bari, Italy
2. Hasan Ali | Department of Mathematical and Computer Sciences, Physical Sciences and Earth Sciences, University of Messina, I-98166 Messina, Italy
3. Giordano Anna | Department of Engineering, University of Messina, I-98166 Messina, Italy
4. Carpentieri Mario | Department of Electrical and Information Engineering, Politecnico di Bari, 70125 Bari, Italy
5. Finocchio Giovanni | Department of Mathematical and Computer Sciences, Physical Sciences and Earth Sciences, University of Messina, I-98166 Messina, Italy