24小时热门版块排行榜    

查看: 835  |  回复: 1

qizhuixia

新虫 (初入文坛)

[求助] 小弟不才,不知道为什么用各种基组算出来的CI的能量(氦)都是一样的,求助 已有1人参与

小弟不才,不知道为什么用各种基组算出来的CI的能量(氦)都是一样的,求助

4-31
$CONTRL SCFTYP=RHF CITYP=GUGA RUNTYP=ENERGY $END
$SYSTEM TIMLIM=1 $END
$BASIS  GBASIS=N31 NGAUSS=4 $END
$GUESS  GUESS=HUCKEL $END
!    look at all state symmetries, by using C1 symmetry
$CIDRT  GROUP=C1 IEXCIT=2 NFZC=0 NDOC=1 NVAL=1 $END
$GUGDIA NSTATE=2 $END
$DATA
He..STO-3G basis
C1
He 2.0
$END
~


6-31G*
$CONTRL SCFTYP=RHF CITYP=GUGA RUNTYP=ENERGY $END
$SYSTEM TIMLIM=1 $END
$BASIS  GBASIS=N31 NGAUSS=6 NDFUNC=1 $END
$GUESS  GUESS=HUCKEL $END
!    look at all state symmetries, by using C1 symmetry
$CIDRT  GROUP=C1 IEXCIT=2 NFZC=0 NDOC=1 NVAL=1 $END
$GUGDIA NSTATE=2 $END
$DATA
He..STO-3G basis
C1
He 2.0
$END
~

算出来的能量都是样的
STATE #    1  ENERGY =      -2.870162139

      CSF      COEF    OCCUPANCY (IGNORING CORE)
      ---      ----    --------- --------- -----
        1    0.997832  20
        3   -0.065749  02

STATE #    2  ENERGY =      -0.948712883

      CSF      COEF    OCCUPANCY (IGNORING CORE)
      ---      ----    --------- --------- -----
        2    0.998465  11
        3    0.054989  02
求助!哪里错了
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
qizhuixia: 金币+30, 很有道理 2019-03-14 23:07:00
这个涉及到基组的定义了 。楼主用的是 4-31G和6-31G 两种基组。这是价层劈裂基组。第一个数字代表非价层电子由多少个gaussian type function来表示。31G代表价层电子是由劈裂成两个基函数 分别由 3个 和 1个gaussian type function组成。 那么对于一般情况这俩基组算出来应该是不同的。但是这算的是He原子,He原子 就没有非价层电子 全都是价层电子。所以无论是4 还是6 都没有区别。只要横杠后面的31G一样 就不会有区别。 可以试试6-311G 算算看,那结果就应该不同。

PS. 你基组里用的NDFUNC=1 也没有用。这个命令是对重原子添加d极化,而所谓重原子是从Li 才开始的。对于 H 和 He加极化 是用 NPFUNC命令
2楼2019-03-14 22:53:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 qizhuixia 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿武理314求调剂 +5 ( ̄~ ̄;) 2026-03-08 6/300 2026-03-15 00:33 by wht0531
[考研] 梁成伟老师课题组欢迎你的加入 +5 一鸭鸭哟 2026-03-14 6/300 2026-03-14 22:29 by lllllcsjsj
[考研] 309求调剂 +4 花与叶@ 2026-03-10 4/200 2026-03-14 21:26 by a不易
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 理学化学调剂 +3 040517FYfy? 2026-03-09 3/150 2026-03-14 03:09 by JourneyLucky
[考研] 0703化学求调剂 +5 很老实人 2026-03-09 5/250 2026-03-14 02:57 by JourneyLucky
[考研] 295复试调剂 +5 简木ChuFront 2026-03-09 5/250 2026-03-14 01:29 by JourneyLucky
[考研] 265求调剂 +9 小木虫085600 2026-03-09 12/600 2026-03-14 01:11 by JourneyLucky
[考研] 311求调剂 +5 牛乳糖的卡卡 2026-03-10 5/250 2026-03-14 00:05 by JourneyLucky
[考研] 复试调剂 +9 Copy267 2026-03-10 9/450 2026-03-13 23:45 by userper
[考研] 0805,333求调剂 +3 112253525 2026-03-10 3/150 2026-03-13 23:42 by JourneyLucky
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
[考研] 293求调剂 +3 世界首富 2026-03-11 3/150 2026-03-13 16:27 by JourneyLucky
[考研] 328化工专硕求调剂 +4 。,。,。,。i 2026-03-12 4/200 2026-03-13 14:44 by JourneyLucky
[考研] 274求调剂0856材料化工 +12 z2839474511 2026-03-11 13/650 2026-03-13 10:39 by peike
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
[考研] 081200计算机科学与技术数一英一306求调剂 +3 intankt 2026-03-08 3/150 2026-03-11 10:50 by dandan413
[考研] 大连大学化学专业研究生调剂 +3 琪久. 2026-03-10 8/400 2026-03-11 10:02 by 琪久.
[考研] 家人们 调剂不迷路 看这里 +8 likeihood 2026-03-09 13/650 2026-03-10 08:09 by likeihood
[考研] 求0856调剂 +5 小力气珂珂 2026-03-08 5/250 2026-03-09 06:42 by houyaoxu
信息提示
请填处理意见