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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】CEP-121G基组收敛，CEP-121G*，SCF迭代就不收敛了，谁知道怎么回事？
162/2返回列表上一页12
查看: 1827  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

gnli

铜虫 (小有名气)
引用回帖:
Originally posted by manson1998 at 2010-09-19 08:49:17:


加上ecp

pbepbe/genecp  优化+频率计算
Na Si O 0
CEP-121G*
****
Ti H 0
CEP-121G
****

谢谢您的建议 我尝试一下
请教一下啊：
gen 和 genecp 不是一样的效果吗？
除非额外定义赝势的时候 才有区别
不知道我理解的对不对
赞一下 回复此楼
争取早日成为合格的菜鸟
11楼2010-09-20 03:51:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家
如果使用ecp基组，genecp = gen + pseudo。ecp基组需要作为价电子部分和芯电子部分，分别输入。
赞一下 回复此楼
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
12楼2010-09-20 12:06:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

manson1998

金虫 (著名写手)

OLED达人
gnli(金币+5): 2010-09-20 19:08:48
引用回帖:
Originally posted by gnli at 2010-09-20 03:51:34:

谢谢您的建议 我尝试一下
请教一下啊：
gen 和 genecp 不是一样的效果吗？
除非额外定义赝势的时候 才有区别
不知道我理解的对不对

如果你不加ecp的话 应该会发现ALPHA电子增加了很多(相比较与加了ecp)。实际上也就是计算了非价电子的内层电子，造成计算量增加以及体系能量的不真实下降。所以一般要加上ecp，在以前的版本里是加上pseudo=read。
赞一下 回复此楼
惺惺相惜情不自禁
13楼2010-09-20 17:43:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gnli

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


如果你不加ecp的话 应该会发现ALPHA电子增加了很多(相比较与加了ecp)。实际上也就是计算了非价电子的内层电子，造成计算量增加以及体系能量的不真实下降。所以一般要加上ecp，在以前的版本里是加上pseudo=read。

谢谢
但是我现在用的机组 cep-121G 和 cep-121G*
并没有涉及到赝势的问题啊
只是不同的元素采取了不同的机组 cep-121G* 机组不是个赝势基组吧？
赞一下 回复此楼
争取早日成为合格的菜鸟
14楼2010-09-20 19:08:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

manson1998

金虫 (著名写手)

OLED达人
引用回帖:
Originally posted by gnli at 2010-09-20 19:08:27:

谢谢
但是我现在用的机组 cep-121G 和 cep-121G*
并没有涉及到赝势的问题啊
只是不同的元素采取了不同的机组 cep-121G* 机组不是个赝势基组吧？

cep121G基组的定义是Stephens/Basch/Krauss 的三重分裂价电子ECP 基组。
我觉得你之前因为是统一基组运算，gaussian默认这样的写法，类似B3LYP/lanl2dz,可能没有什么问题。
一旦混合基组后就必须要加上ecp或者pseudo=read了。欢迎继续讨论
赞一下 回复此楼
惺惺相惜情不自禁
15楼2010-09-20 21:54:16
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 gnli 的主题更新
162/2返回列表上一页12
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭