求助反应力场中的势能面扫描（刚性扫描）

ZZU_SCU · 发表于 Post on 2022-4-5 14:55:29

各位老师好！本人想训练几个反应力场，看是否能用于自己的体系，想做出如下图表。我是这样想的：通过Gaussian的柔性或者刚性扫描得到DFT的势能面，接着提取每一帧结构，保存成lammps识别的data文件，写一个for循环读取每个文件，使用run 0得到体系能量，然后作图。个人感觉这样好麻烦，想问lammps是否可以刚性扫描呢？怎么写这个语句呢？万分感谢

sobereva · 发表于 Post on 2022-4-5 21:38:19

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

ZZU_SCU · 发表于 Post on 2022-4-6 09:16:19

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

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

TinnKuka · 发表于 Post on 2022-6-20 19:30:04

您好，请问您用lammps做出刚性扫描了吗？

ZZU_SCU · 发表于 Post on 2022-6-20 19:51:46

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

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

TinnKuka · 发表于 Post on 2022-6-20 19:58:47

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

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

ZZU_SCU · 发表于 Post on 2022-6-20 20:57:32

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

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

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

[Lammps] 求助反应力场中的势能面扫描（刚性扫描）

浏览过的版块