24小时热门版块排行榜    

查看: 988  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-10-01 21:59 by Equinoxhua
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:52 by Z9YWQ5EAO3qp
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:45 by Z9YWQ5EAO3qp
[有机交流] 同一个分子,一条来自文献,一条来自AI——不告诉你答案,你会选哪条? +4 tianxiaoxian 2026-09-28 7/350 2026-09-30 08:34 by 88817753
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:37 by etmYJ6d2rquH
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:46 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-29 14:30 by etmYJ6d2rquH
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 IDs3scOF0tjC 2026-09-28 6/300 2026-09-29 14:13 by etmYJ6d2rquH
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-29 09:49 by JzYBbHIXWSrW
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
信息提示
请填处理意见