24小时热门版块排行榜    

查看: 1529  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +5 yufeiwaner 2026-08-09 5/250 2026-08-10 08:15 by BlakeReary
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +7 Lanmanbaby 2026-08-09 14/700 2026-08-09 23:33 by Tide man
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 5/250 2026-08-08 21:54 by 13586093586
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 好奇怪的filecode +3 布布和一二 2026-08-08 4/200 2026-08-08 17:46 by zhanghaozhu
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见