24小时热门版块排行榜    

查看: 2065  |  回复: 15

gnli

铜虫 (小有名气)

[交流] 【求助】CEP-121G基组收敛,CEP-121G*,SCF迭代就不收敛了,谁知道怎么回事?

做cluster优化结构和计算光谱
先用 pbepbe/CEP-121G 基组优化 一切正常
然后接着做 pbepbe/gen  优化+频率计算
Na Si O 0
CEP-121G*
****
Ti H 0
CEP-121G
****

SCF 迭代怎么都是129圈不收敛
请教是哪里的问题?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

争取早日成为合格的菜鸟
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

manson1998

金虫 (著名写手)

OLED达人

★ ★
erylingjet(金币+2):感谢交流! 2010-09-17 19:06:07
gnli(金币+2): 2010-09-17 21:35:23
最后是出错了 还是你觉得没算完?
缩小一下步长试试吧
惺惺相惜情不自禁
2楼2010-09-17 17:45:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gnli

铜虫 (小有名气)

引用回帖:
Originally posted by manson1998 at 2010-09-17 17:45:25:
最后是出错了 还是你觉得没算完?
缩小一下步长试试吧

没有出错 就是 121G 的时候每一次迭代都十几圈就收敛了
离子步也就30-50步结构就优化好了
然后接着做 121G、121G* 的时候 SCF每一圈都是129圈没收敛
现在已经跑了好多天了,几十个离子步
我觉得很奇怪
争取早日成为合格的菜鸟
3楼2010-09-17 17:48:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gnli

铜虫 (小有名气)

引用回帖:
Originally posted by manson1998 at 2010-09-17 17:45:25:
最后是出错了 还是你觉得没算完?
缩小一下步长试试吧

这个是 pbepbe/gen
Na Si  CEP -121G
Ti H   CEP-121G*
的迭代结果:

SCF Done:  E(RPBE-PBE) =  -6182.66781085     A.U. after  129 cycles
SCF Done:  E(RPBE-PBE) =  -6186.78468300     A.U. after  129 cycles
SCF Done:  E(RPBE-PBE) =  -6192.78471931     A.U. after  129 cycles
SCF Done:  E(RPBE-PBE) =  -6198.13035441     A.U. after  129 cycles
SCF Done:  E(RPBE-PBE) =  -6205.31909983     A.U. after  129 cycles
SCF Done:  E(RPBE-PBE) =  -6211.92714424     A.U. after  129 cycles
SCF Done:  E(RPBE-PBE) =  -6217.78138323     A.U. after  129 cycles
SCF Done:  E(RPBE-PBE) =  -6222.58569838     A.U. after  129 cycles
SCF Done:  E(RPBE-PBE) =  -6228.00804563     A.U. after  129 cycles
SCF Done:  E(RPBE-PBE) =  -6234.51769181     A.U. after  129 cycles
SCF Done:  E(RPBE-PBE) =  -6241.38936327     A.U. after  129 cycles
SCF Done:  E(RPBE-PBE) =  -6247.68355529     A.U. after  129 cycles
SCF Done:  E(RPBE-PBE) =  -6248.59345416     A.U. after  129 cycles
SCF Done:  E(RPBE-PBE) =  -6249.62935377     A.U. after  129 cycles
SCF Done:  E(RPBE-PBE) =  -6250.65612293     A.U. after  129 cycles
SCF Done:  E(RPBE-PBE) =  -6251.67363038     A.U. after  129 cycles
SCF Done:  E(RPBE-PBE) =  -6252.68249702     A.U. after  129 cycles
SCF Done:  E(RPBE-PBE) =  -6253.52724385     A.U. after  129 cycles
SCF Done:  E(RPBE-PBE) =  -6253.56767553     A.U. after  129 cycles
SCF Done:  E(RPBE-PBE) =  -6256.94645007     A.U. after  129 cycles
SCF Done:  E(RPBE-PBE) =  -6261.05238072     A.U. after  129 cycles

[ Last edited by gnli on 2010-9-17 at 17:51 ]
争取早日成为合格的菜鸟
4楼2010-09-17 17:49:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lihb734

铁杆木虫 (职业作家)

站在计算化学入门的门槛上

★ ★
erylingjet(金币+2):http://www.yydxw.com 2010-09-17 19:06:15
gnli(金币+4): 2010-09-17 21:35:32
方案一:读入pbepbe/CEP-121G 的初始猜测和几何之和重算
方案二:增大SCF的圈数
前途光明,出路难觅!
5楼2010-09-17 18:19:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lihb734

铁杆木虫 (职业作家)

站在计算化学入门的门槛上

★ ★
heyo_123(金币+2):鼓励!! 2010-09-19 09:38:38
erylingjet(金币+2):http://www.yydxw.com 2010-09-17 19:06:15
斑竹竟然打广告啊。
前途光明,出路难觅!
6楼2010-09-17 19:36:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gnli

铜虫 (小有名气)

引用回帖:
Originally posted by lihb734 at 2010-09-17 18:19:40:
方案一:读入pbepbe/CEP-121G 的初始猜测和几何之和重算
方案二:增大SCF的圈数

我是保存了 cep-121G 的chk 文件
然后用 guess=read 和 geom=allcheck 接着算的
cep 121G 和 121G* 基组大小差别不大
所以我一开始以为很快就应该优化结束的
没想到会出现这样的情况
您觉得增大 scf 迭代的次数会有效果吗
感觉像是别的地方也有问题
争取早日成为合格的菜鸟
7楼2010-09-17 20:33:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

没有指定ecp?先比较一下单点能吧
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
8楼2010-09-18 00:31:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gnli

铜虫 (小有名气)

引用回帖:
Originally posted by beefly at 2010-09-18 00:31:12:
没有指定ecp?先比较一下单点能吧

ecp是什么意思啊
争取早日成为合格的菜鸟
9楼2010-09-18 22:34:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

manson1998

金虫 (著名写手)

OLED达人


heyo_123(金币+1):鼓励交流~!~ 2010-09-19 09:38:06
gnli(金币+4): 2010-09-20 03:52:58
引用回帖:
Originally posted by gnli at 2010-09-17 20:33:08:

我是保存了 cep-121G 的chk 文件
然后用 guess=read 和 geom=allcheck 接着算的
cep 121G 和 121G* 基组大小差别不大
所以我一开始以为很快就应该优化结束的
没想到会出现这样的情况
您觉得增大 scf 迭代 ...

加上ecp

pbepbe/genecp  优化+频率计算
Na Si O 0
CEP-121G*
****
Ti H 0
CEP-121G
****
惺惺相惜情不自禁
10楼2010-09-19 08:49:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gnli 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] H口,面上,时间戳7月7日。 +3 plumage5 2026-07-15 3/150 2026-07-15 21:06 by 6543yes
[有机交流] 伯胺上甲基求助 100+3 002tmumu 2026-07-14 3/150 2026-07-15 20:37 by 盛世杂家
[基金申请] 小木虫没落了,除了祈祷帖子,几乎看不到有价值的帖子 +13 botar 2026-07-14 14/700 2026-07-15 19:28 by 寒流32
[基金申请] 会评 +5 无名者登山 2026-07-15 5/250 2026-07-15 18:15 by 杠就是你对
[基金申请] 生生命学部会评 +3 wuke100666 2026-07-15 3/150 2026-07-15 17:11 by kingkocxr
[基金申请] E口会评 +11 布布和一二 2026-07-13 13/650 2026-07-15 16:34 by Kittylucky
[基金申请] 生命口C13面上会评时间 +6 luckinging 2026-07-14 9/450 2026-07-15 14:39 by 雨冰共舞
[教师之家] 内心匮乏 +13 水冰月月野兔 2026-07-12 13/650 2026-07-15 14:21 by jiduquegai
[文学芳草园] 你可曾有过这种感受? +3 Jeanya 2026-07-09 6/300 2026-07-15 11:07 by mbtwt
[论文投稿] 农业机械学报投稿周期 +3 G0DLIN 2026-07-14 7/350 2026-07-15 09:26 by xs74101122
[论文投稿] 跨出版社商投稿 20+4 倪好520 2026-07-13 4/200 2026-07-15 08:35 by bobvan
[基金申请] 没收到消息,再次陪跑 +3 天空星辰fly 2026-07-14 3/150 2026-07-14 19:14 by 6543yes
[基金申请] 不要再数国自然申请书的 filecode 的分隔符个数了 +13 sadpencil 2026-07-12 20/1000 2026-07-14 13:12 by 杠就是你对
[基金申请] 明天E口面上会评 +8 布布和一二 2026-07-12 9/450 2026-07-13 22:59 by QTSorange
[论文投稿] MSER送审了还被拒稿 +4 S778950906 2026-07-08 6/300 2026-07-13 16:22 by tfang
[基金申请] 祈祷青基必中 50+6 hadengry 2026-07-09 16/800 2026-07-12 20:31 by 登山虫虫
[考博] 27届辽宁大学应届毕业生申博 +3 可乐微风大海 2026-07-09 3/150 2026-07-10 21:23 by 贪玩的研究僧
[基金申请] 关于如何从代码看上不上会 +17 且听虎啸 2026-07-09 23/1150 2026-07-10 19:51 by sdfapple719
[有机交流] chemdraw 20+5 ?慕鱼? 2026-07-08 6/300 2026-07-10 15:37 by 紫枫星域
[考研] 在职考研 +6 wu5d 2026-07-08 7/350 2026-07-09 16:40 by 化工小霸王呀
信息提示
请填处理意见