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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【讨论】基组在关键词行输入和用GEN输入的差别
161/2返回列表12下一页
查看: 1491  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

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 ]
回复此楼

» 猜你喜欢
一个人做点好事并不难难的是一辈子做好事!!!
1楼 2009-10-09 10:23:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abbott

金虫 (著名写手)

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

另外 仔细读读说明书， 基组说明部分....
赞一下(11人) 回复此楼
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即可
赞一下(11人) 回复此楼
3楼2009-10-09 12:53:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

a349857460

木虫 (著名写手)
★ ★
lei0736(金币+1,VIP+0):谢谢 10-9 21:19
gongxd325(金币+1,VIP+0):Thanks 10-10 08:33
加pesudo=read
回复此楼
4楼2009-10-09 20:49:15
已阅   回复此楼   关注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 ]
赞一下(10人) 回复此楼
据权威报纸调查，84%的青少年犯罪者有玩电子游戏的经历，故应该禁止经营电子游戏。据我们所知，100%的青少年犯罪者有穿鞋的经历，故制鞋厂应该停业。
5楼2009-10-09 20:51:29
已阅   回复此楼   关注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的回帖

hairan

木虫 (著名写手)

瞌睡木虫
★ ★ ★ ★ ★
gongxd325(金币+2,VIP+0):谢谢你的建议 10-10 10:48
lei0736(金币+3,VIP+0):谢谢 10-10 12:32
引用回帖:
Originally posted by gongxd325 at 2009-10-10 09:00:
谢谢楼上各位，用两种方式输入纯粹为测试目的，因为遇到了问题想通过测试解决。

采用hairan建议的输入方式B3LYP/GENECP计算，得到的基函数数目与B3LYP/CEP-121G输入还是不同，
优化所得结构差别不是很大，但能 ...

这个问题要彻底解决需要你自己去看文献，Gaussian说明书中有CEP系列基组的参考文献，你下载了仔细看看，看到底他们对那些元素使用了赝势，那些没有使用。
一般出现这么大的能量差别，肯定是和赝势基组的使用有关。
另外，为什么要用CEP系列基组，不是很常用，有很多比这好的基组的。
赞一下(11人) 回复此楼
据权威报纸调查，84%的青少年犯罪者有玩电子游戏的经历，故应该禁止经营电子游戏。据我们所知，100%的青少年犯罪者有穿鞋的经历，故制鞋厂应该停业。
7楼2009-10-10 10:28:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gongxd325

至尊木虫 (文坛精英)
由于计算的体系含不同的过渡元素，为了能使用统一的基组，只有CEP, SDD, Lan和UGBS可选，Lan用得普遍些，但SDD和CEP基组更大些，而SDD的膺势对不同元素用的方法似乎不同，UGBS则太大
赞一下 回复此楼
一个人做点好事并不难难的是一辈子做好事!!!
8楼2009-10-10 10:47:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jasminer

铁杆木虫 (著名写手)
★ ★ ★
gongxd325(金币+1,VIP+0):Thanks 10-10 11:22
lei0736(金币+2,VIP+0):谢谢 10-10 12:33
用gen输入的基组视为5d7f纯函数基组，直接输入的基组则各基组有不同的默认值。
如果我没记错的话，cep基组使用的是笛卡尔基组，是否是这个原因造成的能量不同。
引用回帖:
Originally posted by gongxd325 at 2009-10-10 09:00:
谢谢楼上各位，用两种方式输入纯粹为测试目的，因为遇到了问题想通过测试解决。

采用hairan建议的输入方式B3LYP/GENECP计算，得到的基函数数目与B3LYP/CEP-121G输入还是不同，
优化所得结构差别不是很大，但能 ...

赞一下(11人) 回复此楼
9楼2009-10-10 11:17:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hairan

木虫 (著名写手)

瞌睡木虫
★ ★
lei0736(金币+1,VIP+0):谢谢 10-10 15:32
gongxd325(金币+1,VIP+0): 10-10 17:57
引用回帖:
Originally posted by Jasminer at 2009-10-10 11:17:
用gen输入的基组视为5d7f纯函数基组，直接输入的基组则各基组有不同的默认值。
如果我没记错的话，cep基组使用的是笛卡尔基组，是否是这个原因造成的能量不同。

这倒也是一种可能的原因，如果lz计算的分子很大的话。
赞一下(10人) 回复此楼
据权威报纸调查，84%的青少年犯罪者有玩电子游戏的经历，故应该禁止经营电子游戏。据我们所知，100%的青少年犯罪者有穿鞋的经历，故制鞋厂应该停业。
10楼2009-10-10 14:38:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gongxd325 的主题更新
161/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭