Project #3 on Molecular Dynamics Simulations
FCC 金属中的面缺陷
FCC 晶体中,密排面为 {111},它既是滑移面也是共格孪晶面。孪晶关于这个面成镜面对称。
{111} 的另一种面缺陷是层错。层错有两种,本征 (intrinsic) 和非本征 (extrinsic) 。抽出一层
原子形成本征层错,插入一层原子形成非本征层错。需要注意的是,低能量层错都可以由该面上
的剪切 (shearing) 操作得到,比如本征层错就是将某层原子上方所有的原子整体移动
层错是密排面上的原子错排,层错能 (SFE) 是材料的本征属性,
可度量晶面滑移发生的难易。
接下来我们将运用 LAMMPS 计算 FCC 金属 Cu 和 Al 的层错能和孪晶形成能。层错和
孪晶的构型由其他代码生成,LAMMPS 通过
read_data 命令读取构型,进行计算。
例如,计算 Cu 中层错的输入文件 in.isfCu 如下:
boundary p p p
atom_style atomic
read_data isf-Cu
region middle block INF INF INF INF INF 87.659 units box
group middle region middle
pair_style eam/alloy
pair_coeff * * jin_copper_lammps.setfl Cu
min_style sd
minimize 1.0e-8 1.0e-8 10000 10000
compute 3 all pe/atom
compute 4 all ke/atom
compute 5 all coord/atom 3.0
dump 1 middle custom 1 dump.atom id xs ys zs c_3 c_4 c_5
dump_modify 1 format "%d %16.9g %16.9g %16.9g %16.9g %16.9g %g"
run 0
shell mkdir isfCu_dump
shell mv dump.atom isfCu_dump/dump.atom
variable E equal pe
print "--------------------- ISF in Cu, E = $E ------------------"
评论0