| 2 | 1/1 | 返回列表 |
| 查看: 827 | 回復(fù): 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 求助!哪里錯(cuò)了 |
版主 (著名寫手)
|
這個(gè)涉及到基組的定義了 。樓主用的是 4-31G和6-31G 兩種基組。這是價(jià)層劈裂基組。第一個(gè)數(shù)字代表非價(jià)層電子由多少個(gè)gaussian type function來表示。31G代表價(jià)層電子是由劈裂成兩個(gè)基函數(shù) 分別由 3個(gè) 和 1個(gè)gaussian type function組成。 那么對(duì)于一般情況這倆基組算出來應(yīng)該是不同的。但是這算的是He原子,He原子 就沒有非價(jià)層電子 全都是價(jià)層電子。所以無論是4 還是6 都沒有區(qū)別。只要橫杠后面的31G一樣 就不會(huì)有區(qū)別。 可以試試6-311G 算算看,那結(jié)果就應(yīng)該不同。 PS. 你基組里用的NDFUNC=1 也沒有用。這個(gè)命令是對(duì)重原子添加d極化,而所謂重原子是從Li 才開始的。對(duì)于 H 和 He加極化 是用 NPFUNC命令 |
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085602化工求調(diào)劑 +6 | 董boxing 2026-03-10 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 一志愿東北大學(xué)080500,英一數(shù)二343分求調(diào)劑 +8 | Seele_v5 2026-03-04 | 8/400 |
|
|
[考研] 327分求調(diào)劑086 +4 | 西紅柿?小帥 2026-03-09 | 7/350 |
|
|
[考研] 物理學(xué)求調(diào)劑 +4 | sx宋 2026-03-05 | 4/200 |
|
|
[考研] 320求調(diào)劑 +4 | 魏zy 2026-03-08 | 4/200 |
|
|
[考研] 0701-322 求調(diào)劑 +3 | jiliuxian 2026-03-06 | 8/400 |
|
|
[考研] 288求調(diào)劑085600材料與化工 +13 | Daunrin 2026-03-07 | 15/750 |
|
|
[考研] 求調(diào)劑,不管什么專業(yè),我是可塑造的人才一枚,希望遇到知己老師撈撈我 +4 | 13102137290 2026-03-06 | 5/250 |
|
|
[考研] 0307化學(xué)求調(diào)劑 +6 | 0ok0k 2026-03-07 | 6/300 |
|
|
[考研] 一志愿武漢理工085601,初試301分,請(qǐng)問能調(diào)劑到湖北嗎 +3 | 肖yang 2026-03-06 | 3/150 |
|
|
[考研] 085600材料調(diào)劑 總分330 +6 | 池池丶 2026-03-03 | 6/300 |
|
|
[考研] 化學(xué)專業(yè)調(diào)劑 +5 | 好好好1233 2026-03-04 | 6/300 |
|
|
[考研] 274求調(diào)劑 +9 | 一個(gè)學(xué)習(xí)者 2026-03-04 | 9/450 |
|
|
[考研] 材料085601一志愿哈工大317 +4 | 壓迫感行 2026-03-04 | 4/200 |
|
|
[考研] 347求調(diào)劑 +6 | 啊歐歐歐 2026-03-03 | 8/400 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-04 | 3/150 |
|
|
[論文投稿]
100+4
|
Stray2021 2026-03-03 | 4/200 |
|
|
[考研] 376求調(diào)劑 +3 | 王浩然sam 2026-03-04 | 3/150 |
|
|
[考研] 一志愿武漢理工大學(xué)-085602-總分296分-求調(diào)劑 +7 | 紫川葡柚 2026-03-04 | 7/350 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 將未有期 2026-03-04 | 3/150 |
|