Dear All,
I hope you all are doing well. I am simulating TDW in nanowire having two symmetrical square shapes notched at the middle of the wire. When i am applying current density(j) both notch and TDW are moving.To get DW velocity I want the DW position at each point of the nanowire .But here both Notch and DW are moving. Please help me where i am doing wrong.
setGridSize(500, 50, 6)
setCellSize(3e-9, 3e-9, 2e-9)
Msat = 800e3
Aex = 13e-12
Xi = 0.01
alpha = 0.05
m = twodomain(1,0,0, 0,1,0, -1,0,0).transl(-650e-9, 0, 0 )
snapshot(m)
notches := rect(50e-9, 50e-9).Repeat(0, 50*3e-9, 0).Transl(0, 25*3e-9, 0)
setGeom(notches.inverse())
snapshot(geom)
// Remove surface charges from left (mx=1) and right (mx=-1) sides to mimic infinitely long wire. We have to specify the region (0) at the boundaries.
BoundaryRegion := 0
MagLeft := 1
MagRight := -1
ext_rmSurfaceCharge(BoundaryRegion, MagLeft, MagRight)
relax()
ext_centerWall(0) // keep m[0] (m_x) close to zero
// Schedule output
autoSnapshot(m, 50e-12)
tableadd(ext_dwpos) // domain wall position
table autosave(10e-12)
// Run the simulation with current through the sample
pol = 0.56
J = vector(-10e12, 0, 0)
Run(10e-9)