计算化学公社

标题: 用CCSD(T)计算单点能报错 [打印本页]

作者
Author: 快乐如我 时间: 2016-4-21 17:39
标题: 用CCSD(T)计算单点能报错
本帖最后由快乐如我于 2016-4-21 17:40 编辑

请教各位前辈，我已经用bmk/6-31+g(2df,p)将结构都优化好了（含12个原子的C, Si, O, H体系），然后进行单点能计算。
（1）当我用ccsd(t)/aug-cc-pvtz时，显示报错。以前也是用同样的方法算原子数少的分子，就木有这种情况，这是咋回事呢？

（2） aug-cc-pvtz计算时间太长，我换成了Def2TZVP的基组，还是显示报错（L913）...求前辈们指导

作者
Author: sobereva 时间: 2016-4-21 18:24
第一个问题是硬盘空余空间不够。CCSD(T)计算是很占硬盘的，要储存大量的电子积分。没特殊必要的时候绝对不要用弥散，否则慢得多资源消耗得也多得多，还容易有SCF不收敛等问题。建议用cc-pVTZ再试，或者如果硬盘还是腾不出足够的空间储存积分，就甭用CCSD(T)，用双杂化之类。

第二个问题是CCSD迭代不收敛，解决方法：
(, 下载次数 Times of downloads: 60)

作者
Author: 快乐如我 时间: 2016-4-22 02:20

sobereva 发表于 2016-4-21 18:24
第一个问题是硬盘空余空间不够。CCSD(T)计算是很占硬盘的，要储存大量的电子积分。没特殊必要的时候绝对不 ...

我试一下，谢谢老师呀

作者
Author: 18810106930 时间: 2018-11-27 19:24

sobereva 发表于 2016-4-21 18:24
第一个问题是硬盘空余空间不够。CCSD(T)计算是很占硬盘的，要储存大量的电子积分。没特殊必要的时候绝对不 ...

老师，您好！我现在也在做ccsd(T),但是一直L502、L508错误。在网上看到帖子在命令行加了scf=(noincfock,conver=11,maxcyc=1025) iop(5/22=20) 关键词后，可以正常结束。但是这个结果用这些关键词可靠吗，是否具有可比性。期待sob老师的解答！

作者
Author: sobereva 时间: 2018-11-27 21:53

18810106930 发表于 2018-11-27 19:24
老师，您好！我现在也在做ccsd(T),但是一直L502、L508错误。在网上看到帖子在命令行加了scf=(noincfock,c ...

别看网上其它地方乱七八糟的帖子，怎么解决看
解决SCF不收敛问题的方法
http://sobereva.com/61

加maxcyc=1025只会多浪费你9倍的计算资源
iop(5/22=20) 之类乱七八糟关键词别瞎加

用conver=11简直可以说是荒唐

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