oommf軟體的使用以及微磁學的計算？

# 定義π 4*atan(1.0) 或acos(-1.0)
set pi [expr {4*atan(1.0)}]

# 定義μ0 = 4π×10?7 V·s/(A·m)
set mu0 [expr {4*$pi*1e-7}]

# 每一個oommf的mif文件定義一種磁學構造
# 通過定義Oxs_ext對象描述微磁學問題
# &<1&> 容器
#微磁學模擬的幾何空間是通過容器(Atlases)定義的.這個容器可以是一個連續的空間或者幾個不連續的區域.
#以下定義了一個長方體5nm*5nm*200nm
Specify Oxs_BoxAtlas:atlas {
xrange {0 5e-9}
yrange {0 5e-9}
zrange {0 200e-9}
}

# &<2&>網格
#網格是微磁學模擬離散化的一種形式,
#如下建立方體網格2.5nm*2.5nm*2.5nm,每個格子的長寬高和坐標軸平行,並且遍布整個容器，網格邊長需小於交換長度
Specify Oxs_RectangularMesh:mesh {
cellsize {25e-10 25e-10 25e-10}
atlas :atlas
}

# &<3&> 能量
#為了使所有的能量都被包含到微磁學模擬中加以計算,我們需要把所有的能量項都用specify定義,能量項的數目不受限制
# Oxs_UniformExchange 定義標準的6個最近鄰的交換作用,並計算交換能.容器空間中的交換係數相同
Specify Oxs_UniformExchange [subst {
A 31e-12
}]

# Oxs_UniaxialAnisotropy 定義單軸磁晶各項異性能,
Specify Oxs_UniaxialAnisotropy {
K1 1000e3
axis { Oxs_UniformVectorField {
norm 1
vector {0 0 1}
} }
}

# Oxs_UZeeman 均勻磁場的磁場能
# 一個複合列表，每個子列表由7個元素組成：前3個表示該範圍的開始欄位的x，y和z分量，接下來的3個表示x，y和z分量，結束欄位的範圍，
# x,y,z 對應納米線的長寬高
# 最後一個元素指定磁場變化的（線性）步長。
# 三個序列，從0到正的最大，從正的最大到負的最大，從負的最大到正的最大。這是測磁滯回線的通常加磁場的方法
# multiplier 是係數，乘以0.001變T（特斯拉），再除以＄mu0，就是A/m了

Specify Oxs_UZeeman [subst {
multiplier [expr {0.001/$mu0}]
Hrange {
{ 0 0 0 0 20000 50 200 }
{ 0 20000 50 0 -20000 -50 400 }
{ 0 -20000 -500 0 20000 50 400 }
}
}]

# Oxs_Demag 標準的退磁能項
Specify Oxs_Demag {}

# &<4&>演化器
# alpha即α，阻尼係數
# 演化器主要功能是更新步(step)與步(step)之間的磁化狀態.以步進的方式推動微磁學模擬向前演化
# 有兩種演化器,1.利用LLG動力學過程實現的時域演化; 2.通過直接最小化技術鎖定能量平面的最小值

Specify Oxs_RungeKuttaEvolve:evolve {
alpha 0.5
}

# &<5&>驅動
# 驅動則是再整體上協調演化器,讓演化器完成特定的任務和階段
# Oxs_TimeDriver,控制時間演化器
# basename即計算保存時的文件名
Specify Oxs_TimeDriver [subst {
basename Conanorod1Hy
evolver :evolve
stopping_dm_dt 1
mesh :mesh
Ms 12e5
m0 {0.5 0.5 1}