24小时热门版块排行榜    

查看: 1555  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

daixing2011

铜虫 (初入文坛)

[求助] 求助 混合基组问题 很有趣 大家来看下

这问题是混合基组的  可能高人一眼就看破 但小弟新手无奈

算一个体系优化 用的lsda/3-21g 优化好之后 比方说能量 E1=10

后来用混合基组 即 给某些原子用3-21g  某些用3-21g* 我是这样写的

1 2 3 4 0
3-21g
****
5 6 7 8 0
3-21g*
****

运行之后 能算 但是结果和之前只用3-21g算的一模一样 初始构型直接就收敛 能量也完全和E1一样(当然,本次初始构型是之前优化好的,大家懂的)
于是有2个怀疑:1.   3-21g*  也太小了 不能造成影响 所以和纯3-21g结果一样
                           2.   是不是输入文件格式不对 这么写是不是它还在只用3-21g去算,根本没读入我说的3-21g*

为了验证我以上两个疑问 我做了2个实验
实验一:将混合进来的基组加重 变成6-31G*  混合基组写法同上
1 2 3 4 0
3-21g
****
5 6 7 8 0
6-31g*
****
这回算了8步 收敛了 而且能量E2=8 比E1小  看起来很有道理 毕竟基组加大了
(本文写的能量都是随便写的数字 大家别见怪)

实验二:怀疑输入格式有误 改变格式 在第一排****下面加空行 其实我也不知道怎么想的 反正就这么做了  呵呵 格式如下
1 2 3 4 0
3-21g
****

5 6 7 8 0
3-21g*
****
能看到中间的空行吧 ??????
这么写也能算 但是结果真是惊人  好家伙  构型给跑的支离破碎不说  能量超大
(但这个并未收敛
Inaccurate quadrature in CalDSu.
Error termination via Lnk1e in /users/wangzg1/g09/l502.exe at Thu Jul  5 18:49:53 2012
相信大家都懂的)

在这里 我想请教的问题是  我觉得混合基组的写法 还是不加空行是正确的   但是为什么我加了空行  程序也能运行
我想知道 如果加了空行  它这算的到底是什么啊 ?  呵呵
回复此楼
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

daixing2011

铜虫 (初入文坛)

引用回帖:
5楼: Originally posted by beefly at 2012-07-07 00:10:22
在手册的基组部分,注意这句话:
STO-3G and 3-21G accept a * suffix, but this does not actually add any polarization functions.

所以3-21g的极化函数只能自定义

追问下哦  别见怪  那这个意思 是不是说 如果想给3-21g添加极化  用3-21*是没什么用的  那如果我用3-21g(d)呢  这样如何 对了  我是要给C原子加这个   实际上在专门查基组的网站上  对于C原子  只有给出321g  321++g 等等  没有给加极化的
6楼2012-07-09 13:50:29
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

nyzhaoyin

新虫 (正式写手)

但疼  前面的就是对的  能量没变化不恨正常么
2楼2012-07-06 12:04:15
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

gmy1990

荣誉版主 (著名写手)

优秀版主优秀版主

【答案】应助回帖

感谢参与,应助指数 +1
对于你的第一中情况,你可以再做个test:用画好的截个直接用你那两种方法去优化,看能量以及其他参数是什么个情况.这样你就可以找到答案了
加空行,那么后面的就没读,所以优化你也只是对1~4几个原子进行的,5~8就没进行优化,当然分子会给弄的支离破碎,能量也很大.
3楼2012-07-06 16:22:23
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

yanrding

金虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
简单地,你log文件中搜一下primi,看看到底计算用了几个基函数。
稍微复杂一点的,输入#行中加入gfinput,然后在输出log文件中看到底用了哪些基函数。
华北电力大学(北京)数理学院,计算化学,凝聚态物理。欢迎合作!招硕士/博士研究生,招老师。
4楼2012-07-06 22:47:00
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳变了,能看出什么问题? +3 基诺咪客 2026-08-17 3/150 2026-08-17 15:25 by xinshengwu
[基金申请] 感觉是下周放榜了 +3 angus9576 2026-08-17 6/300 2026-08-17 14:32 by 启萌科技
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 6GojgJvkDudM 2026-08-16 6/300 2026-08-17 13:58 by EYK67XJMj64E
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 i7NFEVbjQMM5 2026-08-16 4/200 2026-08-17 12:10 by 4GBAYCdVQoK3
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 6GojgJvkDudM 2026-08-16 4/200 2026-08-17 11:46 by 4GBAYCdVQoK3
[论文投稿] 岩土工程学报什么时候才能终审完呐 25+3 yeager111 2026-08-10 4/200 2026-08-16 20:42 by tfang
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
信息提示
请填处理意见