STT随电流变化的磁化反转

sizeX := 160e-9
sizeY := 80e-9
sizeZ := 5e-9

Nx := 64
Ny := 32
 
setgridsize(Nx, Ny, 1)
setcellsize(sizeX/Nx, sizeY/Ny, sizeZ)
setGeom(ellipse(sizeX, sizeY))

// set up free layer
Msat  = 800e3
Aex   = 13e-12
alpha = 0.01
m     = uniform(1, 0, 0)

// set up spacer layer parameters
lambda       = 1
Pol          = 0.5669
epsilonprime = 0

// set up fixed layer polarization
angle := 20
px := cos(angle * pi/180)
py := sin(angle * pi/180)
fixedlayer = vector(px, py, 0)


area := sizeX*sizeY*pi/4

Jstep:=1.0e-4

tableadd(J)
for Jtot:=0.0;Jtot>=-0.008;Jtot=Jtot-Jstep{
      jc:=Jtot/area
      J=vector(0,0,jc)
      minimize()
      tablesave()
}
for Jtot:=-0.008;Jtot<=0.008;Jtot=Jtot+Jstep{
      jc:=Jtot/area
      J=vector(0,0,jc)
      minimize()
      tablesave()
}
for Jtot:=0.008;Jtot>=0.0;Jtot=Jtot-Jstep{
      jc:=Jtot/area
      J=vector(0,0,jc)
      minimize()
      tablesave()
}