应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 不进行optimize就直接报错收敛失败，求解决
131/2返回列表12下一页
查看: 2699  |  回复: 20
查看全部回帖 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

LinaInverse

木虫 (小有名气)

[求助] 不进行optimize就直接报错收敛失败，求解决

rt，小可做计算有段时间也有些心得了，体系很大，124个原子，元素为C,H,B,N和重元素Ir，方法b3lyp，基组6-31g（d）+lanl2dz，-3电荷，多重度3
试过opt=tight，没用，iop（5/13）=1没试但也不想试
自洽场可以收敛，但是一圈优化都没跑就直接报错收敛失败，求大侠指点
以下是报错语句：
Symmetry not used in FoFCou.
Initial guess <Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 1.0000 <S**2>= 2.0000 S= 1.0000
Defaulting to unpruned grid for atomic number  77.
Defaulting to unpruned grid for atomic number  77.
Defaulting to unpruned grid for atomic number  77.
Defaulting to unpruned grid for atomic number  77.
Defaulting to unpruned grid for atomic number  77.
Defaulting to unpruned grid for atomic number  77.
Defaulting to unpruned grid for atomic number  77.
Defaulting to unpruned grid for atomic number  77.
EnCoef did   100 forward-backward iterations
EnCoef did     1 forward-backward iterations
EnCoef did     5 forward-backward iterations
EnCoef did     2 forward-backward iterations
Restarting incremental Fock formation.
Matrix for removal  4 Erem= -2448.69634785012     Crem= 0.000D+00
Matrix for removal  3 Erem= -2454.71182296237     Crem= 0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UB3LYP) =  -2458.97111905     A.U. after  129 cycles
             Convg  =    0.4215D-04             -V/T =  2.0365
<Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 1.0000 <S**2>= 2.0085 S= 1.0028
<L.S>= 0.000000000000E+00
Annihilation of the first spin contaminant:
S**2 before annihilation     2.0085,   after     2.0000
Convergence failure -- run terminated.
Error termination via Lnk1e in /gpfssan1/home/hyan//g09/l502.exe at Wed Jul  3 19:29:12 2013.
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
秀逗魔导师
1楼 2013-07-04 15:25:22
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

LinaInverse

木虫 (小有名气)
引用回帖:
2楼: Originally posted by Jasminer at 2013-07-04 15:41:32
这个就是自洽场没有收敛。

自洽场不收敛会直接写出SCF Convergence failure的吧
而且这个上面已经写了SCF Done了
赞一下 回复此楼
秀逗魔导师
3楼2013-07-04 17:07:08
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

LinaInverse

木虫 (小有名气)
引用回帖:
4楼: Originally posted by t13340033021 at 2013-07-04 17:13:27
试一下calcfc关键词...

一圈优化都没跑，calcfc关键词能发挥作用吗？
赞一下 回复此楼
秀逗魔导师
5楼2013-07-04 19:17:09
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

LinaInverse

木虫 (小有名气)
引用回帖:
6楼: Originally posted by t13340033021 at 2013-07-04 19:21:09
先解决SCF不收敛
SCF=QC
调整下初始结构
有重金属请使用ECP基组
全电子通常不容易收敛...

主题贴上写了lanl2dz赝势基组了，你猜是用在哪个原子上？
而且问题不是SCF不收敛，如果是scf不收敛报错语句会不一样，最典型的就是经过多少多少圈优化后SCF不收敛这个句子，可是log里没有
个人经验scf=qc对付不收敛得不偿失而且也不是每次都会奏效，iop（5/13=1）对付自洽场不收敛有用但实质是跳过自洽场检验，没解决根本问题

初始结构单线态条件下优化过，正在算频率看会不会有虚频，但我们组做过很多类似配合物，构型问题恐怕不大
赞一下 回复此楼
秀逗魔导师
7楼2013-07-04 22:48:02
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

LinaInverse

木虫 (小有名气)
引用回帖:
8楼: Originally posted by t13340033021 at 2013-07-05 09:01:15
>>>>>>>>>> Convergence criterion not met.
这句话楼主看不懂吗？
帮你出谋划策的人可能会忽视一些信息，但是你自己得上心些

收敛失败和自洽场不收敛是两回事
报错语句我都遇到过，明显不同，这个不是自洽场不收敛的报错语句
赞一下 回复此楼
秀逗魔导师
9楼2013-07-05 10:23:15
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

LinaInverse

木虫 (小有名气)
引用回帖:
8楼: Originally posted by t13340033021 at 2013-07-05 09:01:15
>>>>>>>>>> Convergence criterion not met.
这句话楼主看不懂吗？
帮你出谋划策的人可能会忽视一些信息，但是你自己得上心些

不知道你哪来这样的认为别人怎么样的想法，从一开始用ecp的建议可以看出至少你在应助的时候是不上心的，我自己贴出来的报错信息我没逐句看的可能性你认为大么？
退一步讲即使是你说的问题，SCF=QC也不是个靠谱的解决方法
赞一下 回复此楼
秀逗魔导师
11楼2013-07-05 18:59:03
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

LinaInverse

木虫 (小有名气)
引用回帖:
10楼: Originally posted by yypcnbe at 2013-07-05 13:37:43
你确定是收敛了，这句话的意思只是表示129次后SCF迭代完成了，但没有达到收敛标准，现在才达到scfcon=4,你设置的肯定不止4这么低吧，试试scfcyc=500,或者scf=qc。...

见我之前回复，自洽场不收敛我碰到过不止一次，不是这个报错语句，标志性的就是经过多少圈后自洽场不收敛，而整个log里都没有这句话
而且不论是scf=qc还是scfcyc=500用来解决自洽场不收敛真心不靠谱，自洽场默认的是150圈还是多少我记不确定但是加次数基本没用，我现有资源是8核12G内存尚且嫌qc成本太高，贴个我以前的求助帖，都是有实践基础的
http://muchong.com/bbs/viewthread.php?tid=5157876
赞一下 回复此楼
秀逗魔导师
12楼2013-07-05 19:07:42
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

LinaInverse

木虫 (小有名气)
引用回帖:
13楼: Originally posted by yypcnbe at 2013-07-06 00:06:36
自洽场默认的迭代次数是129次，你看你的那个报错就是129次后没有达到收敛标准，SCF不收敛的话一般加大迭代次数和SCF=qc都可以达到收敛，告诉你另一个命令，就是scf=xqc，只有在迭代次数内没有达到收敛，才会用qc方 ...

就我个人试过的案例，加大圈数和scf=qc或者scf=xqc都是坑货，这三种方法我在一年前的求助帖中就知道并且试用过了
http://muchong.com/bbs/viewthread.php?tid=4007785
个人试过的方法里有用的只有两种：寻找更合理的初始构型和设定或者iop（5/13=1）（极不推荐）
我的版本是g09，不知道是不是你说的这个默认设定，但是加自洽场圈数或者scf=qc，xqc我是绝对不会去试的。检查初始构型和体系设定要靠谱的多
赞一下 回复此楼
秀逗魔导师
14楼2013-07-06 00:20:08
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

LinaInverse

木虫 (小有名气)
引用回帖:
15楼: Originally posted by Jasminer at 2013-07-06 21:51:00
诚然有些回帖没有特别上心，但是lz你这不信那不信一副lz都是对的的样子，那还来求助个p啊~
如你所愿，你就继续去寻找所谓“解决opt收敛失败”的方法吧。

呵呵，已解决，不是以上说的任何一种方法

求助是希望看到穿透性的见解和个人经验的交流，而不是程式化的一个命令语句一个结论，我敢说那些回复加自洽场圈数或者scf=qc的自己都没这么试过几次，而我已经碰到很多次了。如果我觉得求助者实践经历没那么充分，拒绝其意见不应该是很正常也不算不通情理的做法？我的看法肯定不是绝对正确的，那应助者就是么？

不要一副高高在上的样子。虽然这话难听，但这个帖子里的应助者和我曾经求助时碰上的一些大牛，差距还是明显的
赞一下 回复此楼
秀逗魔导师
16楼2013-07-07 00:13:52
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

LinaInverse

木虫 (小有名气)
引用回帖:
17楼: Originally posted by Jasminer at 2013-07-07 12:43:34
呵呵，既然如此，把你的解决方法拿出来把大家彻底打败呗，体现一下“明显的差距”。

说实在的，回帖的人里包括我确实没有一个所谓的“大牛”，你以为为什么？因为你的问题实（qi）践（shi）性（ruo）很（bao）强 ...

如果是交流没问题，但你的语气里没有一丝一毫这样的意味，说白了就是不屑，此等宵小伎俩无需理会，如果还腆着个脸求打的话，答案在我之前的回复里

确实这不是什么高等问题，但是高等问题你玩不来，就是这个问题，以你提出的方法也解决不了，而你在这个问题上体现的思考能力和经验积累，原话奉还
实（qi）践（shi）性（ruo）很（bao）强（le）。
赞一下 回复此楼
秀逗魔导师
18楼2013-07-07 16:30:07
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 LinaInverse 的主题更新
131/2返回列表12下一页
信息提示
关闭
请填处理意见
关闭