lammps实例教程.pdfVIP

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 (a) (b) (c) (d) 图1 (a) Cu 中的 ISF 层错; (b) Al 中的 ISF 层错；(c) Cu 中的孪晶； (d) Al 中的孪晶 层错能和孪晶形成能可由下式计算得到： γ (E =−E ) / A 0 其中，E 为引入层错或孪晶后的体系的能量，E0 为完整晶体体系的能量，A 为层错面或孪晶面 的面积。如下是 EAM 模型计算所得的 Cu 和 Al 的层错能和孪晶形成能： Table 1 Fault energies calculated with LAMMPS