计算化学公社

 找回密码 Forget password
 注册 Register
Views: 14265|回复 Reply: 16
打印 Print 上一主题 Last thread 下一主题 Next thread

[量化理论] 请教算CCSD(T)单点能因T1诊断值太大而能量不准?

[复制链接 Copy URL]

39

帖子

0

威望

651

eV
积分
690

Level 4 (黑子)

请教各位前辈,我用BMK/Def2-TZVP优化一个主族反应物cluster,然后用CCSD(T)计算单点能,因为T1诊断值太大(0.029左右),导致加上单点能后的势能面能量完全偏离了DFT的能量(优化其他产物及过渡态都没有问题,T1值都在0.015左右),加大基组计算CCSD(T)能稍微降低T1值(从0.029降到0.028),但仍然解决不了能量不准的问题。。。
希望有经验的前辈帮帮忙

4133

帖子

4

威望

8907

eV
积分
13120

Level 6 (一方通行)

MOKIT开发者

17#
发表于 Post on yesterday 22:04 | 只看该作者 Only view this author
本帖最后由 zjxitcc 于 2025-8-21 22:16 编辑
嘤嘤嘤 发表于 2025-8-21 21:46
老师,我用的dlpno-ccsd(t)方法做的单点计算,结果显示单点的能量和DFT下的相对能量差了近27kcal/mol。当 ...

(1)先检验UHF波函数稳定性,如果不稳定,则不能做DLPNO-CCSD(T)计算,需要先将波函数优化至稳定。
如何在ORCA中检验波函数稳定性?很简单,输入文件格式为
! UHF 基组 RI加速 格点等等 去掉DLPNO-CC的关键词
%scf
STABPerform true
STABRestartUHFifUnstable true
STABMaxIter 500
STABDTol 1e-5
STABRTol 1e-5
end
这个任务算出来的UHF能量必须与你之前DLPNO-CCSD(T)任务中的UHF能量相等,证明收敛到了同一个解,然后这个任务会给出是否稳定的结论,如果不稳定它会自动优化至稳定。

(2)如果波函数稳定,但DLPNO-CCSD(T)相对能量仍然异常,可以看我前两天对一个自由基体系的详细解答,体系存在两个UHF稳定解,若使用较高能量的UHF解进行后续计算,能量会异常。
“求助ORCA DLPNO-CCSD(t)-F12计算能量异常问题”其中8L的回复
http://bbs.keinsci.com/thread-55176-1-1.html
自动做多参考态计算的程序MOKIT

119

帖子

0

威望

1989

eV
积分
2108

Level 5 (御坂)

16#
发表于 Post on yesterday 21:46 | 只看该作者 Only view this author
sobereva 发表于 2017-6-9 10:22
stable关键词检测

老师,我用的dlpno-ccsd(t)方法做的单点计算,结果显示单点的能量和DFT下的相对能量差了近27kcal/mol。当时普通计算单点是波函数不收敛,后来加了! DLPNO-CCSD(T) TightSCF tightpno kdiis RIJCOSX noautostart miniprint SlowConv noTRAH
! NOPOP NOPRINTMOS
%SCF
CNVShift true
LShift 0.5
ShiftErr 0.1
end这一堆关键词就收敛不报错了,但是能量还是很高。请问这种情况怎么解决啊?这是一个自由基和分子组合的复合物体系。我刚才用ccsd(T1)算了一下T1诊断值是0.252,也还好

6万

帖子

99

威望

5万

eV
积分
120300

管理员

公社社长

15#
发表于 Post on 2018-1-14 23:57:16 | 只看该作者 Only view this author
outing 发表于 2018-1-11 23:01
请教大神问题,假设说,我判断出一个体系应该考虑一下多参考特征,那么这个多参考特征在计算上应该怎么实现 ...

用CASPT2、MRCI之类方法
静态相关强的体系,通常需要用多参考方法才能较准确描述
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

9

帖子

0

威望

1230

eV
积分
1239

Level 4 (黑子)

14#
发表于 Post on 2018-1-11 23:01:28 | 只看该作者 Only view this author
请教大神问题,假设说,我判断出一个体系应该考虑一下多参考特征,那么这个多参考特征在计算上应该怎么实现?    另,这个多参考特征该怎么理解呢?

39

帖子

0

威望

651

eV
积分
690

Level 4 (黑子)

13#
 楼主 Author| 发表于 Post on 2017-6-9 17:33:49 | 只看该作者 Only view this author
liyuanhe211 发表于 2017-6-9 08:02
CCSD(T)与预期不同就怀疑T1的思路感觉就像DFT不好就怀疑自旋污染一样。看对比一下比如顺便产生的HF、MP2等 ...

谢谢老师的回复。
我又检查了顺便产生的HF、MP2等水平差值,同CCSD(T)一样偏差很大,用stable测试波函数稳定性也没有发现问题。。
所以现在我也不知道怎么办了。。不行我就只能用相同DFT方法做单点能校正了。。

6万

帖子

99

威望

5万

eV
积分
120300

管理员

公社社长

12#
发表于 Post on 2017-6-9 11:19:08 | 只看该作者 Only view this author
xujian 发表于 2017-6-9 11:15
请教一下楼主,您是否遇到过CCSD不收敛的问题,不知道您如何解决。我只会用CCSD(T),NOCHECK

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

588

帖子

0

威望

3310

eV
积分
3898

Level 5 (御坂)

11#
发表于 Post on 2017-6-9 11:15:43 | 只看该作者 Only view this author
快乐如我 发表于 2017-6-9 04:14
谢谢您的回复。
我用了很多DFT方法,能量都差不太多,而且跟我实验数据对应很好,用CCSD(T)单点能校正后 ...

请教一下楼主,您是否遇到过CCSD不收敛的问题,不知道您如何解决。我只会用CCSD(T),NOCHECK

6万

帖子

99

威望

5万

eV
积分
120300

管理员

公社社长

10#
发表于 Post on 2017-6-9 10:24:19 | 只看该作者 Only view this author
快乐如我 发表于 2017-6-9 04:14
谢谢您的回复。
我用了很多DFT方法,能量都差不太多,而且跟我实验数据对应很好,用CCSD(T)单点能校正后 ...


多参考方法并不是降低CCSD(T)的T1诊断值,而是对于T1诊断值很高的时候,有必要用多参考方法,以得到靠谱结果。
基于DFT轨道做CCSD(T)倒是让T1诊断值大幅降低,但是结果不中用,见这里的测试
各种后HF方法精度简单横测
http://sobereva.com/378
取消冻核近似,倒是可以让T1诊断值明显降低一些,见这里的对比http://bbs.keinsci.com/forum.php?mod=viewthread&tid=5460
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

5万

eV
积分
120300

管理员

公社社长

9#
发表于 Post on 2017-6-9 10:22:49 | 只看该作者 Only view this author
xujian 发表于 2017-6-9 08:34
李老师,波函数稳定性问题该从哪里看啊

stable关键词检测
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

5万

eV
积分
120300

管理员

公社社长

8#
发表于 Post on 2017-6-9 10:22:25 | 只看该作者 Only view this author
xujian 发表于 2017-6-9 08:44
在这篇文献(JPCA, doi:10.1021/acs.jpca.5b12275)中,有几个过渡态的T1Diag大于0.03,算出来的能量跟高精 ...


不要对T1诊断值过于敏感,此指标过于片面
不代表超过0.02乃至0.03,CCSD(T)结果就不佳,更没法直接说DFT的结果可靠性就差。
这里有一批分子的T1诊断值http://bbs.keinsci.com/forum.php?mod=viewthread&tid=5460
很多超过0.02乃至0.03的,比如四羰基镍,用DFT、CCSD(T)都可以很好描述
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

588

帖子

0

威望

3310

eV
积分
3898

Level 5 (御坂)

7#
发表于 Post on 2017-6-9 08:44:00 | 只看该作者 Only view this author
快乐如我 发表于 2017-6-9 04:14
谢谢您的回复。
我用了很多DFT方法,能量都差不太多,而且跟我实验数据对应很好,用CCSD(T)单点能校正后 ...

在这篇文献(JPCA, doi:10.1021/acs.jpca.5b12275)中,有几个过渡态的T1Diag大于0.03,算出来的能量跟高精度的热力学组合方法W1BD也很接近,比多参考态CASPT2吻合的还更好些,所以T1Diag影响真的有那么大吗?

588

帖子

0

威望

3310

eV
积分
3898

Level 5 (御坂)

6#
发表于 Post on 2017-6-9 08:34:40 | 只看该作者 Only view this author
liyuanhe211 发表于 2017-6-9 08:02
CCSD(T)与预期不同就怀疑T1的思路感觉就像DFT不好就怀疑自旋污染一样。看对比一下比如顺便产生的HF、MP2等 ...

李老师,波函数稳定性问题该从哪里看啊

3097

帖子

29

威望

1万

eV
积分
17098

Level 6 (一方通行)

5#
发表于 Post on 2017-6-9 08:02:31 | 只看该作者 Only view this author
CCSD(T)与预期不同就怀疑T1的思路感觉就像DFT不好就怀疑自旋污染一样。看对比一下比如顺便产生的HF、MP2等水平差值看是否合理。再看看有没有波函数稳定性问题。

44

帖子

0

威望

663

eV
积分
707

Level 4 (黑子)

4#
发表于 Post on 2017-6-9 04:43:20 | 只看该作者 Only view this author
200kJ的误差就不能叫误差了,很可能是错误,估计不是T1diag的问题。用的是同样的基组?

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-22 12:13 , Processed in 0.217763 second(s), 24 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list