版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

qizhuixia

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 91.5
散金: 20
帖子: 16
在線: 12.2小時(shí)
蟲號(hào): 12929110
注冊(cè): 2018-11-22

[求助] 小弟不才，不知道為什么用各種基組算出來(lái)的CI的能量（氦）都是一樣的，求助已有1人參與

小弟不才，不知道為什么用各種基組算出來(lái)的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
~

算出來(lái)的能量都是樣的
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ò)了

回復(fù)此樓

» 猜你喜歡

QE計(jì)算電聲耦合的時(shí)候報(bào)錯(cuò)Error in routine lambda (100)wrong or too many modes 已經(jīng)有2人回復(fù)
cif求助下載已經(jīng)有1人回復(fù)
物理化學(xué)論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有242人回復(fù)
高效微納米儲(chǔ)能材料與產(chǎn)業(yè)化推進(jìn)實(shí)驗(yàn)室蘇州招收碩士研究生（含調(diào)劑）已經(jīng)有5人回復(fù)
求ENDNOTE破解軟件安裝包已經(jīng)有4人回復(fù)
高效微納米儲(chǔ)能材料與產(chǎn)業(yè)化推進(jìn)實(shí)驗(yàn)室蘇州招收碩士研究生（含調(diào)劑）已經(jīng)有0人回復(fù)
高效微納米儲(chǔ)能材料與產(chǎn)業(yè)化推進(jìn)實(shí)驗(yàn)室蘇州招收碩士研究生（含調(diào)劑）已經(jīng)有2人回復(fù)
化學(xué)化工學(xué)院招收化工、化學(xué)、材料等相關(guān)方向研究生（學(xué)碩、專碩都有調(diào)劑名額）已經(jīng)有0人回復(fù)
河海大學(xué) （211.雙一流高校）鈣鈦礦光電實(shí)驗(yàn)室丁勇團(tuán)隊(duì)歡迎有志青年！已經(jīng)有3人回復(fù)
348求調(diào)劑已經(jīng)有7人回復(fù)
黑龍江省原子能研究院化學(xué)課題組招收2026級(jí)化學(xué)專業(yè)或化學(xué)相關(guān)專業(yè)碩士生已經(jīng)有0人回復(fù)

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

paramecium86

版主 (著名寫手)

QC強(qiáng)帖: 5
應(yīng)助: 587 (博士)
貴賓: 0.01
金幣: 26376.6
紅花: 349
帖子: 2025
在線: 1441.9小時(shí)
蟲號(hào): 2466443
注冊(cè): 2013-05-16
性別: GG
專業(yè): 理論和計(jì)算化學(xué)
管轄: 量子化學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
qizhuixia: 金幣+30, 很有道理 2019-03-14 23:07:00

這個(gè)涉及到基組的定義了。樓主用的是 4-31G和6-31G 兩種基組。這是價(jià)層劈裂基組。第一個(gè)數(shù)字代表非價(jià)層電子由多少個(gè)gaussian type function來(lái)表示。31G代表價(jià)層電子是由劈裂成兩個(gè)基函數(shù) 分別由 3個(gè) 和 1個(gè)gaussian type function組成。那么對(duì)于一般情況這倆基組算出來(lái)應(yīng)該是不同的。但是這算的是He原子，He原子就沒(méi)有非價(jià)層電子全都是價(jià)層電子。所以無(wú)論是4 還是6 都沒(méi)有區(qū)別。只要橫杠后面的31G一樣就不會(huì)有區(qū)別。可以試試6-311G 算算看，那結(jié)果就應(yīng)該不同。

PS. 你基組里用的NDFUNC=1 也沒(méi)有用。這個(gè)命令是對(duì)重原子添加d極化,而所謂重原子是從Li 才開始的。對(duì)于 H 和 He加極化是用 NPFUNC命令

贊一下(1人)

回復(fù)此樓

2樓2019-03-14 22:53:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 qizhuixia 的主題更新