24小时热门版块排行榜    

查看: 1499  |  回复: 15
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

gongxd325

至尊木虫 (文坛精英)

[交流] 【讨论】基组在关键词行输入和用GEN输入的差别

最近用B3LYP/GEN计算含CNOAg的体系,发现个问题:

(1)B3LYP/GEN
        C N O 0
        CEP-121G
        ****
        Ag 0
        CEP-121G
        ****
    67 basis functions,   110 primitive gaussians,    70 cartesian basis functions
    Convergence failure -- run terminated.

(2)B3LYP/CEP-121G 能顺利结束
    70 basis functions,   110 primitive gaussians,    70 cartesian basis functions

按道理两种输入方法都是对所有元素采用了CEP-121G基组,但结果居然大不一样:
一个能正常完成另一个不能,而且基函数数目也不等;

两种输入方法究竟有何差别?
用gen做计算的不少,究竟结果的可靠度如何?
gen常出现不收敛的问题(Convergence failure或Inaccurate quadrature in CalDSu),如何解决?

[ Last edited by gongxd325 on 2009-10-9 at 10:28 ]
回复此楼
一个人做点好事并不难难的是一辈子做好事!!!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hairan

木虫 (著名写手)

瞌睡木虫

★ ★ ★ ★ ★ ★
lei0736(金币+4,VIP+0):谢谢 10-9 21:19
gongxd325(金币+2,VIP+0):Thanks 10-10 08:33
第一种方式,你使用了赝势基组却没有读入赝势,当然不能正常结束
估计优化过程中肯定还出现围绕Ag的键长变得异常的短

正确的使用方式是
B3LYP/genecp        或者   B3LYP/gen pseudo=read

        C N O 0
        CEP-121G
        ****
        Ag 0
        CEP-121G
        ****
                                                  !此处有一个空行
        Ag 0
        CEP-121G


第二种方式与上面的写这么多是等效的,所以如果体系内所有的原子都用一种基组,而且是Gaussian自带的基组,那么就不必使用gen相关的关键词

[ Last edited by hairan on 2009-10-9 at 20:56 ]
据权威报纸调查,84%的青少年犯罪者有玩电子游戏的经历,故应该禁止经营电子游戏。据我们所知,100%的青少年犯罪者有穿鞋的经历,故制鞋厂应该停业。
5楼2009-10-09 20:51:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 16 个回答

abbott

金虫 (著名写手)

不要用QQ问我东西

★ ★ ★
lei0736(金币+2,VIP+0):谢谢 10-9 21:19
gongxd325(金币+1,VIP+0):Thanks 10-10 08:33
清仔细看看输出文件的基组信息部分
你会知道有什么不同的....

另外 仔细读读说明书, 基组说明部分....
Chemistry[]==[]Chem[]is[]try!!!
2楼2009-10-09 11:42:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主

★ ★
loovfnd(金币+1,VIP+0):谢谢 10-9 13:17
gongxd325(金币+1,VIP+0):Thanks 10-10 08:33
呵呵,这个问题的通用解决方法:
在命令行中加入gfinput即可
3楼2009-10-09 12:53:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gongxd325

至尊木虫 (文坛精英)

谢谢楼上各位,用两种方式输入纯粹为测试目的,因为遇到了问题想通过测试解决。

采用hairan建议的输入方式B3LYP/GENECP计算,得到的基函数数目与B3LYP/CEP-121G输入还是不同,
优化所得结构差别不是很大,但能量差别非常明显(~100a.u.!),看来两种输入方式不是等效的!

此外,用B3LYP/CEP-121G GFinput输出基组信息,再用B3LYP/GENECP读入,两种方式的基函数数目也不等!

究竟为何?

[ Last edited by gongxd325 on 2009-10-10 at 09:14 ]
一个人做点好事并不难难的是一辈子做好事!!!
6楼2009-10-10 09:00:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见