24小时热门版块排行榜    

查看: 1827  |  回复: 15

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