计算Fe晶格常数

# in.Fe Stefan Bringuier

#设置jump执行入口标志

label lata

#设置晶格常数a的不同取值

variable a index 2.4 2.5 2.6 2.7 2.8 2.9 3.0 3.1 3.2

#模型基本设置

units    metal                       

atom_style  atomic                      

boundary        p p p                       

#以变量a的值创建盒子并填充原子

lattice    bcc $a                      

region    box block 0 5 0 5 0 5      

create_box  1 box                       

create_atoms  1 box                      

mass    1 55.85                       

#设置力场

pair_style eam/fs                          

pair_coeff * * Fe_mm.eam.fs Fe              

#设置邻居列表参数

neighbor  0.3 bin                     

neigh_modify  every 20 delay 0 check no

#输出轨迹文件

dump myDump all atom 10000 dump.${a}

#nve系综

timestep 0.001

fix    1 all nve                      

thermo    500                           

thermo_style custom step pe ke temp vol press 

run  10000                          

#计算单原子能量=总势能/原子个数

variable Ecoh equal pe/count(all)  

print "Lattice parameter (A): $a" 

print "Cohesive Energy (eV/atom):  ${Ecoh}"

#清屏

clear

#获取下一个a的值

next a

#跳转到in.Fe文件lata处重新运行

jump in.Fe lata