24小时热门版块排行榜

>论坛更新日志 (4320)
>文献求助 (323)
>虫友互识 (251)
>基金申请 (228)
>导师招生 (189)
>休闲灌水 (137)
>硕博家园 (136)
>论文投稿 (126)
>考博 (90)
>教师之家 (89)
>考研 (77)
>博后之家 (66)
>论文道贺祈福 (55)
>找工作 (49)
>公派出国 (44)
>绿色求助(高悬赏) (40)

返回列表

qizhuixia

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 91.5
散金: 20
帖子: 16
在线: 12.2小时
虫号: 12929110
注册: 2018-11-22

[求助] 小弟不才，不知道为什么用各种基组算出来的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
求助！哪里错了

回复此楼

» 猜你喜欢

湖南大学刘巧玲课题组2026年博士研究生招生信息已经有28人回复
南方科技大学表面催化反应动力学研究课题组招收2026届博士生已经有1人回复
物理化学论文润色/翻译怎么收费? 已经有141人回复
厦门大学侯旭团队招聘青年教师、博士后及博硕士研究生已经有2人回复
北京邮电大学物理学院老师招收博士生已经有3人回复
哈工大深圳表界面物理化学方向（胶体与界面，光谱探测研究）博士招生已经有0人回复
元素硫的化合价和氧化态已经有1人回复
江苏科技大学2026年招聘启事已经有1人回复
求助：三乙酰丙酮铁（Fe(acac)₃）标准红外(IR)与XRD谱图/数据已经有1人回复
Comsol 6.0 桶内肥料产热产湿过程模拟求助已经有1人回复
International Journal of Nanoscience属于中科院4区吗？已经有0人回复

1楼 2019-03-14 21:58:18

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

paramecium86

版主 (著名写手)

QC强帖: 5
应助: 588 (博士)
贵宾: 0.01
金币: 26684.1
红花: 349
帖子: 2027
在线: 1458.1小时
虫号: 2466443
注册: 2013-05-16
性别: GG
专业: 理论和计算化学
管辖: 量子化学

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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命令

赞一下(1人)