24小时热门版块排行榜    

查看: 1803  |  回复: 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的回帖

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的回帖

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 ]
据权威报纸调查,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系列基组,不是很常用,有很多比这好的基组的。
据权威报纸调查,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输入还是不同,
优化所得结构差别不是很大,但能 ...

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计算的分子很大的话。
据权威报纸调查,84%的青少年犯罪者有玩电子游戏的经历,故应该禁止经营电子游戏。据我们所知,100%的青少年犯罪者有穿鞋的经历,故制鞋厂应该停业。
10楼2009-10-10 14:38:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gongxd325 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 CyTWnZiM0cS6 2026-08-02 3/150 2026-08-02 20:44 by 3VGOm6paSgqv
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 CyTWnZiM0cS6 2026-08-02 5/250 2026-08-02 20:43 by 3VGOm6paSgqv
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 AS6fX1tr1on4 2026-08-02 14/700 2026-08-02 20:11 by 3VGOm6paSgqv
[基金申请] 面上提前没消息,有中的吗 +11 archvillain 2026-08-02 13/650 2026-08-02 20:09 by symmetry
[基金申请] 微信指数没变化,科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 AS6fX1tr1on4 2026-08-02 13/650 2026-08-02 20:04 by 3VGOm6paSgqv
[找工作] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 cBGhfJFuNMca 2026-08-02 13/650 2026-08-02 20:04 by 3VGOm6paSgqv
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 cBGhfJFuNMca 2026-08-02 17/850 2026-08-02 20:03 by 3VGOm6paSgqv
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 cBGhfJFuNMca 2026-08-02 15/750 2026-08-02 19:54 by 3VGOm6paSgqv
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 AS6fX1tr1on4 2026-08-02 8/400 2026-08-02 18:44 by ca3DVAx4OB19
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 AS6fX1tr1on4 2026-08-02 9/450 2026-08-02 18:27 by ca3DVAx4OB19
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 AS6fX1tr1on4 2026-08-02 8/400 2026-08-02 18:24 by ca3DVAx4OB19
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 cBGhfJFuNMca 2026-08-02 20/1000 2026-08-02 18:04 by ca3DVAx4OB19
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 15/750 2026-08-02 16:23 by jHcrHMhbKCFV
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +18 低垂的野花 2026-07-31 24/1200 2026-08-02 16:03 by swallowwh
[基金申请] 2026年国自然面上资助率 +15 布布和一二 2026-07-30 21/1050 2026-08-02 10:56 by iaeyuan
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 10/500 2026-08-02 10:00 by 锦衣卫寒战
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
信息提示
请填处理意见