Seeking Guidance on Reproducing Skyrmion Spin Textures

[rakesh maurya]

Hello Mumax community,

I am trying to reproduce Figure -1 for the three reported cases as shown in the paper. 
Could you please help me understand what part might be missing or incorrect in my .txt input file? By incorporating the necessary corrections, I hope to reproduce the same results.

// ---------- Mesh ----------
a := 0.5e-9       //  lattice parameter
Nx := 288; Ny := 288; Nz := 1
dx := a; dy := a; dz := a
Lx := Nx*dx; Ly := Ny*dy; Lz := Nz*dz
SetMesh(Nx, Ny, Nz, dx, dy, dz, 0,0,0)
SetPBC(1,1,0)

// ---------- Material ----------
J_ex := 1e-3 * 1.602e-19
D_over_J := 0.18
M_spin := 1.0
Aex = J_ex * pow(M_spin,2)/a
Dind = D_over_J * J_ex / a
mu_B := 9.274e-24
Msat = mu_B / pow(a,3)

// ---------- External field ----------
B_over_J := 0.02
B_ext = vector(0,0,B_over_J * J_ex / mu_B)

// ---------- Damping & initial magnetization ----------
alpha = 0.04
m = randomMag()


// ---------- Relax / run ----------
relax()
run(1e-12)
save(m)