CCSD(T,T1DIAG)/cc-pVQZ不收敛

尚艳磊 · 发表于 Post on 2017-12-12 20:26:04

打算用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.
请问是哪里的问题呢

sobereva · 发表于 Post on 2017-12-12 20:29:18

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

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

beefly · 发表于 Post on 2017-12-12 23:15:33

这不是ccsd的问题

尚艳磊 · 发表于 Post on 2017-12-13 09:40:54

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

嗯嗯，谢谢sob老师的答疑

尚艳磊 · 发表于 Post on 2017-12-13 09:41:06

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

谢谢~

尚艳磊 · 发表于 Post on 2017-12-23 14:39:13

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
最后一行是这样的，是不是出错了？
恳请答疑。

zjxitcc · 发表于 Post on 2017-12-23 15:15:55

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

尚艳磊 · 发表于 Post on 2017-12-23 16:24:03

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

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

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

zjxitcc · 发表于 Post on 2017-12-23 16:44:47

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

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

sobereva · 发表于 Post on 2017-12-23 16:54:28

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

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

尚艳磊 · 发表于 Post on 2017-12-23 16:55:03

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

16个核，在集群上计算的

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

sobereva · 发表于 Post on 2017-12-23 16:56:07

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

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

尚艳磊 · 发表于 Post on 2017-12-23 16:59:32

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

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

尚艳磊 · 发表于 Post on 2017-12-23 17:06:46

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

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

sobereva · 发表于 Post on 2017-12-23 17:22:45

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

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

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

[Gaussian/gview] CCSD(T,T1DIAG)/cc-pVQZ不收敛

评分 Rate