计算化学公社

标题: 过渡态计算遇到力很大且不收敛问题 [打印本页]

作者
Author: Garyall 时间: 2022-5-2 10:43
标题: 过渡态计算遇到力很大且不收敛问题
本帖最后由 Garyall 于 2022-5-3 14:45 编辑

各位师兄师姐，本人在计算过渡态时遇到了问题（Pd金属表面C=N双键的加氢）。我把我的操作步骤贴一下，如果有人经历过相似的问题或者看出我问题的，麻烦告诉我，不胜感激！！
1.IS 和FS结构优化，粗算 K点采用111 （图1-3分别表示结构优化时的INCAR， IS/CONTCAR，FS/CONTCAR）
2.新建文件夹，创立IS和FS文件夹，并将优化好的CONTCAR放于相应的文件夹中，并检查相似度为3.394（图4）
3.插入5个点，并将IS的CONTCAR放入00文件夹，把FS的CONTCAR放入06文件夹中（图5）
4.检查距离 nebavoid.pl 每个原子间距是大于1的
5.设置过渡态INCAR（图6） K点取111
6.从第2个离子步开始，2个电子步就跑完能量收敛了，但是力一直不收敛且保持不变（图7）
请各位帮忙看一下哪里出问题了小弟感激不尽！！！

作者
Author: Garyall 时间: 2022-5-2 15:33
有没有遇到过同样的问题

作者
Author: Garyall 时间: 2022-5-2 20:25
是结构优化时精度不够么

作者
Author: Garyall 时间: 2022-5-3 12:05
手动@sob老师 @卡开老师

作者
Author: chands 时间: 2022-5-3 13:54

Garyall 发表于 2022-5-3 12:05
手动@sob老师 @卡开老师

给你一个建议，你把问题说得具体点，什么体系，什么步骤，最好附上输入输出文件，还有截屏图片清晰一点。我作为纯外行看得云里雾里，估计内行不会好到哪里去。你说得具体点，他们说不定顺手就回答了。还有@后面要加用户名全称，否则@不到他们。。。。

作者
Author: Garyall 时间: 2022-5-3 14:41

chands 发表于 2022-5-3 13:54
给你一个建议，你把问题说得具体点，什么体系，什么步骤，最好附上输入输出文件，还有截屏图片清晰一点。 ...

谢谢提醒

作者
Author: Garyall 时间: 2022-5-3 14:49
@sobereva
@卡开发发
两位老师，请求你们的帮助

作者
Author: 卡开发发 时间: 2022-5-3 15:47

Garyall 发表于 2022-5-3 14:49
@sobereva
@卡开发发
两位老师，请求你们的帮助

给不了太多帮助，因为过渡态计算本身就比较考验化学直觉。有这几个问题：
1、吸附的结构可能要小心调整，一开始未解离的分子以及解离后不带H的片段吸附应该是个物理吸附，吸附稳定的构象应该比较多（指的是相对于表面xy位置关系），所以这个过渡态按我理解是比较难找的。
2、表面可能也不够大，相邻晶格之间的分子还有比较强的作用。
3、EDIFF的精度偏低，1E-4对受力计算应该还是有噪声的。
4、IOPT写成IPOT，这里确实可以选个效率稍微高一些的方案。

作者
Author: Garyall 时间: 2022-5-3 18:47

卡开发发发表于 2022-5-3 15:47
给不了太多帮助，因为过渡态计算本身就比较考验化学直觉。有这几个问题：
1、吸附的结构可能要小心调整 ...

谢谢卡开发发老师提的建议。我觉得这几点对我帮助很大，就是少了专业人士的指点。谢谢

作者
Author: 13665491757 时间: 2023-7-6 14:40
请问你这个问题解决了吗？我最近也遇到了这个问题。

作者
Author: 醒醒 时间: 2024-6-13 15:35
请问博主的问题解决了吗，我最近也遇到了相似的问题，过渡态的能量收敛了，但是力一直不收敛

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