LAMMPS计算 SiC 的结合能.pdf

LAMMPS计算 SiC 的结合能 金剑锋 材料学院 东北大学 SiC 结构 （金刚石结构） 碳化硅晶体和金刚石一样 ，是原子晶体 ，空间网状结构 ，如图； 一个 C 与四个 Si 形成正四面体结构,一个晶胞内有 4 个 C 原 子和 4 个 Si 原子。 LAMMPS 内置函数创建 SiC 结构 #原子模型 lattice custom 4.36 a1 1 0 0 a2 0 1 0 a3 0 0 1 basis 0.0 0.0 0.0 basis 0.5 0.5 0.0 basis 0.5 0.0 0.5 basis 0.0 0.5 0.5 basis 0.25 0.25 0.25 basis 0.75 0.25 0.75 basis 0.25 0.75 0.75 basis 0.75 0.75 0.25 region whole block 0 4 0 4 0 4 create_box 2 whole # 两类原子 # basis 1-4为第1类原子 - C ，其他第2类 - Si create_atoms 1 box basis 1 1 basis 2 1 basis 3 1 basis 4 1 basis 5 2 basis 6 2 basis 7 2 basis 8 2 LAMMPS 脚本 - SiC 结合能 units metal boundary p p p atom_style atomic variable a0 equal 4.36 variable N equal 4 lattice custom ${a0} a1 1 0 0 a2 0 1 0 a3 0 0 1 basis 0.0 0.0 0.0 basis 0.5 0.5 0.0 basis 0.5 0.0 0.5 basis 0.0 0.5 0.5 basis 0.25 0.25 0.25 basis 0.75 0.25 0.75 basis 0.75 0.75 0.25 basis 0.25 0.75 0.75 region myreg block 0 $N 0 $N 0 $N units lattice LAMMPS 脚本 - SiC 结合能 create_box 2 myreg # 创建两种类型原子 create_atoms 1 box basis 1 1 basis 2 1 basis 3 1 basis 4 1 basis 5 2 basis 6 2 basis 7 2 basis 8 2 mass 1 12.0107 # C 的原子量 mass 2 28.0855 # Si的原子量 pair_style tersoff pair_coeff * * SiC.tersoff C Si thermo 10 thermo_style custom step temp pe etotal press vol dump 1 all custom 10 SiC.dump id type x y z Tersoff