24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1319  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 327求调剂 +4 拾光任染 2026-04-05 4/200 2026-04-05 20:16 by 南航~万老师
[考研] 367求调剂 +3 谢28 2026-03-30 3/150 2026-04-05 13:27 by huangmoli
[有机交流] 甲醇/二氯 1:15过柱子 5+3 a哎y呦w喂 2026-03-31 3/150 2026-04-05 10:42 by 88817753
[考研] 0703调剂 +11 拾玖壹 2026-04-04 12/600 2026-04-05 10:29 by 果冻大王
[考研] 0854求调剂 +4 assdll 2026-04-04 4/200 2026-04-05 09:44 by zhq0425
[考研] 本科211 分数293请求调剂 +4 莲菜就是藕吧 2026-04-01 4/200 2026-04-04 22:32 by hemengdong
[考研] 材料工程085601数二英一335求调剂 +6 双马尾痞老板2 2026-03-31 6/300 2026-04-04 22:29 by hemengdong
[考研] 320求调剂 +3 一样圆 2026-04-04 3/150 2026-04-04 22:29 by 啵啵啵0119
[考研] 085600调剂 +4 1amJJ 2026-04-02 4/200 2026-04-04 21:53 by hemengdong
[考研] 调剂 +9 19945159693 2026-04-03 10/500 2026-04-04 20:16 by dongzh2009
[考研] 278求调剂 +6 Yy7400 2026-04-03 6/300 2026-04-04 09:53 by zhangdingwa
[考研] 266求调剂 +18 阳阳哇塞 2026-04-01 18/900 2026-04-03 18:38 by zllcz
[考研] 求调剂机会 +5 意染ivy 2026-04-03 5/250 2026-04-03 15:13 by qoooooo614
[考研] 085801 总分275 本科新能源 求调剂 +19 bradoner 2026-04-01 23/1150 2026-04-03 10:07 by linyelide
[考博] 26年申博 +3 staryer 2026-03-30 4/200 2026-04-01 23:21 by ai4pharm
[考研] 070300化学279求调剂 +15 哈哈哈^_^ 2026-03-31 17/850 2026-04-01 21:37 by 给你你注意休息
[考研] 311求调剂 +10 李芷新1 2026-03-31 10/500 2026-04-01 14:38 by chenqifeng666
[考研] 085601 329分调剂 +6 yzsa12 2026-03-31 6/300 2026-03-31 15:23 by yanflower7133
[考研] 323分 食品与营养调剂 +3 嘿ooo 2026-03-31 3/150 2026-03-31 09:38 by longlotian
[考研] 297 地理学070500 复试求调剂 +3 小圆圈圈ooo 2026-03-30 3/150 2026-03-30 21:05 by 余震yz
信息提示
请填处理意见