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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】CEP-121G基组收敛，CEP-121G*，SCF迭代就不收敛了，谁知道怎么回事？
162/2返回列表上一页12
查看: 2065  |  回复: 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
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 生命口面上会评结束了吗大家有知道结果吗？现在还没有结果正常不？ +3 碧水00 2026-07-16 5/250 2026-07-16 07:39 by jaasxulu
[基金申请] 生生命学部会评 +4 wuke100666 2026-07-15 4/200 2026-07-16 00:12 by 啦啦小妮子
[有机交流] 伯胺上甲基求助 100+4 002tmumu 2026-07-14 4/200 2026-07-15 22:01 by 学员oN9JQD
[基金申请] H口，面上，时间戳7月7日。 +3 plumage5 2026-07-15 3/150 2026-07-15 21:06 by 6543yes
[基金申请] 小木虫没落了，除了祈祷帖子，几乎看不到有价值的帖子 +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 杠就是你对
[基金申请] E口会评 +11 布布和一二 2026-07-13 13/650 2026-07-15 16:34 by Kittylucky
[基金申请] E口青C会评结束了吗？ 40+3 小小书虫c 2026-07-14 9/450 2026-07-15 16:04 by 安静的知了
[基金申请] 生命口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
[基金申请] 国自然面上D口祈祷 +4 monkeynana 2026-07-14 4/200 2026-07-15 09:14 by 白水煮
[论文投稿] 跨出版社商投稿 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
[基金申请] 祈祷青基必中 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
信息提示
关闭
请填处理意见
关闭