create_atoms type style args keyword values ... 举例如下: # 在整个模拟区域中添加类型为1的原子(即添加Al原子) create_atoms 1 region whole # 模拟环境中的原子信息 mass 1 26.981 # Al mass 2 55.845 # Fe 9、mass mass命令为某一种或几种类型的原子设置质量。 10、pair_style 与 pair_coeff 设置...
Last command: pair_coeff 3 1 test.txt CN 用了pair_coeff指定每个原子对之间的力,结果报错 Illegal pair_coeff command,明明没有违反语法,而且在自己电脑上跑的也没问题,一放超算上就报错。 这里在于,pair_coeff i j ffield label这个命令对原子i与原子j之间指定受力,i必须≤j,要按顺序来,所以应该改为 p...
region box block 0 4 0 4 0 4 # x,y,z 各方向上的晶胞重复单元数,也即区域大小 create_box 1 box # 将上述区域指定为模拟的盒子create_atoms 1 box # 将原子按晶格填满盒子pair_style eam # 选取 Cu 的 EAM 势作为模型pair_coeff * * Cu_u3.eam # EAM 势文件名称 run 0 # 运行 0 步,仅为...
region box block 0 4 0 4 0 4 # x,y,z 各方向上的晶胞重复单元数,也即区域大小 create_box 1 box # 将上述区域指定为模拟的盒子create_atoms 1 box # 将原子按晶格填满盒子pair_style eam # 选取 Cu 的 EAM 势作为模型pair_coeff * * Cu_u3.eam # EAM 势文件名称 run 0 # 运行 0 步,仅为...
这里,我们虽然已经将部分原子的类别做了转换,但是并没有对新设定的原子种类进行定义。因此,需要使用mass命令赋予新生成的原子以质量。 赋予原子质量 还有一点需要注意的是,之前的算例中,都只涉及单一一种原子,具有相同的势函数。但是,在本例中,设计Fe,Ni,Cr,Co,Cu五种原子的五个势函数。因此,需要在pair_coeff中...
mass 1 64 mass 2 27 pair_style hybrid meam/c meam/c lj/cut 10 pair_coeff * * meam/c 1 library.Cu.meam Cu Cu.meam Cu NULL pair_coeff * * meam/c 2 library-Al.meam Al Al.meam NULL Al pair_coeff 1 2 lj/cut 0.4 2.47 #Cu-Al ...
之后，设置与之相对于的参数，分别使用pair_coeff、bond_coeff、angle_coeff、dihedral_coeff、improper_...
improper_coeff, improper_style, kspace_modify, kspace_ style, pair_coeff, pair_modify, pair_style, pair_write,special_bonds Settings: communicate, dipole, group, mass, min_ modify, min_style, neigh_modify, neighbor, reset_timestep,
mass 5 63 #Cu #设置势函数 pair_style eam/alloy #选用eam合金势 pair_coeff * * FeNiCrCoCu-with-ZBL.eam.alloy Fe Ni Cr Co Cu #选择文献提供的Fe Ni Cr Co Cu势函数[1] min_style cg #选用共轭梯度法对模拟盒子进行结构优化 minimize 1.0e-6 1.0e-7 10000 10000 ...
mass 5 63 #Cu #设置势函数 pair_style eam/alloy pair_coeff * * FeNiCrCoCu-with-ZBL.eam.alloy Fe Ni Cr Co Cu #cg和sd交替进行能量最小化 min_style cg minimize 1.0e-6 1.0e-7 10000 10000 run 0 min_style sd minimize 1.0e-6 1.0e-7 10000 10000 run 0 min_style...