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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【讨论】基组在关键词行输入和用GEN输入的差别
北京石油化工学院2026年研究生招生接收调剂公告
161/2返回列表12下一页
查看: 1546  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 328求调剂 +7 嗯滴的基本都 2026-03-27 7/350 2026-03-28 04:19 by fmesaito
[考研] 0856材料化工调剂 总分330 +10 zhubinhao 2026-03-27 10/500 2026-03-28 03:34 by fmesaito
[考研] 张芳铭-中国农业大学-环境工程专硕-298 +3 手机用户 2026-03-26 3/150 2026-03-27 21:39 by lbsjt
[考研] 材料292调剂 +12 橘颂思美人 2026-03-23 12/600 2026-03-27 15:44 by caszguilin
[考研] 279 分 求调剂 +4 睡个好觉_16 2026-03-24 4/200 2026-03-27 15:05 by 醉在风里
[考研] 复试调剂，一志愿南农083200食品科学与工程 +5 XQTJZ 2026-03-26 5/250 2026-03-27 14:49 by 狂炫麦当当
[考研] 08开头275求调剂 +4 拉谁不重要 2026-03-26 4/200 2026-03-27 14:12 by Delta2012
[材料工程] 一志愿C9材料与化工专业总分300求调剂 +7 曼111 2026-03-24 8/400 2026-03-27 13:53 by 杨杨杨紫
[考研] 298调剂 +3 jiyingjie123 2026-03-27 3/150 2026-03-27 11:57 by wxiongid
[考研] 0703化学一志愿南京师范大学303求调剂 +3 zzffylgg 2026-03-24 3/150 2026-03-27 10:42 by shangxh
[考研] 071000生物学求调剂，初试成绩343 +6 小小甜面团 2026-03-25 6/300 2026-03-26 23:01 by 不吃魚的貓
[考研] 一志愿北化求调剂 +3 Jsman 2026-03-22 3/150 2026-03-26 21:06 by ajpv风雷
[考研] 321求调剂 +6 wasdssaa 2026-03-26 6/300 2026-03-26 20:57 by sanrepian
[考研] 085601求调剂总分293英一数二 +4 钢铁大炮 2026-03-24 4/200 2026-03-26 16:28 by dick_runner
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研经历 +6 cccchenso 2026-03-23 6/300 2026-03-25 22:25 by 544594351
[考研] 求调剂 +3 QiMing7 2026-03-25 3/150 2026-03-25 21:13 by 给你你注意休息
[考研] 284求调剂 +15 Zhao anqi 2026-03-22 15/750 2026-03-25 12:51 by wht0531
[考研] 求调剂 +6 研研，接电话 2026-03-24 7/350 2026-03-24 17:01 by barlinike
[考研] 291求调剂 +3 HanBeiNingZC 2026-03-24 3/150 2026-03-24 16:34 by barlinike
[考研] 336化工调剂 +4 王大坦1 2026-03-23 5/250 2026-03-23 18:32 by allen-yin
信息提示
关闭
请填处理意见
关闭