计算化学公社

标题: 求助不计算分子内LJ的LAMMPS力场设置 [打印本页]

作者
Author: 刘桂林 时间: 2023-8-29 16:02
标题: 求助不计算分子内LJ的LAMMPS力场设置
请问各位前辈，
LAMPPS设置立场时，分子内部包含键、角、二面角以及coul力，不计算lj。但分子间相互包含lj/coul。
这样的力场参数应该如何设置？
我找到了这个命令，但好像该命令也把coul排除了。
neigh_modify exclude molecule/intra rigid
感谢各位大佬了~
Journal ofCatalysis 227 (2004) 227–241

作者
Author: Graphite 时间: 2023-8-29 22:44
这想法有点奇怪，文献想表达的也有些迷惑。
如果一定想这样做，目前唯一的办法是把每个原子都定义成不同类型，假如你有1000个原子，就定义1000个原子类型。
虽然这样感觉很粗犷暴力，读取文件也可能会慢点，但理论上不会造成性能损失。这几千行pair_coeff相当于排列组合、复制粘贴，所以可以写个脚本自动生成。
pair_coeff的写法本身也可以用通配符，能省些行数

作者
Author: 刘桂林 时间: 2023-8-30 08:41

Graphite 发表于 2023-8-29 22:44
这想法有点奇怪，文献想表达的也有些迷惑。
如果一定想这样做，目前唯一的办法是把每个原子都定义成不同类 ...

哦，这方式确实很暴力，但可行，
如果 LAMMPS 原生代码并不提供相关的设置选项，
那可能是文章本身的公式的问题，或者我理解有误了，
我再去看看相关文献验证一下，
谢谢您了~

欢迎光临计算化学公社 (http://bbs.keinsci.com/)