MP2单点能量异常高

xueduomei · 发表于 Post on 2019-10-29 20:02:11

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

yflchx · 发表于 Post on 2019-10-29 20:30:05

设置%nproc=36的依据是什么？

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

xueduomei · 发表于 Post on 2019-10-29 20:59:40

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

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

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

yflchx · 发表于 Post on 2019-10-29 21:26:39

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

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

xueduomei · 发表于 Post on 2019-10-29 21:54:42

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

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

xueduomei · 发表于 Post on 2019-10-29 21:56:07

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

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

yflchx · 发表于 Post on 2019-10-29 21:59:02

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

遇到过这种情况，

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

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

xueduomei · 发表于 Post on 2019-10-29 22:08:49

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

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

谢谢您了

sobereva · 发表于 Post on 2019-10-30 03:57:30

这个问题跟核数一般不会有什么直接关系
当前这个体系严重不适合用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)才可能是这个数量级。

xueduomei · 发表于 Post on 2019-10-30 09:26:34

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

谢谢社长

，您的帖子我看了

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

yflchx · 发表于 Post on 2019-10-30 16:25:03

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

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

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

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

sobereva · 发表于 Post on 2019-10-31 05:56:16

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

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

xueduomei · 发表于 Post on 2019-10-31 09:09:34

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

好的，谢谢社长

yflchx · 发表于 Post on 2019-10-31 10:37:55

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

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

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

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

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

xueduomei · 发表于 Post on 2019-11-1 21:30:58

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

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

谢谢您

，辛苦

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

[新手求助] MP2单点能量异常高

浏览过的版块