lammps建模时,我们使用create_atoms创建原子,有些情况下,需要删除一些原子,用到的命令为:delete_atoms。

delete_atoms语法:

delete_atoms style args keyword value ...

下面介绍delete_atoms常见的用法

1.删除组内或区域内原子

整组或者整个区域删除原子,直接使用group或者region关键字即可,删除原子后原子ID已经不连续,使用compress yes可对ID序号进行压缩,也就是重新排序,产生连续的原子ID。

#删除edge组内所有原子
delete_atoms group edge
#删除sphere区域内所有原子,并且不压缩ID
delete_atoms region sphere compress no

2. 删除重叠原子

在不同区域内填充原子,有可能会造成两区域边界部分原子重叠,产生重叠的原子之间受力过大容易飞出box,出现“lost atoms”错误。

因此,在建模时可删除重叠原子,用的命令为:

delete_atoms overlap args = cutoff group1-ID group2-ID

group1-ID、group2-ID为两个原子组,cutoff确定原子间距,间距在cutoff内并且分别处于两个组内的原子之一(group1-ID组)会被删除。

如:

delete_atoms overlap 0.3 group1 group2

3. 随机删除原子

按照一定比例随机删除原子,比如随机产生空位,语法为:

delete_atoms porosity args = region-ID fraction seed

fraction为删除原子比例,seed为随机数种子。

delete_atoms porosity cube 0.1 482793

4. 删除分子

分子含有多个原子,如果仅删除一个原子,如H2O仅删除一个H或者一个O原子,在模拟的时候可能会报错。

这时候就需要删除整个H2O分子,加上关键字mol即可。

delete_atoms  overlap 0.5 Water Cu  mol yes

以上为delete_atoms的主要用法,更高级的使用说明请查询官网手册。

更多lammps教程,关注微信公众号:lammps加油站