昨天根据一篇文献的描述构建了一个Ni的111表面模型,采用的是LAMMPS内置的命令建模,但是建出来的模型是这样的:
与文献里的256个原子模型的描述少了两个原子,只有254个,in文件的内容如下:
##初始化
units metal
boundary p p p
atom_style atomic
##定义并构建金属Ni模型
lattice fcc 3.52 orient x -1 1 0 orient y -1 -1 2 orient z 1 1 1
region box block 0 4 0 3 0 2
create_box 1 box
region Ni block 0 4 0 3 0 1
create_atoms 1 region Ni
##最后把模型输出来
mass 1 58.69
write_data Ni2.lmp