| 4 | 1/1 | 返回列表 |
| 查看: 1028 | 回復(fù): 3 | |||
qzhaosdu金蟲(chóng) (著名寫(xiě)手)
|
[交流]
【討論】關(guān)于G03里guess=only的一個(gè)疑惑。
|
|
G03里GUESS關(guān)鍵詞里說(shuō)到guess=alter的用二重態(tài)NH2舉例。 首先,用了guess=only得出初始猜測(cè),說(shuō)明書(shū)里是這樣的 Initial guess orbital symmetries. Alpha Orbitals: Occupied (A1) (A1) (B2) (B1) (A1) Virtual (A1) (B2) Beta Orbitals: Occupied (A1) (A1) (B2) (A1) Virtual (B1) (A1) (B2) 但是我算了一下卻是 Orbital symmetries: Alpha Orbitals: Occupied (A1) (A1) (B2) (A1) (B1) Virtual (A1) (B2) Beta Orbitals: Occupied (A1) (A1) (B2) (A1) Virtual (B1) (A1) (B2) 也就是說(shuō)說(shuō)明書(shū)里的占據(jù)軌道里4和5是反著的。 但是用了OPT UHF去做的時(shí)候,程序在算第一步SCF的時(shí)候卻給的是和說(shuō)明書(shū)一樣的結(jié)果。 正常來(lái)說(shuō),Alpha和Beta的對(duì)稱(chēng)性應(yīng)該是一一對(duì)應(yīng)相同的吧?(我也不知道這樣說(shuō)對(duì)不對(duì)?) ************************** 另外我剛看了下G09的說(shuō)明書(shū),GUESS里給的是正常的結(jié)果,也就是說(shuō)G03里寫(xiě)錯(cuò)了? ************************** 另外還有,我算了一個(gè)Ni+...C2H2的C2V結(jié)構(gòu),guess=only給的猜測(cè)是2-B2,但是b3lyp/6-31+g*算的時(shí)候給的初始猜測(cè)卻是2-B1. ************************** 所以說(shuō),G03里這個(gè)guess=only這個(gè)關(guān)鍵詞是不是有問(wèn)題?(起碼第一個(gè)SCF是有錯(cuò)誤的) [ Last edited by qzhaosdu on 2010-12-9 at 17:50 ] |
» 搶金幣啦!回帖就可以得到:
+1/469
+1/99
+1/87
+1/81
+1/76
+1/70
+1/67
+1/66
+2/50
+1/37
+1/36
+1/33
+1/30
+1/17
+1/17
+5/15
+1/5
+1/2
+1/1
+1/1
金蟲(chóng) (正式寫(xiě)手)
|
第一個(gè)問(wèn)題:不同版本的程序所用的初始猜測(cè)方法可能有所不同(即使方法相同也可能參數(shù)或算法有所不同)。不同的方法給出的初始猜測(cè)不一定相同,但只要方法合適,優(yōu)化后的結(jié)果應(yīng)該是一致的。 第二個(gè)問(wèn)題:“正常來(lái)說(shuō)占據(jù)軌道和非占據(jù)軌道的對(duì)稱(chēng)性應(yīng)該是一一對(duì)應(yīng)相同的吧?”這一理解是錯(cuò)誤的!首先因?yàn)檎紦?jù)軌道的數(shù)目與空軌道的數(shù)目一般是不同的,故不可能建立一一對(duì)應(yīng)關(guān)系。其次,即使二者的數(shù)目相同,也沒(méi)有任何物理或幾何原因能在二者之間建立這種關(guān)系。 第三個(gè)問(wèn)題:有些開(kāi)殼層體系特別是自由基,初始猜測(cè)的對(duì)稱(chēng)性常常是與基組有關(guān)的。原因在于你使用了單Slater行列式近似,而這種近似對(duì)于有些開(kāi)殼層體系并不是好的近似。 |
鐵桿木蟲(chóng) (著名寫(xiě)手)
新蟲(chóng) (正式寫(xiě)手)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 286求調(diào)劑 +3 | lemonzzn 2026-03-16 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 6/300 |
|
|
[考研] 材料專(zhuān)碩306英一數(shù)二 +9 | z1z2z3879 2026-03-16 | 12/600 |
|
|
[考研] 08工科 320總分 求調(diào)劑 +4 | 梨花珞晚風(fēng) 2026-03-17 | 4/200 |
|
|
[考研] 材料工程專(zhuān)碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 283求調(diào)劑 +3 | 聽(tīng)風(fēng)就是雨; 2026-03-16 | 3/150 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷,論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢(mèng)微 2026-03-15 | 3/150 |
|
|
[考研] 288求調(diào)劑 +4 | 奇點(diǎn)0314 2026-03-14 | 4/200 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
[考研] (081700)化學(xué)工程與技術(shù)-298分求調(diào)劑 +12 | 11啦啦啦 2026-03-11 | 35/1750 |
|
|
[考研] 315求調(diào)劑 +9 | 小羊小羊_ 2026-03-11 | 10/500 |
|
|
[考研] 310求調(diào)劑 +3 | 【上上簽】 2026-03-11 | 3/150 |
|
|
[考研] 土木第一志愿276求調(diào)劑,科研和技能十分豐富,求新興方向的導(dǎo)師收留 +3 | 土木小天才 2026-03-12 | 3/150 |
|
|
[考研] 085600材料與化工 309分請(qǐng)求調(diào)劑 +7 | dtdxzxx 2026-03-12 | 8/400 |
|
|
[考博] 福州大學(xué)楊黃浩課題組招收2026年專(zhuān)業(yè)學(xué)位博士研究生,2026.03.20截止 +3 | Xiangyu_ou 2026-03-12 | 3/150 |
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|