lammps实例3.pdf
需积分: 50 191 浏览量
2020-09-27
16:18:33
上传
评论
收藏 328KB PDF 举报
Project #3 on Molecular Dynamics Simulations
FCC 金属中的面缺陷
FCC 晶体中，密排面为 {111}，它既是滑移面也是共格孪晶面。孪晶关于这个面成镜面对称。
{111} 的另一种面缺陷是层错。层错有两种，本征 (intrinsic) 和非本征 (extrinsic) 。抽出一层
原子形成本征层错，插入一层原子形成非本征层错。需要注意的是，低能量层错都可以由该面上
的剪切 (shearing) 操作得到，比如本征层错就是将某层原子上方所有的原子整体移动
1
112 .
6
<>
层错是密排面上的原子错排，层错能 (SFE) 是材料的本征属性，
可度量晶面滑移发生的难易。
接下来我们将运用 LAMMPS 计算 FCC 金属 Cu 和 Al 的层错能和孪晶形成能。层错和
孪晶的构型由其他代码生成，LAMMPS 通过
read_data 命令读取构型，进行计算。
例如，计算 Cu 中层错的输入文件 in.isfCu 如下：
units metal
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 ------------------"
资源评论
