计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 分子模拟 (Molecular Modeling) › 求助不计算分子内LJ的LAMMPS力场设置

求助不计算分子内LJ的LAMMPS力场设置

查看数: 1044 | 评论数: 2 | 收藏 Add to favorites 0

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

刘桂林

发布时间: 2023-8-29 16:02

正文摘要:

请问各位前辈， LAMPPS设置立场时，分子内部包含键、角、二面角以及coul力，不计算lj。但分子间相互包含lj/coul。这样的力场参数应该如何设置？我找到了这个命令，但好像该命令也把coul排除了。 neigh_modify ...

回复 Reply

刘桂林 发表于 Post on 2023-8-30 08:41:21

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

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

Graphite 发表于 Post on 2023-8-29 22:44:26

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register