计算化学公社

标题: CCSD(T,T1DIAG)/cc-pVQZ不收敛 [打印本页]

作者
Author: 尚艳磊 时间: 2017-12-12 20:26
标题: CCSD(T,T1DIAG)/cc-pVQZ不收敛
打算用CCSD来外推CBS，前几个都算好了，这个Q基组，老是不收敛错误信息是这样：
Restarting incremental Fock formation.
CoulSu:  requested number of processors reduced to: 6 ShMem 1 Linda.
CoulSu:  requested number of processors reduced to: 6 ShMem 1 Linda.
CoulSu:  requested number of processors reduced to: 6 ShMem 1 Linda.
CoulSu:  requested number of processors reduced to: 6 ShMem 1 Linda.
CoulSu:  requested number of processors reduced to: 6 ShMem 1 Linda.
CoulSu:  requested number of processors reduced to: 6 ShMem 1 Linda.
CoulSu:  requested number of processors reduced to: 6 ShMem 1 Linda.
CoulSu:  requested number of processors reduced to: 6 ShMem 1 Linda.
CoulSu:  requested number of processors reduced to: 6 ShMem 1 Linda.
CoulSu:  requested number of processors reduced to: 6 ShMem 1 Linda.
CoulSu:  requested number of processors reduced to: 6 ShMem 1 Linda.
CoulSu:  requested number of processors reduced to: 6 ShMem 1 Linda.
CoulSu:  requested number of processors reduced to: 6 ShMem 1 Linda.
CoulSu:  requested number of processors reduced to: 6 ShMem 1 Linda.
CoulSu:  requested number of processors reduced to: 6 ShMem 1 Linda.
CoulSu:  requested number of processors reduced to: 6 ShMem 1 Linda.
CoulSu:  requested number of processors reduced to: 6 ShMem 1 Linda.
CoulSu:  requested number of processors reduced to: 6 ShMem 1 Linda.
CoulSu:  requested number of processors reduced to: 6 ShMem 1 Linda.
CoulSu:  requested number of processors reduced to: 6 ShMem 1 Linda.
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UHF) =  -358.262239769    A.U. after  129 cycles
         NFock=128  Conv=0.75D-04    -V/T= 2.0015
<Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 0.5000 <S**2>= 0.7741 S= 0.5120
<L.S>= 0.000000000000E+00
Annihilation of the first spin contaminant:
S**2 before annihilation    0.7741, after    0.7504
Convergence failure -- run terminated.
Error termination via Lnk1e in /home3/shijc/g09/l502.exe at Wed Dec 13 02:00:17 2017.
Job cpu time:    0 days 12 hours 43 minutes 29.9 seconds.
File lengths (MBytes):  RWF= 345 Int=    0 D2E=    0 Chk=    7 Scr=    2

总是出现下面这样的语句：
Restarting incremental Fock formation.
CoulSu:  requested number of processors reduced to: 6 ShMem 1 Linda.
请问是哪里的问题呢

作者
Author: sobereva 时间: 2017-12-12 20:29
还是老生常谈的问题
解决SCF不收敛问题的方法
http://sobereva.com/61
建议优先尝试用cc-pVTZ收敛的波函数当初猜。

默认情况每隔一定圈数就会做完整的Fock矩阵生成，并且重启Incremental Fock过程。你没用#P，当前碰到难收敛情况，迭代圈数非常多了，所以这个提示也出来很多遍。

作者
Author: beefly 时间: 2017-12-12 23:15
这不是ccsd的问题

作者
Author: 尚艳磊 时间: 2017-12-13 09:40

sobereva 发表于 2017-12-12 20:29
还是老生常谈的问题
解决SCF不收敛问题的方法
http://sobereva.com/61

嗯嗯，谢谢sob老师的答疑

作者
Author: 尚艳磊 时间: 2017-12-13 09:41

beefly 发表于 2017-12-12 23:15
这不是ccsd的问题

谢谢~

作者
Author: 尚艳磊 时间: 2017-12-23 14:39

sobereva 发表于 2017-12-12 20:29
还是老生常谈的问题
解决SCF不收敛问题的方法
http://sobereva.com/61

sob老师，算了快1个月的任务，log文件很久没更新了，但是rwf文件还在更新。
我看了输出文件：
NORM(A)= 0.11902784D+01
T1 Diagnostic =    0.02133096
Dominant configurations:
***********************
Spin Case       I J A B       Value
ABAB       24 23 26 27    0.105255D+00
Largest amplitude= 1.05D-01
最后一行是这样的，是不是出错了？
恳请答疑。

作者
Author: zjxitcc 时间: 2017-12-23 15:15
没问题，不是报错。但是从你这计算时间来看，体系太大了，选错程序了，应该用ORCA的DLPNO-CCSD(T)来做，牺牲一点点计算精度，换取节省的大量时间。用高斯的CCSD(T)算这么大体系是作死，因为你这才刚算完CCSD，才要开始(T)部分的计算，这后面指不定还要15天，30天还是多久，更可能直接硬盘写满、内存不够，或者停个电一切回到解放前。
做大计算前，要估算计算量，超过一个月的任务，通常是不会有人做的。

作者
Author: 尚艳磊 时间: 2017-12-23 16:24

zjxitcc 发表于 2017-12-23 15:15
没问题，不是报错。但是从你这计算时间来看，体系太大了，选错程序了，应该用ORCA的DLPNO-CCSD(T)来做，牺 ...

10个原子，想做一个Benchmark的。。。。。

谢谢前辈，刚才简单了解了一下ORCA，我试一下ORCA这个软件，多谢前辈！

作者
Author: zjxitcc 时间: 2017-12-23 16:44

尚艳磊发表于 2017-12-23 16:24
10个原子，想做一个Benchmark的。。。。。
谢谢前辈，刚才简单了解了一下ORCA，我试一下ORCA这个软 ...

10个原子不算大。应该最多半个月能做完的····你是不是机器资源（核数、内存）太少了哈哈哈，不能用拖拉机来算啊

作者
Author: sobereva 时间: 2017-12-23 16:54

尚艳磊发表于 2017-12-23 16:24
10个原子，想做一个Benchmark的。。。。。
谢谢前辈，刚才简单了解了一下ORCA，我试一下ORCA这个软 ...

看你用什么基组，具体什么体系。诸如10个原子的有机体系和10个全是重原子构成的体系，耗时差多了去了。用cc-pVTZ和QZ的耗时也是数个数量级的差距。
如果没经验，瞎用CCSD(T)算的话，几个月出不来一个有用的数据都有可能

作者
Author: 尚艳磊 时间: 2017-12-23 16:55

zjxitcc 发表于 2017-12-23 16:44
10个原子不算大。应该最多半个月能做完的····你是不是机器资源（核数、内存）太少了哈哈哈，不能用拖 ...

16个核，在集群上计算的

我看大家都说的挺快的，所以我觉得我的怎么这么慢呢

作者
Author: sobereva 时间: 2017-12-23 16:56

尚艳磊发表于 2017-12-23 16:55
16个核，在集群上计算的
我看大家都说的挺快的，所以我觉得我的怎么这么慢呢

和你的计算条件、具体设定、具体体系估计根本没可比性
CCSD(T)/cc-pVTZ算有机分子10个原子倒是挺快

作者
Author: 尚艳磊 时间: 2017-12-23 16:59

sobereva 发表于 2017-12-23 16:56
和你的计算条件、具体设定、具体体系估计根本没可比性
CCSD(T)/cc-pVTZ算有机分子10个原子倒是挺快

嗯，我的是10个原子都是CHON的，CCSD(T)/cc-pVTZ倒是挺快的（我用window跑的也很快），就是不知道CCSD(T)/cc-pVQZ（集群上）怎么就很慢了

作者
Author: 尚艳磊 时间: 2017-12-23 17:06

sobereva 发表于 2017-12-23 16:56
和你的计算条件、具体设定、具体体系估计根本没可比性
CCSD(T)/cc-pVTZ算有机分子10个原子倒是挺快

sob老师，我算的是双分子的反应，因为在算之前查了一下网上大家说的情况，10个个原子的有机体系，最慢也就1个月左右多。可是已经跑了一个月了还是没算完。我用的1个节点16个核。
想请教下，会不会有其他什么可能导致的原因会算的这么慢?主要是所有的CCSD(T)/cc-pVTZ倒确实挺快额。

作者
Author: sobereva 时间: 2017-12-23 17:22

尚艳磊发表于 2017-12-23 16:59
嗯，我的是10个原子都是CHON的，CCSD(T)/cc-pVTZ倒是挺快的（我用window跑的也很快），就是不知道CCSD(T) ...

QZ和TZ的耗时有天壤之别...耗时简直就像电动车和坦克的差别
死活非要用QZ，用ORCA或Molpro，还要再快点，可以再利用F12

作者
Author: 尚艳磊 时间: 2017-12-23 18:25

sobereva 发表于 2017-12-23 17:22
QZ和TZ的耗时有天壤之别...耗时简直就像电动车和坦克的差别
死活非要用QZ，用ORCA或Molpro，还要再快 ...

这就很尴尬了。。。。
谢谢sob老师答疑，我看看试一下ORCA。
多谢了，提前祝您，圣诞快乐！

作者
Author: youxx189 时间: 2023-9-26 10:46

zjxitcc 发表于 2017-12-23 15:15
没问题，不是报错。但是从你这计算时间来看，体系太大了，选错程序了，应该用ORCA的DLPNO-CCSD(T)来做，牺 ...

您好  在算ccsd(T)时，log文件显示
Spin Case       I J A B       Value
AA          73       76          -0.382844D+00
AA          74       76          0.226944D+00
ABAB       73 73 76 76    -0.128267D+00
Largest amplitude= 3.83D-01
然后任务掉了  也没有报错之类的  请问这是什么原因啊

作者
Author: zjxitcc 时间: 2023-9-26 11:04

youxx189 发表于 2023-9-26 10:46
您好  在算ccsd(T)时，log文件显示
Spin Case       I J A B       Value
AA       ...

极大概率是内存不足，或硬盘被写满了，从这两个角度检查即可。注意，如果机器有250GB内存，你只写了100GB，这显然是输入文件不合理，设置的内存不足；如果机器内存只有200GB，你却写了220GB，这显然超过了机器上限，也不行。

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