Syntax
//geometry
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 // Saturation magnetization (A/m)
Aex = 13e-12 // Exchange stiffness (J/m)
alpha = 0.01 // Landau-Lifshitz damping constant
m = uniform(1, 0, 0)
// set up spacer layer parameters
lambda = 1 // Slonczewski Lambda parameter
Pol = 0.5669 // Electrical current polarization
epsilonprime = 0 // Slonczewski secondary STT term
// set up fixed layer polarization
angle := 20
px := cos(angle * pi/180)
py := sin(angle * pi/180)
fixedlayer = vector(px, py, 0) // Slonczewski fixed layer polarization
// send current
Jtot := -0.008 // total current in A
area := sizeX/2*sizeY/2*pi // R1 x R2 x pi [m2]
jc := Jtot / area // current density in A/m2
J = vector(0, 0, jc)
// schedule output & run
tableautosave(10e-12)
AutoSnapshot(m, 100e-12)
SnapshotFormat = "png"
run(1e-9)
Results
| 0 (ns) |
0.1 |
0.2 |
0.3 |
0.4 |
|
|
|
|
|
| 0.5 |
0.6 |
0.7 |
0.8 |
0.9 |
|
|
|
|
|
