请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 3113|回复: 10

[新手求助] 用CCSD(T)优化结构老是报错L114

[复制链接]

39

帖子

0

威望

565

eV
积分
604

Level 4 (黑子)

发表于 2016-10-14 17:07:28 | 显示全部楼层 |阅读模式
请教各位前辈,用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

2万

帖子

25

威望

2万

eV
积分
45687

管理员

公社社长

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

评分

参与人数 1eV +1 收起 理由
jackie + 1 赞!

查看全部评分

北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

39

帖子

0

威望

565

eV
积分
604

Level 4 (黑子)

 楼主| 发表于 2016-10-14 23:19:38 | 显示全部楼层
sobereva 发表于 2016-10-14 18:21
CCSD(T)本来就不适合作为结构优化目的,比CCSD优化结构没什么改进,本身单点就比CCSD耗时高很多,再加上没 ...

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

2万

帖子

25

威望

2万

eV
积分
45687

管理员

公社社长

发表于 2016-10-15 00:56:08 | 显示全部楼层
快乐如我 发表于 2016-10-14 23:19
谢谢老师的回复,我本来不用CCSD(T)来优化的,是因为这个cluster用不同的DFT方法得出来的稳定结构不一样 ...


用CCSD/cc-pVTZ,或者降一点到def2-TZVP,对这种体系绝对优化得动。用CCSD(T)结合6-311G*耗时多还不准
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

39

帖子

0

威望

565

eV
积分
604

Level 4 (黑子)

 楼主| 发表于 2016-10-15 05:32:53 | 显示全部楼层
sobereva 发表于 2016-10-15 00:56
用CCSD/cc-pVTZ,或者降一点到def2-TZVP,对这种体系绝对优化得动。用CCSD(T)结合6-311G*耗时多还不准

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

46

帖子

0

威望

632

eV
积分
678

Level 4 (黑子)

发表于 2018-3-14 09:15:09 | 显示全部楼层
sobereva 发表于 2016-10-14 18:21
CCSD(T)本来就不适合作为结构优化目的,比CCSD优化结构没什么改进,本身单点就比CCSD耗时高很多,再加上没 ...

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

2万

帖子

25

威望

2万

eV
积分
45687

管理员

公社社长

发表于 2018-3-14 11:34:44 | 显示全部楼层
清微 发表于 2018-3-14 09:15
请问老师,CCSD(T)没有解析梯度,只有数值梯度,那么CCSD(T)可以计算频率么?如果可以,数值梯度下的频率 ...

可以算频率,但显然超级慢,而且全数值计算肯定会有一定数值噪音
CFOUR倒是有解析Hessian
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

46

帖子

0

威望

632

eV
积分
678

Level 4 (黑子)

发表于 2018-3-14 12:29:49 | 显示全部楼层
在高斯下有没有做过对比,freq=analytic(默认)和freq=numerical的区别?比如偶极矩,极化张量,振动频率和自由能等。

2万

帖子

25

威望

2万

eV
积分
45687

管理员

公社社长

发表于 2018-3-14 21:25:12 | 显示全部楼层
清微 发表于 2018-3-14 12:29
在高斯下有没有做过对比,freq=analytic(默认)和freq=numerical的区别?比如偶极矩,极化张量,振动频率 ...

just找个小体系去compare
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

32

帖子

0

威望

172

eV
积分
204

Level 3 能力者

发表于 2019-6-11 18:33:10 | 显示全部楼层
本帖最后由 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”, 是否还有其它解决办法?

2万

帖子

25

威望

2万

eV
积分
45687

管理员

公社社长

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

没细节没法说
实在找不出问题就尝试用别的程序,特别是有CCSD(T)解析梯度的程序诸如PSI4、CFOUR、Molpro
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2019-9-16 09:00 , Processed in 0.163821 second(s), 24 queries .

快速回复 返回顶部 返回列表