Cu间隙原子形成能

#模型基本设置

units   metal

dimension   3

boundary    p  p    p      

atom_style  atomic

#建模

lattice fcc 3.615

region  box block -4 4 -4 4 -4 4

create_box  2 box

lattice fcc 3.615

create_atoms    1 box

#力场参数设置

pair_style  eam/alloy

pair_coeff * * FeCuNi.eam.alloy Cu Cu

#计算势能

compute eng all pe/atom 

compute eatoms all reduce sum c_eng

#设置输出

thermo 10

thermo_style custom step pe c_eatoms 

#能量最小化

min_style cg

minimize 1e-15 1e-15 5000 5000

#计算原子数量

variable N equal count(all)

variable No equal $N

#计算初始能量E0

variable E equal "c_eatoms"

variable Eo equal $E

#产生间隙原子

create_atoms 2 single 0 -1.8075 0 units box

write_data  cu.data

#设置输出

thermo 10

thermo_style custom step pe lx ly lz press pxx pyy pzz c_eatoms 

#再一次能量最小化

min_style cg

minimize 1e-15 1e-15 5000 5000

write_data  cu2.data

#计算插入间隙原子后的能量Ef

variable Ef equal "c_eatoms"

#计算间隙原子形成能

variable Ei equal (${Ef}-((${No}+1)/${No})*${Eo})

#输出间隙原子形成能

print "Interstitial formation energy = ${Ei}"