24小时热门版块排行榜    

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

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

beefly

专家顾问 (职业作家)

地沟油冶炼专家

在手册的基组部分,注意这句话:
STO-3G and 3-21G accept a * suffix, but this does not actually add any polarization functions.

所以3-21g的极化函数只能自定义
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
5楼2012-07-07 00:10:22
已阅   关注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的回帖

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的极化函数只能自定义

还有就是怎么送金币啊  虽然不多 一点心意吧   我好想就10个金币 哈哈
7楼2012-07-09 14:06:58
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 daixing2011 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 23jxep3nCNZb 2026-09-14 5/250 2026-09-18 11:28 by HjTo14L7GlnE
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-18 10:30 by yS4CkeWDhD6V
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 s3fFTmArrBt6 2026-09-14 6/300 2026-09-18 10:28 by yS4CkeWDhD6V
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 C79jjtjAKjEn 2026-09-15 3/150 2026-09-18 03:42 by cITaGg3p5edr
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-17 22:39 by cuPZTDXS3VOd
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-17 15:29 by 0ngsZVSCTS7z
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-17 14:27 by ObVzqQQh1rrE
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 0pYnUiPDfdnk 2026-09-14 4/200 2026-09-17 14:15 by ObVzqQQh1rrE
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 40ms4Wlo4umh 2026-09-14 5/250 2026-09-17 14:05 by ObVzqQQh1rrE
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 23jxep3nCNZb 2026-09-14 5/250 2026-09-17 13:29 by ObVzqQQh1rrE
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-17 13:15 by ObVzqQQh1rrE
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-17 13:05 by ObVzqQQh1rrE
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-17 12:39 by ObVzqQQh1rrE
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-17 05:38 by k06BKNblrGRH
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-17 03:14 by Tql5LQhh5rLK
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 6F5UbRU2I5hL 2026-09-14 6/300 2026-09-16 17:27 by i3KBRD19ERkG
[考博] 上海工程技术大学激光智能制造课题组|2027级博士研究生招生公告 +7 水士口 2026-09-16 8/400 2026-09-16 13:47 by kkklens
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +8 两三岁ss 2026-09-14 8/400 2026-09-16 09:46 by 水士口
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-16 03:14 by DgNGHc3h5tPl
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
信息提示
请填处理意见