计算化学公社

标题: 求助反应力场中的势能面扫描（刚性扫描） [打印本页]

作者
Author: ZZU_SCU 时间: 2022-4-5 14:55
标题: 求助反应力场中的势能面扫描（刚性扫描）
各位老师好！本人想训练几个反应力场，看是否能用于自己的体系，想做出如下图表。我是这样想的：通过Gaussian的柔性或者刚性扫描得到DFT的势能面，接着提取每一帧结构，保存成lammps识别的data文件，写一个for循环读取每个文件，使用run 0得到体系能量，然后作图。个人感觉这样好麻烦，想问lammps是否可以刚性扫描呢？怎么写这个语句呢？万分感谢 (, 下载次数 Times of downloads: 11)

作者
Author: sobereva 时间: 2022-4-5 21:38
我不用lammps，不过想提醒的是，图中的DFT的势能曲线是错的，8成是一直当成闭壳层算，而没考虑超过不稳定点时要当对称破缺方式计算，见下文里关于断键的扫描说明
详谈使用Gaussian做势能面扫描
http://sobereva.com/474（http://bbs.keinsci.com/thread-12660-1-1.html）

作者
Author: ZZU_SCU 时间: 2022-4-6 09:16

sobereva 发表于 2022-4-5 21:38
我不用lammps，不过想提醒的是，图中的DFT的势能曲线是错的，8成是一直当成闭壳层算，而没考虑超过不稳定点 ...

谢谢卢老师的提醒，我会使用对称破缺的方式计算。

作者
Author: TinnKuka 时间: 2022-6-20 19:30
您好，请问您用lammps做出刚性扫描了吗？

作者
Author: ZZU_SCU 时间: 2022-6-20 19:51

TinnKuka 发表于 2022-6-20 19:30
您好，请问您用lammps做出刚性扫描了吗？

使用lammps分别计算每个构型的能量，然后连接成曲线就可以了。

作者
Author: TinnKuka 时间: 2022-6-20 19:58

ZZU_SCU 发表于 2022-6-20 19:51
使用lammps分别计算每个构型的能量，然后连接成曲线就可以了。

好的谢谢回复~ 我也正在学习训练反应力场，但还是起步状态。请问你的力场现在怎么样啦？

作者
Author: ZZU_SCU 时间: 2022-6-20 20:57

TinnKuka 发表于 2022-6-20 19:58
好的谢谢回复~ 我也正在学习训练反应力场，但还是起步状态。请问你的力场现在怎么样啦？

我不训练力场了，我只是验证一下已有的力场能不能用

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