Stoner–Wohlfarth model

link: https://en.wikipedia.org/wiki/Stoner%E2%80%93Wohlfarth_model

When the size of the studied system is below the exchange length, the magnetization of the system is uniform and thus can be described using a single 3d vector. In this situation, the demagnetization field can be calculated by simply multiplying the demagnetization tensor and the magnetization. Therefore, the demagnetization energy is equivalent to an effective anisotropy. Assuming that the external field $\overrightarrow{H}=(H,0,0)$, easy axis $\hat{u}=(\cos \theta ,\sin \theta ,0)$ and unit magnetization vector $\hat{m}=(\cos \varphi ,\sin \varphi ,0)$, we have

$$\begin{array}{rl}E& =-K(\overrightarrow{m}\cdot \hat{u}{)}^2-{\mu }_0{M}_{s}H\cos \varphi \\ & =-\frac{K}{2}[1+\cos (2(\theta -\varphi ))+4h\cos \varphi ]\end{array}$$

where $h=H/H_k$ and $H_k=2K/({\mu }_0{M}_s)$. In the equilibrium state, the first derivative of the energy with respect to the magnetization direction is zero, i.e.,

$$\frac{\partial E}{\partial \varphi }=K[-\sin (2(\theta -\varphi ))+2h\sin \varphi ]=0$$

In principle, there will be a solution of $\varphi$ for each given $h$ and $\theta$. The hysteresis loops can be constructed by plotting $\cos \varphi$ as a function of $h$. The switching field can be obtained by extra setting the second derivative of the energy with respect to the magnetization direction to zero, i.e.,

$$\frac{{\partial }^{2}E}{\partial {\varphi }^2}=2K[\cos (2(\theta -\varphi ))+h\cos \varphi ]=0$$

The obtained switching field is

$$h_s=\frac{{(1-t^2+t^4)}^{1/2}}{1+t^2}$$

where $t={\tan }^{1/3}\theta$. Specifically, $h_s=1/2$ if $\theta =\pi /4$. In this example, we use MicroMagnetic to demostrate this result.

We chose the system to be a cubic sample, so the demagnetization tensor is $N_x=N_y=N_z=1/3$. That is, the demagnetization itself does not contribute to the effective anisotropy. So in this simulation, we have ignored the demagnetization field. The MicroMagnetic script is shown below:

using MicroMagnetic

We create a mesh for a cubic geometry 4nm x 4nm x 4nm

mesh = FDMesh(; nx=4, ny=4, nz=4, dx=1e-9, dy=1e-9, dz=1e-9);

We create a simulation instance with the LLG driver and the BS23 integrator.

sim = Sim(mesh, driver="LLG", integrator="BS23", name="SW")

sim.driver.alpha = 0.5
sim.driver.integrator.tol = 1e-7

set_Ms(sim, 1e6)
init_m0(sim, (-1,1,0))

add_exch(sim, 1.3e-11)
add_anis(sim, 5e4, axis=(1, 1, 0))
add_zeeman(sim, (0,0,0))

[ Info: MicroSim (FD) has been created.
[ Info: Exchange has been added.
[ Info: Uniaxial Anisotropy has been added.
[ Info: Static Zeeman has been added.

We use the hysteresis function to simulate the hysteresis loops.

Hs = [i*mT for i=-100:5:100]
hysteresis(sim, Hs, direction=(1,0,0), full_loop=false, stopping_dmdt=0.05, output="vts")

[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=439, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=87, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=88, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=90, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=91, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=93, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=94, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=96, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=98, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=99, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=101, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=103, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=105, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=107, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=109, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=111, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=113, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=115, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=117, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=120, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=122, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=125, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=129, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=133, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=138, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=144, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=150, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=158, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=168, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=180, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=782, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=102, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=100, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=98, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=96, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=95, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=93, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=92, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=90, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=88, Done!
[ Info: Running Driver : MicroMagnetic.LLG{Float64}.
[ Info: max_dmdt is less than stopping_dmdt=0.05 @steps=87, Done!

For the used anisotropy $K_u=5\times {10}^4$ A/m${}^3$, the expected switch field is $H_c=(1/2)H_K=39788.7$ A/m. We plot the hysteresis loops using the plot_ts function

using CairoMakie
fig = plot_ts("SW_llg.txt", x_key="Hx", ["m_x", "m_y"], x_unit=1/mT, xlabel="H (mT)", ylabel="m", mirror_loop=true);