计算化学公社

标题: MP2单点能量异常高 [打印本页]

作者
Author: xueduomei 时间: 2019-10-29 20:02
标题: MP2单点能量异常高
请教一下，B3LYP/6-311g(3df)优化能量zero-point Energies=-20108.381208,用MP2做单点能量HF=-20088.5959644\MP2=-65306.8319173，其他同类型体系单点计算都正常，就这一个结果异常，请问一下这是什么原因呢？谢谢各位老师。输出文件压缩在附件中

作者
Author: yflchx 时间: 2019-10-29 20:30
设置%nproc=36的依据是什么？

一个节点36核？减少%nproc试试。不要用Linda跨节点并行试试（如果用了）。

作者
Author: xueduomei 时间: 2019-10-29 20:59

yflchx 发表于 2019-10-29 20:30
设置%nproc=36的依据是什么？

一个节点36核？减少%nproc试试。不要用Linda跨节点并行试试（如果用了）。 ...

您好，我这个是用36核的小服务器算的

作者
Author: yflchx 时间: 2019-10-29 21:26

xueduomei 发表于 2019-10-29 20:59
您好，我这个是用36核的小服务器算的

减少核数试试；或者换个服务器，核数少一点的，比如12核。

作者
Author: xueduomei 时间: 2019-10-29 21:54

yflchx 发表于 2019-10-29 21:26
减少核数试试；或者换个服务器，核数少一点的，比如12核。

您好，我想问核数会影响MP2单点异常吗？

作者
Author: xueduomei 时间: 2019-10-29 21:56

xueduomei 发表于 2019-10-29 21:54
您好，我想问核数会影响MP2单点异常吗？

我这个体系900多个基函数，核数少的机器我们组内存也小 g_write

作者
Author: yflchx 时间: 2019-10-29 21:59

xueduomei 发表于 2019-10-29 21:54
您好，我想问核数会影响MP2单点异常吗？

遇到过这种情况，

但是对你这个具体的问题，我并不能确定问题所在，所以建议你试试。

我把你的输入文件提交了一下，看看能否算对。

作者
Author: xueduomei 时间: 2019-10-29 22:08

yflchx 发表于 2019-10-29 21:59
遇到过这种情况，

但是对你这个具体的问题，我并不能确定问题所在，所以建议你试试。

谢谢您了

作者
Author: sobereva 时间: 2019-10-30 03:57
这个问题跟核数一般不会有什么直接关系
当前这个体系严重不适合用MP2，一看结构就知道静态相关肯定很强（因为成键诡异），MP2相关能严重离谱并不意外

改成DFT，快得多而且结果好得多。对这类体系用MP2，就算结果看似没异常，发表文章的时候也很容易被内行审稿人喷，远不如继续用优化时候用的B3LYP

注意看此文了解MP2如今的地位和仅存的价值
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272（http://bbs.keinsci.com/thread-536-1-1.html）

此外，6-311g(3df)这种基组也严重不适合给后HF计算用，用cc-pVTZ耗时相仿佛，而结果好得多得多，仔细看
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）

zero-point Energies=-20108.381208这个描述明显是错误的，ZPE绝对不可能这么大，U(0)才可能是这个数量级。

作者
Author: xueduomei 时间: 2019-10-30 09:26

sobereva 发表于 2019-10-30 03:57
这个问题跟核数一般不会有什么直接关系
当前这个体系严重不适合用MP2，一看结构就知道静态相关肯定很强（ ...

谢谢社长

，您的帖子我看了

。
还有一点疑问想问老师：
1. 我的优化其实用B3LYP/6-311G**就可以，优化我没必要用6-311G(3df)基组？
2. 我的需要改成DFT-D3(BJ)改善综合性？
3. 我的单点考虑到双杂化泛函计算资源有限，如果方法继续用B3LYP，级别用B3LYP/def2TZVPP会更好点？

作者
Author: yflchx 时间: 2019-10-30 16:25

sobereva 发表于 2019-10-30 03:57
这个问题跟核数一般不会有什么直接关系
当前这个体系严重不适合用MP2，一看结构就知道静态相关肯定很强（ ...

http://bbs.keinsci.com/thread-10113-1-1.html

简单分子也会出这种问题。调整核数可能就解决。

当然，MP2已经被各种泛函超越，前提是泛函选对。

作者
Author: sobereva 时间: 2019-10-31 05:56

xueduomei 发表于 2019-10-30 09:26
谢谢社长，您的帖子我看了。
还有一点疑问想问老师：
1. 我的优化其实用B3LYP/6-311G** ...

1 是
2 已经不是改善“综合性”的层面了，绝大部分问题，特别是包括你当前的体系，合适的泛函都吊打MP2
3 你的体系用ORCA跑双杂化毫无压力。只能用Gaussian算普通泛函的话，算单点用B3LYP/def2TZVPP可以，但最好加上D3

作者
Author: xueduomei 时间: 2019-10-31 09:09

sobereva 发表于 2019-10-31 05:56
1 是
2 已经不是改善“综合性”的层面了，绝大部分问题，特别是包括你当前的体系，合适的泛函都吊打MP2
...

好的，谢谢社长

作者
Author: yflchx 时间: 2019-10-31 10:37

xueduomei 发表于 2019-10-29 22:08
谢谢您了

State=1-A\HF=-20088.5959645\MP2=-20091.5071129\

输出文件附上（涉及账号、节点的内容我删除了）。内存没那么大，所以MP2后加了SemiDirect。

另外，使用的版本也有差距，所以这个技术问题到底出在哪里并不十分明朗。

鉴于你要放弃MP2，这个问题也就不重要了。倒是可以看看用双杂化是不是还出这问题。

作者
Author: xueduomei 时间: 2019-11-1 21:30

yflchx 发表于 2019-10-31 10:37
State=1-A\HF=-20088.5959645\MP2=-20091.5071129\

输出文件附上（涉及账号、节点的内容我删除了）。 ...

谢谢您

，辛苦

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