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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】CEP-121G基组收敛，CEP-121G*，SCF迭代就不收敛了，谁知道怎么回事？
162/2返回列表上一页12
查看: 1906  |  回复: 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
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料专硕306英一数二 +3 z1z2z3879 2026-03-16 3/150 2026-03-16 09:19 by Demonsssss
[考研] 290求调剂 +5 孔志浩 2026-03-12 10/500 2026-03-16 09:01 by 余晖&
[考研] 085600材料与化工 求调剂 +12 enenenhui 2026-03-13 13/650 2026-03-16 08:30 by Linda Hu
[考研] 化学调剂0703 +7 啊我我的 2026-03-11 7/350 2026-03-15 23:03 by 凌千颂111
[考研] 调剂 +8 调剂的考研学生 2026-03-09 8/400 2026-03-15 22:14 by Winj1e
[考博] 欢迎申博同学联系 +3 天道酬勤2026686 2026-03-10 7/350 2026-03-15 19:03 by 天道酬勤2026686
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 308求调剂 +4 是Lupa啊 2026-03-09 4/200 2026-03-14 02:06 by tranquil_ya
[考研] 328，0703考生求调剂，一志愿为东北师范大学 +4 观素律 2026-03-09 5/250 2026-03-14 01:24 by JourneyLucky
[考研] 一志愿郑大070303，338分，求调剂 +4 dadawaf 2026-03-10 5/250 2026-03-14 01:20 by lsw010101
[考研] 四川大学085601材料工程专硕 初试294求调剂 +4 祝我们好在冬天 2026-03-11 4/200 2026-03-13 21:39 by peike
[考研] 315求调剂 +9 小羊小羊_ 2026-03-11 10/500 2026-03-13 21:13 by SXNU李老师
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 307求调剂 +5 超级伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 0703一志愿211 285分求调剂 +4 ly3471z 2026-03-13 4/200 2026-03-13 13:00 by JourneyLucky
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
信息提示
关闭
请填处理意见
关闭