24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1318  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +7 张zic 2026-04-05 8/400 2026-04-05 22:57 by Hdyxbekcb
[考研] 332求调剂 +14 小小孟... 2026-04-05 15/750 2026-04-05 22:27 by hemengdong
[考研] 考研生物与医药调剂 +8 铁憨憨123425 2026-03-31 8/400 2026-04-05 20:28 by lys0704
[考研] 材料专硕 调剂 +14 CXN123456 2026-04-03 14/700 2026-04-05 17:18 by Hdyxbekcb
[考研] 282求调剂 +3 aaa车辆 2026-04-01 3/150 2026-04-05 17:03 by yulian1987
[考研] 一志愿上海海洋大学083200食品学硕,求调剂,接受其他专业083200 +4 what张 2026-04-04 5/250 2026-04-05 14:07 by chw1980_0
[考研] 346分的生物与医药08600求调剂 +5 常雨阳上岸 2026-04-05 6/300 2026-04-05 13:42 by imissbao
[考研] 一志愿江南大学085501机械工程专硕326分,本科佳木斯大学 +5 顾若浮生 2026-04-03 9/450 2026-04-05 09:57 by 1753564080
[考研] 材料调剂 +15 一样YWY 2026-04-01 15/750 2026-04-04 22:23 by hemengdong
[考研] 一志愿武理材料工程302调剂环化或化工 +19 Doleres 2026-03-31 20/1000 2026-04-04 16:44 by 啊俊!
[考研] 288求调剂 一志愿哈工大 材料与化工 +39 洛神哥哥 2026-03-31 41/2050 2026-04-03 21:51 by qlm5820
[考研] 338求调剂 +4 zzz,,r 2026-04-03 4/200 2026-04-03 16:39 by lijunpoly
[考研] 求调剂 +9 akdhjs 2026-03-31 11/550 2026-04-03 13:32 by akdhjs
[考研] 一志愿中国科学院大学265求调剂 +9 恬淡ye 2026-03-31 10/500 2026-04-03 11:10 by txp1986
[考研] 专硕 351 086100 也是考的材科基 本科也是材料 +8 202451007219 2026-04-02 8/400 2026-04-03 09:50 by 蓝云思雨
[考研] 07生物学求调剂 一志愿同济大学359分 +3 LAMC. 2026-03-30 3/150 2026-04-02 10:26 by 18828373951
[考研] 292求调剂 +17 木虫er12138 2026-04-01 17/850 2026-04-01 21:37 by 七度不信任
[考研] 285求调剂 +7 AZMK 2026-03-30 13/650 2026-04-01 17:00 by 七度不信任
[考研] 材料专业调剂 +5 啦啦啦哭 2026-03-31 6/300 2026-04-01 16:48 by JourneyLucky
[考研] 254材料与化工求调剂 +3 翰冬林楠 2026-03-30 4/200 2026-03-31 17:53 by yishunmin
信息提示
请填处理意见