| 查看: 841 | 回复: 1 | ||
[求助]
小弟不才，不知道为什么用各种基组算出来的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 求助！哪里错了 |
回复此楼» 猜你喜欢
能源环境新材料研究团队 招收 新能源、冶金、材料、化工、环境 等工科类研究生
已经有0人回复
-
武汉纺织大学教育部重点实验室调剂
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有283人回复
-
理论与计算化学方向调剂——初试388分（数学144）
已经有1人回复
-
能源环境新材料研究团队 研究生调剂招生
已经有0人回复
-
齐齐哈尔大学材料科学与工程学院有调剂名额
已经有0人回复
-
苏州科技大学 物理学专业 接收调剂生
已经有0人回复
-
硕士研究生调剂
已经有0人回复
-
苏州科技大学 物理学专业 接收调剂生
已经有0人回复
-
苏州科技大学 物理学专业 接收调剂生
已经有0人回复
-
高效微纳米储能材料与产业化推进实验室 苏州 招收硕士研究生（含调剂）
已经有0人回复
paramecium86
版主 (著名写手)
- QC强帖: 5
- 应助: 588 (博士)
- 贵宾: 0.01
- 金币: 26424.6
- 红花: 349
- 帖子: 2026
- 在线: 1442.8小时
- 虫号: 2466443
- 注册: 2013-05-16
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
qizhuixia: 金币+30, 很有道理 2019-03-14 23:07:00
感谢参与，应助指数 +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
