计算化学公社

 找回密码 Forget password
 注册 Register
Views: 18218|回复 Reply: 16
打印 Print 上一主题 Last thread 下一主题 Next thread

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

[复制链接 Copy URL]

39

帖子

0

威望

651

eV
积分
690

Level 4 (黑子)

请教各位前辈,用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

6万

帖子

99

威望

5万

eV
积分
120205

管理员

公社社长

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

评分 Rate

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

查看全部评分 View all ratings

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

39

帖子

0

威望

651

eV
积分
690

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2016-10-14 23:19:38 | 只看该作者 Only view this author
sobereva 发表于 2016-10-14 18:21
CCSD(T)本来就不适合作为结构优化目的,比CCSD优化结构没什么改进,本身单点就比CCSD耗时高很多,再加上没 ...

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

6万

帖子

99

威望

5万

eV
积分
120205

管理员

公社社长

4#
发表于 Post on 2016-10-15 00:56:08 | 只看该作者 Only view this author
快乐如我 发表于 2016-10-14 23:19
谢谢老师的回复,我本来不用CCSD(T)来优化的,是因为这个cluster用不同的DFT方法得出来的稳定结构不一样 ...


用CCSD/cc-pVTZ,或者降一点到def2-TZVP,对这种体系绝对优化得动。用CCSD(T)结合6-311G*耗时多还不准
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

39

帖子

0

威望

651

eV
积分
690

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2016-10-15 05:32:53 | 只看该作者 Only view this author
sobereva 发表于 2016-10-15 00:56
用CCSD/cc-pVTZ,或者降一点到def2-TZVP,对这种体系绝对优化得动。用CCSD(T)结合6-311G*耗时多还不准

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

46

帖子

0

威望

835

eV
积分
881

Level 4 (黑子)

6#
发表于 Post on 2018-3-14 09:15:09 | 只看该作者 Only view this author
sobereva 发表于 2016-10-14 18:21
CCSD(T)本来就不适合作为结构优化目的,比CCSD优化结构没什么改进,本身单点就比CCSD耗时高很多,再加上没 ...

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

6万

帖子

99

威望

5万

eV
积分
120205

管理员

公社社长

7#
发表于 Post on 2018-3-14 11:34:44 | 只看该作者 Only view this author
清微 发表于 2018-3-14 09:15
请问老师,CCSD(T)没有解析梯度,只有数值梯度,那么CCSD(T)可以计算频率么?如果可以,数值梯度下的频率 ...

可以算频率,但显然超级慢,而且全数值计算肯定会有一定数值噪音
CFOUR倒是有解析Hessian
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

46

帖子

0

威望

835

eV
积分
881

Level 4 (黑子)

8#
发表于 Post on 2018-3-14 12:29:49 | 只看该作者 Only view this author
在高斯下有没有做过对比,freq=analytic(默认)和freq=numerical的区别?比如偶极矩,极化张量,振动频率和自由能等。

6万

帖子

99

威望

5万

eV
积分
120205

管理员

公社社长

9#
发表于 Post on 2018-3-14 21:25:12 | 只看该作者 Only view this author
清微 发表于 2018-3-14 12:29
在高斯下有没有做过对比,freq=analytic(默认)和freq=numerical的区别?比如偶极矩,极化张量,振动频率 ...

just找个小体系去compare
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

35

帖子

0

威望

192

eV
积分
228

Level 3 能力者

10#
发表于 Post on 2019-6-11 18:33:10 | 只看该作者 Only view this author
本帖最后由 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”, 是否还有其它解决办法?

6万

帖子

99

威望

5万

eV
积分
120205

管理员

公社社长

11#
发表于 Post on 2019-6-12 20:01:49 | 只看该作者 Only view this author
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程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

8

帖子

0

威望

63

eV
积分
71

Level 2 能力者

12#
发表于 Post on 2021-12-22 10:38:39 | 只看该作者 Only view this author
sobereva 发表于 2016-10-14 18:21
CCSD(T)本来就不适合作为结构优化目的,比CCSD优化结构没什么改进,本身单点就比CCSD耗时高很多,再加上没 ...

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

4114

帖子

4

威望

8875

eV
积分
13069

Level 6 (一方通行)

MOKIT开发者

13#
发表于 Post on 2021-12-22 10:50:10 | 只看该作者 Only view this author
a35078337 发表于 2021-12-22 10:38
请问sobereva老师,审稿人建议用CCSD(T)做优化,请问有相关文献说明CCSD优化和CCSD(T)做优化并无相应差 ...

大概率是你输入文件写的不对,没写清楚要优化的变量。
自动做多参考态计算的程序MOKIT

6万

帖子

99

威望

5万

eV
积分
120205

管理员

公社社长

14#
发表于 Post on 2021-12-22 11:56:02 | 只看该作者 Only view this author
a35078337 发表于 2021-12-22 10:38
请问sobereva老师,审稿人建议用CCSD(T)做优化,请问有相关文献说明CCSD优化和CCSD(T)做优化并无相应差 ...

这属于基本常识

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


北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

8

帖子

0

威望

63

eV
积分
71

Level 2 能力者

15#
发表于 Post on 2021-12-22 13:58:52 | 只看该作者 Only view this author
sobereva 发表于 2021-12-22 11:56
这属于基本常识

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

好的,谢谢sobereva老师

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-18 07:11 , Processed in 0.166212 second(s), 24 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list