1.create_box 这个关键字看似简单实则让人难受。
①在写这个关键字的时候切记:一一对应原则,即:
bond/types 对应 extra/bond/per/atom
angle/types 对应 extra/angle/per/atom
dihedral/types 对应 extra/dihedral/per/atom
improper/types 对应 extra/improper/per/atom
当第一列关键字被用到的时候,第二列紧随其后加进来,不然就会报错:Molecule auto special bond generation overflow (../molecule.cpp:1142)
②必写关键字:extra/improper/per/atom
该关键字必须写,否则报错:同上。该关键字系数必须 >= 所建粒子(分子)中原子的个数
2.create_atoms 该关键字没什么太难的,注意一点:
seed = random # seed 即正整数,并不是让选择random/seed (这个错误让我尴尬)。
lammps--create_atoms and create_box
原文:https://www.cnblogs.com/FY-successful/p/13275904.html