24小时热门版块排行榜    

查看: 1820  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见