neel relaxation of 8 nM spherical particle code example

[darthboy2]

Seems to work. This is the region 2 demag z component at the center of the sphereSee anything wrong in the code? 

Are the mumax parameters directly usable in this equation?

Can I use ku1 as K in other words?

I ask this because I think I should see more switching than I do for a 10 nM particle (not shown).

Is the cell size small enough? The dt?

// Set meshgridsand geometry
mygrid:=16
mycellsize:=1e-9
mydiameter:=8e-9
d:=mydiameter

SetGridsize(mygrid,mygrid,mygrid)
SetCellsize(mycellsize,mycellsize,mycellsize)

EdgeSmooth = 8

SetGeom(Ellipsoid(mydiameter,mydiameter,mydiameter))

defregion(1,Ellipsoid(mydiameter, mydiameter, mydiameter))
defregioncell(2,mygrid/2,mygrid/2,mygrid/2)

//saveas(regions,"regions")

//tableAdd(m.Region(1))  
//tableAdd(m.Region(2))  
tableAdd(B_demag.Region(2))  

// Defineeffectivesaturation magnetizationand anisotropy
Mdef:=4.76e5*(1-2*0.68e-9/d) *(1-2*0.68e-9/d) * (1-2*0.68e-9/d)
Kdef:=1.04e4 + (6/d) * 3.9e-5
// Set magnetic parameters

//Msat=Mdef
Msat.SetRegion(1, Mdef)
Msat.SetRegion(2, Mdef)

//Aex=1.3e-11
Aex.SetRegion(1, 1.3e-11)
Aex.SetRegion(2, 1.3e-11)

//Alpha=0.3
Alpha.SetRegion(1, .3)
Alpha.SetRegion(2, .3)

//Ku1=Kdef
Ku1.SetRegion(1, Kdef)
Ku1.SetRegion(2, Kdef)

myanisU:=vector(0, 0, 1)
anisU.SetRegion(1, myanisU)
anisU.SetRegion(2, myanisU)

// Set initialmagnetization and run
//m=RandomMag()
//saveas(m,"m_random")
//saveas(B_demag.region(1),"demag_random")
//tablesave()

//B_ext = vector(0,0,0)
m=uniform(0, 0, 1)
//relax( )
tablesave()
save(m)
//saveas(m,"m_post-relax")
//saveas(B_demag.region(1),"demag_post_relax")
//exit()

Temp=300
fixdt=1e-14

//B_ext = vector(0.0,0.0,0.0)

tableautosave(100e-12)

autosave(m,100e-12)

run(10.0e-9)