The Heistracher et al. standard problem tests domain wall pinning at the interface of a two-phase magnetic rod. It is sensitive to discontinuities in the exchange constant $A$, uniaxial anisotropy $K$, and spontaneous polarisation $J_s$ between phases, making it a targeted test for correct interface handling.
P. Heistracher, C. Abert, F. Bruckner, T. Schrefl, D. Suess, J. Magn. Magn. Mater. 548, 168875 (2022).
Problem setup#
A cuboid rod (80 nm × 1 nm × 1 nm) is divided into two equal halves:
- Phase I ($x < 0$): soft or hard material parameters
- Phase II ($x \geq 0$): always hard material
| Parameter | Hard | Soft |
|---|---|---|
| $A$ (J/m) | $1.0 \times 10^{-11}$ | $2.5 \times 10^{-12}$ |
| $K$ (J/m³) | $1.0 \times 10^{6}$ | $1.0 \times 10^{5}$ |
| $J_s$ (T) | 1.00 | 0.25 |
The uniaxial easy axis is along $x$, damping $\alpha = 1$, and demagnetisation is neglected. A linearly increasing external field $\mu_0 H_x(t) = r \cdot t$ with $r = 2 \times 10^7$ T/s is applied. The initial magnetisation creates a domain wall at the interface.
The pinning field $\mu_0 H_p$ is defined as the field at which $\langle m_x \rangle$ exceeds 0.999.
Analytical pinning field#
$$H_p = \frac{2 K^{\mathrm{II}}}{J_s^{\mathrm{II}}} \cdot \frac{1 - \varepsilon_K \varepsilon_A}{\left(1 + \sqrt{\varepsilon_J \varepsilon_A}\right)^2}$$
where $\varepsilon_A = A^{\mathrm{I}}/A^{\mathrm{II}}$, $\varepsilon_K = K^{\mathrm{I}}/K^{\mathrm{II}}$, $\varepsilon_J = J_s^{\mathrm{I}}/J_s^{\mathrm{II}}$, and $H_p$ is in A/m (multiply by $\mu_0$ for Tesla). See Eq. (28) in Heistracher et al.
Test cases#
| Case | Description | $\mu_0 H_p^{\mathrm{an}}$ (T) |
|---|---|---|
| akj | Jump in A, K, Js | 1.568 |
| ak | Jump in A, K | 1.089 |
| aj | Jump in A, Js | 1.206 |
| a | Jump in A only | 0.838 |
| kj | Jump in K, Js | 1.005 |
| k | Jump in K only | 0.565 |
| j | Jump in Js only | 0.000 |
| none | No jump | 0.000 |
Solvers#
| Solver | Type |
|---|---|
| magnum.af | Finite-difference (GPU) |
| OOMMF | Finite-difference (CPU) |
| mumax3.11 | Finite-difference (GPU) |
All solvers use 80 cubic cells with 1 nm edge length.
Numerical pinning fields#
| Case | $\mu_0 H_p^{\mathrm{an}}$ (T) | magnum.af (T) | OOMMF (T) | mumax3.11 (T) |
|---|---|---|---|---|
| akj | 1.568 | 1.585 | 1.585 | 1.585 |
| a | 0.838 | 0.868 | 0.868 | 0.868 |
| aj | 1.206 | 1.256 | 1.256 | 1.256 |
Case akj: jump in A, K, Js#

Case a: jump in A only#

Case aj: jump in A and Js#
