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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】CEP-121G基组收敛，CEP-121G*，SCF迭代就不收敛了，谁知道怎么回事？
52/1返回列表
查看: 1872  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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圈不收敛
请教是哪里的问题？
回复此楼

» 猜你喜欢
争取早日成为合格的菜鸟
1楼 2010-09-17 17:36:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gnli

铜虫 (小有名气)
引用回帖:
Originally posted by manson1998 at 2010-09-20 21:54:16:

cep121G基组的定义是Stephens/Basch/Krauss 的三重分裂价电子ECP 基组。
我觉得你之前因为是统一基组运算，gaussian默认这样的写法，类似B3LYP/lanl2dz,可能没有什么问题。
一旦混合基组后就必须要加上ecp或者 ...

我测试了如下写法，但是问题仍然存在，SCF 迭代还是不收敛，计算非常耗时
#P pbepbe/genecp opt freq  iop(5/13=1)  guess=read geom=allcheck

Na Si O 0
CEP-121G*
****
Ti H 0
CEP-121G
****

Ti H 0
CEP-121G

也测试了 ：
#P pbepbe/cep-121G*  opt freq  iop(5/13=1)  guess=read geom=allcheck
出现同样的SCF 问题

真是崩溃啊 文献上明明写的就是这么算的 为啥就不收敛呢？
请你帮忙看一下这篇文献中关于算法的描述。难道我理解错了？
JPCC（2007， 111，  1776-1782）method 部分 第二段
非常感谢 我们可以继续讨论吗？

[ Last edited by gnli on 2010-9-22 at 17:48 ]
赞一下 回复此楼
争取早日成为合格的菜鸟
16楼2010-09-22 17:44:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 16 个回答

manson1998

金虫 (著名写手)

OLED达人
★ ★
erylingjet(金币+2):感谢交流！ 2010-09-17 19:06:07
gnli(金币+2): 2010-09-17 21:35:23
最后是出错了 还是你觉得没算完？
缩小一下步长试试吧
赞一下(1人) 回复此楼
惺惺相惜情不自禁
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的回帖
查看全部 16 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭