计算化学公社

标题: 用CCSD(T)优化结构老是报错L114 [打印本页]

作者
Author: 快乐如我 时间: 2016-10-14 17:07
标题: 用CCSD(T)优化结构老是报错L114
请教各位前辈，用CCSD(T)优化一个结构，老是提示出现如下错误，不知道如何解决?
输入文件和报错信息如下：
%chk=Si2O4.chk
%nprocshared=4
%mem=8000MB
# opt=maxcycle=1000 UCCSD(T)/6-311g(d)
scf=(qc,maxcycle=600)

Molecule Name

1 2
O
O                1          B1
Si                1          B2 2          A1
Si                1          B3 3          A2 2          D1 0
O                4          B4 1          A3 3          D2 0
O                3          B5 1          A4 4          D3 0

B1          2.35300000
B2          1.65796366
B3          1.65775229
B4          1.54210000
B5          1.54180000
A1          44.79717960
A2          89.58700152
A3          134.78982192
A4          134.79717960
D1          0.00000000
D2          180.00000000
D3          180.00000000
报错信息：
HESSIAN DOES NOT HAVE THE DESIRED LOCAL STRUCTURE
TAKING SIMPLE RFO STEP
SEARCHING FOR LAMDA THAT MINIMIZES ALONG ALL MODES

*****************************************
*** UNABLE TO DETERMINE LAMDA IN FmD114 **
*****************************************

Error termination via Lnk1e in /afs/math/software/gaussian/g09.D01/@sys/l114.exe

作者
Author: sobereva 时间: 2016-10-14 18:21
CCSD(T)本来就不适合作为结构优化目的，比CCSD优化结构没什么改进，本身单点就比CCSD耗时高很多，再加上没有解析梯度耗时颇高。
而且，这种档次的基组和CCSD(T)完全不搭，即便是对于优化来说。
maxcyc设这么大根本无意义
SCF若无特殊必要甭用qc，qc本身还有收敛到非最稳定波函数的可能

作者
Author: 快乐如我 时间: 2016-10-14 23:19

sobereva 发表于 2016-10-14 18:21
CCSD(T)本来就不适合作为结构优化目的，比CCSD优化结构没什么改进，本身单点就比CCSD耗时高很多，再加上没 ...

谢谢老师的回复，我本来不用CCSD(T)来优化的，是因为这个cluster用不同的DFT方法得出来的稳定结构不一样，所以我想对比一下CCSD(T)的结果看究竟哪个结构能量更低。基组我最初选的是Def2-TZVP，就是因为太耗时了，所以先把基组改小了，但是不论用哪个基组，都会报上面相同的错误。。希望得到老师的指点。

作者
Author: sobereva 时间: 2016-10-15 00:56

快乐如我发表于 2016-10-14 23:19
谢谢老师的回复，我本来不用CCSD(T)来优化的，是因为这个cluster用不同的DFT方法得出来的稳定结构不一样 ...

用CCSD/cc-pVTZ，或者降一点到def2-TZVP，对这种体系绝对优化得动。用CCSD(T)结合6-311G*耗时多还不准

作者
Author: 快乐如我 时间: 2016-10-15 05:32

sobereva 发表于 2016-10-15 00:56
用CCSD/cc-pVTZ，或者降一点到def2-TZVP，对这种体系绝对优化得动。用CCSD(T)结合6-311G*耗时多还不准

好的，我试一下，非常非常感谢老师的指导。

作者
Author: 清微 时间: 2018-3-14 09:15

sobereva 发表于 2016-10-14 18:21
CCSD(T)本来就不适合作为结构优化目的，比CCSD优化结构没什么改进，本身单点就比CCSD耗时高很多，再加上没 ...

请问老师，CCSD(T)没有解析梯度，只有数值梯度，那么CCSD(T)可以计算频率么？如果可以，数值梯度下的频率和解析梯度下的频率相差不大么？

作者
Author: sobereva 时间: 2018-3-14 11:34

清微发表于 2018-3-14 09:15
请问老师，CCSD(T)没有解析梯度，只有数值梯度，那么CCSD(T)可以计算频率么？如果可以，数值梯度下的频率 ...

可以算频率，但显然超级慢，而且全数值计算肯定会有一定数值噪音
CFOUR倒是有解析Hessian

作者
Author: 清微 时间: 2018-3-14 12:29
在高斯下有没有做过对比，freq=analytic（默认）和freq=numerical的区别？比如偶极矩，极化张量，振动频率和自由能等。

作者
Author: sobereva 时间: 2018-3-14 21:25

清微发表于 2018-3-14 12:29
在高斯下有没有做过对比，freq=analytic（默认）和freq=numerical的区别？比如偶极矩，极化张量，振动频率 ...

just找个小体系去compare

作者
Author: paitel 时间: 2019-6-11 18:33
本帖最后由 paitel 于 2019-6-11 18:35 编辑

sobereva 发表于 2016-10-15 00:56
用CCSD/cc-pVTZ，或者降一点到def2-TZVP，对这种体系绝对优化得动。用CCSD(T)结合6-311G*耗时多还不准

用 CCSD(T)/cc-pVTZ 优化还是会报同样的错 “UNABLE TO DETERMINE LAMDA IN FmD114”, 是否还有其它解决办法？

作者
Author: sobereva 时间: 2019-6-12 20:01

paitel 发表于 2019-6-11 18:33
用 CCSD(T)/cc-pVTZ 优化还是会报同样的错 “UNABLE TO DETERMINE LAMDA IN FmD114”, 是否还有其它解决 ...

没细节没法说
实在找不出问题就尝试用别的程序，特别是有CCSD(T)解析梯度的程序诸如PSI4、CFOUR、Molpro

作者
Author: a35078337 时间: 2021-12-22 10:38

sobereva 发表于 2016-10-14 18:21
CCSD(T)本来就不适合作为结构优化目的，比CCSD优化结构没什么改进，本身单点就比CCSD耗时高很多，再加上没 ...

请问sobereva老师，审稿人建议用CCSD（T）做优化，请问有相关文献说明CCSD优化和CCSD(T)做优化并无相应差别吗？CCSD（T）优化老是L114，难以解决。

作者
Author: zjxitcc 时间: 2021-12-22 10:50

a35078337 发表于 2021-12-22 10:38
请问sobereva老师，审稿人建议用CCSD（T）做优化，请问有相关文献说明CCSD优化和CCSD(T)做优化并无相应差 ...

大概率是你输入文件写的不对，没写清楚要优化的变量。

作者
Author: sobereva 时间: 2021-12-22 11:56

a35078337 发表于 2021-12-22 10:38
请问sobereva老师，审稿人建议用CCSD（T）做优化，请问有相关文献说明CCSD优化和CCSD(T)做优化并无相应差 ...

这属于基本常识

以后提问直接说完整报错，说不清楚就上传输出文件，别就问个L几几几，这是最忌讳的提问方式。

(, 下载次数 Times of downloads: 46)

作者
Author: a35078337 时间: 2021-12-22 13:58

sobereva 发表于 2021-12-22 11:56
这属于基本常识

以后提问直接说完整报错，说不清楚就上传输出文件，别就问个L几几几，这是最忌讳的提 ...

好的，谢谢sobereva老师

作者
Author: a35078337 时间: 2021-12-22 13:59

zjxitcc 发表于 2021-12-22 10:50
大概率是你输入文件写的不对，没写清楚要优化的变量。

好的，谢谢您~，我再研究研究

作者
Author: danni511 时间: 2022-5-4 13:06
请问sob老师，我分别用PBE-D3，TPSS-D3和CCSD，基组def2-tzvp优化过渡金属团簇并计算了频率，对于同一个输入文件，PBE和TPSS泛函方法得到的频率都有一个200左右的虚频，而CCSD计算后却没有虚频，请问老师是CCSD计算的准确，PBE和TPSS泛函方法计算的不准确吗，谢谢老师

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